Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Approximation 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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Approximation of the Frequency-Amplitude Curve Using the Homotopy Analysis Method

The Frequency-Amplitude (F-A) curve on power system oscillation under a large disturbance characterizes how a natural oscillation mode transitions to nonlinear oscillations with growing amplitudes and decaying frequencies. The existing formulation of the F-A curve is derived by solving elliptical integrals on oscillation of a single-machine-infinite-bus equivalent about the targeted oscillation mode. The formula is in a form of infinite series and needs to sum a large number of terms for satisfactory accuracy. This paper introduces an explicit, approximate expression obtained from the Homotopy Analysis Method on the F-A curve. The proposed F-A curve expression is derived from an SMIB system and verified on the IEEE 3-machine 9-bus system to show how the oscillation frequency of a dominant mode varies with oscillation amplitude under large disturbances.

frequency-amplitude (F-A) curve↗
DOE OSTI2017

Annotations on the virtual element method for second-order elliptic problems

This document contains working annotations on the Virtual Element Method (VEM) for the approximate solution of diffusion problems with variable coefficients. To read this document you are assumed to have familiarity with concepts from the numerical discretization of Partial Differential Equations (PDEs) and, in particular, the Finite Element Method (FEM). This document is not an introduction to the FEM, for which many textbooks (also free on the internet) are available. Eventually, this document is intended to evolve into a tutorial introduction to the VEM (but this is really a long-term goal).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Magnetic Properties of La-Doped ZnO(0001)–Zn Polar Surface with and without Vacancies: a First-Principle Study

The effects of La doping with and without vacancies on the magnetism of a ZnO(0001)–Zn polar surface were studied using the generalized gradient approximation + U (GGA + U) method of density functional theory to determine the magnetism source of La-doped ZnO. Results revealed that La doping with O vacancy did not cause ZnO ferromagnetism, but La-doped ZnO with Zn vacancy was ferromagnetic. The net magnetic moment was primarily related to the unpaired electron amount of the O 2p state and slightly modified by orbital and anomalous magnetic moments. The anomalous magnetic moment originated from the radiating magnetism moment via quantum electrodynamics and was confirmed to exist via theoretical calculations of magnetism. Formation energy showed that La doping easily generated Zn vacancies in ZnO, but magnetic quenching occurred when the ratio of La ions to Zn vacancies was 2:1. The unpaired O 2p electrons induced by Zn vacancies were the magnetism source of La-doped ZnO. Moreover, the Curie temperature of La-doped ZnO with Zn vacancy reached room temperature. Results showed that the ZnO(0001)–Zn polar surface can be utilized as the surface of ZnO-based diluted magnetic semiconductors because of its unipolarity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Optimal, centralized dynamic curbside parking space zoning

In this paper we formulate a dynamic mixed integer program for optimally zoning curbside parking spaces subject to transportation policy-inspired constraints and regularization terms. First, we illustrate how given some objective of curb zoning valuation as a function of zone type (paid parking, bus stop, etc.), dynamically rezoning involves unrolling this optimization program over a fixed time horizon. Second, we implement two different solution methods given an example curb zoning valuation. In the first method, we solve long horizon dynamic zoning problems via approximate dynamic programming. In the second method, we employ Dantzig-Wolfe decomposition to break-up the mixed-integer program into a master problem and several sub-problems that can be solved in parallel. This speeds up the computational solve-time of the MIP considerably. We present simulation results and comparisons of the different employed techniques on vehicle arrival-rate data obtained for a neighborhood in downtown Seattle, Washington, USA.

Nazir, Mohammad Nawaf↗
DOE OSTI2015

TU-F-CAMPUS-I-01: Investigation of the Effective Dose From Bolus Tracking Acquisitions at Different Anatomical Locations in the Chest for CT

