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 1,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Superconvergence results of Legendre spectral projection methods for Volterra integral equations of second kind

In this paper, Legendre spectral projection methods are applied for the Volterra integral equations of second kind with a smooth kernel. We prove that the approximate solutions of the Legendre Galerkin and Legendre collocation methods converge to the exact solution with the order O(n{sup −r}) in L{sup 2}-norm and order O(n{sup −r+(1/2)}) in infinity norm, and the iterated Legendre Galerkin solution converges with the order O(n{sup −2r}) in both L{sup 2}-norm and infinity norm, whereas the iterated Legendre collocation solution converges with the order O(n{sup −r}) in both L{sup 2}-norm and infinity norm, n being the highest degree of Legendre polynomials employed in the approximation and r being the smoothness of the kernels. We have also considered multi-Galerkin method and its iterated version, and prove that the iterated multi-Galerkin solution converges with the order O(n{sup −3r}) in both infinity and L{sup 2} norm. Numerical examples are given to illustrate the theoretical results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Neural network approximation of Bayesian models for the inference of ion and electron temperature profiles at W7-X

Here, we describe a method for training a neural network (NN) to approximate the full model Bayesian inference of plasma profiles from x-ray imaging diagnostic measurements. The modeling is carried out within the Minerva Bayesian modeling framework where models are defined as a set of assumptions, prior beliefs on parameter values and physics knowledge. The goal is to use NNs for fast ion and electron temperature profile inversion from measured image data. The NN is trained solely on artificial data generated by sampling from the joint distribution of the free parameters and model predictions. The training is carried out in such a way that the mapping learned by the network constitutes an approximation of the full model Bayesian inference. The analysis is carried out on images constituted of 20 × 195 pixels corresponding to binned lines of sight and spectral channels, respectively. Through the full model inference, it is possible to infer electron and ion temperature profiles as well as impurity density profiles. When the network is utilized for the inference of the temperature profiles, the analysis time can be reduced down to a few tens of microseconds for a single time point, which is a drastic improvement if compared to the ≈4 h long Bayesian inference. The procedure developed for the generation of the training set does not rely on diagnostic-specific features, and therefore it is in principle applicable to any other model developed within the Minerva framework. The trained NN has been tested on data collected during the first operational campaign at W7-X, and compared to the full model Bayesian inference results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Parametrically designed surface topography on CAD models of additively manufactured lattice structures for improved design validation

Additively Manufactured (AM) lattice structures of the laser powder bed fusion (LPBF) process generally have different surface geometries depending on the overhang angle and location. This means that the design validation stage is often challenging, considering it is based on an ideal Computer Aided Design (CAD) model that is not truly representative, since the AM part will be different in terms of dimensional accuracy and surface finish. Previous studies have relied on the design of surface textures that are independent of the overhang angle, or techniques based on directly using X-ray Computed Tomography (XCT) data. In this paper, a new technique for designing surface texture on the CAD of lattices has been investigated and correlated with areal surface roughness parameters. After extracting areal surface parameters from the XCT data of the lattice, the method allowed for the design of this surface texture on the part CAD taking in consideration if its up skin or down skin, without using XCT data as the main input. By applying this method, it is possible to obtain a lattice CAD model with a designed surface texture and geometry that is more representative of the actual AM lattice. The mean deviation between the CAD model with the designed surface and the XCT was a third of the one between the XCT and the initial CAD. The proposed method allows for designing and replicating LPBF AM surfaces on the CAD of a lattice, taking into consideration the dimensional deviation caused by AM surfaces, especially on overhangs. Assuming that an LPBF AM process is stable and produces approximately the same AM surface, this method can be used to predict the geometry of a lattice, providing a cost-effective and efficient methodology that minimizes the time needed for design validation.

42 ENGINEERING↗
DOE OSTI2016

Nonlinear bioheat transfer models and multi-objective numerical optimization of the cryosurgery operations

Numerical computation of the three dimensional problem of the freezing interface propagation during the cryosurgery coupled with the multi-objective optimization methods is used in order to improve the efficiency and safety of the cryosurgery operations performing. Prostate cancer treatment and cutaneous cryosurgery are considered. The heat transfer in soft tissue during the thermal exposure to low temperature is described by the Pennes bioheat model and is coupled with an enthalpy method for blurred phase change computations. The finite volume method combined with the control volume approximation of the heat fluxes is applied for the cryosurgery numerical modeling on the tumor tissue of a quite arbitrary shape. The flux relaxation approach is used for the stability improvement of the explicit finite difference schemes. The method of the additional heating elements mounting is studied as an approach to control the cellular necrosis front propagation. Whereas the undestucted tumor tissue and destucted healthy tissue volumes are considered as objective functions, the locations of additional heating elements in cutaneous cryosurgery and cryotips in prostate cancer cryotreatment are considered as objective variables in multi-objective problem. The quasi-gradient method is proposed for the searching of the Pareto front segments as the multi-objective optimization problem solutions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Covariance matrices for variance-suppressed simulations

