Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “consistency”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Thermodynamically consistent algorithms for models of incompressible multiphase polymer solutions with a variable mobility

Here we present a general strategy for developing structure and property preserving numerical algorithms for thermodynamically consistent models of incompressible multiphase polymer solutions with a variable mobility. We first present a formalism to derive thermodynamically consistent, incompressible, multiphase polymer models. Then, we develop the general strategy, known as the supplementary variable method, to devise thermodynamically consistent numerical approximations to the models. We illustrate the numerical strategy using newly developed models of incompressible diblock copolymer solutions coupled with an electric and a magnetic field, respectively. Mesh refinement is conducted to verify convergence rates of the developed schemes. Some numerical examples are given to exhibit underlying dynamics absent from and driven by the external fields, respectively, highlighting differences between models with the variable and constant mobilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A thermodynamically consistent discretization of 1D thermal-fluid models using their metriplectic 4-bracket structure

Thermodynamically consistent models in continuum physics, i.e. models which satisfy the first and second laws of thermodynamics, may be expressed using the metriplectic formalism. In this work, we leverage the structures underlying this modeling formalism to preserve thermodynamic consistency in discretizations of a fluid model. The procedure relies (1) on ensuring that the spatial semi-discretization retains certain symmetries and degeneracies of the Poisson and metriplectic 4-brackets, and (2) on the use of an appropriate energy conserving time-stepping method. Here, the minimally simple yet nontrivial example of a one-dimensional thermal-fluid model is treated. It is found that preservation of the requisite symmetries and degeneracies of the 4-bracket is relatively simple to ensure in Galerkin spatial discretizations, suggesting a path forward for thermodynamically consistent discretizations of more complex fluid models using more specialized Galerkin methods.

Hamiltonian structure↗
DOE OSTI2021

A consistent and conservative model and its scheme for N -phase- M -component incompressible flows

Here, we propose a consistent and conservative model for multiphase and multicomponent incompressible flows, where there can be arbitrary numbers of phases and components. Each phase has a background fluid called the pure phase, each pair of phases is immiscible, and components are dissolvable in some specific phases. The model is developed based on the multiphase Phase-Field model including the contact angle boundary condition, the diffuse domain approach, and the analyses on the proposed consistency conditions for multiphase and multicomponent flows. The model conserves the mass of individual pure phases, the amount of each component in its dissolvable region, and thus the mass of the fluid mixture, and the momentum of the flow. It ensures that no fictitious phases or components can be generated and that the summation of the volume fractions from the Phase-Field model is unity everywhere so that there is no local void or overfilling. It satisfies a physical energy law and it is Galilean invariant. A corresponding numerical scheme is developed for the proposed model, whose formal accuracy is 2nd-order in both time and space. It is shown to be consistent and conservative and its solution is demonstrated to preserve the Galilean invariance and energy law. Numerical tests indicate that the proposed model and scheme are effective and robust to study various challenging multiphase and multicomponent flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Statistical dynamics of classical systems: A self-consistent field approach

We develop a self-consistent field theory for particle dynamics by extremizing the functional integral representation of a microscopic Langevin equation with respect to the collective fields. Although our approach is general, here we formulate it in the context of polymer dynamics to highlight satisfying formal analogies with equilibrium self-consistent field theory. An exact treatment of the dynamics of a single chain in a mean force field emerges naturally via a functional Smoluchowski equation, while the time-dependent monomer density and mean force field are determined self-consistently. As a simple initial demonstration of the theory, leaving an application to polymer dynamics for future work, we examine the dynamics of trapped interacting Brownian particles. For binary particle mixtures, we observe the kinetics of phase separation.

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

Frank–Kasper Phases of Diblock Copolymer Melts: Self-Consistent Field Results of Two Commonly Used Models

