Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite difference”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Average-atom model with Siegert states

In plasmas, electronic states can be well-localized bound states or itinerant free states, or something in between. In self-consistent treatments of plasma electronic structure such as the average-atom model, all states must be accurately resolved in order to achieve a converged numerical solution. Furthermore, this is a challenging numerical and algorithmic problem in large part due to the continuum of free states which is relatively expensive and difficult to resolve accurately. Siegert states are an appealing alternative. They form a complete eigenbasis with a purely discrete spectrum while still being equivalent to a representation in terms of the usual bound states and free states. However, many of their properties are unintuitive, and it is not obvious that they are suitable for self-consistent plasma electronic structure calculations. Here it is demonstrated that Siegert states can be used to accurately solve an average-atom model and offer advantages over the traditional finite-difference approach, including a concrete physical picture of pressure ionization and continuum resonances.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Integrability technique for fluid flow induced deformation of a boundary hair

The deformation of a dense carpet of hair due to Stokes flow in a channel can be described by a nonlinear integrodifferential equation for the shape of a single hair, which possesses several solutions for a given choice of parameters. Although it was posed in a previous study and it bears a resemblance to the pendulum problem from mechanics, this equation has not been analytically solved until now. Despite the presence on an integral with a nonlinear functional dependence on the dependent variable, the system is integrable. We compare the analytically obtained solution to a finite-difference numerical approach, identify the physically realizable solution branch, and briefly study the solution structure through a conserved energylike quantity. Time-dependent fluid-structure interactions are a rich and complex subject to investigate, and we argue that the solution discussed herein can be used as a basis for understanding these systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Mass and momentum transport in the Tilted Rocket Rig experiment

Here, high-fidelity large eddy simulations of the inclined Rayleigh-Taylor instability in the Tilted Rocket Rig experimental configuration [Smeeton and Youngs, AWE Report No. O 35/87 (1987)] are performed using a tenth-order compact finite difference code. These simulations are analyzed for spatial distributions of turbulent kinetic energy, turbulent mass flux velocity, species mass fraction flux, species mass fraction variance, and Favre-averaged Reynolds stresses at two time instances, t = 37 ms and t = 52 ms. Additionally, the vertical distribution of the components of the unclosed budget equations over the center of the domain for these quantities are examined. The dominant terms of these budget equations are further decomposed to examine the principal contributions to these terms in each axis. Notably, the principal contribution to the horizontal turbulent mass flux velocity budget is found to be from a term which is commonly neglected in many Reynolds-averaged Navier-Stokes models.

42 ENGINEERING↗
DOE OSTI2023

Phonon-Induced Localization of Excitons in Molecular Crystals from First Principles

The spatial extent of excitons in molecular systems underpins their photophysics and utility for optoelectronic applications. Phonons are reported to lead to both exciton localization and delocalization. However, a microscopic understanding of phonon-induced (de)localization is lacking, in particular, how localized states form, the role of specific vibrations, and the relative importance of quantum and thermal nuclear fluctuations. Here, we present a first-principles study of these phenomena in solid pentacene, a prototypical molecular crystal, capturing the formation of bound excitons, exciton-phonon coupling to all orders, and phonon anharmonicity, using density functional theory, the ab initio GW-Bethe-Salpeter equation approach, finite-difference, and path integral techniques. In this study, we find that for pentacene zero-point nuclear motion causes uniformly strong localization, with thermal motion providing additional localization only for Wannier-Mott-like excitons. Anharmonic effects drive temperature-dependent localization, and, while such effects prevent the emergence of highly delocalized excitons, we explore the conditions under which these might be realized

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Solving the electronic structure problem for over 100000 atoms in real space

Using a real-space high-order finite-difference approach, we investigate the electronic structure of large spherical silicon nanoclusters. Within Kohn-Sham density functional theory and using pseudopotentials, we report the self-consistent field convergence of a system with over 100000 atoms: a Si 107,641 ⁢H 9,084 nanocluster with a diameter of 16 nm. Our approach uses Chebyshev-filtered subspace iteration to speed up the convergence of the eigenspace, and blockwise Hilbert space-filling curves to speed up sparse matrix-vector multiplications, all of which are implemented in the parsec code. For the largest system, we utilized 2048 nodes (114 688 cores) on the Frontera machine in the Texas Advanced Computing Center. Our quantitative analysis of the electronic structure shows how it gradually approaches its bulk counterpart as a function of nanocluster size. The band gap is enlarged due to quantum confinement in nanoclusters, but decreases as the system size increases, as expected. In conclusion, our work serves as a proof of concept for the capacity of the real-space approach in efficiently parallelizing very large calculations using high-performance computer platforms, which can straightforwardly be replicated in other systems with more than 10 5 atoms.

