Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Formal Methods”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

A semi-Lagrangian transport method for kinetic problems with application to dense-to-dilute polydisperse reacting spray flows

For sprays, as described by a kinetic disperse phase model strongly coupled to the Navier–Stokes equations, the resolution strategy is constrained by accuracy objectives, robustness needs, and the computing architecture. In order to leverage the good properties of the Eulerian formalism, we introduce a deterministic particle-based numerical method to solve transport in physical space, which is simple to adapt to the many types of closures and moment systems. The method is inspired by the semi-Lagrangian schemes, developed for Gas Dynamics. We show how semi-Lagrangian formulations are relevant for a disperse phase far from equilibrium and where the particle–particle coupling barely influences the transport; i.e., when particle pressure is negligible. The particle behavior is indeed close to free streaming. The new method uses the assumption of parcel transport and avoids to compute fluxes and their limiters, which makes it robust. It is a deterministic resolution method so that it does not require efforts on statistical convergence, noise control, or post-processing. All couplings are done among data under the form of Eulerian fields, which allows one to use efficient algorithms and to anticipate the computational load. This makes the method both accurate and efficient in the context of parallel computing. After a complete verification of the new transport method on various academic test cases, we demonstrate the overall strategy's ability to solve a strongly-coupled liquid jet with fine spatial resolution and we apply it to the case of high-fidelity Large Eddy Simulation of a dense spray flow. A fuel spray is simulated after atomization at Diesel engine combustion chamber conditions. The large, parallel, strongly coupled computation proves the efficiency of the method for dense, polydisperse, reacting spray flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Simulation-Based Valuation of Transactive Energy Systems

Transactive Energy (TE) has been recognized as a promising technique for integrating responsive loads and distributed energy resources as well as advancing grid modernization. To help the industry better understand the value of TE and compare different TE schemes in a systematic and transparent manner, a comprehensive simulation-based TE valuation method is developed. The method has the following salient features: 1) it formally defines the valuation scenarios, use cases, baseline and valuation metrics; 2) an open-source simulation platform for transactive energy systems has been developed by integrating transmission, distribution and building simulators, and plugin TE and non-TE agents through the Framework for Network Co-Simulation (FNCS); 3) transparency and flexibility of the valuation is enhanced through separation of simulation and valuation, base valuation metrics and final valuation metrics. In conclusion, a valuation example based on the Smart Grid Interoperability Panel (SGIP) Use Case 1 is provided to demonstrate the developed TE simulation program and the valuation method.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2016

On Improving Analytical Models of Cosmic Reionization for Matching Numerical Simulations

The methods for studying the epoch of cosmic reionization vary from full radiative transfer simulations to purely analytical models. While numerical approaches are computationally expensive and are not suitable for generating many mock catalogs, analytical methods are based on assumptions and approximations. We explore the interconnection between both methods. First, we ask how the analytical framework of excursion set formalism can be used for statistical analysis of numerical simulations and visual representation of the morphology of ionization fronts. Second, we explore the methods of training the analytical model on a given numerical simulation. We present a new code which emerged from this study. Its main application is to match the analytical model with a numerical simulation. Then, it allows one to generate mock reionization catalogs with volumes exceeding the original simulation quickly and computationally inexpensively, meanwhile reproducing large scale statistical properties. These mock catalogs are particularly useful for CMB polarization and 21cm experiments, where large volumes are required to simulate the observed signal.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

ON IMPROVING ANALYTICAL MODELS OF COSMIC REIONIZATION FOR MATCHING NUMERICAL SIMULATION

The methods for studying the epoch of cosmic reionization vary from full radiative transfer simulations to purely analytical models. While numerical approaches are computationally expensive and are not suitable for generating many mock catalogs, analytical methods are based on assumptions and approximations. We explore the interconnection between both methods. First, we ask how the analytical framework of excursion set formalism can be used for statistical analysis of numerical simulations and visual representation of the morphology of ionization fronts. Second, we explore the methods of training the analytical model on a given numerical simulation. We present a new code which emerged from this study. Its main application is to match the analytical model with a numerical simulation. Then, it allows one to generate mock reionization catalogs with volumes exceeding the original simulation quickly and computationally inexpensively, meanwhile reproducing large-scale statistical properties. These mock catalogs are particularly useful for cosmic microwave background polarization and 21 cm experiments, where large volumes are required to simulate the observed signal.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

The many-body Wigner Monte Carlo method for time-dependent ab-initio quantum simulations