We constructed phase diagrams of conformationally asymmetric diblock copolymer A-B melts using the polymer self-consistent field (SCF) calculations of both the dissipative particle dynamics chain (DPDC) model (i.e., compressible melts of discrete Gaussian chains with the DPD non-bonded potential) and the “standard” model (i.e., incompressible melts of continuous Gaussian chains with the Dirac δ-function non-bonded potential) in the χN-ε plane, where χN and ε characterize, respectively, the repulsion and conformational asymmetry between the A and B blocks, at the A-block volume fraction f = 0.2 and 0.3. Consistent with previous SCF calculations of the “standard” model, σ and A15 are the only stable Frank–Kasper (FK) phases among the five FK (i.e., σ, A15, C14, C15 and Z) phases considered. The stability of σ and A15 is due to their delicate balance between the energetic and entropic contributions to the Helmholtz free energy per chain of the system, which, within our parameter range, increases in the order of σ/A15, Z, and C14/C15. While in general the SCF phase diagrams of these two models are qualitatively consistent, A15 is not stable for the DPDC model at the copolymer chain length N = 10 and f = 0.3; any differences in the SCF phase diagrams are solely due to the differences between these two models.

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

Jacobian consistency of a one-parametric class of smoothing Fischer–Burmeister functions for SOCCP

The Jacobian consistency of smoothing functions plays an important role for achieving the rapid convergence of Newton methods or Newton-like methods with an appropriate parameter control. In this paper, we study the properties, derive the computable formula for the Jacobian matrix and prove the Jacobian consistency of a one-parametric class of smoothing Fischer–Burmeister functions for second-order cone complementarity problems proposed by Tang et al. (Comput Appl Math 33:655–669, 2014). Then we apply its Jacobian consistency to a smoothing Newton method with the appropriate parameter control presented by Chen et al. (Math Comput 67:519–540, 1998), and show the global convergence and local quadratic convergence of the algorithm for solving the SOCCP under rather weak assumptions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Some comments on thermodynamic consistency for equilibrium mixture equations of state

We investigate sufficient conditions for thermodynamic consistency for equilibrium mixtures. Such models assume that the mass fraction average of the material component equations of state, when closed by a suitable equilibrium condition, provide a composite equation of state for the mixture. Here, we show that the two common equilibrium models of component pressure/temperature equilibrium and volume/temperature equilibrium (Dalton, 1808) define thermodynamically consistent mixture equations of state and that other equilibrium conditions can be thermodynamically consistent provided appropriate values are used for the mixture specific entropy and pressure.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

A consistent and conservative scheme for MHD flows with complex boundaries on an unstructured Cartesian adaptive system

The numerical simulation of Magnetohydrodynamics (MHD) flows with complex boundaries has been a topic of great interest in the development of a fusion reactor blanket for the difficulty to accurately simulate the Hartmann layers and side layers along arbitrary geometries. An adaptive version of a consistent and conservative scheme has been developed for simulating the MHD flows. Besides, the present study forms the first attempt to apply the cut-cell approach for irregular wall-bounded MHD flows, which is more flexible and conveniently implemented under adaptive mesh refinement (AMR) technique. It employs a Volume-of-Fluid (VOF) approach to represent the fluid–conducting wall interface that makes it possible to solve the fluid–solid coupling magnetic problems, emphasizing at how electric field solver is implemented when conductivity is discontinuous in cut-cell. For the irregular cut-cells, the conservative interpolation technique is applied to calculate the Lorentz force at cell-center. On the other hand, it will be shown how consistent and conservative scheme is implemented on fine/coarse mesh boundaries when using AMR technique. Then, the applied numerical schemes are validated by five test simulations and excellent agreement was obtained for all the cases considered, simultaneously showed good consistency and conservative properties.

42 ENGINEERING↗
DOE OSTI2016

Evaluating the consistency of location of the most severe acute skin reaction and highest skin dose measured by thermoluminescent dosimeter during radiotherapy for breast cancer

