Search NASA⌕ Search

SEARCH · Search NASA

Results for “generalized finite differences”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A multi-dimensional Child–Langmuir law for any diode geometry

While prior theoretical studies of multi-dimensional space-charge limited current (SCLC) assumed emission from a small patch on infinite electrodes, none have considered emission from an entire finite electrode. In this paper, we apply variational calculus (VC) and conformal mapping, which have previously been used to derive analytic solutions for SCLC density (SCLCD) for nonplanar one-dimensional geometries, to obtain mathematical relationships for any multi-dimensional macroscopic diode with finite cathode and anode. We first derive a universal mathematical relationship between space-charge limited potential and vacuum potential for any diode and apply this technique to determine SCLCD for an eccentric spherical diode. We then apply VC and the Schwartz–Christoffel transformation to derive an exact equation for SCLCD in a general two-dimensional planar geometry with emission from a finite emitter. Particle-in-cell simulations using VSim agreed within 4%–13% for a range of ratios of emitter width to gap distance using the thinnest electrodes practical for the memory constraints of our hardware, with the difference partially attributed to the theory's assumption of infinitesimally thin electrodes. After generalizing this approach to determine SCLCD for any orthogonal diode as a function of only the vacuum capacitance and vacuum potential, we derive an analytical formulation of the three-dimensional Child–Langmuir law for finite parallel rectangular and disk geometries. These results demonstrate the utility for calculating SCLCD for any diode geometry using vacuum capacitance and vacuum potential, which are readily obtainable for many diode geometries, to guide experiment and simulation development.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Doubly self-consistent field theory of grafted polymers under simple shear in steady state

We present a generalization of the numerical self-consistent mean-field theory of polymers to the case of grafted polymers under simple shear. The general theoretical framework is presented, and then applied to three different chain models: rods, Gaussian chains, and finitely extensible nonlinear elastic (FENE) chains. The approach is self-consistent at two levels. First, for any flow field, the polymer density profile and effective potential are calculated self-consistently in a manner similar to the usual self-consistent field theory of polymers, except that the calculation is inherently two-dimensional even for a laterally homogeneous system. Second, through the use of a modified Brinkman equation, the flow field and the polymer profile are made self-consistent with respect to each other. For all chain models, we find that reasonable levels of shear cause the chains to tilt, but it has very little effect on the overall thickness of the polymer layer, causing a small decrease for rods, and an increase of no more than a few percent for the Gaussian and FENE chains. Using the FENE model, we also probe the individual bond lengths, bond correlations, and bond angles along the chains, the effects of the shear on them, and the solvent and bonded stress profiles. We find that the approximations needed within the theory for the Brinkman equation affect the bonded stress, but none of the other quantities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2011

Frequency-Selective Near-Field Radiative Heat Transfer between Photonic Crystal Slabs: A Computational Approach for Arbitrary Geometries and Materials

Here, we demonstrate the possibility of achieving enhanced frequency-selective near-field radiative heat transfer between patterned (photonic-crystal) slabs at designable frequencies and separations, exploiting a general numerical approach for computing heat transfer in arbitrary geometries and materials based on the finite-difference time-domain method. Our simulations reveal a tradeoff between selectivity and near-field enhancement as the slab-slab separation decreases, with the patterned heat transfer eventually reducing to the unpatterned result multiplied by a fill factor (described by a standard proximity approximation). We also find that heat transfer can be further enhanced at selective frequencies when the slabs are brought into a glide-symmetric configuration, a consequence of the degeneracies associated with the nonsymmorphic symmetry group.

36 MATERIALS SCIENCE↗
DOE OSTI2015

BiGlobal Secondary Instabilities of Roughness-Induced Transient Growth

Roughness-induced transient disturbances distort laminar boundary layers and can lead to premature transition if the roughness is of sufficient height and shape. If transition does occur, it does so through a high-frequency "secondary" instability of the roughness wake. A BiGlobal secondaryinstability code has been developed to solve the stability problem downstream of roughness. The code implements multiple solvers for incompressible/compressible, inviscid/viscous, and temporal/spatial stability. High-order finite differences, sparse matrix linear algebra, LU decomposition, and Arnoldi iterations are used to solve the resulting generalized eigenvalue problem. The code has been used to evaluate the stability properties of experimental roughness-induced transient growth on a flat-plate at low-speed and a slightly blunted cone at M∞= 6. The resulting mode shapes, growth rates, and spatial N-factors show excellent agreement with experimental data. The viscous mode shapes show the best agreement with the experimental disturbance profiles. Yet, execution time can be reduced by an order-of-magnitude by using the inviscid solvers with only modest differences in growth rates.

