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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Finite field-dependent symmetries in perturbative quantum gravity

In this paper we discuss the absolutely anticommuting nilpotent symmetries for perturbative quantum gravity in general curved spacetime in linear and non-linear gauges. Further, we analyze the finite field-dependent BRST (FFBRST) transformation for perturbative quantum gravity in general curved spacetime. The FFBRST transformation changes the gauge-fixing and ghost parts of the perturbative quantum gravity within functional integration. However, the operation of such symmetry transformation on the generating functional of perturbative quantum gravity does not affect the theory on physical ground. The FFBRST transformation with appropriate choices of finite BRST parameter connects non-linear Curci–Ferrari and Landau gauges of perturbative quantum gravity. The validity of the results is also established at quantum level using Batalin–Vilkovisky (BV) formulation. -- Highlights: •The perturbative quantum gravity is treated as gauge theory. •BRST and anti-BRST transformations are developed in linear and non-linear gauges. •BRST transformation is generalized by making it finite and field dependent. •Connection between linear and non-linear gauges is established. •Using BV formulation the results are established at quantum level also.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Excited-state quantum phase transitions in systems with two degrees of freedom: II. Finite-size effects

This article extends our previous analysis Stránský et al. (2014) of Excited-State Quantum Phase Transitions (ESQPTs) in systems of dimension two. We focus on the oscillatory component of the quantum state density in connection with ESQPT structures accompanying a first-order ground-state transition. It is shown that a separable (integrable) system can develop rather strong finite-size precursors of ESQPT expressed as singularities in the oscillatory component of the state density. The singularities originate in effectively 1-dimensional dynamics and in some cases appear in multiple replicas with increasing excitation energy. Using a specific model example, we demonstrate that these precursors are rather resistant to proliferation of chaotic dynamics. - Highlights: • Oscillatory components of state density and spectral flow studied near ESQPTs. • Enhanced finite-size precursors of ESQPT caused by fully/partly separable dynamics. • These precursors appear due to criticality of a subsystem with lower dimension. • Separability-induced finite-size effects disappear in case of fully chaotic dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

A finite-time exponent for random Ehrenfest gas

We consider the motion of a system of free particles moving on a plane with regular hard polygonal scatterers arranged in a random manner. Calling this the Ehrenfest gas, which is known to have a zero Lyapunov exponent, we propose a finite-time exponent to characterize its dynamics. As the number of sides of the polygon goes to infinity, when polygon tends to a circle, we recover the usual Lyapunov exponent for the Lorentz gas from the exponent proposed here. To obtain this result, we generalize the reflection law of a beam of rays incident on a polygonal scatterer in a way that the formula for the circular scatterer is recovered in the limit of infinite number of vertices. Thus, chaos emerges from pseudochaos in an appropriate limit. - Highlights: • We present a finite-time exponent for particles moving in a plane containing polygonal scatterers. • The exponent found recovers the Lyapunov exponent in the limit of the polygon becoming a circle. • Our findings unify pseudointegrable and chaotic scattering via a generalized collision rule. • Stretch and fold:shuffle and cut :: Lyapunov:finite-time exponent :: fluid:granular mixing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Tomograms for open quantum systems: In(finite) dimensional optical and spin systems

Tomograms are obtained as probability distributions and are used to reconstruct a quantum state from experimentally measured values. We study the evolution of tomograms for different quantum systems, both finite and infinite dimensional. In realistic experimental conditions, quantum states are exposed to the ambient environment and hence subject to effects like decoherence and dissipation, which are dealt with here, consistently, using the formalism of open quantum systems. This is extremely relevant from the perspective of experimental implementation and issues related to state reconstruction in quantum computation and communication. These considerations are also expected to affect the quasiprobability distribution obtained from experimentally generated tomograms and nonclassicality observed from them. -- Highlights: •Tomograms are constructed for open quantum systems. •Finite and infinite dimensional quantum systems are studied. •Finite dimensional systems (phase states, single & two qubit spin states) are studied. •A dissipative harmonic oscillator is considered as an infinite dimensional system. •Both pure dephasing as well as dissipation effects are studied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Quantum dilute droplets of dipolar bosons at finite temperature