ABSTRACT Cosmological N-body simulations provide numerical predictions of the structure of the Universe against which to compare data from ongoing and future surveys, but the growing volume of the Universe mapped by surveys requires correspondingly lower statistical uncertainties in simulations, usually achieved by increasing simulation sizes at the expense of computational power. It was recently proposed to reduce simulation variance without incurring additional computational costs by adopting fixed-amplitude initial conditions. This method has been demonstrated not to introduce bias in various statistics, including the two-point statistics of galaxy samples typically used for extracting cosmological parameters from galaxy redshift survey data, but requires us to revisit current methods for estimating covariance matrices of clustering statistics for simulations. In this work, we find that it is not trivial to construct covariance matrices analytically for fixed-amplitude simulations, but we demonstrate that ezmock (Effective Zel’dovich approximation mock catalogue), the most efficient method for constructing mock catalogues with accurate two- and three-point statistics, provides reasonable covariance matrix estimates for such simulations. We further examine how the variance suppression obtained by amplitude-fixing depends on three-point clustering, small-scale clustering, and galaxy bias, and propose intuitive explanations for the effects we observe based on the ezmock bias model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Path Integral Monte Carlo Simulations of Iron Plasmas (Final Technical Report)

This documents is the final technical report for our grant entitled "Path Integral Monte Carlo Simulations of Iron Plasmas" that focused on developing path integral Monte Carlo (PIMC) computer simulations. This techniques will be developed to study plasmas composed of heavier elements including iron and other third row elements. Equations of state (EOS) and transport properties will be derived in the regime of warm dense matter (WDM) and dense plasmas where existing first-principles methods cannot be applied. While standard density functional theory (DFT) has been used to accurately predict the structure of many solids and liquids up to temperatures on the order of 100,000 K, this method is not applicable at much higher temperature because the number of partially occupied electronic orbitals reaches intractably large numbers or the use of finite-temperature free energy functionals in orbital-free DFT introduces an uncontrolled approximation. Here we focus on PIMC methods that become more and more efficient with increasing temperatures and still include all electronic correlation effects. In this approach, electronic excitations increase the efficiency rather than reduce it. While it had commonly been assumed this method could only be applied to elements without core electrons, we showed that PIMC with free-particle nodes works well for first-row elements (PRL 108 (2012) 115502). Most recently, we extended the applicability range of all-electron PIMC to second-row elements by adopting localized nodal surfaces (PRL 115 (2015) 176403). To simulate third-row elements efficiently under WDM conditions, we propose a new method to remove core electrons by introducing pseudo-nodes. We explain our approach step by step and present preliminary results. We focus our method development on getting PIMC simulations of iron to work because of its fundamental importance for WDM and astrophysics. Then we move on to krypton and copper-doped beryllium, a ICF ablator material. We plan to continue working on key second-row material such as Na, Mg, MgO, Al, silica, and silicon-doped plastic ablators. Our collaborators at LLNL, will use our PIMC EOS data both as comparisons to existing semi-empirical, EOS-generating schemes, and as input for continuum radiation hydrodynamics simulations. We will establish an efficient pipeline from PIMC to macroscopic continuum studies of materials response. An emphasis will be placed on benchmarking such methods for plasmas of heavy elements at the very high temperatures (~100 eV) and low densities that are generated when Hohlraum radiation heats the ablator material in indirect drive laser experiments. Results from changes to the EOS will be of immeasurable importance to the designers at the National Ignition Facility (NIF) and at other facilities. Starting with our EOS of Cu-doped Be, our second collaborator at LLE, will perform real-time simulations of laser fusion experiments at the Omega laser and at the NIF to determine how sensitive the compression path depends on the ablator EOS. Since our collaborator also has experience in performing orbital-free DFT calculations, we propose to compare predictions from this method with PIMC results. In joint publications, we plan to analyze the accuracy of different free-energy functionals in order to understand why existing orbital-free DFT calculations do not predict compression peaks along the shock Hugoniot curve that we see with PIMC. The peaks are caused by the ionization of various electron shells. Their accurate characterization is important to compare with experimental results. We will break new ground by developing PIMC techniques that can simulate iron and all other third row elements in the plasma and WDM regimes. We introduce the concept of pseudo-nodes for the efficient treatment of the core-electrons. The EOS and transport properties will be derived and published online in the form of a new WDM database. Our PIMC EOS calculations will benchmark and possibly replace semi-analytical EOS tables like QEOS or SESAME, which will impact the hydrocode simulation community and will affect the design of NIF targets. Our PIMC results will help to improve the accuracy of orbital-free DFT simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Microscopic optical potentials for calcium isotopes