BiGlobal↗
DOE OSTI2014

Calculation of the extinction cross section and lifetime of a gold nanoparticle using FDTD simulations

The electromagnetic theory of light explains the behavior of light in most of the domains quite accurately. The problem arises when the exact solution of the Maxwell's equation is not present, in case of objects with arbitrary geometry. To find the extinction cross-section and lifetime of the gold nanoparticle, the software FDTD solutions 8.6 by Lumerical is employed. The extinction cross-sections and lifetimes of Gold nanospheres of different sizes and arrangements are studied using pulse lengths of the order of femtoseconds. The decay constant and other properties are compared. Further, the lifetimes are calculated using frequency and time domain calculations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Diagonal ordering operation technique applied to Morse oscillator

We generalize the technique called as the integration within a normally ordered product (IWOP) of operators referring to the creation and annihilation operators of the harmonic oscillator coherent states to a new operatorial approach, i.e. the diagonal ordering operation technique (DOOT) about the calculations connected with the normally ordered product of generalized creation and annihilation operators that generate the generalized hypergeometric coherent states. We apply this technique to the coherent states of the Morse oscillator including the mixed (thermal) state case and get the well-known results achieved by other methods in the corresponding coherent state representation. Also, in the last section we construct the coherent states for the continuous dynamics of the Morse oscillator by using two new methods: the discrete–continuous limit, respectively by solving a finite difference equation. Finally, we construct the coherent states corresponding to the whole Morse spectrum (discrete plus continuous) and demonstrate their properties according the Klauder’s prescriptions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Performance and Self-Consistency of the Generalized Dielectric Dependent Hybrid Functional

Here, we analyze the performance of the recently proposed screened exchange constant functional (SX) on the GW100 test set, and we discuss results obtained at different levels of self-consistency. The SX functional is a generalization of dielectric dependent hybrid functionals to finite systems; it is nonempirical and depends on the average screening of the exchange interaction. We compare results for ionization potentials obtained with SX to those of CCSD(T) calculations and experiments, and we find excellent agreement, on par with recent state of the art methods based on many body perturbation theory. Applying SX perturbatively to correct PBE eigenvalues yields improved results in most cases, except for ionic molecules, for which wave function self-consistency is instead crucial. Calculations where wave functions and the screened exchange constant (α SX ) are determined self-consistently, and those where α SX is fixed to the value determined within PBE, yield results of comparable accuracy. Perturbative G 0 W 0 corrections of eigenvalues obtained with self-consistent αSX are small on average, for all molecules in the GW100 test set.

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

An Electrical Resistance Diagnostic for Conductivity Monitoring in Laser Powder Bed Fusion

With the growing interest in metal additive manufacturing using laser powder bed fusion (LPBF), there is a need for advanced in-situ nondestructive evaluation (NDE) methods that can dynamically monitor manufacturing process-related variations, that can be used as a feedback mechanism to further improve the manufacturing process, leading to parts with improved microstructural properties and mechanical properties. Current NDE techniques either lack sensitivity beyond build layer, are costly or time-consuming, or are not compatible for in-situ integration. In this research, we develop an electrical resistance diagnostic for in-situ monitoring of powder fused regions during laser powder bed fusion printing. The technique relies on injecting current into the build plate and detecting voltage differences from conductive variations during printing using a simple, cheap four-point electrode array directly connected to the build plate. A computational model will be utilized to determine sensitivities of the approach, and preliminary experiments will be performed during the printing process to test the overall approach.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Single-cone finite difference scheme for the (2+1)D Dirac von Neumann equation

An explicit finite difference scheme is presented for the von Neumann equation for (2+1)D Dirac fermions. It is founded upon a staggered space–time grid which ensures a single-cone energy dispersion and performs the time-derivative in one sweep using a three-step leap-frog procedure. It enables a space–time-resolved numerical treatment of the mixed-state dynamics of Dirac fermions within the effective single-particle density matrix formalism. Energy-momentum dispersion, stability and convergence properties are derived. Elementary numerical tests to demonstrate stability properties use parameters which pertain to topological insulator surface states. A method for the simulation of charge injection from an electric contact is presented and tested numerically. Potential extensions of the scheme to a Dirac–Lindblad equation, real-space–time Green's function formulations, and higher-order finite-difference schemes are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

Transient analysis of printed lines using finite-difference time-domain method

