Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical 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,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

BPS fivebrane stars. Part II. Fluctuations

We investigate quantum fluctuations of metric components in coherent 1/2-BPS bound states of n1 fundamental strings and n 5 NS5-branes. The leading order contribution in an expansion in 1/(n 1 n 5 ) is calculated via a combination of analytical and numerical methods. We find that the fluctuations are small away from a tiny distance from the source, comparable to the 6d Planck scale. Comparing this result with an analysis in the literature of fluctuations in the maximally mixed state, we conclude that the large fluctuations previously found for the latter are statistical rather than quantum in nature, and that perturbative string theory provides an accurate description of these backgrounds.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Continuum-mediated self-interacting dark matter

Dark matter may self-interact through a continuum of low-mass states. This happens if dark matter couples to a strongly-coupled nearly-conformal hidden sector. This type of theory is holographically described by brane-localized dark matter interacting with bulk fields in a slice of 5D anti-de Sitter space. The long-range potential in this scenario depends on a non-integer power of the spatial separation, in contrast to the Yukawa potential generated by the exchange of a single 4D mediator. The resulting self-interaction cross section scales like a non-integer power of velocity. We identify the Born, classical and resonant regimes and investigate them using state-of-the-art numerical methods. We demonstrate the viability of our continuum-mediated framework to address the astrophysical small-scale structure anomalies. Investigating the continuum-mediated Sommerfeld enhancement, we demonstrate that a pattern of resonances can occur depending on the non-integer power. We conclude that continuum mediators introduce novel power-law scalings which open new possibilities for dark matter self-interaction phenomenology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Complexity growth in integrable and chaotic models

We use the SYK family of models with N Majorana fermions to study the complexity of time evolution, formulated as the shortest geodesic length on the unitary group manifold between the identity and the time evolution operator, in free, integrable, and chaotic systems. Initially, the shortest geodesic follows the time evolution trajectory, and hence complexity grows linearly in time. We study how this linear growth is eventually truncated by the appearance and accumulation of conjugate points, which signal the presence of shorter geodesics intersecting the time evolution trajectory. By explicitly locating such “shortcuts” through analytical and numerical methods, we demonstrate that: (a) in the free theory, time evolution encounters conjugate points at a polynomial time; consequently complexity growth truncates at O($\sqrt{N}$), and we find an explicit operator which “fast-forwards” the free N-fermion time evolution with this complexity, (b) in a class of interacting integrable theories, the complexity is upper bounded by O(poly(N)), and (c) in chaotic theories, we argue that conjugate points do not occur until exponential times O(e N ), after which it becomes possible to find infinitesimally nearby geodesics which approximate the time evolution operator. Finally, we explore the notion of eigenstate complexity in free, integrable, and chaotic models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

De Sitter decays to infinity

Bubbles of nothing are a class of vacuum decay processes present in some theories with compactified extra dimensions. We investigate the existence and properties of bubbles of nothing in models where the scalar pseudomoduli controlling the size of the extra dimensions are stabilized at positive vacuum energy, which is a necessary feature of any realistic model. We map the construction of bubbles of nothing to a four-dimensional Coleman-De Luccia problem and establish necessary conditions on the asymptotic behavior of the scalar potential for the existence of suitable solutions. We perform detailed analyses in the context of five-dimensional theories with metastable dS 4 x S 1 vacua, using analytic approximations and numerical methods to calculate the decay rate. We find that bubbles of nothing sometimes exist in potentials with no ordinary Coleman-De Luccia decay process, and that in the examples we study, when both processes exist, the bubble of nothing decay rate is typically faster. Our methods can be generalized to other stabilizing potentials and internal manifolds.

classical theories of gravity↗
DOE OSTI2022

Nonperturbative bounds on scattering of massive scalar particles in d ≥ 2