0-dimensional systems↗
DOE OSTI2025

Systematic study of the validity of the eikonal model including uncertainties

Nuclear reactions at intermediate beam energies are often interpreted using the eikonal model. In the analysis of complex reaction probes, where few-body reaction methods are needed, the eikonal method may be used as an efficient way for describing the fragment-target reaction process. In this work, we perform a systematic study to test the validity of the eikonal approximation for nucleon-nucleus reactions. We also quantify uncertainties due to the nucleon optical potential on reaction observables. We inspect the validity of the eikonal model and its semiclassical correction by comparing it to exact solutions (obtained from solving the optical-model equation with a finite-differences method) for a wide range of reactions. We also study the effect of relativistic corrections, both kinematic and dynamic, by effectively incorporating the relativistic effects at intermediate energies. The uncertainties from a Bayesian global optical potential (KDUQ) are propagated to the observables of interest. Our study includes neutron and proton reactions on 27 Al , 40 Ca , 90 Zr , and 208 Pb , for a wide range of energies 𝐸 lab = 0–400 MeV. We calculate neutron-total cross sections (elastic and reactions) as well as proton-absorption cross sections as a function of beam energy, using the eikonal model, the eikonal model with a semiclassical correction, and the exact solution. Here, we also compute angular distributions for the methods above. Our results show that for the proton-absorption cross section, the eikonal model can be used down to around 60 MeV and the semiclassical correction extends its use to 30 MeV. However, the validity of the eikonal model for the neutron-total cross section only goes down to ≈120 MeV, a range extended to ≈ 50 MeV when using the semiclassical correction. We find the semiclassical correction to the eikonal model to be less effective in describing the angular distributions. The 1⁢𝜎 uncertainty intervals on the observables we studied is less than 5% for most of the energies considered, but increases rapidly for higher energies, namely energies outside the range of KDUQ (𝐸 lab > 200MeV).

Cluster models↗
DOE OSTI2024

Rotordynamic Analysis and Comparative Study of High-Speed Outer Rotor Permanent Magnet Motor Designs

High–power density electric motor designs are a requirement in aerospace and automotive applications. Outer rotor permanent magnet motors can offer high power density but have mechanical challenges such as structural stability and rotodynamic issues. In this work, a rotodynamic study was performed for two outer rotor permanent magnet motor designs. The first design was a cantilever design in which the rotor was suspended at one end, supported by four bearings; in the second design, the rotor was simply supported by two bearings in each end. Two different finite element method–based approaches, solid rotor and beam rotor methods, were used to extract the critical speed.

Barua, Himel↗
DOE OSTI2022

Analysis of First-Order Gratings in Silicon Photonic Waveguides

A simple thin film effective index analysis for first order gratings in Si photonic waveguides is shown to provide highly accurate results for reflected and transmitted power spectrums as long as the waveguide remains single mode and non-radiating. A cover layer can be added to the grating region of a Si photonic waveguide to increase the strength of the grating, modify transition losses from the input waveguide to the grating waveguide region, and/or modify the width of the reflectivity spectrum. For a given grating period, the peak reflection and spectral width of the reflectivity decrease as the duty cycle is decreased or increased from ~50%. For both radiating and multimode structures, the coupling between all modes, power radiated towards the superstrate (upwards), power radiated downwards (substrate) and transmitted power analyzed by Floquet-Bloch, Eigenmode Expansion and Finite Difference Time Domain methods show excellent agreement. Coupling coefficients calculated using analytic formulas are shown to be accurate only for shallow grating depths.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Development of an Optical Library for Coevaporated CdSe x Te 1– x

The conversion efficiency of CdTe solar cells may be improved by bandgap engineering, i.e., changing the bandgap value through the addition of Se in the absorber. The Se alloying enables a short-circuit current density improvement, as it leads to a bandgap energy value decrease. Furthermore, it has been associated with increased minority carrier lifetimes, assuring high open-circuit voltage values. An Se gradient profile control can further optimize the solar cell performance. Thus, an optical model baseline of the CdSe x Te 1–x (CST) compound was developed. Spectroscopic ellipsometry measurements were conducted to accurately extract the optical constants of ten CST layers deposited through coevaporation with x varying from 0 to 1. Using the measured dielectric function spectra from the discrete CST layers with varying x, and considering the composition-induced shift in the critical point energies, an energy-shift model was employed to develop the accurate optical library for the CST compound for any x value to provide data for future modeling and optimization. Furthermore, the library accuracy was validated through optical simulations of the quantum efficiency of a graded CST solar cell using the finite-difference time-domain method by replicating the Se profile in the absorber layer measured through secondary ion mass spectrometry.