The aim of ab-initio approaches is the simulation of many-body quantum systems from the first principles of quantum mechanics. These methods are traditionally based on the many-body Schrödinger equation which represents an incredible mathematical challenge. In this paper, we introduce the many-body Wigner Monte Carlo method in the context of distinguishable particles and in the absence of spin-dependent effects. Despite these restrictions, the method has several advantages. First of all, the Wigner formalism is intuitive, as it is based on the concept of a quasi-distribution function. Secondly, the Monte Carlo numerical approach allows scalability on parallel machines that is practically unachievable by means of other techniques based on finite difference or finite element methods. Finally, this method allows time-dependent ab-initio simulations of strongly correlated quantum systems. In order to validate our many-body Wigner Monte Carlo method, as a case study we simulate a relatively simple system consisting of two particles in several different situations. We first start from two non-interacting free Gaussian wave packets. We, then, proceed with the inclusion of an external potential barrier, and we conclude by simulating two entangled (i.e. correlated) particles. The results show how, in the case of negligible spin-dependent effects, the many-body Wigner Monte Carlo method provides an efficient and reliable tool to study the time-dependent evolution of quantum systems composed of distinguishable particles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Sparse maps—A systematic infrastructure for reduced-scaling electronic structure methods. I. An efficient and simple linear scaling local MP2 method that uses an intermediate basis of pair natural orbitals

In this work, a systematic infrastructure is described that formalizes concepts implicit in previous work and greatly simplifies computer implementation of reduced-scaling electronic structure methods. The key concept is sparse representation of tensors using chains of sparse maps between two index sets. Sparse map representation can be viewed as a generalization of compressed sparse row, a common representation of a sparse matrix, to tensor data. By combining few elementary operations on sparse maps (inversion, chaining, intersection, etc.), complex algorithms can be developed, illustrated here by a linear-scaling transformation of three-center Coulomb integrals based on our compact code library that implements sparse maps and operations on them. The sparsity of the three-center integrals arises from spatial locality of the basis functions and domain density fitting approximation. A novel feature of our approach is the use of differential overlap integrals computed in linear-scaling fashion for screening products of basis functions. Finally, a robust linear scaling domain based local pair natural orbital second-order Möller-Plesset (DLPNO-MP2) method is described based on the sparse map infrastructure that only depends on a minimal number of cutoff parameters that can be systematically tightened to approach 100% of the canonical MP2 correlation energy. With default truncation thresholds, DLPNO-MP2 recovers more than 99.9% of the canonical resolution of the identity MP2 (RI-MP2) energy while still showing a very early crossover with respect to the computational effort. Based on extensive benchmark calculations, relative energies are reproduced with an error of typically <0.2 kcal/mol. The efficiency of the local MP2 (LMP2) method can be drastically improved by carrying out the LMP2 iterations in a basis of pair natural orbitals. While the present work focuses on local electron correlation, it is of much broader applicability to computation with sparse tensors in quantum chemistry and beyond.

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

Noiseless Vlasov–Poisson simulations with linearly transformed particles

We introduce a deterministic discrete-particle simulation approach, the Linearly-Transformed Particle-In-Cell (LTPIC) method, that employs linear deformations of the particles to reduce the noise traditionally associated with particle schemes. Formally, transforming the particles is justified by local first order expansions of the characteristic flow in phase space. In practice the method amounts of using deformation matrices within the particle shape functions; these matrices are updated via local evaluations of the forward numerical flow. Because it is necessary to periodically remap the particles on a regular grid to avoid excessively deforming their shapes, the method can be seen as a development of Denavit's Forward Semi-Lagrangian (FSL) scheme (Denavit, 1972 [8]). However, it has recently been established (Campos Pinto, 2012 [20]) that the underlying Linearly-Transformed Particle scheme converges for abstract transport problems, with no need to remap the particles; deforming the particles can thus be seen as a way to significantly lower the remapping frequency needed in the FSL schemes, and hence the associated numerical diffusion. To couple the method with electrostatic field solvers, two specific charge deposition schemes are examined, and their performance compared with that of the standard deposition method. Finally, numerical 1d1v simulations involving benchmark test cases and halo formation in an initially mismatched thermal sheet beam demonstrate some advantages of our LTPIC scheme over the classical PIC and FSL methods. Benchmarked test cases also indicate that, for numerical choices involving similar computational effort, the LTPIC method is capable of accuracy comparable to or exceeding that of state-of-the-art, high-resolution Vlasov schemes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

A feasibility study of dynamic adaptive radiotherapy for nonsmall cell lung cancer

