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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Numerical simulation of scattering of acoustic waves by inelastic bodies using hypersingular boundary integral equation

A numerical method for solving a problem of diffraction of acoustic waves by system of solid and thin objects based on the reduction the problem to a boundary integral equation in which the integral is understood in the sense of finite Hadamard value is proposed. To solve this equation we applied piecewise constant approximations and collocation methods numerical scheme. The difference between the constructed scheme and earlier known is in obtaining approximate analytical expressions to appearing system of linear equations coefficients by separating the main part of the kernel integral operator. The proposed numerical scheme is tested on the solution of the model problem of diffraction of an acoustic wave by inelastic sphere.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Flux tube spectra from approximate integrability at low energies

We provide a detailed introduction to a method we recently proposed for calculating the spectrum of excitations of effective strings such as QCD flux tubes. The method relies on the approximate integrability of the low-energy effective theory describing the flux tube excitations and is based on the thermodynamic Bethe ansatz. The approximate integrability is a consequence of the Lorentz symmetry of QCD. For excited states, the convergence of the thermodynamic Bethe ansatz technique is significantly better than that of the traditional perturbative approach. We apply the new technique to the lattice spectra for fundamental flux tubes in gluodynamics in D = 3 + 1 and D = 2 + 1, and to k-strings in gluodynamics in D = 2 + 1. We identify a massive pseudoscalar resonance on the worldsheet of the confining strings in SU(3) gluodynamics in D = 3 + 1, and massive scalar resonances on the worldsheet of k = 2.3 strings in SU(6) gluodynamics in D = 2 + 1.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Analysis of Electromagnetic Scattering From a Dihedral Cut Into a Sphere Using Vector Spherical Harmonics

Proper calibration or verification of the calibration of fully - polarimetric radar requires more than a simple specular target, such as a conducting sphere or corner reflector. The dihedral reflector is a useful radar - scattering target for calibrating the polarization response of a radar system, because it possesses the ability to change the direction of the scattered electric - field vector in a special way, while a sphere or trihedral target does not. However, the scattered electromagnetic field of a sphere can be computed with arbitrary precision using the well - known Mie series, which is obtained by expanding the scattered field as a series of vector spherical harmonics. In contrast, the computation of the scattered electromagnetic field of dihedral reflectors typically relies on numerical methods or various approximate solutions of Maxwell's equations. Precision calibration of wide - bandwidth polarimetric radar has not been adequately achieved for some applications using this approach. This report describes a method for computing the scattered electromagnetic field of a specially shaped dihedral, derived through the application of a series expansion in vector spherical harmonics, with the goal of achieving the same precision that is available for the sphere.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Distributed and communication-efficient solutions to linear equations with special sparse structure

In this paper we report two distributed and communication-efficient algorithms based on the multi-agent system are proposed to solve a system of linear equations with the Laplacian sparse system matrix. One algorithm is based on the gradient descent method in optimization. In this algorithm, the agents only share partial information instead of all of their collective state vectors to save significant communication. The other algorithm is obtained by approximating Newton’s method for a faster convergence rate. Although it requires twice as much communication as the first one, it is still communication-efficient given the low dimension of the information shared among agents. The convergence at a linear rate is proved for both algorithms, and a comprehensive comparison of their convergence rate, communication burden, and computation costs is also performed. The proposed algorithms can be applied to various systems to solve those problems that can be modeled as a system of linear equations with a Laplacian sparse system matrix. Simulation results with the electric power system illustrate their effectiveness.

42 ENGINEERING↗
DOE OSTI2013

Rhombohedrally Distorted γ-Au 5– x Zn 8+ y Phases in the Au–Zn System