We construct nucleonic microscopic optical potentials by combining the Green's function approach with the coupled-cluster method for 40 Ca and 48 Ca. For the computation of the ground states of 40 Ca and 48 Ca, we use the coupled-cluster method in the singles-and-doubles approximation, while for the A = ±1 nuclei we use particle-attached and particle-removed equation-of-motion method truncated at two-particle–one-hole and one-particle–two-hole excitations, respectively. Our calculations are based on the chiral nucleon-nucleon and three-nucleon interaction NNLO sat , which reproduces the charge radii of 40 Ca and 48 Ca, and the chiral nucleon-nucleon interaction NNLO opt . In all cases considered here, we observe that the overall form of the neutron scattering cross section is reproduced for both interactions, but the imaginary part of the potential, which reflects the loss of flux in the elastic channel, is negligible. Furthermore, the latter points to neglected many-body correlations that would appear beyond the coupled-cluster truncation level considered in this work. We show that, by artificially increasing the parameter η in the Green's function, practical results can be further improved.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2017

A compatible high-order meshless method for the Stokes equations with applications to suspension flows

A stable numerical solution of the steady Stokes problem requires compatibility between the choice of velocity and pressure approximation that has traditionally proven problematic for meshless methods. Here in this work, we present a discretization that couples a staggered scheme for pressure approximation with a divergence-free velocity reconstruction to obtain an adaptive, high-order, finite difference-like discretization that can be efficiently solved with conventional algebraic multigrid techniques. We use analytic benchmarks to demonstrate equal-order convergence for both velocity and pressure when solving problems with curvilinear geometries. In order to study problems in dense suspensions, we couple the solution for the flow to the equations of motion for freely suspended particles in an implicit monolithic scheme. The combination of high-order accuracy with fully-implicit schemes allows the accurate resolution of stiff lubrication forces directly from the solution of the Stokes problem without the need to introduce sub-grid lubrication models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Bayesian prior construction for uncertainty quantification in first-principles statistical mechanics

First-principles statistical mechanics enables the prediction of thermodynamic and kinetic properties of materials, but is computationally expensive. Many approaches require surrogate models to calculate energies within Monte Carlo or molecular dynamics simulations. Inexpensive surrogates such as cluster expansions enable otherwise intractable calculations by interpolating data from higher accuracy methods, such as Density Functional Theory (DFT). Surrogate models introduce uncertainty into downstream calculations, in addition to any uncertainty inherent to DFT calculations. Bayesian frameworks address this by quantifying uncertainty and incorporating expert knowledge through priors. However, constructing effective priors remains challenging. This work introduces and describes practical strategies for building Bayesian cluster expansions, focusing on basis truncation, hyperparameter selection, and ground state replication. We analyze multiple basis truncation schemes, compare cross-validation to the evidence-approximation for hyperparameter optimization, and provide methods to find and enforce ground-state-preserving models through priors. Additionally, we compare the uncertainties between different approximations to DFT (LDA, PBE, SCAN) against the uncertainty introduced with the use of cluster expansion surrogate models. These approaches are demonstrated on the BCC Li x Mg 1-x and Li x Al 1-x alloys, which are both of interest for solid-state Li batteries. Our results provide guidelines for constructing and utilizing Bayesian cluster expansions, thereby improving the transparency of materials modeling. Furthermore, the approaches and insights developed in this work can be transferred to a wide range of cluster expansion surrogate models, including the atomic cluster expansion and related machine-learned interatomic potential architectures.

Alloy theory↗
DOE OSTI2014

Electronic excitation spectra of molecules in solution calculated using the symmetry-adapted cluster-configuration interaction method in the polarizable continuum model with perturbative approach