Purpose: The final state of the tumor at the end of a radiotherapy course is dependent on the doses given in each fraction during the treatment course. This study investigates the feasibility of using dynamic adaptive radiotherapy (DART) in treating lung cancers assuming CBCT is available to observe midtreatment tumor states. DART adapts treatment plans using a dynamic programming technique to consider the expected changes of the tumor in the optimization process. Methods: DART is constructed using a stochastic control formalism framework. It minimizes the total expected number of tumor cells at the end of a treatment course, which is equivalent to maximizing tumor control probability, subject to the uncertainty inherent in the tumor response. This formulation allows for nonstationary dose distributions as well as nonstationary fractional doses as needed to achieve a series of optimal plans that are conformal to the tumor over time, i.e., spatiotemporally optimal plans. Sixteen phantom cases with various sizes and locations of tumors and organs-at-risk (OAR) were generated using in-house software. Each case was planned with DART and conventional IMRT prescribing 60 Gy in 30 fractions. The observations of the change in the tumor volume over a treatment course were simulated using a two-level cell population model. Monte Carlo simulations of the treatment course for each case were run to account for uncertainty in the tumor response. The same OAR dose constraints were applied for both methods. The frequency of replanning was varied between 1, 2, 5 (weekly), and 29 times (daily). The final average tumor dose and OAR doses have been compared to quantify the potential dosimetric benefits of DART. Results: The average tumor max, min, mean, and D95 doses using DART relative to these using conventional IMRT were 124.0%–125.2%, 102.1%–114.7%, 113.7%–123.4%, and 102.0%–115.9% (range dependent on the frequency of replanning). The average relative maximum doses for the cord and esophagus, mean doses for the heart and lungs, and D05 for the unspecified tissue resulting 84%–102.4%, 99.8%–106.9%, 66.9%–85.6%, 58.2%–78.8%, and 85.2%–94.0%, respectively. Conclusions: It is feasible to apply DART to the treatment of NSCLC using CBCT to observe the midtreatment tumor state. Potential increases in the tumor dose and reductions in the OAR dose, particularly for parallel OARs with mean or dose–volume constraints, could be achieved using DART compared to nonadaptive IMRT.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Limitations of the TG-43 formalism for skin high-dose-rate brachytherapy dose calculations

Purpose: In skin high-dose-rate (HDR) brachytherapy, sources are located outside, in contact with, or implanted at some depth below the skin surface. Most treatment planning systems use the TG-43 formalism, which is based on single-source dose superposition within an infinite water medium without accounting for the true geometry in which conditions for scattered radiation are altered by the presence of air. The purpose of this study is to evaluate the dosimetric limitations of the TG-43 formalism in HDR skin brachytherapy and the potential clinical impact. Methods: Dose rate distributions of typical configurations used in skin brachytherapy were obtained: a 5 cm × 5 cm superficial mould; a source inside a catheter located at the skin surface with and without backscatter bolus; and a typical interstitial implant consisting of an HDR source in a catheter located at a depth of 0.5 cm. Commercially available HDR{sup 60}Co and {sup 192}Ir sources and a hypothetical {sup 169}Yb source were considered. The Geant4 Monte Carlo radiation transport code was used to estimate dose rate distributions for the configurations considered. These results were then compared to those obtained with the TG-43 dose calculation formalism. In particular, the influence of adding bolus material over the implant was studied. Results: For a 5 cm × 5 cm{sup 192}Ir superficial mould and 0.5 cm prescription depth, dose differences in comparison to the TG-43 method were about −3%. When the source was positioned at the skin surface, dose differences were smaller than −1% for {sup 60}Co and {sup 192}Ir, yet −3% for {sup 169}Yb. For the interstitial implant, dose differences at the skin surface were −7% for {sup 60}Co, −0.6% for {sup 192}Ir, and −2.5% for {sup 169}Yb. Conclusions: This study indicates the following: (i) for the superficial mould, no bolus is needed; (ii) when the source is in contact with the skin surface, no bolus is needed for either {sup 60}Co and {sup 192}Ir. For lower energy radionuclides like {sup 169}Yb, bolus may be needed; and (iii) for the interstitial case, at least a 0.1 cm bolus is advised for {sup 60}Co to avoid underdosing superficial target layers. For {sup 192}Ir and {sup 169}Yb, no bolus is needed. For those cases where no bolus is needed, its use might be detrimental as the lack of radiation scatter may be beneficial to the patient, although the 2% tolerance for dose calculation accuracy recommended in the AAPM TG-56 report is not fulfilled.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2017