The region of the Au–Zn phase diagram encompassing γ-brass-type phases has been studied experimentally from 45 to 85 atom % Zn. The γ phases were obtained directly from the pure elements by heating to 680 °C in evacuated silica tubes, followed by annealing at 300 °C. Powder X-ray and single-crystal diffraction studies show that γ-“Au5Zn8” phases adopt a rhombohedrally distorted Cr5Al8 structure type rather than the cubic Cu5Zn8 type. The refined compositions from two single crystals extracted from the Zn- and Au-rich loadings are Au4.27(3)Zn8.26(3)γ0.47 (I) and Au4.58(3)Zn8.12(3)γ0.3 (II), respectively (γ = vacancy). These (I and II) refinements indicated both nonstatistical mixing of Au and Zn atoms as well as partially ordered vacancy distributions. The structures of these γ phases were solved in the acentric space group R3m (No. 160, Z = 6), and the observed lattice parameters from powder patterns were found to be a = 13.1029(6) and 13.1345(8) Å and c = 8.0410(4) and 8.1103(6) Å for crystals I and II, respectively. According to single-crystal refinements, the vacancies were found on the outer tetrahedron (OT) and octahedron (OH) of the 26-atom cluster. Single-crystal structural refinement clearly showed that the vacancy content per unit cell increases with increasing Zn, or valence-electron concentration. Electronic structure calculations, using the tight-binding linear muffin-tin orbital method with the atomic-sphere approximation (TB-LMTO-ASA) method, indicated the presence of a well-pronounced pseudogap at the Fermi level for “Au5Zn8” as the representative composition, an outcome that is consistent with the Hume–Rothery interpretation of γ brass.

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

Quasielastic least-squares reverse time migration of PS reflections

Conventional analysis of amplitude variation with offset for elastic PS reflections is based on analytical reflection coefficients in a layered medium, and wave-equation-based PS migration is mainly used to produce a structural image. To overcome this problem, we have developed a least-squares reverse time migration (LSRTM) method for elastic PS reflections based on a quasielastic wave equation. Here, the quasielastic wave equation can accurately model PS reflections with elastic amplitudes under the first-order Born approximation. Our LSRTM method inverts for perturbations of the S-wave velocity and density by minimizing the L 2 norm of the difference between recorded and predicted PS reflections modeled using a quasielastic wave equation. We refer to our new method as quasielastic LSRTM of PS reflections. Numerical tests on synthetic and field data indicate that our method can properly handle the amplitudes of elastic PS reflections and provides an accurate estimate of the perturbations of S-wave velocity and density. Extending the method to the 3D case is not straightforward and might require incorporating certain PS data processing techniques into the inversion itself.

58 GEOSCIENCES↗
DOE OSTI2018

Finite-Difference Methods for Fractional Differential Equations of Order 1/2

In this work, we study approximations of solutions of fractional differential equations of order 1/2. We present a new method of approximation and obtain the order of convergence. The presentation is given within the abstract framework of a semidiscrete approximation scheme, which includes finite-element methods, finite-difference schemes, and projection methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

QMCPACK: Advances in the development, efficiency, and application of auxiliary field and real-space variational and diffusion quantum Monte Carlo

We review recent advances in the capabilities of the open source ab initio Quantum Monte Carlo (QMC) package QMCPACK and the workflow tool Nexus used for greater efficiency and reproducibility. The auxiliary field QMC (AFQMC) implementation has been greatly expanded to include k-point symmetries, tensor-hypercontraction, and accelerated graphical processing unit (GPU) support. These scaling and memory reductions greatly increase the number of orbitals that can practically be included in AFQMC calculations, increasing the accuracy. Advances in real space methods include techniques for accurate computation of bandgaps and for systematically improving the nodal surface of ground state wavefunctions. Results of these calculations can be used to validate application of more approximate electronic structure methods, including GW and density functional based techniques. To provide an improved foundation for these calculations, we utilize a new set of correlation-consistent effective core potentials (pseudopotentials) that are more accurate than previous sets; these can also be applied in quantum-chemical and other many-body applications, not only QMC. These advances increase the efficiency, accuracy, and range of properties that can be studied in both molecules and materials with QMC and QMCPACK.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

On modeling subgrid-scale macro-structures in narrow twisted channels

Porosity-based subgrid topography models often fail to capture the effects of subgrid-scale topographic features in the interior of coarse grid cells. Existing approaches that modify bottom roughness or a drag coefficient are inadequate for macro-structures (large emergent or submerged obstacles) in subgrid-scale narrow twisted channels. Such structures partially block the cross-sectional area and provide enhanced topographic dissipation – effects that are not well represented by a drag coefficient that scales on a coarse-grid cell-averaged velocity and the cell volume. The relative alignment between mesh and flow further complicates this problem as it makes the subgrid model sensitive to mesh design. In the present study, three new approaches for simulating subgrid-scale macro-structures in narrow channels are proposed. The interior partial-blocking effect of structures is modeled as reduction of grid face-area. The sheltering of flow volumes around obstacles, which leads to topographic dissipation, is modeled by reducing the cell volume in the momentum equation (only). A mesh-shift procedure is designed to optimize mesh alignment for identifiable subgrid features. Combining the three subgrid methods improves the approximation of surface elevation and in-channel flow rate with a coarse-grid model. Tests are conducted for channelized flow using both synthetic domains and real marsh topography. Finally, the new methods reduce the overall mesh dependency of the subgrid model and provides stronger physical connection between effects of macro-structures and their geometry at coarse grid scales.

