Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Greens function”

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,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

A Kinematic Rupture Model Generator Using Irikura's Recipe

In this project we developed GEN_SRF4 a computer program for generating kinematic rupture models, compatible with the SRF format, using Irikura and Miyake (2011) asperity-­based earthquake rupture model (IM2011, hereafter). IM2011, also known as Irkura’s recipe, has been widely used to model and simulate ground motion from earthquakes in Japan. An essential part of the method is its kinematic rupture generation technique, which is based on a deterministic rupture asperity modeling approach. The source model simplicity and efficiency of IM2011 at reproducing ground motion from earthquakes recorded in Japan makes it attractive to developers and users of the Southern California Earthquake Center Broadband Platform (SCEC BB platform). Besides writing the code the objective of our study was to test the transportability of IM2011 to broadband simulation methods used by the SCEC BB platform. Here we test it using the Graves and Pitarka (2010) method, implemented in the platform. We performed broadband (0.1- -10 Hz) ground motion simulations for a M6.7 scenario earthquake using rupture models produced with both GEN_SRF4 and rupture generator of Graves and Pitarka (2016), (GP2016 hereafter). In the simulations we used the same Green’s functions, and same high frequency approach for calculating the low-­frequency and high-­frequency parts of ground motion, respectively.

58 GEOSCIENCES↗
DOE OSTI2017

Manipulating Light with Transition Metal Clusters, Organic Dyes, and Metal Organic Frameworks

The primary goals of our research program is to develop and apply state-of-the-art first-principles methods to predict electronic and optical properties of three systems of significant scientific and technological interest: transition metal clusters, organic dyes, and metal-organic frameworks. These systems offer great opportunities to manipulate light for a wide ranging list of energy-related scientific problems and applications. During this grant period, we focused our investigations on the development, implementation, and benchmarking of many-body Green’s function methods (GW approximation and the Bethe-Salpeter equation) to examine excited-state properties of transition metal/transition-metal-oxide clusters and organic molecules that comprise the building blocks of dyes and metal-organic frameworks.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Analytical Solutions of Transient Drift-Diffusion in P-N Junction Pixel Sensors

Radiation detection in applications ranging from high energy physics to medical imaging rely on solid state detectors, often hybrid pixel detectors with (1) reverse biased p-n junction pixel sensors and (2) readout ASICs, attached by flip chip-bonding. Transient signals characteristics are important in, e.g., matching ASIC and sensor design, modeling and optimizing detector parameters and describing timing and charge sharing properties. Currently analytical forms of transient signals are available for only a few limited cases (e.g., drift or diffusion) or for the steady state (which is not relevant for high energy radiation detection). Tools are available for (relatively slow) numerical evaluation of the transient charge transport. We present here the first analytical solutions of partial differential equations de scribing drift-diffusion-recombination charge transport in planar p-n junction sensors in a variety of conditions: (1) undepleted, (2) fully depleted, (3) taking into account the gradual velocity saturation, and (4) overdepleted. We deduce the Green’s functions which can be applied to any detection problem through simple convolution with the initial conditions. We compare the analytical solutions with Monte Carlo simulations and industry standard simulations (Synopsys Sentaurus), demonstrating good agreement. Using the analytical equations enables fast modeling of the influence of various detector parameters on tracking, imaging and timing performance, describing performance and enabling optimizations for different applications. Finally, we illustrate this model with applications in 3D+T (x,y,z,time) photon tracking and 4D+T (x,y,θ,φ,time) relativistic charged particle tracking.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

The use of atmospheric prediction models to invert infrasound for linear-equivalent time domain moment tensors: Source Physics Experiment, Phase 1

