Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Finite momentum Cooper pairing in three-dimensional topological insulator Josephson junctions

Abstract Unconventional superconductivity arising from the interplay between strong spin–orbit coupling and magnetism is an intensive area of research. One form of unconventional superconductivity arises when Cooper pairs subjected to a magnetic exchange coupling acquire a finite momentum. Here, we report on a signature of finite momentum Cooper pairing in the three-dimensional topological insulator Bi 2 Se 3 . We apply in-plane and out-of-plane magnetic fields to proximity-coupled Bi 2 Se 3 and find that the in-plane field creates a spatially oscillating superconducting order parameter in the junction as evidenced by the emergence of an anomalous Fraunhofer pattern. We describe how the anomalous Fraunhofer patterns evolve for different device parameters, and we use this to understand the microscopic origin of the oscillating order parameter. The agreement between the experimental data and simulations shows that the finite momentum pairing originates from the coexistence of the Zeeman effect and Aharonov–Bohm flux.

36 MATERIALS SCIENCE↗
DOE OSTI2020

PRISMS-PF: A general framework for phase-field modeling with a matrix-free finite element method

Abstract A new phase-field modeling framework with an emphasis on performance, flexibility, and ease of use is presented. Foremost among the strategies employed to fulfill these objectives are the use of a matrix-free finite element method and a modular, application-centric code structure. This approach is implemented in the new open-source PRISMS-PF framework. Its performance is enabled by the combination of a matrix-free variant of the finite element method with adaptive mesh refinement, explicit time integration, and multilevel parallelism. Benchmark testing with a particle growth problem shows PRISMS-PF with adaptive mesh refinement and higher-order elements to be up to 12 times faster than a finite difference code employing a second-order-accurate spatial discretization and first-order-accurate explicit time integration. Furthermore, for a two-dimensional solidification benchmark problem, the performance of PRISMS-PF meets or exceeds that of phase-field frameworks that focus on implicit/semi-implicit time stepping, even though the benchmark problem’s small computational size reduces the scalability advantage of explicit time-integration schemes. PRISMS-PF supports an arbitrary number of coupled governing equations. The code structure simplifies the modification of these governing equations by separating their definition from the implementation of the numerical methods used to solve them. As part of its modular design, the framework includes functionality for nucleation and polycrystalline systems available in any application to further broaden the phenomena that can be used to study. The versatility of this approach is demonstrated with examples from several common types of phase-field simulations, including coarsening subsequent to spinodal decomposition, solidification, precipitation, grain growth, and corrosion.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Domain decomposition and multiscale mortar mixed finite element methods for linear elasticity with weak stress symmetry

Two non-overlapping domain decomposition methods are presented for the mixed finite element formulation of linear elasticity with weakly enforced stress symmetry. The methods utilize either displacement or normal stress Lagrange multiplier to impose interface continuity of normal stress or displacement, respectively. By eliminating the interior subdomain variables, the global problem is reduced to an interface problem, which is then solved by an iterative procedure. The condition number of the resulting algebraic interface problem is analyzed for both methods. A multiscale mortar mixed finite element method for the problem of interest on non-matching multiblock grids is also studied. It uses a coarse scale mortar finite element space on the non-matching interfaces to approximate the trace of the displacement and impose weakly the continuity of normal stress. A priori error analysis is performed. It is shown that, with appropriate choice of the mortar space, optimal convergence on the fine scale is obtained for the stress, displacement, and rotation, as well as some superconvergence for the displacement. Computational results are presented in confirmation of the theory of all proposed methods.

Mathematics↗
DOE OSTI2011

Wigner Crystallization in Spin-Polarized Coupled Electron Quantum Layers: Finite Width Effects

We study the effect of finite width on Wigner crystallization in a spin-polarized coupled electron layers. Correlation among carriers are treated beyond the static mean-field theories by using the quantum version of Singwi, Tosi, Land and Sjoelander (qSTLS) theory over a wide range of layer parameters viz. density parameter r{sub sl} and interlayer spacing d. Interestingly we find that the inclusion of finite width of layers lowered the critical density for the onset of Wigner crystal (WC) as compared to the similar recent studies of spin-polarized coupled electron layers without finite width effects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2012

Study of the deconfinement phase transition in a finite volume with massive particles: Hydrodynamics of the system near the transition

We study the finite size and finite mass effects for the thermal deconfinement phase transition in Quantum Chromodynamics (QCD), using a simple model of coexistence of hadronic (H) gas and quark-gluon plasma (QGP) phases in a finite volume. We consider the equations of state of the two phases with the QGP containing two massless u and d quarks and massive s quarks, and a hadronic gas of massive pions, and we probe the system near the transition. For this, we examine the behavior of the most important hydrodynamical quantities describing the system, at a vanishing chemical potential ({mu}= 0), with temperature and energy density.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2012