42 ENGINEERING↗
DOE OSTI2021

Numerically exact generalized Green's function cluster expansions for electron-phonon problems

We generalize the family of approximate momentum average methods to formulate a numerically exact, convergent hierarchy of equations whose solution provides an efficient algorithm to compute the Green's function of a particle dressed by bosons suitable in the entire parameter regime. We use this approach to extract ground-state properties and spectral functions. Our approximation-free framework, dubbed the generalized Green's function cluster expansion (GGCE), allows access to exact numerical results in the extreme adiabatic limit, where many standard methods struggle or completely fail. We showcase the performance of the method, specializing three important models of charge-boson coupling in solids and molecular complexes: the molecular Holstein model, which describes coupling between charge density and local distortions, the Peierls model, which describes modulation of charge hopping due to intersite distortions, and a more complex Holstein + Peierls system with couplings to two different phonon modes, paradigmatic of charge-lattice interactions in organic crystals. Furthermore, the GGCE serves as an efficient approach that can be systematically extended to different physical scenarios, thus providing a tool to model the frequency dependence of dressed particles in realistic settings.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Dependence of the structure and dynamics of liquid silicon on the choice of density functional approximation

In addition to its technological relevance, silicon poses a challenge for first principles simulations because it undergoes a semiconductor-to-metal transition upon melting. Moreover, the resulting metallic liquid contains a mixture of metallic and covalent bonding. This coexistence of fundamentally different interactions is difficult to describe within approximate density functional methods, which oftentimes cannot accurately describe these two extremes simultaneously. Here, we report an investigation of the structure, dynamics, and thermodynamics of liquid silicon using ab initio molecular dynamics simulations with three density functional approximations: the local density approximation, the Perdew-Burke-Ernzerhof generalized gradient approximation, and the strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation. We demonstrate that SCAN describes this liquid with better accuracy than the other often-used functionals because it can simultaneously capture covalent and metallic bonding with similar high accuracy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Electron impact ionization in dense plasmas

Here, the distorted-wave with exchange (DWE) method is employed to evaluate electron impact ionization cross sections in dense electron–ion plasmas. Bound and continuum electron wave functions are obtained from a non-relativistic average-atom code. Plots of DWE cross sections are presented for 1s electrons in Li and Be plasmas and 2p electrons in Na and Mg plasmas. For each of these elements, cross sections are evaluated at metallic density in a range of temperatures from 10 to 100 eV and at their respective melting points. Resonances in the cross sections appear near the incident energy threshold at high temperatures. The origin of these resonances is discussed. In general, the distorted wave (DW) method without exchange is found to be a good approximation to the DWE method for electron impact ionization calculations. In the resonance region, however, exchange effects are found to be very important and cannot be neglected.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

Accuracy and convergence properties of the fixed-stress iterative solution of two-way coupled poromechanics

Here, the paper deals with the numerical solution of Biot's equations of coupled consolidation obtained by a mixed formulation combining continuous Galerkin finite–element and multipoint flux approximation finite–volume methods. The solution algorithm relies on the recently developed fixed–stress solution scheme, in which first the flow problem and then the mechanical one are addressed iteratively. We show that the algorithm can be interpreted as a particular block triangular preconditioning strategy applied within a Richardson iteration. The key component to the success of the preconditioner is the sparse approximation to the Schur complement based on a pressure space mass matrix scaled by a weighting factor that depends element–wise on the inverse of a suitable bulk modulus. The accuracy of the method is assessed, making use of well–known analytical solutions from the literature. Numerical results demonstrate robustness and low computational cost of the fixed–stress scheme in accurately capturing the two–way coupling between deformation and pressure.

58 GEOSCIENCES↗
DOE OSTI2015

Dispersion characteristics of planar grating with arbitrary grooves for terahertz Smith-Purcell radiation