A perturbative approximation of the state specific polarizable continuum model (PCM) symmetry-adapted cluster-configuration interaction (SAC-CI) method is proposed for efficient calculations of the electronic excitations and absorption spectra of molecules in solutions. This first-order PCM SAC-CI method considers the solvent effects on the energies of excited states up to the first-order with using the zeroth-order wavefunctions. This method can avoid the costly iterative procedure of the self-consistent reaction field calculations. The first-order PCM SAC-CI calculations well reproduce the results obtained by the iterative method for various types of excitations of molecules in polar and nonpolar solvents. The first-order contribution is significant for the excitation energies. The results obtained by the zeroth-order PCM SAC-CI, which considers the fixed ground-state reaction field for the excited-state calculations, are deviated from the results by the iterative method about 0.1 eV, and the zeroth-order PCM SAC-CI cannot predict even the direction of solvent shifts in n-hexane for many cases. The first-order PCM SAC-CI is applied to studying the solvatochromisms of (2,2{sup ′}-bipyridine)tetracarbonyltungsten [W(CO){sub 4}(bpy), bpy = 2,2{sup ′}-bipyridine] and bis(pentacarbonyltungsten)pyrazine [(OC){sub 5}W(pyz)W(CO){sub 5}, pyz = pyrazine]. The SAC-CI calculations reveal the detailed character of the excited states and the mechanisms of solvent shifts. The energies of metal to ligand charge transfer states are significantly sensitive to solvents. The first-order PCM SAC-CI well reproduces the observed absorption spectra of the tungsten carbonyl complexes in several solvents.

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

Optimization problems governed by systems of PDEs with uncertainties

This paper reviews current theoretical and numerical approaches to optimization problems governed by partial differential equations (PDEs) that depend on random variables or random fields. Such problems arise in many engineering, science, economics and societal decision-making tasks. This paper focuses on problems in which the governing PDEs are parametrized by the random variables/fields, and the decisions are made at the beginning and are not revised once uncertainty is revealed. Examples of such problems are presented to motivate the topic of this paper, and to illustrate the impact of different ways to model uncertainty in the formulations of the optimization problem and their impact on the solution. A linear–quadratic elliptic optimal control problem is used to provide a detailed discussion of the set-up for the risk-neutral optimization problem formulation, study the existence and characterization of its solution, and survey numerical methods for computing it. Different ways to model uncertainty in the PDE-constrained optimization problem are surveyed in an abstract setting, including risk measures, distributionally robust optimization formulations, probabilistic functions and chance constraints, and stochastic orders. Furthermore, approximation-based optimization approaches and stochastic methods for the solution of the large-scale PDE-constrained optimization problems under uncertainty are described. Some possible future research directions are outlined.