We invert far field infrasound data for the equivalent seismo-acoustic time domain moment tensor to assess the effects of variable atmospheric models as well as to quantify the relative contributions of two presumed source phenomena. The infrasound data was produced by a series of underground chemical explosions that were conducted during the Source Physics Experiment, (SPE) which was originally designed to study explosion-generated seismo-acoustic signal phenomena. The goal of the work presented herein is two-fold: the first goal is to investigate the sensitivity of the estimated time domain moment tensors to variability of the estimated atmospheric model. The second goal is to determine the relative contribution of two possible source mechanisms to the observed infrasonic wave field. Rather than using actual atmospheric observations to estimate the necessary atmospheric Green's functions, we build a series of atmospheric models that rely on publicly available, regional atmospheric observations and the assumption that the acoustic energy results from a linear combination of an underground isotropic explosion and surface spall. The atmospheric observations are summarized and interpolated onto a 3D grid to produce a model of sound speed at the time of the experiment. For each of four SPE acoustic datasets that we invert, we produced a suite of three atmospheric models, based on ten years of regional meteorological observations: an average model, which averages the atmospheric conditions for ten years prior to each SPE event, as well as two extrema models. We find that the inversion yields relatively repeatable results for the estimated spall source. Conversely, the estimated isotropic explosion source is highly variable. This suggests that the majority of the observed acoustic energy is produced by the spall source and/or our modeling of the elastic energy propagation, and it's subsequent conversion to acoustic energy via linear elastic-to-acoustic coupling at the free surface, is too simplistic.

58 GEOSCIENCES↗
DOE OSTI2018

Building and Running TDAAPS Models: WRF Postdictions

This document serves to guide a researcher through the process of running the Weather Research and Forecasting (WRF) model and incorporating observations into coarse resolution reanalysis products to model atmospheric conditions at high (50 m) resolution. This documentation is specific to WRF and the WRF Preprocessing System (WPS) version 3.8.1 and the Objective Analysis (OBSGRID) code released on April 8, 2016. Output from WRF serves as an input into the Time-Domain Atmospheric Acoustic Propagation Suite (TDAAPS) which performs staggered-grid finite difference modeling of the acoustic velocity pressure system to produce Green's functions through these atmospheric models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Quantum Monte Carlo Calculations in Nuclear Theory

To a remarkable extent, atomic nuclei can be described as collections of point-like particles whose dynamics are dictated by a nonrelativistic Hamiltonian involving two- and three-nucleon potentials.Ab initio approaches are aimed at solving the many-body Schr$ \overline{o} $dinger equation associated with the nuclear Hamiltonian, which is a highly non-trivial problem because of the nonperturbative nature and the strong spin-isospin dependence of the nuclear forces. Within this framework it is possible to disentangle the theoretical uncertainty coming from modeling the nuclear potential and currents from that due to the approximations that are usually inherent in many-body techniques.Our many-body method of choice is quantum Monte Carlo (QMC), in particular Green’s function Monte Carlo (GFMC), which allows solving the nuclear Schrodinger equation with the required 1%accuracy level for both the ground- and the low-lying excited states of A ≤ 12 (A is the number of nucleons) nuclei. Since 2000 we have been achieving excellent agreement with experiment for nuclei with increasing number of nucleons by using the Argonne V18 two-nucleon interaction and the Illinois three-nucleon interactions (AV18+IL7). During this period many other groups have been developing chiral effective field theory (χEFT) potentials; these are generally nonlocal and notsuitable for QMC methods. We and collaborators have recently developed a local χEFT potential including the Δ-isobar degrees of freedom which can be used with QMC.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Quantification of MagLIF Morphology using the Mallat Scattering Transformation