Purpose: Stationary table acquisitions (Bolus tracking) in X-ray Computed Tomography (CT) can Result in dose length products (DLP) comparable to spiral scans. It is today unclear whether or not the effective dose (E) for Bolus Tracking can be approximated using target region specific conversion factors (E/DLP). The purpose of this study was to investigate how E depends on the anatomical location of the Bolus Tracking in relation to Chest CT scans with the same DLP. Methods: Effective doses were approximated for the ICRP 110 adult Reference Male (AM) and adult Reference Female (FM) computational voxel phantoms using software for CT dose approximations (pre-simulated MC data). The effective dose was first approximated for a Chest CT scan using spiral technique and a CTDIvol (32 cm) of 6 mGy. The effective dose from the spiral scan was then compared to E approximated for contiguous Bolus Tracking acquisitions (1 cm separation), with a total collimation of 1 cm, over different locations of the chest of the voxel phantoms. The number of rotations used for the Bolus Tracking acquisitions was adjusted to yield the same DLP (32 cm) as the spiral scan. Results: Depending on the anatomical location of the Bolus Tracking, E ranged by factors of 1.3 to 6.8 for the AM phantom and 1.4 to 3.3 for the AF phantom, compared to the effective dose of the spiral scans. The greatest E for the Bolus Tracking acquisitions was observed for anatomical locations coinciding with breast tissue. This can be expected as breast tissue has a high tissue weighting factor in the calculation of E. Conclusion: For Chest CT scans, the effective dose from Bolus Tracking is highly dependent on the anatomical location where the scan is administered and will not always accurately be represented using target region specific conversion factors.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-G-IeP3-04: Effective Dose Measurements in Fast Kvp Switch Dual Energy Computed Tomography

Purpose: The objective of this study was two-fold: (a) to test a new approach to approximating organ dose by using the effective energy of the combined 80kV/140kV beam in dual-energy (DE) computed tomography (CT), and (b) to derive the effective dose (ED) in the abdomen-pelvis protocol in DECT. Methods: A commercial dual energy CT scanner was employed using a fast-kV switch abdomen/pelvis protocol alternating between 80 kV and 140 kV. MOSFET detectors were used for organ dose measurements. First, an experimental validation of the dose equivalency between MOSFET and ion chamber (as a gold standard) was performed using a CTDI phantom. Second, the ED of DECT scans was measured using MOSFET detectors and an anthropomorphic phantom. For ED calculations, an abdomen/pelvis scan was used using ICRP 103 tissue weighting factors; ED was also computed using the AAPM Dose Length Product (DLP) method and compared to the MOSFET value. Results: The effective energy was determined as 42.9 kV under the combined beam from half-value layer (HVL) measurement. ED for the dual-energy scan was calculated as 16.49 ± 0.04 mSv by the MOSFET method and 14.62 mSv by the DLP method. Conclusion: Tissue dose in the center of the CTDI body phantom was 1.71 ± 0.01 cGy (ion chamber) and 1.71 ± 0.06 (MOSFET) respectively; this validated the use of effective energy method for organ dose estimation. ED from the abdomen-pelvis scan was calculated as 16.49 ± 0.04 mSv by MOSFET and 14.62 mSv by the DLP method; this suggests that the DLP method provides a reasonable approximation to the ED.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Structural transformation features in titanium crystallite under mechanical loading

The features of defect structure generation and development in titanium crystallites were studied on the basis of molecular dynamics method. Interatomic interaction was described using many-body potentials calculated in the approximation of the embedded atom method. It is shown that local structural changes begin to occur in the crystallite at the achievement of a threshold strain value, which is accompanied by a dramatic decrease of potential energy. The features of the formation of local structural changes, which are precursors of classical stacking faults, were studied.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Taylor approximation variance reduction for approximation errors in PDE-constrained Bayesian inverse problems

In numerous applications, surrogate models are used as a replacement for accurate parameter-to-observable mappings when solving large-scale inverse problems governed by partial differential equations (PDEs). The surrogate model may be a computationally cheaper alternative to the accurate parameter-to-observable mappings and/or may ignore additional unknowns or sources of uncertainty. The Bayesian approximation error (BAE) approach provides a means to account for the induced uncertainties and approximation errors, i.e. the errors between the accurate parameter-to-observable mapping and the surrogate. The statistics of these errors are, however, in general unknown a priori, and are thus calculated using Monte Carlo sampling. Although the sampling is typically carried out offline, i.e. before considering the data, the process can still represent a computational bottleneck. In this work, we develop a scalable computational approach for reducing the costs associated with the sampling stage of the BAE approach. Specifically, we consider the Taylor expansion of the accurate and surrogate forward models with respect to the uncertain parameter fields either as a control variate for variance reduction or as a means to directly and efficiently approximate the mean and covariance of the approximation errors. We propose efficient methods for evaluating the expressions for the mean and covariance of the Taylor approximations based on linear(-ized) PDE solves. Furthermore, the proposed approach is independent of the dimension of the uncertain parameter, depending instead on the intrinsic dimension of the data, ensuring scalability to high-dimensional problems. The potential benefits of the proposed approach are demonstrated for two high-dimensional inverse problems governed by PDE examples, namely for the estimation of a distributed Robin boundary coefficient in a linear diffusion problem, and for a coefficient estimation problem governed by a nonlinear diffusion problem.