We systematically study the properties of dipolar Bose gases with two- and three-body contact interactions at finite temperature using the Hartree–Fock–Bogoliubov–Popov approximation. In uniform case, we obtain an exciting new extension of the seminal Lee–Huang–Yang corrected equation of state that depends explicitly on the thermal fluctuations and on the coupling constant of the three-body interaction. We investigate, on the other hand, the effects of thermal fluctuations on the occurrence and stability of a droplet state in a Bose–Einstein condensate with strong dipole–dipole interactions. We find that at finite temperature, the droplet phase appears as a narrow peak surrounded by a broader thermal halo. We show that the number of particles inside the droplet decays with increasing temperature. - Highlights: • Investigation of a dipolar Bose gas with two- and three-body interactions. • Derivation of a generalized Lee–Huang–Yang corrected equation. • Studying effects of thermal fluctuations on the nucleation of droplets in a dilute dipolar BEC. • Solving the finite-temperature GP equation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Quantitative coherence witness for finite dimensional states

We define the stringent coherence witness as an observable whose mean value vanishes for all incoherent states but nonzero for some coherent states. Such witnesses are proved to exist for any finite-dimension states. Not only is the witness efficient in testing whether a state is coherent, but also its mean value can quantitatively reveal the amount of coherence. For an unknown state, the modulus of the mean value of a normalized witness provides a tight lower bound to the l{sub 1}-norm of coherence. When we have some previous knowledge of a state, the optimal witness which has the maximal mean value is derived. It is proved that for any finite dimension state, the mean value of the optimal witness, which we call the witnessed coherence, equals the l{sub 1}-norm of coherence. In the case that the witness is fixed and the incoherent operations are allowed, the maximal mean value can reach the witnessed coherence if and only if certain relations between the fixed witness and the initial state are satisfied. Our results provide a way to directly measure the coherence in arbitrary finite dimension states and an operational interpretation of the l{sub 1}-norm of coherence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

A new numerically stable sequential algorithm for coupled finite-strain elastoplastic geomechanics and flow

In this work, we propose a new numerically stable sequential implicit method for coupled flow and finite-strain multiplicative elastoplastic geomechanics. We find from stability analysis that the sequential method that solves the flow problem first by fixing the first Piola-Kirchhoff total stress, solving the geomechanics (solid deformation) problem at the next step, is unconditionally stable (contractive and B-stable). In this sense, this method named the fixed first Piola-Kirchhoff stress method is an extension of the fixed stress method in coupled flow and infinitesimal geomechanics. We also study the fixed second Piola-Kirchhoff stress method, comparing it with the fixed first Piola-Kirchhoff stress method, because the constitutive relations are formulated by the second Piola-Kirchhoff total stress, although fixing the second Piola-Kirchhoff total stress field does not provide theoretical unconditional stability. In space discretization, we use the finite element method for the geomechanics problem with the total Lagrangian approach, while employing the finite volume method for the flow problem. Geometrical nonlinearity from the total Lagrangian approach results in full-tensor permeability even if the initial permeability is isotropic. To consider the full-tensor permeability accurately, we employ the multipoint flux approximation in solving the flow problem. In time discretization, the backward Euler method is used. Then, we perform numerical experiments of the two fixed-stress sequential methods with various scenarios, and find superiority of their numerical stability. Specifically, the two fixed-stress sequential methods provide almost identical results for elasticity. For J2 plasticity, the two sequential methods yield slightly different numerical results although the trends of the results are similar. For the Drucker-Prager plasticity, we find that the results from the two methods are almost identical. For all the test cases, the two sequential methods are numerically stable.

58 GEOSCIENCES↗
DOE OSTI2021

Finite-rate entrainment effects on nitrogen oxide (NO x ) emissions in staged combustors