The morphology of the stagnated plasma resulting from Magnetized Liner Inertial Fusion (MagLIF) is measured by imaging the self-emission x-rays coming from the multi-keV plasma, and the evolution of the imploding liner is measured by radiographs. Equivalent diagnostic response can be derived from integrated rad-MHD simulations from programs such as Hydra and Gorgon. There have been only limited quantitative ways to compare the image morphology, that is the texture, of simulations and experiments. We have developed a metric of image morphology based on the Mallat Scattering Transformation (MST), a transformation that has proved to be effective at distinguishing textures, sounds, and written characters. This metric has demonstrated excellent performance in classifying ensembles of synthetic stagnation images. We use this metric to quantitatively compare simulations to experimental images, cross experimental images, and to estimate the parameters of the images with uncertainty via a linear regression of the synthetic images to the parameter used to generate them. This coordinate space has proved very adept at doing a sophisticated relative back-ground subtraction in the MST space. This was needed to compare the experimental self emission images to the rad-MHD simulation images. We have also developed theory that connects the transformation to the causal dynamics of physical systems. This has been done from the classical kinetic perspective and from the field theory perspective, where the MST is the generalized Green's function, or S-matrix of the field theory in the scale basis. From both perspectives the first order MST is the current state of the system, and the second order MST are the transition rates from one state to another. An efficient, GPU accelerated, Python implementation of the MST was developed. Future applications are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Acoustic Source Parameters and Yield Estimation for SPE Phase I and Forensics Surface Events

Acoustic signals recorded for explosive events are sensitive to pressure changes in the source region and can provide valuable information about source characteristics. However, acoustic wave propagation in the atmosphere is significantly affected by weather conditions (e.g., temperature, pressure, and wind), and the atmospheric propagation effects should be considered properly to obtain accurate source information. Acoustic waveform inversion method has been increasing useful to determine quantitative source parameters by incorporating accurate Green’s function into the inversion and by effectively removing propagation path effects. This paper reviews recent progress in the acoustic waveform inversion and its application to explosion yield/emplacement (e.g., height-of-burst or depth-of-burial) estimation. The inversion method is applied to acoustic signals recorded during Source Physics Experiment Phase I and Forensics Surface Event conducted at the Nevada National Security Site (NNSS). Inverted acoustic source strengths are related to the explosion yield and emplacement, and their estimation uncertainties are explored by the comparison with the ground-truth data.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2020

Atomistic Studies of Individual Impurities and Impurity Complexes in III-V Semiconductors. Final report

The central hypothesis of our investigation is that reduced symmetries significantly influence the properties of individual impurities and impurity complexes in semiconductors. To disentangle the properties of the local environments of these impurities and impurity complexes from the behavior of an ensemble of impurities, at U. Iowa we performed theoretical calculations of specific local environments and developed the theoretical technology to calculate spin and charge transport through those local environments, in which the contact is a localized probe (scanning tunneling microscopy tip). These calculations were designed to be brought into comparison with experimental measurements made by Prof. Gupta’s group at Ohio State University. The principal environment that we focused on is the (110) surface of GaAs and other III-V semiconductors, because pristine surfaces can be revealed through a natural cleavage plane. We studied how the p-d hybridization of transition-metal dopants is influenced by the energetic location and the local environment’s symmetry-breaking of the d levels, and how the 4f-shell states of rare earth impurities couple to the surroundings. We calculated long-distance properties using tight-binding (linear combination of atomic orbitals, LCAO) Green’s function methods. The spin character of the dopants is often coherent on the time-scale of electronic transport through the region, and so correlations of the spin character of a dopant with that of a spinpolarized STM tip, or another dopant in the vicinity, can influence the measured conductivity from STM. We constructed a robust theory for these effects that can be used to describe quantitatively the measurements made by spin-polarized STM.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Multiscale Approach to Fast ModSim for Laser Processing of Metals for Future Nuclear Deterrence Environments

Predicting performance of parts produced using laser-metal processing remains an out- standing challenge. While many computational models exist, they are generally too computationally expensive to simulate the build of an engineering-scale part. This work develops a reduced order thermal model of a laser-metal system using analytical Green's function solutions to the linear heat equation, representing a step towards achieving a full part performance prediction in an "overnight" time frame. The developed model is able to calculate a thermal history for an example problem 72 times faster than a traditional FEM method. The model parameters are calibrated using a non-linear solution and microstructures and residual stresses calculated and compared to a non-linear case. The calibrated model shows promising agreement with a non-linear solution.

36 MATERIALS SCIENCE↗
DOE OSTI2018