In this paper, a novel method of getting the dispersion relations in planar grating with arbitrary grooves for terahertz Smith-Purcell radiation is investigated analytically. The continuous profile of the groove is approximately replaced by a series of rectangular steps. By making use of field matches method and the continuity of transverse admittance, the universal dispersion equation for grating with arbitrarily shaped grooves is derived. By solving the dispersion equation in presence of electron beam, the growth rate is obtained directly and the dependence on beam parameters is analyzed. Comparisons of the dispersion characteristics among some special groove shapes have been made by numerical calculation. The results show that the rectangular-step approximation method provides a novel approach to obtain the universal dispersion relation for grating with arbitrary grooves for Smith-Purcell radiation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

A modified Monte Carlo 'local importance function transform' method

The Local Importance Function Transform (LIFT) method uses an approximation of the contribution transport problem to bias a forward Monte-Carlo (MC) source-detector simulation [1-3]. Local (cell-based) biasing parameters are calculated from an inexpensive deterministic adjoint solution and used to modify the physics of the forward transport simulation. In this research, we have developed a new expression for the LIFT biasing parameter, which depends on a cell-average adjoint current to scalar flux (J{sup *}/{phi}{sup *}) ratio. This biasing parameter differs significantly from the original expression, which uses adjoint cell-edge scalar fluxes to construct a finite difference estimate of the flux derivative; the resulting biasing parameters exhibit spikes in magnitude at material discontinuities, causing the original LIFT method to lose efficiency in problems with high spatial heterogeneity. The new J{sup *}/{phi}{sup *} expression, while more expensive to obtain, generates biasing parameters that vary smoothly across the spatial domain. The result is an improvement in simulation efficiency. A representative test problem has been developed and analyzed to demonstrate the advantage of the updated biasing parameter expression with regards to solution figure of merit (FOM). For reference, the two variants of the LIFT method are compared to a similar variance reduction method developed by Depinay [4, 5], as well as MC with deterministic adjoint weight windows (WW). (authors)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Electron-ion collision-frequency for x-ray Thomson scattering in dense plasmas

Two methods are presented to calculate the electron-ion collision-frequency in dense plasmas using an average-atom model. The first one is based on the Kubo-Greenwood approach. The second one uses the Born and Lenard-Balescu approximations. The two methods are used to calculate x-ray Thomson scattering spectra. Illustrations are shown for dense beryllium and aluminum plasmas. Comparisons with experiment are presented in the case of an x-ray Thomson scattering spectrum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

A Fresh Math Perspective Opens New Possibilities for Computational Chemistry

By reformulating the TDDFT problem as a matrix function approximation, making use of a special transformation and taking advantage of the underlying symmetry with respect to a non-Euclidean metric, Yang and his colleagues were able to apply the Lanczos algorithm and a Kernal Polynomial Method (KPM) to approximate the absorption spectrum of several molecules. Both of these algorithms require relatively low-memory compared to non-symmetrical alternatives, which is the key to the computational savings.

Vu, Linda↗
DOE OSTI2024

A data-driven sensor placement approach for detecting voltage violations in distribution systems

Stochastic fluctuations in power injections from distributed energy resources (DERs) combined with load variability can cause constraint violations (e.g., exceeded voltage limits) in electric distribution systems. To monitor grid operations, sensors are placed to measure important quantities such as the voltage magnitudes. Here, in this paper, we consider a sensor placement problem which seeks to identify locations for installing sensors that can capture all possible violations of voltage magnitude limits. We formulate a bilevel optimization problem that minimizes the number of sensors and avoids false sensor alarms in the upper level while ensuring detection of any voltage violations in the lower level. This problem is challenging due to the nonlinearity of the power flow equations and the presence of binary variables. Accordingly, we employ recently developed conservative linear approximations of the power flow equations that overestimate or underestimate the voltage magnitudes. By replacing the nonlinear power flow equations with conservative linear approximations, we can ensure that the resulting sensor locations and thresholds are sufficient to identify any constraint violations. Additionally, we apply various problem reformulations to significantly improve computational tractability while simultaneously ensuring an appropriate placement of sensors. Lastly, we improve the quality of the results via an approximate gradient descent method that adjusts the sensor thresholds. We demonstrate the effectiveness of our proposed method for several test cases, including a system with multiple switching configurations.

24 POWER TRANSMISSION AND DISTRIBUTION↗