Quantum state discrimination bounds for finite sample size

In the problem of quantum state discrimination, one has to determine by measurements the state of a quantum system, based on the a priori side information that the true state is one of the two given and completely known states, {rho} or {sigma}. In general, it is not possible to decide the identity of the true state with certainty, and the optimal measurement strategy depends on whether the two possible errors (mistaking {rho} for {sigma}, or the other way around) are treated as of equal importance or not. Results on the quantum Chernoff and Hoeffding bounds and the quantum Stein's lemma show that, if several copies of the system are available then the optimal error probabilities decay exponentially in the number of copies, and the decay rate is given by a certain statistical distance between {rho} and {sigma} (the Chernoff distance, the Hoeffding distances, and the relative entropy, respectively). While these results provide a complete solution to the asymptotic problem, they are not completely satisfying from a practical point of view. Indeed, in realistic scenarios one has access only to finitely many copies of a system, and therefore it is desirable to have bounds on the error probabilities for finite sample size. In this paper we provide finite-size bounds on the so-called Stein errors, the Chernoff errors, the Hoeffding errors, and the mixed error probabilities related to the Chernoff and the Hoeffding errors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

Exact analysis of particle dynamics in combined field of finite duration laser pulse and static axial magnetic field

Dynamics of a charged particle is studied in the field of a relativistically intense linearly polarized finite duration laser pulse in the presence of a static axial magnetic field. For a finite duration laser pulse whose temporal shape is defined by Gaussian profile, exact analytical expressions are derived for the particle trajectory, momentum, and energy as function of laser phase. From the solutions, it is shown that, unlike for the monochromatic plane wave case, resonant phase locking time between the particle and laser pulse is finite. The net energy transferred to the particle does not increase monotonically but tends to saturate. It is further shown that appropriate tuning of cyclotron frequency of the particle with the characteristic frequency in the pulse spectrum can lead to the generation of accelerated particles with variable energies in MeV-TeV range.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Viscous potential flow analysis of electrified miscible finitely conducting fluid through porous media

In this work, a viscous potential flow analysis is used to investigate capillary surface waves between two horizontal finite fluid layers. The two layers have finite conductivities and admit mass and heat transfer. A general dispersion relation is derived. The presence of finite conductivities together with the dielectric permeabilities makes the horizontal electric field play a dual role in the stability criterion. The phenomenon of negative viscosity is observed. A new growth rate parameter, depending on the kinematical viscosity of the lower fluid layer, is found and has a stabilizing effect on the unstable modes. The growth rates and neutral stability curve are given and applied to air-water interface. The effects of various parameters are discussed for the Kelvin-Helmholtz and the Rayleigh-Taylor instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Finite beta effects on low- and high-frequency magnetosonic waves in a two-ion-species plasma

A magnetosonic wave propagating perpendicular to a magnetic field in a two-ion-species plasma has two branches, high-frequency and low-frequency modes. The finite beta effects on these modes are analyzed theoretically on the basis of the three-fluid model with finite ion and electron pressures. First, it is shown that the Korteweg-de Vries (KdV) equation for the low-frequency mode is valid for amplitudes ε<ε{sub max}, where the upper limit of the amplitude ε{sub max} is given as a function of β (β is the ratio of the kinetic and magnetic energy densities), the density ratio, and the cyclotron frequency ratio of two ion species. Next, the linear dispersion relation and KdV equation for the high-frequency mode are derived, including β as a factor. In addition, the theory for heavy ion acceleration by the high-frequency mode pulse and the pulse damping due to this energy transfer in a finite beta plasma are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Finite Mach number spherical shock wave, application to shock ignition

A converging and diverging spherical shock wave with a finite initial Mach number M{sub s0} is described by using a perturbative approach over a small parameter M{sub s}{sup −2}. The zeroth order solution is the Guderley's self-similar solution. The first order correction to this solution accounts for the effects of the shock strength. Whereas it was constant in the Guderley's asymptotic solution, the amplification factor of the finite amplitude shock Λ(t)∝dU{sub s}/dR{sub s} now varies in time. The coefficients present in its series form are iteratively calculated so that the solution does not undergo any singular behavior apart from the position of the shock. The analytical form of the corrected solution in the vicinity of singular points provides a better physical understanding of the finite shock Mach number effects. The correction affects mainly the flow density and the pressure after the shock rebound. In application to the shock ignition scheme, it is shown that the ignition criterion is modified by more than 20% if the fuel pressure prior to the final shock is taken into account. A good agreement is obtained with hydrodynamic simulations using a Lagrangian code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Finite-temperature second-order many-body perturbation and Hartree–Fock theories for one-dimensional solids: An application to Peierls and charge-density-wave transitions in conjugated polymers