Axially staged combustors offer the potential for achieving low NO x emissions from gas turbines at the elevated temperatures required to achieve at least 65% combined cycle efficiencies. While there is evidence that premixing between the main burner products and secondary stream reduces NO x , it is important to characterize the specific requirements needed to achieve low emissions at gas turbine conditions. This study examines the sensitivity of NO x to finite-rate, large-scale entrainment of the main and secondary streams under the simplifying assumption of infinitely-fast small-scale mixing. Essentially, this allows us to isolate the effects of limited large-scale entrainment rates by using a homogeneous/uniform condition for the entrained and reacting gases. We use a reduced-order reactor network model to examine a generic staged-combustor whose inputs are physical time scales that embody the finite-rate entrainment characteristics. With simulations conducted over a large parameter space at typical operating conditions (25 atm, 650 K air), the results show that, even when small-scale mixing is infinitely fast and the reaction zone is uniform, entrainment significantly affects NO x emissions due to its influence on the equivalence ratio — and thus the temperature and time — at which the entrained mixture ignites. Fuel-air staging is shown to be vital to NO x reduction in most practical cases where entrainment times exceed roughly 1 ms and the secondary stream finishes entraining before the main burner products. A constrained NO x minimization indicates that using the leanest possible main burner and re-routing as much air as possible to the secondary stage is key to low NO x with finite-rate entrainment; for example, less than 10 ppm NO x can be achieved if the secondary stage contains 20% of the total mass flow.

03 NATURAL GAS↗
DOE OSTI2012

Efficient gradient estimation using finite differencing and likelihood ratios for kinetic Monte Carlo simulations

While many optimization and control methods for stochastic processes require gradient information from the process of interest, obtaining gradient information from experiments is prohibitively expensive and time-consuming. As a result, such information is often obtained from stochastic process simulations. Computing gradients efficiently and accurately from stochastic simulations is challenging, especially for simulations involving computationally expensive models with significant inherent noise. In this work, we analyze and characterize the applicability of two gradient estimation methods for kinetic Monte Carlo simulations: finite differencing and likelihood ratio. We developed a systematic method for choosing an optimal perturbation size for finite differencing and discuss, for both methods, important implementation issues such as scaling with respect to the number of elements in the gradient vector. Through a series of numerical experiments, the methods were compared across different time and size regimes to characterize the precision and accuracy associated with each method. We determined that the likelihood ratio method is appropriate for estimating gradients at short (transient) times or for systems with small population sizes, whereas finite differencing is better-suited for gradient estimation at long times (steady state) or for systems with large population sizes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Finite-volume scheme for anisotropic diffusion

In this paper, we apply a special finite-volume scheme, limited to smooth temperature distributions and Cartesian grids, to test the importance of connectivity of the finite volumes. The area of application is nuclear fusion plasma with field line aligned temperature gradients and extreme anisotropy. We apply the scheme to the anisotropic heat-conduction equation, and compare its results with those of existing finite-volume schemes for anisotropic diffusion. Also, we introduce a general model adaptation of the steady diffusion equation for extremely anisotropic diffusion problems with closed field lines.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Higher-order finite-difference formulation of periodic Orbital-free Density Functional Theory

We present a real-space formulation and higher-order finite-difference implementation of periodic Orbital-free Density Functional Theory (OF-DFT). Specifically, utilizing a local reformulation of the electrostatic and kernel terms, we develop a generalized framework for performing OF-DFT simulations with different variants of the electronic kinetic energy. In particular, we propose a self-consistent field (SCF) type fixed-point method for calculations involving linear-response kinetic energy functionals. In this framework, evaluation of both the electronic ground-state and forces on the nuclei are amenable to computations that scale linearly with the number of atoms. We develop a parallel implementation of this formulation using the finite-difference discretization. We demonstrate that higher-order finite-differences can achieve relatively large convergence rates with respect to mesh-size in both the energies and forces. Additionally, we establish that the fixed-point iteration converges rapidly, and that it can be further accelerated using extrapolation techniques like Anderson's mixing. We validate the accuracy of the results by comparing the energies and forces with plane-wave methods for selected examples, including the vacancy formation energy in Aluminum. Overall, the suitability of the proposed formulation for scalable high performance computing makes it an attractive choice for large-scale OF-DFT calculations consisting of thousands of atoms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Treating network junctions in finite volume solution of transient gas flow models