14 SOLAR ENERGY↗
DOE OSTI2021

A Semi-Analytical Solution Approach for Solving Constant-Coefficient First-Order Partial Differential Equations

Simulation and control of many dynamic systems involve solving partial differential equations (PDE). This letter proposes a semi-analytical solution (SAS) approach for fast and high-quality solution of first-order PDEs. The region of interest of the studied PDE is divided into a grid, and an SAS is derived for each grid cell in the form of the multivariate polynomials, of which the coefficients are identified using initial value and boundary value conditions. The solutions are solved in a “time-stepping” manner, i.e. within one time step, the coefficients of the SAS are identified and the initial value of the next time step is evaluated. This approach achieves a significantly larger grid cell than the widely used finite difference method, and thus enhances the computational efficiency significantly. Furthermore, the simulation result on the natural gas pipeline model demonstrates the advantages of SAS in accuracy and computational efficiency.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Advanced Quantum Poisson Solver in the NISQ era

The Poisson equation has many applications across the broad areas of science and engineering. Most quantum algorithms for the Poisson solver presented so far, either suffer from lack of accuracy and/or are limited to very small sizes of the problem, and thus have no practical usage. Here we present an advanced quantum algorithm for solving the Poisson equation with high accuracy and dynamically tunable problem size. After converting the Poisson equation to the linear systems through the finite difference method, we adopt the Harrow-Hassidim-Lloyd (HHL) algorithm as the basic framework. Particularly, in this work we present an advanced circuit that ensures the accuracy of the solution by implementing non-truncated eigenvalues through eigenvalue amplification as well as by increasing the accuracy of the controlled rotation angular coefficients, which are the critical factors in the HHL algorithm. We show that our algorithm not only increases the accuracy of the solutions, but also composes more practical and scalable circuits by dynamically controlling problem size in the NISQ devices. We present both simulated and experimental results, and discuss the sources of errors. Finally, we conclude that overall results on the quantum hardware are dominated by the error in the CNOT gates.

Robson, Walter↗
DOE OSTI2020

Multi-Energy Microgrid Planning Considering Heat Flow Dynamics

In this paper, we present a multi-energy microgrid optimal planning method, considering the intra-hour dynamics of the heating system as constraints of the energy dispatch, and consequently of the sizing of the Combined Heat and Power units. The fixed-mass flow rate method and finite difference technique are applied to linearize the dynamic constraints. Moreover, to determine the values of the mass flow rates, a critical scenario representing the slowest heat dynamic process is put forward. Finally, a realistic microgrid case study is presented to illustrate the effect of considering the heat flow dynamics in the overall solution.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Solving Seismic Wave Equations on Variable Velocity Models With Fourier Neural Operator

Here, in the study of subsurface seismic imaging, solving the acoustic wave equation is a pivotal component in existing models. The advancement of deep learning (DL) enables solving partial differential equations (PDEs), including wave equations, by applying neural networks to identify the mapping between the inputs and the solution. This approach can be faster than traditional numerical methods when numerous instances are to be solved. Previous works that concentrate on solving the wave equation by neural networks consider either a single velocity model or multiple simple velocity models, which is restricted in practice. Instead, inspired by the idea of operator learning, this work leverages the Fourier neural operator (FNO) to effectively learn the frequency domain seismic wavefields under the context of variable velocity models. We also propose a new framework paralleled FNO (PFNO) for efficiently training the FNO-based solver given multiple source locations and frequencies. Numerical experiments demonstrate the high accuracy of both FNO and PFNO with complicated velocity models in the OpenFWI datasets. Furthermore, the cross-dataset generalization test verifies that PFNO adapts to out-of-distribution velocity models. Finally, PFNO admits higher computational efficiency on large-scale testing datasets than the traditional finite-difference method. The aforementioned advantages endow the FNO-based solver with the potential to build powerful models for research on seismic waves.

58 GEOSCIENCES↗
DOE OSTI2023

Mechanism of structural colors in binary mixtures of nanoparticle-based supraballs

Inspired by structural colors in avian species, various synthetic strategies have been developed to produce noniridescent, saturated colors using nanoparticle assemblies. Nanoparticle mixtures varying in particle chemistry and size have additional emergent properties that affect the color produced. For complex multicomponent systems, understanding the assembled structure and a robust optical modeling tool can empower scientists to identify structure-color relationships and fabricate designer materials with tailored color. Here, we demonstrate how we can reconstruct the assembled structure from small-angle scattering measurements using the computational reverse-engineering analysis for scattering experiments method and use the reconstructed structure in finite-difference time-domain calculations to predict color. We successfully, quantitatively predict experimentally observed color in mixtures containing strongly absorbing nanoparticles and demonstrate the influence of a single layer of segregated nanoparticles on color produced. The versatile computational approach that we present is useful for engineering synthetic materials with desired colors without laborious trial-and-error experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Rehydrated glass embayments record the cooling of a Yellowstone ignimbrite