Finite-temperature extensions of ab initio Gaussian-basis-set spin-restricted Hartree–Fock (HF) and second-order many-body perturbation (MP2) theories are implemented for infinitely extended, periodic, one-dimensional solids and applied to the Peierls and charge-density-wave (CDW) transitions in polyyne and all-trans polyacetylene. The HF theory predicts insulating CDW ground states for both systems in their equidistant structures at low temperatures. In the same structures, they turn metallic at high temperatures. Starting from the “dimerized” low-temperature equilibrium structures, the systems need even higher temperatures to undergo a Peierls transition, which is accompanied by geometric as well as electronic distortions from dimerized to non-dimerized forms. The conventional finite-temperature MP2 theory shows a sign of divergence in any phase at any nonzero temperature and is useless. The renormalized finite-temperature MP2 (MP2R) theory is divergent only near metallic electronic structures, but is well behaved elsewhere. MP2R also predicts CDW and Peierls transitions occurring at two different temperatures. The effect of electron correlation is primarily to lower the Peierls transition temperature.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

Systematic study of finite-size effects in quantum Monte Carlo calculations of real metallic systems

We present a systematic and comprehensive study of finite-size effects in diffusion quantum Monte Carlo calculations of metals. Several previously introduced schemes for correcting finite-size errors are compared for accuracy and efficiency, and practical improvements are introduced. In particular, we test a simple but efficient method of finite-size correction based on an accurate combination of twist averaging and density functional theory. Our diffusion quantum Monte Carlo results for lithium and aluminum, as examples of metallic systems, demonstrate excellent agreement between all of the approaches considered.

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

Mean force on a finite-sized rigid particle, droplet, or bubble in a viscous compressible medium