We conducted this prospective study to evaluate whether the location of the most severe acute skin reaction matches the highest skin dose measured by thermoluminescent dosimeter (TLD) during adjuvant radiotherapy (RT) for patients with breast cancer after breast conservative surgery. To determine whether TLD measurement can reflect the location of the most severe acute skin reaction, 80 consecutive patients were enrolled in this prospective study. We divided the irradiated field into breast, axillary, inframammary fold, and areola/nipple areas. In 1 treatment session when obvious skin reaction occurred, we placed the TLD chips onto the 4 areas and measured the skin dose. We determined whether the highest measured skin dose area is consistent with the location of the most severe skin reaction. The McNemar test revealed that the clinical skin reaction and TLD measurement are more consistent when the most severe skin reaction occurred at the axillary area, and the p = 0.0108. On the contrary, TLD measurement of skin dose is less likely consistent with clinical observation when the most severe skin reaction occurred at the inframammary fold, breast, and areola/nipple areas (all the p > 0.05). Considering the common site of severe skin reaction over the axillary area, TLD measurement may be an appropriate way to predict skin reaction during RT.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2023

Probing the consistency of cosmological contours for supernova cosmology

As the scale of cosmological surveys increases, so does the complexity in the analyses. This complexity can often make it difficult to derive the underlying principles, necessitating statistically rigorous testing to ensure the results of an analysis are consistent and reasonable. This is particularly important in multi-probe cosmological analyses like those used in the Dark Energy Survey (DES) and the upcoming Legacy Survey of Space and Time, where accurate uncertainties are vital. In this paper, we present a statistically rigorous method to test the consistency of contours produced in these analyses and apply this method to the Pippin cosmological pipeline used for type Ia supernova cosmology with the DES. We make use of the Neyman construction, a frequentist methodology that leverages extensive simulations to calculate confidence intervals, to perform this consistency check. A true Neyman construction is too computationally expensive for supernova cosmology, so we develop a method for approximating a Neyman construction with far fewer simulations. We find that for a simulated dataset, the 68% contour reported by the Pippin pipeline and the 68% confidence region produced by our approximate Neyman construction differ by less than a percent near the input cosmology; however, they show more significant differences far from the input cosmology, with a maximal difference of 0.05 in $Ω$ M and 0.07 in w. In conclusion, this divergence is most impactful for analyses of cosmological tensions, but its impact is mitigated when combining supernovae with other cross-cutting cosmological probes, such as the cosmic microwave background.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Lowering the Scaling of Self-Consistent Field Methods by Combining Tensor Hypercontraction and a Density Difference Ansatz

We present the tensor hypercontraction difference self-consistent field (SCF) method, an approach that reduces the formal computational scaling of traditional naive self-consistent field methods from 𝑂(𝑁 4 ) to 𝑂(𝑁 3 ) with system size 𝑁. The scaling reduction is achieved by developing a new technique for constructing the tensor hypercontraction decomposition based on the fundamental approximation made in density fitting. Combining this scheme with the difference self-consistent field methodology, we achieve a method that enables 𝑂(𝑁 3 ) scaling SCF calculations with only 𝑂(𝑁 2 ) storage requirements. In conclusion, our proof-of-concept numerical tests demonstrate robust performance with errors in total energies below 8 × 10 –4 E h and with sub 1 kcal/mol errors for relative energies.

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

A Coverage Self-Consistent Microkinetic Model for Vapor-Phase Formic Acid Decomposition over Pd/C Catalysts