First-principles Green-Kubo method for thermal conductivity calculations

We present a first-principles approach to calculate the phonon thermal conductivity based on the Green-Kubo formalism. In this approach, the density functional theory energy is distributed to each atom, and a two-step method in the molecular dynamics is introduced to avoid the atomic position R wrapping problem in a periodic system when the heat current is calculated. We show that this first-principles Green-Kubo approach is particularly suitable for disordered systems like amorphous and liquid, where the thermal conductivities are small due to strong phonon scattering but difficult to be calculated using anharmonic interaction energy. We have applied our method to liquid Ar, liquid Si, and amorphous Si. In this work, the calculated thermal conductivities agree well with previous theoretical and experimental results. We have also compared our method to previous works combining first-principles simulations with the Green-Kubo formalism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Quantum mechanics of a constrained particle on an ellipsoid: Bein formalism and Geometric momentum

In this work we apply the Dirac method in order to obtain the classical relations for a particle on an ellipsoid. We also determine the quantum mechanical form of these relations by using Dirac quantization. Then by considering the canonical commutation relations between the position and momentum operators in terms of curved coordinates, we try to propose the suitable representations for momentum operator that satisfy the obtained commutators between position and momentum in Euclidean space. We see that our representations for momentum operators are the same as geometric one.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A Methodology for the Analysis of Water Oxidation Electrocatalysts in the Absence of Limiting Current that Avoids the Pitfalls of Existing Methods

Water oxidation is an important reaction studied as a way to generate electrons from water, to promote water splitting and the formation of green hydrogen. When using electrodes to drive homogeneous water oxidation catalysis, cyclic voltammograms are analyzed to provide catalytic rate constants. There are two main methods, foot-of-the-wave analysis (FOWA) and limiting current analysis. FOWA relies on approximations inherent to analyzing water oxidation catalysis, such as determining the formal potential of the catalytic intermediate, E 0 cat . Limiting current methods are the optimal way to analyze catalyst performance but rely on observable limiting current, which is virtually never seen in water oxidation. To avoid those issues, a method is proposed for analyzing nonideal cyclic voltammetry waveshapes in water oxidation: by analyzing rate data across a large range of potentials, an optimal potential, E 0 cat , can be obtained, where catalytic current, i cat , is nearly independent of scan rate and has a linear dependency on buffer concentration. Here, the method is applied to four homogeneous water oxidation catalysts with prior extensive electrochemical elucidation, all of which lack an ideal, purely kinetic waveshape in cyclic voltammetry. Application of the method avoids the biases of the other methods cited for the kinetic analyses of water oxidation catalysts.

14 SOLAR ENERGY↗
DOE OSTI2014

Rashba-Zeeman-effect-induced spin filtering energy windows in a quantum wire

We perform a numerical study on the spin-resolved transport in a quantum wire (QW) under the modulation of both Rashba spin-orbit coupling (SOC) and a perpendicular magnetic field by using the developed Usuki transfer-matrix method in combination with the Landauer-Büttiker formalism. Wide spin filtering energy windows can be achieved in this system for unpolarized spin injection. In addition, both the width of energy window and the magnitude of spin conductance within these energy windows can be tuned by varying Rashba SOC strength, which can be apprehended by analyzing the energy dispersions and spin-polarized density distributions inside the QW, respectively. Further study also demonstrates that these Rashba-SOC-controlled spin filtering energy windows show a strong robustness against disorders. These findings may not only benefit to further understand the spin-dependent transport properties of a QW in the presence of external fields but also provide a theoretical instruction to design a spin filter device.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Spin projection with double hybrid density functional theory

A spin projected double-hybrid density functional theory is presented that accounts for different scaling of opposite and same spin terms in the second order correction. This method is applied to three dissociation reactions which in the unprojected formalism exhibit significant spin contamination with higher spin states. This gives rise to a distorted potential surface and can lead to poor geometries and energies. The projected method presented is shown to improve the description of the potential over unprojected double hybrid density functional theory. Comparison is made with the reference states of the two double hybrid functionals considered here (B2PLYP and mPW2PLYP) in which the projected potential surface is degraded by an imbalance in the description of dynamic and static correlation.

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

Polyad quantum numbers and multiple resonances in anharmonic vibrational studies of polyatomic molecules