Comprehensive studies of ultra-wideband pulses and electromagnetic coupling on printed coupled lines have been performed using full-wave 3D finite-difference time-domain analysis. Effects of unequal phase velocities of coupled modes, coupling between line traces, and the frequency dispersion on the waveform fidelity and crosstalk have been investigated in detail. To discriminate the contributions of different mechanisms into pulse evolution, single and coupled microstrip lines without (ϵ r = 1) and with (ϵ r > 1) dielectric substrates have been examined. To consistently compare the performance of the coupled lines with substrates of different permittivities and transients of different characteristic times, a generic metric similar to the electrical wavelength has been introduced. The features of pulse propagation on coupled lines with layered and pedestal substrates and on the irregular traces have been explored. Finally, physical interpretations of the simulation results are discussed in the paper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

Space-charge-based electrostatic plasma confinement involving relaxed plasma species

A numerical study is reported on the equilibrium properties of a surface-emitted or edge-confined non-drifting plasma. A self-consistent finite-differences evaluation of the electrostatic potential is carried out for a non-neutral plasma that follows a Boltzmann density distribution. The non-neutral plasma generates an electrostatic potential that has an extremum at the geometric center. Poisson's equation is solved for different ratios of the non-neutral plasma size to the edge Debye length. The profiles of the electrostatic potential and the plasma density are presented for different values of that ratio. A second plasma species is then introduced for two-plasma-species confinement studies, with one species confined by the space charge of the other, while each species follows a Boltzmann density distribution. An equilibrium in which a neutral region forms is found. An equilibrium is also found in which the two species have equal temperatures and charge states.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Discontinuous Skeletal Gradient Discretisation methods on polytopal meshes

Here, in this work we develop arbitrary-order Discontinuous Skeletal Gradient Discretisations (DSGD) on general polytopal meshes. Discontinuous Skeletal refers to the fact that the globally coupled unknowns are broken polynomials on the mesh skeleton. The key ingredient is a high-order gradient reconstruction composed of two terms: (i) a consistent contribution obtained mimicking an integration by parts formula inside each element and (ii) a stabilising term for which sufficient design conditions are provided. An example of stabilisation that satisfies the design conditions is proposed based on a local lifting of high-order residuals on a Raviart–Thomas–Nédélec subspace. We prove that the novel DSGDs satisfy coercivity, consistency, limit-conformity, and compactness requirements that ensure convergence for a variety of elliptic and parabolic problems. Lastly, links with Hybrid High-Order, non-conforming Mimetic Finite Difference and non-conforming Virtual Element methods are also studied. Numerical examples complete the exposition.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Communication: Partial linearized density matrix dynamics for dissipative, non-adiabatic quantum evolution

An approach for treating dissipative, non-adiabatic quantum dynamics in general model systems at finite temperature based on linearizing the density matrix evolution in the forward-backward path difference for the environment degrees of freedom is presented. We demonstrate that the approach can capture both short time coherent quantum dynamics and long time thermal equilibration in an application to excitation energy transfer in a model photosynthetic light harvesting complex. Results are also presented for some nonadiabatic scattering models which indicate that, even though the method is based on a ''mean trajectory'' like scheme, it can accurately capture electronic population branching through multiple avoided crossing regions and that the approach offers a robust and reliable way to treat quantum dynamical phenomena in a wide range of condensed phase applications.

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

Magnetic reconnection in the presence of externally driven and self-generated turbulence

Magnetic reconnection is an important process that violates flux freezing and induces change of magnetic field topology in conducting fluids and, as a consequence, converts magnetic field energy into particle energy. It is thought to be operative in laboratory, heliophysical, and astrophysical plasmas. These environments exhibit wide variations in collisionality, ranging from collisionless in the Earth's magnetosphere to highly collisional in molecular clouds. A common feature among these plasmas is, however, the presence of turbulence. We review the present understanding of the effects of turbulence on the reconnection rate, discussing both how strong pre-existing turbulence modifies Sweet-Parker reconnection and how turbulence may develop as a result of reconnection itself. In steady state, reconnection rate is proportional to the aspect ratio of the diffusion region. Thus, two general MHD classes of models for fast reconnection have been proposed, differing on whether they keep the aspect ratio finite by increasing the width due to turbulent broadening or shortening the length of the diffusion layer due to plasmoid instability. One of the consequences of the plasmoid instability model is the possibility that the current sheet thins down to collisionless scales where kinetic effects become dominant. As a result, kinetic effects may be of importance for many astrophysical applications which were considered to be in the realm of MHD. Whether pre-existing turbulence can significantly modify the transition to the kinetic regime is not currently known. Although most studies of turbulent reconnection have been based on MHD, recent advances in kinetic simulations are enabling 3D studies of turbulence and reconnection in the collisionless regime. A summary of these recent works, highlighting similarities and differences with the MHD models of turbulent reconnection, as well as comparison with in situ observations in the magnetosphere and in the solar wind, are presented. The paper concludes with a list of important open questions and suggestions for future work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