We study two-to-two scattering amplitudes of a scalar particle of mass m. For simplicity, we assume the presence of Z 2 symmetry and that the particle is Z 2 odd. We consider two classes of amplitudes: the fully nonperturbative ones and effective field theory (EFT) ones with a cut-off scale M. Using the primal numerical method which allows us to impose full non-linear unitarity, we construct novel bounds on various observables in 2 ≤ d ≤ 4 space-time dimensions for both classes of amplitudes. We show that our bounds are much stronger than the ones obtained by using linearized unitarity or positivity only. We discuss applications of our bounds to constraining EFTs. Finally, we compare our bounds to the amplitude in Φ 4 theory computed perturbatively at weak coupling, and find that they saturate the bounds.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Compact feature-aware Hermite-style high-order surface reconstruction

High-order surface reconstruction is now a commonly used technique for mesh generation and adaptation, geometric and physical modeling, and high-order numerical methods for solving partial differential equations (PDEs). However, surface reconstruction from a relatively coarse mesh remained a challenging problem, especially for surfaces with sharp features. Here in this paper, we introduce a new method to address this challenge by improving the previous state of the art, including continuous moving frames (CMF) and weighted averaging of local fittings (WALF) (Eng Comput 28 (2012)), in two aspects. First, we significantly improve the robustness of reconstruction from coarse meshes using a Hermite-style least squares approximation to incorporate normals of the surface and tangents of the feature curves. Second, we ensure both G 0 continuity and high-order accuracy of the reconstruction near sharp features by generating parametric surface elements with an iterative feature-aware parameterization. We present the theoretical framework of our method and compare it against point-based methods in terms of accuracy and stability. We demonstrate the use of the proposed technique in generating high-order meshes for finite element methods, and show that it enables nearly identical solutions as using the meshes generated from exact geometry, while allowing additional flexibility.

42 ENGINEERING↗
DOE OSTI2014

Implementation of a thermomechanical model for the simulation of selective laser melting

Selective laser melting (SLM) is an additive manufacturing process in which multiple, successive layers of metal powders are heated via laser in order to build a part. Modeling of SLM requires consideration of both heat transfer and solid mechanics. The present work describes continuum modeling of SLM as envisioned for eventual support of part-scale modeling of this fabrication process to determine end-state information such as residual stresses and distortion. The determination of the evolving temperatures is dependent on the material, the state of the material (powder or solid), the specified heating, and the configuration. Similarly, the current configuration is dependent on the temperatures, the powder-solid state, and the constitutive models. A multi-physics numerical formulation is required to solve such problems. This article describes the problem formulation, numerical method, and constitutive parameters necessary to solve such a problem. Additionally, various verification and example problems are simulated in the parallel, multi-physics finite element code Diablo, and the results presented herein.

42 ENGINEERING↗
DOE OSTI2017

An algorithm for continuum modeling of rocks with multiple embedded nonlinearly-compliant joints [Continuum modeling of elasto-plastic media with multiple embedded nonlinearly-compliant joints]

Here, we present a numerical method for modeling the mechanical effects of nonlinearly-compliant joints in elasto-plastic media. The method uses a series of strain-rate and stress update algorithms to determine joint closure, slip, and solid stress within computational cells containing multiple “embedded” joints. This work facilitates efficient modeling of nonlinear wave propagation in large spatial domains containing a large number of joints that affect bulk mechanical properties. We implement the method within the massively parallel Lagrangian code GEODYN-L and provide verification and examples. We highlight the ability of our algorithms to capture joint interactions and multiple weakness planes within individual computational cells, as well as its computational efficiency. We also discuss the motivation for developing the proposed technique: to simulate large-scale wave propagation during the Source Physics Experiments (SPE), a series of underground explosions conducted at the Nevada National Security Site (NNSS).

42 ENGINEERING↗
DOE OSTI2019

Multidimensional sum-up rounding for integer programming in optimal experimental design

Here we present a numerical method for approximating the solution of convex integer programs stemming from optimal experimental design. The statistical setup consists of a Bayesian framework for linear inverse problems for which the direct relationship is described by a discretized integral equation. Specifically, we aim to find the optimal sensor placement from a set of candidate locations where data are collected with measurement error. The convex objective function is a measure of the uncertainty, described here by the trace or log-determinant of the posterior covariance matrix, for the discretized linear inverse problem solution. The resulting convex integer program is relaxed, producing a lower bound. An upper bound is obtained by extending the sum-up rounding approach to multiple dimensions. For this extension, we analyze its accuracy as a function of the discretization mesh size for a rectangular domain. We show asymptotic optimality of the integer solution defining the upper bound for different experimental design criteria (A- and D-optimal), by proving the convergence to zero of the gap between the upper and lower bounds as the mesh size goes to zero. The technique is illustrated on a two-dimensional gravity surveying problem for both A-optimal and D-optimal sensor placement where our designs yield better results compared with a thresholding rounding approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Modeling of Gas Dynamic Processes in the Interelectrode Gap of a Multielement Thermionic Electricity Generating Channel by Solving the Kinetic Equation