A finite volume scheme for the numerical solution of a non-isothermal non-adiabatic compressible flow model for gas transportation networks on non-flat topography is introduced. Unlike standard Euler equations, the model takes into account wall friction, variable height and heat transfer between the pipe and the environment which are source terms. The case of one single pipe was considered in a previous reference by the authors, , where a finite volume method with upwind discretization of the flux and source terms has been proposed in order to get a well-balanced scheme. The main goal of the present paper is to go a step further by considering a network of pipes. The main issue is the treatment of junctions for which container-like 2D finite volumes are introduced. The couplings between pipes (1D) and containers (2D) are carefully described and the conservation properties are analyzed. Numerical tests including real gas networks are solved showing the performance of the proposed methodology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

A high-order multiscale finite-element method for time-domain acoustic-wave modeling

Accurate and efficient wave equation modeling is vital for many applications in such as acoustics, electromagnetics, and seismology. However, solving the wave equation in large-scale and highly heterogeneous models is usually computationally expensive because the computational cost is directly proportional to the number of grids in the model. We develop a novel high-order multiscale finite-element method to reduce the computational cost of time-domain acoustic-wave equation numerical modeling by solving the wave equation on a coarse mesh based on the multiscale finite-element theory. In contrast to existing multiscale finite-element methods that use only first-order multiscale basis functions, our new method constructs high-order multiscale basis functions from local elliptic problems which are closely related to the Gauss–Lobatto–Legendre quadrature points in a coarse element. Essentially, these basis functions are not only determined by the order of Legendre polynomials, but also by local medium properties, and therefore can effectively convey the fine-scale information to the coarse-scale solution with high-order accuracy. Numerical tests show that our method can significantly reduce the computation time while maintain high accuracy for wave equation modeling in highly heterogeneous media by solving the corresponding discrete system only on the coarse mesh with the new high-order multiscale basis functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

A single crystal plasticity finite element formulation with embedded deformation twins

Deformation twinning is an important plastic deformation mechanism in some polycrystalline metals such as titanium and magnesium. In this paper, we present a novel crystal plasticity finite element framework that accounts for deformation twinning explicitly, in addition to crystallographic slip. Within this computational framework, deformation twins are treated as weak discontinuities embedded within individual finite elements, such that a jump in the velocity gradient field is introduced (via the discretized gradient operator) between the twinned and untwinned crystalline regions, taking into account compatibility and traction continuity conditions at the interface between these two regions. The deformation gradient is multiplicatively split into elastic and plastic parts in the untwinned region, as is customary in finite-deformation crystal plasticity formulations. A different multiplicative decomposition of the deformation gradient into elastic, plastic (slip), and twinning parts is adopted in the twinned region, allowing deformation twinning to be accounted for as an additional mode of plastic deformation. Here, a stochastic model is used to predict twin nucleation at grain boundaries, and the evolution of the length and thickness of the twinned region under subsequent deformation is taken into account.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Extension of the Kohn-Sham formulation of density functional theory to finite temperature

Based on Mermin's extension of the Hohenberg and Kohn theorems to non-zero temperature, the Kohn-Sham formulation of density functional theory (KS-DFT) is generalized to finite temperature. Here, we show that present formulations are inconsistent with Mermin's functional containing expressions, in particular describing the Coulomb energy, that defy derivation and are even in violation of rules of logical inference. More; current methodology is in violation of fundamental laws of both quantum and classical mechanics. Based on this feature, we demonstrate the impossibility of extending the KS formalism to finite temperature through the self-consistent solutions of the single-particle Schrödinger equation of T>0. Guided by the form of Mermin's functional that depends on the eigenstates of a Hamiltonian, determined at T>0 we base our extension of KS-DFT on the determination of the excited states of a non-interacting system at the zero of temperature. The resulting formulation is consistent with that of Mermin constructing the free energy at T>0 in terms of the excited states of a non-interacting Hamiltonian (system) that, within the KS formalism, are described by Slater determinants. To determine the excited states at T=0 use is made of the extension of the Hohenberg and Kohn theorems to excited states presented in previous work applied here to a non-interacting collection of replicas of a non-interacting N-particle system, whose ground state density is taken to match that of K non-interacting replicas of an interacting N-particle system at T>0. The formalism allows for an ever denser population of the excitation spectrum of a Hamiltonian, within the KS approximation. The form of the auxiliary potential, (Kohn-Sham potential), is formally identical to that in the ground state formalism with the contribution of the Coulomb energy provided by the derivative of the Coulomb energy in all excited states taken into account. Once the excited states are determined, the minimum of the free energy within the KS formalism follows immediately in the form of Mermin's functional, but with the exact excited states in that functional represented by Slater determinants obtained through self-consistency conditions at the zero of temperature. Lastly, it is emphasized that, in departure from all existing formulations, no self-consistency conditions are implemented at finite T; as we show, in fact, such formulations are rigorously blocked.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Simulation of crack induced nonlinear elasticity using the combined finite-discrete element method