An iterative approach utilizing density functional theory (DFT, PW91-GGA)-informed mean-field microkinetic models and reaction kinetics experiments is used to determine the reaction mechanism and the active site for formic acid (HCOOH, FA) decomposition over a Pd/C catalyst. Models parametrized using DFT energetics on clean Pd(100) and Pd(111) required large corrections to the DFT energetics for capturing our experimental data. Further, both Pd(111) and Pd(100) models predicted a high coverage of adsorbed CO (CO*), inconsistent with the assumption of a clean surface at which the rate parameters for these models were calculated. To better represent the active site under reaction conditions and explicitly account for the presence of CO*, subsequent microkinetic models were formulated using DFT energetics that were calculated on partially (5/9 ML) CO*-covered Pd (111) and (100) facets. Upon parameter adjustment, the resultant 5/9 ML CO*-covered Pd(100) model, although consistent in terms of CO* coverage, was unable to capture the dehydration path measured in the experiments and was, therefore, deemed not to offer an accurate representation of the active site for FA decomposition over Pd/C. In contrast, a partially CO*-covered Pd(111) model was better at representing the catalytic active site, as in addition to being consistent in terms of CO* coverages, it required small adjustments of the DFT parameters to accurately capture the experimental data set (both dehydrogenation and dehydration). Our results suggest that the reaction occurs via the spectroscopically elusive carboxyl (COOH*) intermediate and that spectator CO*-assisted decomposition pathways play an important role under typical experimental conditions. In addition, our study highlights the importance of striving for coverage self-consistent microkinetic models and for including spectator-assisted mechanisms in order to develop an improved picture of the active site under reaction conditions.

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

Correlation consistent basis sets for the atoms In–Xe

In this work, the correlation consistent family of Gaussian basis sets has been expanded to include all-electron basis sets for In–Xe. The methodology for developing these basis sets is described, and several examples of the performance and utility of the new sets have been provided. Dissociation energies and bond lengths for both homonuclear and heteronuclear diatomics demonstrate the systematic convergence behavior with respect to increasing basis set quality expected by the family of correlation consistent basis sets in describing molecular properties. Comparison with recently developed correlation consistent sets designed for use with the Douglas-Kroll Hamiltonian is provided.

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

The breaks and the hidden components in the power-law spectra of synchrotron radiation of the self-consistent current structures

Widespread use of a broken-power-law description of the spectra of synchrotron emission of various plasma objects requires an analysis of origin and a proper interpretation of spectral components. We show that, for a self-consistent magnetic configuration in a collisionless plasma, these components may be angle-dependent according to an anisotropic particle momentum distribution and may have no counterparts in a particle energy distribution. That has never been studied analytically and is in contrast to a usual model of synchrotron radiation, assuming an external magnetic field and a particle ensemble with isotropic momentum distribution. We demonstrate that for the wide intervals of observation angle the power-law spectra and, in particular, the positions and number of spectral breaks may be essentially different for the cases of the self-consistent and not-self-consistent magnetic fields in current structures responsible for the synchrotron radiation of the ensembles of relativistic particles with the multi-power-law energy distributions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Repartitioned Brillouin-Wigner perturbation theory with a size-consistent second-order correlation energy

Second-order Møller-Plesset perturbation theory (MP2) often breaks down catastrophically in small-gap systems, leaving much to be desired in its performance for myriad chemical applications such as noncovalent interactions, thermochemistry, and dative bonding in transition metal complexes. This divergence problem has reignited interest in Brillouin-Wigner perturbation theory (BWPT), which is regular at all orders but lacks size consistency and extensivity, severely limiting its application to chemistry. In this work, we propose an alternative partitioning of the Hamiltonian that leads to a regular BWPT perturbation series that, through the second order, is size-extensive, size-consistent (provided its Hartree–Fock reference is also), and orbital invariant. Here, our second-order size-consistent Brillouin-Wigner (BW-s2) approach can describe the exact dissociation limit of H 2 in a minimal basis set, regardless of the spin polarization of the reference orbitals. More broadly, we find that BW-s2 offers improvements relative to MP2 for covalent bond breaking, noncovalent interaction energies, and metal/organic reaction energies, although rivaling coupled-cluster with single and double substitutions for thermochemical properties.

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

Polarization consistent dielectric screening in polarizable continuum model calculations of solvation energies