In the theory of anharmonic vibrations of a polyatomic molecule, mixing the zero-order vibrational states due to cubic, quartic and higher-order terms in the potential energy expansion leads to the appearance of more-or-less isolated blocks of states (also called polyads), connected through multiple resonances. Such polyads of states can be characterized by a common secondary integer quantum number. This polyad quantum number is defined as a linear combination of the zero-order vibrational quantum numbers, attributed to normal modes, multiplied by non-negative integer polyad coefficients, which are subject to definition for any particular molecule. According to Kellman's method [J. Chem. Phys. 93, 6630 (1990)], the corresponding formalism can be conveniently described using vector algebra. In the present work, a systematic consideration of polyad quantum numbers is given in the framework of the canonical Van Vleck perturbation theory (CVPT) and its numerical-analytic operator implementation for reducing the Hamiltonian to the quasi-diagonal form, earlier developed by the authors. It is shown that CVPT provides a convenient method for the systematic identification of essential resonances and the definition of a polyad quantum number. The method presented is generally suitable for molecules of significant size and complexity, as illustrated by several examples of molecules up to six atoms. The polyad quantum number technique is very useful for assembling comprehensive basis sets for the matrix representation of the Hamiltonian after removal of all non-resonance terms by CVPT. In addition, the classification of anharmonic energy levels according to their polyad quantum numbers provides an additional means for the interpretation of observed vibrational spectra.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Wave function continuity and the diagonal Born-Oppenheimer correction at conical intersections

We demonstrate that though exact in principle, the expansion of the total molecular wave function as a sum over adiabatic Born-Oppenheimer (BO) vibronic states makes inclusion of the second-derivative nonadiabatic energy term near conical intersections practically problematic. In order to construct a well-behaved molecular wave function that has density at a conical intersection, the individual BO vibronic states in the summation must be discontinuous. When the second-derivative nonadiabatic terms are added to the Hamiltonian, singularities in the diagonal BO corrections (DBOCs) of the individual BO states arise from these discontinuities. In contrast to the well-known singularities in the first-derivative couplings at conical intersections, these singularities are non-integrable, resulting in undefined DBOC matrix elements. Though these singularities suggest that the exact molecular wave function may not have density at the conical intersection point, there is no physical basis for this constraint. Instead, the singularities are artifacts of the chosen basis of discontinuous functions. We also demonstrate that continuity of the total molecular wave function does not require continuity of the individual adiabatic nuclear wave functions. We classify nonadiabatic molecular dynamics methods according to the constraints placed on wave function continuity and analyze their formal properties. Based on our analysis, it is recommended that the DBOC be neglected when employing mixed quantum-classical methods and certain approximate quantum dynamical methods in the adiabatic representation.

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

Bismuth Detection in Alkaline Electrolyte via Anodic Stripping Voltammetry for Battery Separator Evaluation

Abstract Anodic stripping voltammetry (ASV) has been widely used for the detection of several heavy metal ions in neutral and acidic solution, in many cases employing electrodes and/or solutions incorporating Bi. In this work we demonstrate that Bi(OH) 4 − ion concentration can be measured in highly alkaline 8.5 M KOH solution using ASV. The addition of Pb in similar concentrations to the Bi(OH) 4 − being measured is shown to improve both the sensitivity and precision of the method. When the Pb additive is employed, a formal limit of detection of 8.5 ppb is achieved, compared to 17.3 ppb when the Pb additive is not used. Due to the use of Bi additives in alkaline battery chemistries, it follows that separators which limit Bi(OH) 4 − diffusion into the bulk electrolyte and away from the electrodes are of interest. For this purpose, we utilize ASV to determine Bi(OH) 4 − diffusion rates through Celgard 3501, cellophane 350P00, and Nafion 211. Bi(OH) 4 − crossover rates, as determined by ASV, are shown to be repeatable and consistent with expectations from the known separator structure.

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

The photon content of the neutron

In this work, we complete our CT18qed study with the neutron’s photon parton distribution function (PDF), which is essential for the nucleus scattering phenomenology. Two methods, CT18lux and CT18qed, based on the LUXqed formalism and the DGLAP evolution, respectively, to determine the neutron’s photon PDF have been presented. Various low-Q 2 non-perturbative variations have been carefully examined, which are treated as additional uncertainties on top of those induced by quark and gluon PDFs. The impacts of the momentum sum rule as well as isospin symmetry violation have been explored and turned out to be negligible. A detailed comparison with other neutron’s photon PDF sets has been performed, which shows a great improvement in the precision and a reasonable uncertainty estimation. Finally, two phenomenological implications are demonstrated with photon-initiated processes: neutrino-nucleus W-boson production, which is important for the near-future TeV–PeV neutrino observations, and the axion-like particle production at a high-energy muon beam-dump experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