Bayesian approximation error↗
DOE OSTI2018

Generalized Jacobi–Galerkin method for nonlinear fractional differential algebraic equations

In this paper, we provide an approximate approach based on the Galerkin method to solve a class of nonlinear fractional differential algebraic equations. The fractional derivative operator in the Caputo sense is utilized and the generalized Jacobi functions are employed as trial functions. The existence and uniqueness theorem as well as the asymptotic behavior of the exact solution are provided. It is shown that some derivatives of the solutions typically have singularity at origin dependence on the order of the fractional derivative. The influence of the perturbed data on the exact solutions along with the convergence analysis of the proposed scheme is also established. Some illustrative examples provided to demonstrate that this novel scheme is computationally efficient and accurate.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Structural, electronic and elastic properties of REIr{sub 2} (RE=La and Ce) Laves phase compounds

REIr{sub 2} (RE = La and Ce) Laves phase intermetallic compounds were investigated with respect to their structural, electronic and elastic properties using full potential linearized augmented plane wave (FP-LAPW) method within generalized gradient approximation (GGA) as implemented in WIEN2k code. The ground state properties such as lattice constants (a{sub 0}), bulk modulus (B), pressure derivative of bulk modulus (B′) and density of state at Fermi level N(E{sub F}) have been obtained by optimization method. The electronic structure (BS, TDOS and PDOS) reveals that these Laves phase compounds are metallic in nature. The calculated elastic constants indicate that these compounds are mechanically stable at ambient pressure and found to be ductile in nature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

High-throughput, one-step screening, cloning and expression based on the lethality of DpnI in Escherichia coli

Highlights: • The extreme toxicity of the DpnI from Streptococcus pneumoniae in bacteria containing Gm6ATC, the product of Dam methylation. • One-step screening, cloning and expression based on the lethality of DpnI in Escherichia coli. • The efficiency of using this method for cloning was approximately 100%. • Rapid one-step selection of soluble proteins. The manipulation of recombinant DNA has been an integral step in molecular biology to date. A number of strategies have been developed over the years, as traditional cloning methods are time consuming, have high backgrounds and low efficiency and are often limited by the number of suitable restriction sites available. Here, we constructed a series of new positive-selection-based cloning vectors that overcome most of the above mentioned drawbacks and can be applied in both eukaryotic and prokaryotic systems. This strategy is based on the extreme toxicity of DpnI in wild-type E. coli and the inactivation of this lethality by the introduction of target gene within multiple cloning sites. There are no rapid approaches for identifying soluble proteins for high-throughput screening. In this study, we combined this highly efficient cloning strategy with rapid identification of soluble proteins to construct vectors with multiple fusion tags, such as MBP, GST, CBD, NusA, and Sumo, to generate enzymes with potential diagnostic, industrial or therapeutic applications. Thus, this versatile positive-selection-based technology is appropriate for routine cloning, DNA library construction, and high-throughput screening for the expression of proteins of interest.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Vibrational coherence and energy transfer in two-dimensional spectra with the optimized mean-trajectory approximation

The optimized mean-trajectory (OMT) approximation is a semiclassical method for computing vibrational response functions from action-quantized classical trajectories connected by discrete transitions that represent radiation-matter interactions. Here, we extend the OMT to include additional vibrational coherence and energy transfer processes. This generalized approximation is applied to a pair of anharmonic chromophores coupled to a bath. The resulting 2D spectra are shown to reflect coherence transfer between normal modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

The arbitrary-order virtual element method for linear elastodynamics models. Convergence, stability and dispersion-dissipation analysis.