The Relative Importance of Assumed Infrasound Mechanisms on the Linear Inversion of Infrasound Time Series at the Source Physics Experiment

We invert far field infrasound data for the equivalent seismo-acoustic time domain moment tensor to assess the relative importance of two assumed seismoacoustic source mechanisms. The infrasound data were produced by a four of the underground chemical explosions that were conducted during the Source Physics Experiment (SPE). For each SPE event that we invert, we produce three set of atmospheric Green's functions: an average model based on ten years of atmospheric data, as well as two extrema models designed to maximize the variability of atmospheric conditions for the given time-of-day and day-of-year for each SPE event. To parameterize the inversion, we assume that the source of infrasonic energy results from the linear combination of explosion-induced surface spall and linear seismic-to-elastic mode conversion at the Earth's free surface. We find that the inversion yields relatively repeatable results for the estimated spall source whereas the estimated isotropic explosion source is highly variable. This suggests that the majority of the observed acoustic energy is produced by the spall source and/or our modeling of the elastic energy propagation, and data are subsequent conversion to acoustic energy via linear elastic-to-acoustic coupling at the surface, is too simplistic.

58 GEOSCIENCES↗
DOE OSTI2018

Building and Running TDAAPS Models: WRF Predictions

This document serves to guide a researcher through the process of predicting atmospheric conditions in a region of interest utilizing the Weather Research and Forecasting (WRF) model. This documentation is specific to WRF and WRF Preprocessing System (WPS) version 3.8.1. WRF is an atmospheric prediction system designed for meteorological research and numerical atmospheric prediction. In WRF, simulations may be generated utilizing real data or idealized atmospheric conditions. Output from WRF serves as input into the Time-Domain Atmospheric Acoustic Propagation Suite (TDAAPS) which performs staggered-grid finite difference modeling of the acoustic velocity pressure system to produce Green's functions through these atmospheric models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Nuclear Theory from First Principles to Forefront Experiments (Final Report for DE-SC0018638)

The Lee research group is a part of the Nuclear Lattice Effective Field Theory Collaboration, which has developed and performed ab initio lattice simulations of nuclear structure, scattering, reactions, and many-body systems. The other senior members of the collaboration include Ulf-G. Meißner at Bonn/Julich, Evgeny Epelbaum and Hermann Krebs at Bochum, Timo Lahde and Thomas Luu at Julich, and Gautam Rupak at Mississippi State. Our letter “Ab initio alpha-alpha scattering” was featured in a Nature News and Views article. Another letter “Nuclear binding near a quantum phase transition” was highlighted in a Viewpoint article in the online APS journal Physics as well as a news article in the IOP publication Physics World (September 21, 2016). Our letter “Eigenvector continuation with subspace learning” was also highlighted a Synopsis article in Physics. There are many promising ab initio approaches being used to calculate the properties of few-and many-nucleon systems. This includes the no-core shell model, symmetry-adapted no-core shell model quantum Monte Carlo, auxiliary-field diffusion Monte Carlo, self-consistent Green’s functions, many-body perturbation theory, in-medium similarity renormalization group, and coupled cluster methods. Lattice effective field theory is another ab initio approach which combines the framework of effective field theory with lattice Monte Carlo methods to allow favorable scaling from few- to many-body systems. Perhaps the most important aspect of lattice effective field theory is that its strengths and weaknesses are orthogonal to that of other ab initio methods. For example, lattice effective field theory has little difficulty in probing cluster structures and collectivity, while such features are much more difficult using other methods. Furthermore it can be used to compute superfluid condensate fractions as well as the phase diagram of strongly matter and the density and temperature dependence of clustering. Lattice effective field theory was first used in simulations of infinite nuclear matter and infinite neutron matter with pions and without pions. In addition to the efforts by our collaboration, there have been recent efforts by other groups as well.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

On the fractional Laplacian of variable order