Heinkenschloss, Matthias [Rice Univ., Houston, TX ↗
DOE OSTI2016

SU-F-T-661: Dependence of Gold Nano Particles Cluster Morphology On Dose Enhancement of Photon Radiation Therapy Apply for Radiation Biology Effect

Purpose: Injected gold nano particles (GNPs) to a body for dose enhancement are known to form in the tumorcell cluster morphology. We investigated the dependence of dose enhancement on the morphology characteristic with an approximated morphology model by using Monte Carlo simulations. Methods: For MC simulation, TOPAS version 2.0P-03 was used. GNP cluster morphology was approximated as a body center cubic(BCC) model by placing 8 GNPs at the corner and one at the center of cube with length from 2.59 µm to 0.25 µm located in a 4 µm length water filled cube phantom. 4 µm length square shaped beams of poly-energetic 50, 260 kVp photons were irradiated to the water filled cube phantom with 100 nm diameter GNPs in it. Dose enhancement ratio(DER) was computed as a function of distance from the surface of the GNP at the cube center for 18 cubes geometries. For scoring particles, 10 nm width of concentric shell shaped detector was constructed up to 100 nm from the center. Total dose in a sphere of 100 nm radius of detector were normalized to 2.59 µm length cube morphology. To verified biological effect of BCC model applied to cell survival curve fitting. Results: DER increase as the distance of the GNPs reduces. DER was largest for 0.25 µm length cube. Dependence of GNP distance DER increment was 1.73, 1.60 for 50 kVp, 260 kVp photons, respectively. Also, Using BCC model applied to cell survival curve was well prediction. Conclusion: DER with GNPs was larger when they are closely packed in the phantom. Therefore, better therapeutic effects can be expected with close-packed GNPs. This research was supported by the NRF funded by the Ministry of Science, ICT & Future Planning (2012M3A9B6055201 and 2012R1A1A2042414), Samsung Medical Center grant[GFO1130081].

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Quantum breathers in Heisenberg ferromagnetic chains with Dzyaloshinsky-Moriya interaction

We present an analytical study on quantum breathers in one-dimensional ferromagnetic XXZ chains with Dzyaloshinsky-Moriya interaction by means of the time-dependent Hartree approximation and the semidiscrete multiple-scale method. The stationary localized single-boson wave functions are obtained and these analytical solutions are checked by numerical simulations. With such stationary localized single-boson wave functions, we construct quantum breather states. Furthermore, the role of the Dzyaloshinsky-Moriya interaction is discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Ab Initio Investigation of Structural Stability and Electronic and Magnetic Properties of the Half-Heusler Alloys: MTiSb (M = Fe, Co, and Ni)

The structural, electronic, and magnetic properties of the half-Heusler MTiSb (M = Fe, Co, Ni) compounds have been examined using the first-principles calculations with the full-potential linear augmented plane wave (FP-LAPW) method within the density functional theory. The electronic structures and magnetic properties of MTiSb (M = Fe, Co, Ni) are studied. It is found that the calculated lattice parameter is in good agreement with the theoretical values. Using the general gradient approximation (GGA), we observe that for the three compounds, the α phase is more stable than β and γ phases. CoTiSb is a semiconductor, NiTiSb is a metal, while the FeTiSb compound is a half-metallic, and the total magnetic moments are 0.00085, − 0.26074, and 0.99353 μ B, respectively, which is in agreement with the Slater–Pauling rule, Furthermore, the origin for the appearance of the half-metallic band gap in the FeTiSb compound was also discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Systematically improvable multiscale solver for correlated electron systems

The development of numerical methods capable of simulating realistic materials with strongly correlated electrons, with controllable errors, is a central challenge in quantum many-body physics. Here we describe a framework for a general multiscale method based on embedding a self-energy of a strongly correlated subsystem into a self-energy generated by a method able to treat large weakly correlated systems approximately. As an example, we present the embedding of an exact diagonalization self-energy into a self-energy generated from self-consistent second-order perturbation theory. Using a quantum impurity model, generated from a cluster dynamical mean field approximation to the two-dimensional Hubbard model, as a benchmark, we illustrate that our method allows us to obtain accurate results at a fraction of the cost of typical Monte Carlo calculations. We test the method in multiple regimes of interaction strengths and dopings of the model. Here, the general embedding framework we present avoids difficulties such as double counting corrections, frequency-dependent interactions, or vertex functions. As it is solely formulated at the level of the single-particle Green's function, it provides a promising route for the simulation of realistic materials that are currently difficult to study with other methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Early results of microwave transmission experiments through an overly dense rectangular plasma sheet with microparticle injection

These experiments utilize a linear hollow cathode to create a dense, rectangular plasma sheet to simulate the plasma layer surrounding vehicles traveling at hypersonic velocities within the Earth's atmosphere. Injection of fine dielectric microparticles significantly reduces the electron density and therefore lowers the electron plasma frequency by binding a significant portion of the bulk free electrons to the relatively massive microparticles. Measurements show that microwave transmission through this previously overly dense, impenetrable plasma layer increases with the injection of alumina microparticles approximately 60 μm in diameter. This method of electron depletion is a potential means of mitigating the radio communications blackout experienced by hypersonic vehicles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A Feynman-Kac based numerical method for the exit time probability of a class of transport problems

The exit time probability, which gives the likelihood that an initial condition leaves a prescribed region of the phase space of a dynamical system at, or before, a given time, is arguably one of the most natural and important transport problems. In this work, we present an accurate and efficient numerical method for computing this probability for systems described by non-autonomous (time-dependent) stochastic differential equations (SDEs) or their equivalent Fokker-Planck partial differential equations. The method is based on the direct approximation of the Feynman-Kac formula that establishes a link between the adjoint Fokker-Planck equation and the forward SDE. The Feynman-Kac formula is approximated using the Gauss-Hermite quadrature rules and piecewise cubic Hermite interpolating polynomials, and a GPU accelerated matrix representation is used to compute the entire time evolution of the exit time probability using a single pass of the algorithm. The method is unconditionally stable, exhibits second order convergence in space, first order convergence in time, and it is straightforward to parallelize. Applications are presented to the advection diffusion of a passive tracer in a fluid flow exhibiting chaotic advection, and to the runaway acceleration of electrons in a plasma in the presence of an electric field, collisions, and radiation damping. Benchmarks against analytical solutions as well as comparisons with explicit and implicit finite difference standard methods for the adjoint Fokker-Planck equation are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