This work is devoted to the application of a numerical method of solving Boltzmann’s equation for modeling the behavior of radionuclides (Kr, Xe, Rb, Sr, Cs, Ba) in the cavity of the interelectrode gap of a multielement electricity generating channel. Modeling methods were developed and software system implemented in the course of this work. Calculations were performed for two structural arrangements: with uni- and bilateral extraction of radionuclides into the vacuum-cesium system. Data on the pressure and flow distributions were obtained. The krypton and xenon pressure near the collector as functions of their flow and the cesium gas pressure at the gap exits were determined computationally.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2016

Sandia fracture challenge 2: Sandia California's modeling approach

The second Sandia Fracture Challenge illustrates that predicting the ductile fracture of Ti-6Al-4V subjected to moderate and elevated rates of loading requires thermomechanical coupling, elasto-thermo-poro-viscoplastic constitutive models with the physics of anisotropy and regularized numerical methods for crack initiation and propagation. We detail our initial approach with an emphasis on iterative calibration and systematically increasing complexity to accommodate anisotropy in the context of an isotropic material model. Blind predictions illustrate strengths and weaknesses of our initial approach. We then revisit our findings to illustrate the importance of including anisotropy in the failure process. Furthermore, mesh-independent solutions of continuum damage models having both isotropic and anisotropic yields surfaces are obtained through nonlocality and localization elements.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Fermions in Two Dimensions: Scattering and Many-Body Properties

Ultracold atomic Fermi gases in two dimensions (2D) are an increasingly popular topic of research. The interaction strength between spin-up and spin-down particles in two-component Fermi gases can be tuned in experiments, allowing for a strongly interacting regime where the gas properties are yet to be fully understood. We have probed this regime for 2D Fermi gases by performing T = 0 ab initio diffusion Monte Carlo calculations. The many-body dynamics are largely dependent on the two-body interactions; therefore, we start with an in-depth look at scattering theory in 2D. We show the partial-wave expansion and its relation to the scattering length and effective range. Then, we discuss our numerical methods for determining these scattering parameters. Here, we close out this discussion by illustrating the details of bound states in 2D. Transitioning to the many-body system, we also use variationally optimized wave functions to calculate ground-state properties of the gas over a range of interaction strengths. We show results for the energy per particle and parametrize an equation of state. We then proceed to determine the chemical potential for the strongly interacting gas.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Strain Field Analysis in Electronic Components by ESPI: Bad Thermal Contact and Damage Evaluation

Designing electronic components and devices is a complex task that must include also considerations about thermo-mechanical behavior of the component itself. This aspect, in fact, can affect the reliability of the device as a consequence of the thermal stresses that are introduced in the component while it is working. Stress concentration can arise, in fact, due to different thermal expansion coefficients of the materials involved. Determining the thermo-mechanical response by numerically methods is not simple and not always possible so that experimental methods are preferable. In this paper, specifically, the electronic speckle pattern interferometry technique was adopted to measure full-field strain of the component during exercise. By using an optical technique it is possible to get information about the behavior of the component without being in contact with the component itself so that the dissipation coefficient is not altered. The system was able to detect gradients of deformation in the component itself connected to the different distribution of internal current in the component. Moreover, it was observed, that the presence of bad thermal contact between the package and the heat sink can be revealed because it results in a different thermos-mechanical behavior of the component. Finally, a critical defect was introduced in the component and it was observed the way the presence of damage affects the in-plane displacement of the component. Very huge differences were observed leading to the consideration that this kind of approach could be conveniently adopted also as a damage detection tool.

42 ENGINEERING↗
DOE OSTI2018