We present a novel definition of variable-order fractional Laplacian on R n based on a natural generalization of the standard Riesz potential. Our definition holds for values of the fractional parameter spanning the entire open set (0, n/2). We then discuss some properties of the fractional Poisson’s equation involving this operator and we compute the corresponding Green’s function, for which we provide some instructive examples for specific problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

An Overview of Gemma FY20 Verification Activities

Gemma verification activities for FY20 can be divided into three categories: the development of specialized quadrature rules, initial progress towards the development of manufactured solutions for code verification, and automated code-verification testing. In the method-of-moments implementation of the electric-field integral equation, the presence of a Green’s function in the four-dimensional integrals yields singularities in the integrand when two elements are nearby. To address these challenges, we have developed quadrature rules to integrate the functions through which the singularities can be characterized. Code verification is necessary to develop confidence in the implementation of the numerical methods in Gemma. Therefore, we have begun investigating the use of manufactured solutions to more thoroughly verify Gemma. Manufactured solutions provide greater flexibility for testing aspects of the code; however, the aforementioned singularities provide challenges, and existing work is limited in rigor and quantity. Finally, we have implemented automated code-verification testing using the VVTest framework to automate the mesh refinement and execution of a Gemma simulation to generate mesh convergence data. This infrastructure computes the observed order of accuracy from these data and compares it with the theoretical order of accuracy to either develop confidence in the implementation of the numerical methods or detect coding errors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

3D orthorhombic earth model effects on seismic source characterization

Most earth materials are anisotropic with regard to seismic wave-speeds, especially materials such as shales, or where oriented fractures are present. However, the base assumption for many numerical simulations is to treat earth materials as isotropic media. This is done for simplicity, the apparent weakness of anisotropy in the far field, and the lack of well-characterized anisotropic material properties for input into numerical simulations. One approach for addressing the higher complexity of actual geologic regions is to model the material as an orthorhombic medium. We have developed an explicit time-domain, finite-difference (FD) algorithm for simulating three-dimensional (3D) elastic wave propagation in a heterogeneous orthorhombic medium. The objective of this research is to investigate the errors and biases that result from modeling a non-isotropic medium as an isotropic medium. This is done by computing “observed data” by using synthetic, anisotropic simulations with the assumption of an orthorhombic, anisotropic earth model. Green’s functions for an assumed isotropic earth model are computed and then used an inversion designed to estimate moment tensors with the “observed” data. One specific area of interest is how shear waves, which are introduced in an anisotropic model even for an isotropic explosion, affect the characterization of seismic sources when isotropic earth assumptions are made. This work is done in support of the modeling component of the Source Physics Experiment (SPE), a series of underground chemical explosions at the Nevada National Security Site (NNSS).

58 GEOSCIENCES↗
DOE OSTI2021

An Overview of Gemma FY2021 Verification Activities

Though the method-of-moments implementation of the electric-field integral equation plays an important role in computational electromagnetics, it provides many code-verification challenges due to the different sources of numerical error and their possible interactions. Matters are further complicated by singular integrals, which arise from the presence of a Green's function. In this report, we document our research to address these issues, as well as its implementation and testing in Gemma.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

An Infinite Domain 3D Poisson Solver Based on the Barnes-Hut Algorithm

We present a domain decomposition method for the solution of the 3D Poisson equation with infinite domain boundary conditions. The method is based on an application of the Barnes-Hut tree particle scheme adapted to gridded data. Long range interactions are computed using the first two terms in the Cartesian multipole expansion of Green’s function convoluted with the charge while short range computations are performed using Hockney’s domain doubling algorithm. A standard domain decomposition strategy requires O(N 2 ) applications of Hockney’s algorithm, where N is the number of subdomains that intersect that charge support, while in the present approach only O(Nlog 2 N) such computations suffice. The discretization scheme employed is a sixth order Mehrstellen approximation of the 3D Laplace opera tor. The method exhibits satisfactory accuracy at a substantially reduced computational cost compared to the full domain decomposition Hockney’s algorithm.

97 MATHEMATICS AND COMPUTING↗