In this paper, a force formulation to compute the axial acoustic mean second-order force on finite-sized compressible and rigid particles is presented. The flow inside and outside the spherical inclusion is considered viscous and compressible. Other than for volumetric pulsations of the bubble/droplet, the sphericity of the inclusion is maintained (taken to be unity). A far-field derivation approach has been used to compute the force due to standing and traveling waves; and the force is expressed as a multipole expansion (infinite series). In case of a bubble and a rigid particle, there exist three length scales that govern the mean second-order force: mean radius of the spherical inclusion (R 0 ), wavelength of the incoming acoustics (λ), and the momentum diffusion thickness of the ambient fluid (δ o ). While R 0 and λ are arbitrary, we assume the viscous length scale is negligibly small compared to the acoustic wavelength. In case of a droplet, however, the following additional parameters (inside to outside fluid ratios) also play a role: density ratio ($\sim\atop{p}$), viscosity ratio ($\sim\atop{μ}$, and speed of sound ratio ($\sim\atop{c}$). The force expression yields the correct behavior in several limiting cases considered: (i) inviscid bubble and droplets with R 0 /λ « 1, (ii) inviscid bubbles with finite R 0 /λ, and (iii) finite size rigid immovable particles. In general, while the monopole alone is sufficient to capture the force for small bubbles, higher-order terms are found to be important when R 0 /λ ≥ 0.02. In addition to reporting similar behavior for droplets, we study the effect arising from $\sim\atop{p}$, $\sim\atop{μ}$, $\sim\atop{c}$, and δ o /R 0 on the mean second-order force.

42 ENGINEERING↗
DOE OSTI2016

Extending the density functional embedding theory to finite temperature and an efficient iterative method for solving for embedding potentials

A key element in the density functional embedding theory (DFET) is the embedding potential. We discuss two major issues related to the embedding potential: (1) its non-uniqueness and (2) the numerical difficulty for solving for it, especially for the spin-polarized systems. To resolve the first issue, we extend DFET to finite temperature: all quantities, such as the subsystem densities and the total system’s density, are calculated at a finite temperature. This is a physical extension since materials work at finite temperatures. We show that the embedding potential is strictly unique at T > 0. To resolve the second issue, we introduce an efficient iterative embedding potential solver. We discuss how to relax the magnetic moments in subsystems and how to equilibrate the chemical potentials across subsystems. The solver is robust and efficient for several non-trivial examples, in all of which good quality spin-polarized embedding potentials were obtained. We also demonstrate the solver on an extended periodic system: iron body-centered cubic (110) surface, which is related to the modeling of the heterogeneous catalysis involving iron, such as the Fischer-Tropsch and the Haber processes. This work would make it efficient and accurate to perform embedding simulations of some challenging material problems, such as the heterogeneous catalysis and the defects of complicated spin configurations in electronic materials.

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

Insights on finite size effects in ab initio study of CO adsorption and dissociation on Fe 110 surface

Adsorption and dissociation of hydrocarbons on metallic surfaces represent crucial steps on the path to carburization, eventually leading to dusting corrosion. While adsorption of CO molecules on Fe surface is a barrier-less exothermic process, this is not the case for the dissociation of CO into C and O adatoms and the diffusion of C beneath the surface that are found to be associated with large energy barriers. In practice, these barriers can be affected by numerous factors that combine to favour the CO-Fe reaction such as the abundance of CO and other hydrocarbons as well as the presence of structural defects. From a numerical point of view, studying these factors is challenging and a step-by-step approach is necessary to assess, in particular, the influence of the finite box size on the reaction parameters for adsorption and dissociation of CO on metal surfaces. Here, we use density functional theory (DFT) total energy calculations with the climbing-image nudged elastic band method to estimate the adsorption energies and dissociation barriers for different CO coverages with surface supercells of different sizes. We further compute the effect of periodic boundary condition for DFT calculations and find that the contribution from van der Waals interaction in the computation of adsorption parameters is important as they contribute to correcting the finite-size error in small systems. The dissociation process involves carbon insertion into the Fe surface causing a lattice deformation that requires a larger surface system for unrestricted relaxation. We show that, in the larger surface systems associated with dilute CO-coverages, C-insertion is energetically more favourable, leading to a significant decrease in the dissociation barrier. This observation suggests that a large surface system with dilute coverage is necessary for all similar metal-hydrocarbon reactions in order to study their fundamental electronic mechanisms, as an isolated phenomenon, free from finite-size effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Nudged elastic band method for solid-solid transition under finite deformation

Solid-state nudged elastic band (SSNEB) methods can be used for finding solid-solid transition paths when solids are subjected to external stress fields. However, previous SSNEB methods may lead to inaccurate barriers and deviated reaction paths for transitions under stress and finite deformation due to an inaccurate evaluation of the external work contributions in enthalpies. In this paper, a finite deformation nudged elastic band (FD-NEB) method is formulated for finding transition paths of solids under finite deformation. Applications of FD-NEB to a phase transition of silicon from the diamond phase to the β-tin phase under uniaxial compression are presented. Here, the results are compared with those from the generalized solid-state nudged elastic band method.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Real-time density-matrix coupled-cluster approach for closed and open systems at finite temperature

We extend the coupled-cluster method to correlated quantum dynamics of both closed and open systems at finite temperatures using the thermofield formalism. The approach expresses the time-dependent density matrix in an exponential ansatz and describes time-evolution along the Keldysh path contour. A distinct advantage of the approach is exact trace-preservation as a function of time, ensuring conservation of probability and particle number. Furthermore, the method avoids the computation of correlated bra-states, simplifying the computational implementation. We develop the method in a thermal quasiparticle representation, which allows seamless connection to the projection method and diagrammatic techniques of the traditional coupled-cluster formalism. For comparison, we also apply the thermofield framework to the density-matrix renormalization-group method to obtain reference results for closed and open systems at finite temperature. We test the singles and doubles approximation to the density-matrix coupled-cluster method on the correlated electronic dynamics of the single-impurity Anderson model, demonstrating that the new method successfully captures the correlated dynamics of both closed systems at finite temperature and driven-dissipative open systems. This encouraging performance motivates future applications to nonequilibrium quantum many-body dynamics in realistic systems.

Shushkov, Philip (ORCID:0000000266136117)↗
DOE OSTI2021

Conservation laws in coupled cluster dynamics at finite temperature

We extend the finite-temperature Keldysh non-equilibrium coupled cluster theory (Keldysh-CC) [A. F. White and G. K.-L. Chan, J. Chem. Theory Comput. 15, 6137–6253 (2019)] to include a time-dependent orbital basis. When chosen to minimize the action, such a basis restores local and global conservation laws (Ehrenfest’s theorem) for all one-particle properties while remaining energy conserving for time independent Hamiltonians. We present the time-dependent Keldysh orbital-optimized coupled cluster doubles method in analogy with the formalism for zero-temperature dynamics, extended to finite temperatures through the time-dependent action on the Keldysh contour. To demonstrate the conservation property and understand the numerical performance of the method, we apply it to several problems of nonequilibrium finite-temperature dynamics: a 1D Hubbard model with a time-dependent Peierls phase, laser driving of molecular H2, driven dynamics in warm-dense silicon, and transport in the single impurity Anderson model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