Thickness-Dependent Current Density and Flux Distribution in Thin Current-Carrying Superconducting Films Exposed to a Magnetic Field

With the rapid development of high-temperature superconductors, it is of great significance to get the precise current density and flux distribution within thin high-temperature superconducting films subjected to a transport current and an applied magnetic field. The transport current distribution and flux density in thin high-temperature superconducting films are calculated by a numerical method based on the Kim model and exponential model in this paper. The influences of transport current, applied magnetic field, width, and thickness of a superconducting film on the current distribution are discussed. The results reveal that the thickness has a significant effect on the critical current density of superconducting films.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Interaction of Harmonic Longitudinal Shear Waves with V-Shaped Inclusions

We solve the problem of determination of the stressed state in the vicinity of a tunnel rigid inclusion whose cross section consists of two segments originating out from the same point. The inclusion is placed in an infinite elastic medium with propagating plane harmonic longitudinal shear waves. The problem is reduced to the solution of a system of two singular integral equations with fixed singularities. For the approximate solution of this system, we apply a numerical method that takes into account the true asymptotics of the unknown functions and is based on the use of special quadrature formulas for the evaluation of singular integrals.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Design of all-optical XOR and XNOR logic gates based on Fano resonance in plasmonic ring resonators

In this paper, compact all-optical XOR and XNOR gates based on different configurations of metal–insulator–metal plasmonic ring resonators (PRRs) have been proposed. Square and octagon-shaped rings have been used in this case. The logic gates have been simulated using the two-dimensional finite-difference time-domain numerical method, with a conventional perfectly matched layer as the absorbing boundary condition of the area under simulation. The phenomenon of Fano resonance is employed by the proposed gates to excite ON/OFF states. As a result, a large value of contrast ratio (C.R.) is obtained. The results show that the values of C.R. of XNOR and XOR gates for square-shaped PRR are 22.66 and 22.9 dB, respectively and the values of C.R. of XNOR and XOR gates for octagon-shaped PRR are 23.01 and 23.52 dB, respectively. Also, it is shown that the octagon-shaped gates have higher transmission ratios than other proposed configurations. The proposed optical logic gates can be used as key components in optical communications and for designing compact plasmonic devices.

42 ENGINEERING↗
DOE OSTI2019

Analyzing honeycomb photonic crystal waveguides by Dirichlet-to-Neumann maps

As an analog to graphene, honeycomb photonic crystals (PhCs) have attracted a great deal of interest in recent years. The additional degrees of freedom in a honeycomb PhC are useful in designing and optimizing waveguides, can be used to control the symmetry of the system and to realize novel optical devices. In this paper, we present an efficient numerical method for analyzing two-dimensional honeycomb PhC waveguides. Our method is based on a special supercell that covers one period of the honeycomb PhC waveguide and the so-called Dirichlet-to-Neumann map of the supercell. The method gives rise to a linear eigenvalue value problem with relatively small matrices, and it is validated by comparing with previous works. As an application of the method, we calculate edge modes for a few different PhCs. For a honeycomb PhC with dielectric rods in air, we found an edge mode at the zigzag edge, without using special modified rods at the edge. For a honeycomb PhC with air holes, we found a new edge mode when the PhC is bounded by a perfectly electric conductor.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Nanoscale plasmonic detector of wave intensity difference and uni-directional waveguide

In this paper, a new nanoscale plasmonic differential detector of wave intensities is proposed. This structure can be useful for calibration applications and for finding out if a wave intensity is high enough to pass a threshold level or not. Next, a nanoscale plasmonic uni-directional waveguide is presented. The distinction of this waveguide with similar cases is that the one-way wave-guiding is provided by adding a properly located controlling source, without employing any non-reciprocal material. This feature can be used in many applications that only uni-directional wave propagation is desirable. Both structures have been designed in metal–insulator–metal structure including cavity resonators. The operations are based on the coupling of the wave from input waveguides to the cavity resonators and decoupling to output waveguides in the cavity resonance wavelengths. The simulation results of performance of the structures, as elements applicable to photonic integrated circuits, have been obtained using the numerical method of the finite difference time domain.

77 NANOSCIENCE AND NANOTECHNOLOGY↗