An efficient method to estimate the probability density of seismic Green's functions

We present a computationally efficient method to approximate the probability distribution of seismic Green's functions given the uncertainty of an Earth model. The method is based on the Karhunen-Loève (KL) theorem and an approximation of the Green's function (or seismogram) covariance. Using Monte Carlo (MC) simulations as a control case, we demonstrate that our KL-based method can accurately reproduce a probability distribution of seismograms that results from an uncertain Earth model for a MC-derived seismogram covariance. We then describe a method to estimate the covariance of the seismograms resulting from those Earth models that is not based on MC simulations. We use the estimated Green's function covariance in conjunction with our KL-based method to produce a Green's function probability distribution, and compare that distribution to a Green's function probability distribution produced using a MC finite difference method. We find that the Green's function probability distribution approximated using our KL-based method generally mimics that produced using the MC simulations, especially for direct-arriving body waves. However the accuracy of the KL-based method generally decreases for later times in the simulated Green's function distribution.

58 GEOSCIENCES↗
DOE OSTI2012

The COMET method in 3-D hexagonal geometry

The hybrid stochastic-deterministic coarse mesh radiation transport (COMET) method developed at Georgia Tech now solves reactor core problems in 3-D hexagonal geometry. In this paper, the method is used to solve three preliminary test problems designed to challenge the method with steep flux gradients, high leakage, and strong asymmetry and heterogeneity in the core. The test problems are composed of blocks taken from a high temperature test reactor benchmark problem. As the method is still in development, these problems and their results are strictly preliminary. Results are compared to whole core Monte Carlo reference solutions in order to verify the method. Relative errors are on the order of 50 pcm in core eigenvalue, and mean relative error in pin fission density calculations is less than 1% in these difficult test cores. The method requires the one-time pre-computation of a response expansion coefficient library, which may be compiled in a comparable amount of time to a single whole core Monte Carlo calculation. After the library has been computed, COMET may solve any number of core configurations on the order of an hour, representing a significant gain in efficiency over other methods for whole core transport calculations. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2013

Coupling of kinetic Monte Carlo simulations of surface reactions to transport in a fluid for heterogeneous catalytic reactor modeling

We have developed a method to couple kinetic Monte Carlo simulations of surface reactions at a molecular scale to transport equations at a macroscopic scale. This method is applicable to steady state reactors. We use a finite difference upwinding scheme and a gap-tooth scheme to efficiently use a limited amount of kinetic Monte Carlo simulations. In general the stochastic kinetic Monte Carlo results do not obey mass conservation so that unphysical accumulation of mass could occur in the reactor. We have developed a method to perform mass balance corrections that is based on a stoichiometry matrix and a least-squares problem that is reduced to a non-singular set of linear equations that is applicable to any surface catalyzed reaction. The implementation of these methods is validated by comparing numerical results of a reactor simulation with a unimolecular reaction to an analytical solution. Furthermore, the method is applied to two reaction mechanisms. The first is the ZGB model for CO oxidation in which inevitable poisoning of the catalyst limits the performance of the reactor. The second is a model for the oxidation of NO on a Pt(111) surface, which becomes active due to lateral interaction at high coverages of oxygen. This reaction model is based on ab initio density functional theory calculations from literature.

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

A projection gradient method for computing ground state of spin-2 Bose–Einstein condensates

In this paper, a projection gradient method is presented for computing ground state of spin-2 Bose–Einstein condensates (BEC). We first propose the general projection gradient method for solving energy functional minimization problem under multiple constraints, in which the energy functional takes real functions as independent variables. We next extend the method to solve a similar problem, where the energy functional now takes complex functions as independent variables. We finally employ the method into finding the ground state of spin-2 BEC. The key of our method is: by constructing continuous gradient flows (CGFs), the ground state of spin-2 BEC can be computed as the steady state solution of such CGFs. We discretized the CGFs by a conservative finite difference method along with a proper way to deal with the nonlinear terms. We show that the numerical discretization is normalization and magnetization conservative and energy diminishing. Numerical results of the ground state and their energy of spin-2 BEC are reported to demonstrate the effectiveness of the numerical method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