A polarization consistent framework, where dielectric screening is affected consistently in polarizable continuum model (PCM) calculations, is employed for the study of solvation energies. The computational framework combines a screened range-separated-hybrid functional (SRSH) with PCM calculations, SRSH-PCM, where dielectric screening is imposed in both PCM self-consistent reaction field (SCRF) iterations and the electronic structure Hamiltonian. We begin by demonstrating the impact of modifying the Hamiltonian to include such dielectric screening in SCRF iterations by considering the solutions of electrostatically embedded Hartree–Fock (HF) exact exchange equations. Long-range screened HF-PCM calculations are shown to capture properly the linear dependence of gap energy of frontier orbitals on the inverse of the dielectric constant, whereas unscreened HF-PCM orbital energies are fallaciously semi-constant with respect to the dielectric constant and, therefore, inconsistent with the ionization energy gaps. Similar trends affect density functional theory (DFT) calculations that aim to achieve predictive quality. Importantly, the dielectric screened calculations are shown to significantly affect DFT- and HF PCM-based solvation energies, where screened solvation energies are smaller compared to the unscreened values. Importantly, SRSH-PCM, therefore, appears to reduce the tendency of DFT-PCM to overestimate solvation energies, where we find the effect to increase with the dielectric constant and the polarity of the molecular solute, trends that enhance the quality of DFT-PCM calculations of solvation energy. Understanding the relationship of dielectric screening in the Hamiltonian and DFT-PCM calculations can ultimately benefit on-going efforts for the design of predictive and parameter free descriptions of solvation energies.

Chemistry↗
DOE OSTI2025

Thermodynamic consistency and fluctuations in mesoscopic stochastic simulations of reactive gas mixtures

It is essential that mesoscopic simulations of reactive systems reproduce the correct statistical distributions at thermodynamic equilibrium. By considering a compressible fluctuating hydrodynamics (FHD) simulation method of ideal gas mixtures undergoing reversible reactions described by the chemical Langevin equations, we show that thermodynamic consistency in reaction rates and the use of instantaneous temperatures for the evaluation of reaction rates is required for fluctuations for the overall system to be correct. We then formulate the required properties of a thermodynamically consistent reaction (TCR) model. As noted in the literature, while reactions are often discussed in terms of forward and reverse rates, these rates should not be modeled independently because they must be compatible with thermodynamic equilibrium for the system. Using a simple TCR model where each chemical species has constant heat capacity, we derive the explicit condition that the forward and reverse reaction rate constants must satisfy in order for the system to be thermodynamically consistent. We perform equilibrium and non-equilibrium simulations of ideal gas mixtures undergoing a reversible dimerization reaction to measure the fluctuational behavior of the system numerically. We confirm that FHD simulations with the TCR model give the correct static structure factor of equilibrium fluctuations. For the statistically steady simulation of a gas mixture between two isothermal walls with different temperatures, we show using the TCR model that the temperature variance agrees with the corresponding thermodynamic-equilibrium temperature variance in the interior of the system, whereas noticeable deviations are present in regions near walls, where chemistry is far from equilibrium.

Polimeno, Matteo [University of California, Merced↗
DOE OSTI2020

Improved surrogates in inertial confinement fusion with manifold and cycle consistencies

Neural networks have become the method of choice in surrogate modeling because of their ability to characterize arbitrary, high-dimensional functions in a data-driven fashion. This paper advocates for the training of surrogates that are 1) consistent with the physical manifold, resulting in physically meaningful predictions, and 2) cyclically consistent with a jointly trained inverse model; i.e., backmapping predictions through the inverse results in the original input parameters. We find that these two consistencies lead to surrogates that are superior in terms of predictive performance, are more resilient to sampling artifacts, and tend to be more data efficient. Using inertial confinement fusion (ICF) as a test-bed problem, we model a one-dimensional semianalytic numerical simulator and demonstrate the effectiveness of our approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