We design the conforming virtual element method for the numerical approximation of the two dimensional elastodynamics problem. We prove stability and convergence of the semi-discrete approximation and derive optimal error estimates under $\textit{h}$-refinement in both the energy and the $L^2$ norms, and optimal error estimates under $\textit{p}$-refinement in the energy norm. The performance of the proposed virtual element method is assessed on a set of different computational meshes, including non-convex cells up to order four in the h-refinement setting. Exponential convergence is also experimentally observed under p-refinement. Finally, we present a dispersion-dissipation analysis for both the semi-discrete and fully-discrete schemes, showing that polygonal meshes behave as classical simplicial/quadrilateral grids in terms of dispersion-dissipation properties.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Rigorous Computation of Diffusion Coefficients for Expanding Maps

For real analytic expanding interval maps, a novel method is given for rigorously approximating the diffusion coefficient of real analytic observables. As a theoretical algorithm, our approximation scheme is shown to give quadratic exponential convergence to the diffusion coefficient. The method for converting this rapid convergence into explicit high precision rigorous bounds is illustrated in the setting of Lanford’s map x↦2x+(1/2)x(1−x) (mod 1).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Density-functional theory of material design: fundamentals and applications-I

This article is part-I of a review of density-functional theory (DFT) that is the most widely used method for calculating electronic structure of materials. The accuracy and ease of numerical implementation of DFT methods has resulted in its extensive use for materials design and discovery and has thus ushered in the new field of computational material science. In this article, we start with an introduction to Schrödinger equation and methods of its solutions. After presenting exact results for some well-known systems, difficulties encountered in solving the equation for interacting electrons are described. How these difficulties are handled using the variational principle for the energy to obtain approximate solutions of the Schrödinger equation is discussed. The resulting Hartree and Hartree–Fock theories are presented along with results they give for atomic and solid-state systems. We then describe Thomas–Fermi theory and its extensions which were the initial attempts to formulate many-electron problem in terms of electronic density of a system. Having described these theories, we introduce modern DFT by discussing Hohenberg–Kohn theorems that form its foundations. We then go on to discuss Kohn–Sham (KS) formulation of DFT in its exact form. Next, local density approximation (LDA) is introduced and solutions of KS equation for some representative systems, obtained using the LDA, are presented. We end part-I of the review describing the contents of part-II.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Assessing the distinguishable cluster approximation based on the triple bond-breaking in the nitrogen molecule

Obtaining the correct potential energy curves for the dissociation of multiple bonds is a challenging problem for ab initio methods which are affected by the choice of a spin-restricted reference function. Coupled cluster (CC) methods such as CCSD (coupled cluster singles and doubles model) and CCSD(T) (CCSD + perturbative triples) correctly predict the geometry and properties at equilibrium but the process of bond dissociation, particularly when more than one bond is simultaneously broken, is much more complicated. New modifications of CC theory suggest that the deleterious role of the reference function can be diminished, provided a particular subset of terms is retained in the CC equations. The Distinguishable Cluster (DC) approach of Kats and Manby [J. Chem. Phys. 139, 021102 (2013)], seemingly overcomes the deficiencies for some bond-dissociation problems and might be of use in quasi-degenerate situations in general. DC along with other approximate coupled cluster methods such as ACCD (approximate coupled cluster doubles), ACP-D45, ACP-D14, 2CC, and pCCSD(α, β) (all defined in text) falls under a category of methods that are basically obtained by the deletion of some quadratic terms in the double excitation amplitude equation for CCD/CCSD (coupled cluster doubles model/coupled cluster singles and doubles model). Here these approximate methods, particularly those based on the DC approach, are studied in detail for the nitrogen molecule bond-breaking. The N{sub 2} problem is further addressed with conventional single reference methods but based on spatial symmetry-broken restricted Hartree–Fock (HF) solutions to assess the use of these references for correlated calculations in the situation where CC methods using fully symmetry adapted SCF solutions fail. The distinguishable cluster method is generalized: 1) to different orbitals for different spins (unrestricted HF based DCD and DCSD), 2) by adding triples correction perturbatively (DCSD(T)) and iteratively (DCSDT-n), and 3) via an excited state approximation through the equation of motion (EOM) approach (EOM-DCD, EOM-DCSD). The EOM-CC method is used to identify lower-energy CC solutions to overcome singularities in the CC potential energy curves. It is also shown that UHF based CC and DC methods behave very similarly in bond-breaking of N{sub 2}, and that using spatially broken but spin preserving SCF references makes the CCSD solutions better than those for DCSD.

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