Numerical simulation of nonlinear elastic wave propagation in solids with cracks is indispensable for decoding the complicated mechanisms associated with the nonlinear ultrasonic techniques in NonDestructive Testing (NDT). Here, we introduce a two-dimensional implementation of the combined finitediscrete element method (FDEM), which merges the finite element method (FEM) and the discrete element method (DEM), to explicitly simulate the crack induced nonlinear elasticity in solids with both horizontal and inclined cracks. In the FDEM model, the solid is discretized into finite elements to capture the wave propagation in the bulk material, and the finite elements along the two sides of the crack also behave as discrete elements to track the normal and tangential interactions between crack surfaces. The simulation results show that for cracked models, nonlinear elasticity is generated only when the excitation amplitude is large enough to trigger the contact between crack surfaces, and the nonlinear behavior is very sensitive to the crack surface contact. The simulations reveal the influence of normal and tangential contact on the nonlinear elasticity generation. Moreover, the results demonstrate the capabilities of FDEM for decoding the causality of nonlinear elasticity in cracked solid and its potential to assist in Non-Destructive Testing (NDT).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Finite Strain Analysis of Shear and Compressional Wave Velocities

Published shear-modulus measurements for a wide variety of materials (Ar, Xe, H 2 , He, NaCl, H 2 O-VII, MgO, stishovite, bridgmanite) show that the Eulerian (spatial) description of energy vs. strain fits both finite- and infinitesimal-strain (e.g., wave velocity) elasticity data under high pressure. The Eulerian (spatial) formulations do so better than the Lagrangian (material) finite-strain description, with differences of 1% to 60% in both P- and S -wave velocities for these materials at the pressures of Earth's mantle. The findings are significant in extending to shear previous findings that compressional (pressure-volume and bulk-modulus) measurements are also best fit using the spatial formulation. Our assessment empirically documents that a self-consistent Eulerian finite-strain equation of state offers a reliable means of describing the thermodynamic and elastic properties of planetary interiors.

58 GEOSCIENCES↗
DOE OSTI2022

Finite Element Modeling of Diffusion in Fractured Porous Media by Using Hierarchical Material Properties

Abstract Fractured porous media challenge modeling approaches due to high computational costs and excessive mesh refinement imposed by the extreme scale variability of fractures and the heterogeneity of the surrounding porous rock. To overcome such difficulties, we utilize the hierarchical finite element method ( Hi ‐FEM) that has been developed previously to simulate the electrical potential distribution in complex geologic environments. The method employs the hierarchical basis functions in classical finite element analysis to enable representation of material properties on each dimensional component of a given 3D unstructured finite element, thereby inherently allowing for interactions at the boundary between fracture and a host rock. In this study, we extend its application to transient fluid flow and heat conduction in the Laplace domain. Time‐domain flow solutions are obtained by numerical inverse Laplace transform. We evaluate the accuracy of the method using different flow models and demonstrate its robustness for large‐scale, rock mass models featuring complex fracture networks. Moreover, for the computation of nodal Darcian velocity fields in fractured porous media where the fractures are represented as 2D features, a new approach that employs the Yeh's Galerkin model for both volume and facet elements is proposed. Results show that Hi ‐FEM can produce accurate flow solutions for fractured porous media without any need of coupling or transfer mechanism while still being computationally economical and numerically robust, even for large‐scale simulations.

Beskardes, G. D.↗