Abstract Hydration fronts penetrate 50–135 μm into glassy rhyolite embayments hosted in quartz crystals from the Mesa Falls Tuff in the Yellowstone Plateau volcanic field. The hydration fronts occur as steep enrichments that reach 2.4 ± 0.6 wt% H2O at the embayment opening, representing much higher values than interior concentrations of 0.9 ± 0.2 wt% H2O. Molecular water accounts for most of the water enrichment. Water speciation indicates the hydration fronts comprise absorbed meteoric water that modified the original magmatic composition of the rhyolitic glass. We used finite difference diffusion models to demonstrate that glass rehydration was likely produced over a few decades as the ignimbrite cooled. Such temperatures and time scales are consistent with rare firsthand observations of decadal hydrothermal systems associated with cooling ignimbrites at Mount Pinatubo (Philippines) and the Valley of Ten Thousand Smokes (Alaska).

Geology↗
DOE OSTI2021

Computing Edge States without Hard Truncation

We present a numerical method which accurately computes the discrete spectrum and associated bound states of semi-infinite Hamiltonians which model electronic “edge” states localized at boundaries of one- and two-dimensional crystalline materials. The problem is nontrivial since arbitrarily large finite “hard” (Dirichlet) truncations of the Hamiltonian in the infinite bulk direction tend to produce spurious bound states partially supported at the truncation. Our method, which overcomes this difficulty, is to compute the Green's function of the semi-infinite Hamiltonian by imposing an appropriate boundary condition in the bulk direction; then, the spectral data is recovered via Riesz projection. We demonstrate our method's effectiveness by studies of edge states at a graphene zig-zag edge in the presence of defects modeled both by a discrete tight-binding model and a continuum PDE model under finite difference discretization. Our method may also be used to study states localized at domain wall-type edges in one- and two-dimensional materials where the edge Hamiltonian is infinite in both directions; we demonstrate this for the case of a tight-binding model of distinct honeycomb structures joined along a zig-zag edge. Here, we expect our method to be useful for designing novel devices based on precise wave-guiding by edge states.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Improved Water Sound Synthesis using Coupled Bubbles

We introduce a practical framework for synthesizing bubble-based water sounds that captures the rich inter-bubble coupling effects responsible for low-frequency acoustic emissions from bubble clouds. We propose coupled-bubble oscillator models with regularized singularities, and techniques to reduce the computational cost of time stepping with dense, time-varying mass matrices. Airborne acoustic emissions are estimated using finite-difference time-domain (FDTD) methods. We propose a simple, analytical surface-acceleration model, and a sample-and-hold GPU wavesolver that is simple and faster than prior CPU wavesolvers. Sound synthesis results are demonstrated using bubbly flows from incompressible, two-phase simulations, as well as procedurally generated examples using single-phase FLIP fluid animations. Our results demonstrate sound simulations with hundreds of thousands of bubbles, and perceptually significant frequency transformations with fuller low-frequency content.

Computer Science↗
DOE OSTI2025

Tensorized Interior Radiative Heat Transfer for a Scalable and Calibrated Building Energy Simulator

Building energy simulation is a critical tool for developing and testing advanced control strategies, such as Reinforcement Learning (RL), to provide demand flexibility and affordable energy costs. The recently introduced Smart Buildings Control Suite (sbsim) provides a lightweight, scalable, and data-calibrated simulation environment based on a 2D finite-difference model. However, the initial model primarily focused on conductive and convective heat transfer, neglecting the significant impact of long-wave radiative heat exchange between interior surfaces. This paper presents a significant extension to the sbsim framework by incorporating a physically-grounded model for interior radiative heat transfer. Our primary contribution is the development and integration of a fully tensorized radiative heat transfer module, which preserves the computational efficiency and scalability of the original simulator. This was achieved by developing a pipeline for view factor calculation, including an algorithm to identify directly seeing surfaces within complex floor plans, and formulating the net radiation equations for efficient execution on modern hardware accelerators. We validate the numerical accuracy of our tensorized implementation by comparing its results against a traditional iterative approach, demonstrating identical outcomes. This enhancement increases the physical fidelity of sbsim, enabling more accurate training of RL agents for building energy optimization.

Ham, Sang woo↗