Demonstration of Full Core Reactor Depletion with MPACT

MPACT is a three-dimensional (3-D) whole core transport code capable of generating sub-pin level power distributions. This feature is accomplished by obtaining the integral transport solutions to the heterogeneous reactor problem in which the actual detailed geometrical configuration of fuel components such as the pellet and cladding is modelled explicitly.. The cross section data needed for the neutron transport calculation are obtained directly from a multigroup microscopic cross section library similar to those used in lattice physics codes. Hence MPACT involves neither a priori homogenization nor group condensation for the core spatial solution. The integral transport solution is obtained by means of the method of characteristics (MOC) that employs discrete ray tracing. Since the direct application of the 3-D MOC capability in MPACT to 3-D core configuration requires considerable amounts of memory and computing time for practical reactor applications, an alternative approximate 3-D solution method was implemented in MPACT based on a 2D-1D approach which employs planar MOC solutions in the framework of the 3-D coarse mesh finite difference (CMFD) formulation. The axial coupling is resolved by one-dimensional (1-D) diffusion solutions and the planar and axial problems are coupled through the transverse leakage. The use of a lower order 1-D solution in the axial direction is justified by the fact that most heterogeneity in the core occurs in the radial direction rather than the axial. However, it is possible to use one of the higher order methods in MPACT such as the 1-D SN kernel if improved solution accuracy is desired. The CMFD formulation which originally served as an efficient nodal formulation is used in MPACT to accelerate the whole core transport calculation. The basic mesh in the CMFD formulation is a pin cell which is much coarser than the flat source regions defined for MOC calculations. The concept of dynamic homogenization of group constants for the pin cells provides the basis of the CMFD formulation in whole core transport calculation. The intra-cell flux distribution determined from the MOC calculation is used to generate the homogenized cell constants while the MOC cell surface averaged currents are used to determine the radial nodal coupling coefficients. This dynamic implementation of the equivalence theory formulation ensures the same transport solution is obtained with CMFD as the one obtained using unaccelerated MOC calculations alone. The ability of MPACT to predict the Hot, Zero Power (HZP) condition of a reactor was reported in CASL-U-2014-0045-000. The purpose of this report is to demonstrate the capability of MPACT to simulate the operation of a reactor at Hot, Full Power (HFP) conditions during normal operation throughout a full operating cycle. The time scale for a reactor burnup cycle is typically one to two years and therefore short time phenomena such as the effect of delayed neutrons can be neglected. Instead the time dependence can be treated quasi-statically using the depletion algorithm presented in CASL-U-2013-0276-001. The full power operation requires thermal hydraulic feedback to be present during the entire simulation which will utilize the work in CASL-U-2014-0051-000. However, due to the significant computational burden of the coupled physics simulation, a simplified thermal hydraulic model is used for the analysis reported here which is principally intended to demonstrate the full core neutronics depletion. In addition to the thermal hydraulic feedback, the effect of Boron 10 depletion should be modeled in cases where major boron additions are not performed for extended periods of time. Each of these capabilities will be described in the following sections.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

A Ranked-Orbital Approach to Select Active Spaces for High-Throughput Multireference Computation

The past decade has seen a great increase in the application of high-throughput computation to a variety of important problems in chemistry. However, one area which has been resistant to the high-throughput approach is multireference wave function methods, in large part due to the technicalities of setting up these calculations and in particular the not always intuitive challenge of active space selection. As we look toward a future of applying high-throughput computation to all areas of chemistry, it is important to prepare these methods for large-scale automation. Here, we propose a ranked-orbital approach to select active spaces with the goal of standardizing multireference methods for high-throughput computation. This method allows for the meaningful comparison of different active space selection schemes and orbital localizations, and we demonstrate the utility of this approach across 1120 multireference calculations for the excitation energies of small molecules. Our results reveal that it is helpful to distinguish the method used to generate orbitals from the method of ranking orbitals in terms of importance for the active space. Additionally, we propose our own orbital ranking scheme that estimates the importance of an orbital for the active space through a pair-interaction framework from orbital energies and features of the Hartree–Fock exchange matrix. Here, we call this new scheme the “approximate pair coefficient” (APC) method and we show that it performs quite well for the test systems presented.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