Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical solution”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Verification of an energy-conserving semi-implicit electrostatic particle-in-cell scheme for modeling high-density plasma at scale

A verification study of a semi-implicit energy-conserving electrostatic particle-in-cell algorithm is presented. The algorithm relaxes the time-step and mesh-size constraints that require resolution of the plasma period and Debye length associated with traditional explicit momentum-conserving particle-in-cell algorithms. Physical implications and applicability of using the semi-implicit scheme for modeling high-density plasmas are discussed. Where possible, numerical results are compared against analytical solutions. The simulation results indicate that the algorithm is stable at time steps larger than twice the inverse plasma frequency and cell sizes larger than the Debye length. It is found that the algorithm gives adequate results, provided that the distribution function and the spatiotemporal scales dictating the physics of the problem are resolved. As such, the algorithm may provide a robust method for kinetic modeling of high-density plasmas at scale.

Cyclotron resonance↗
DOE OSTI2021

Symmetrically tuned large-volume conic shell-cavities for axion searches

In an earlier paper [1], a new class of thin-shell cavities were proposed to evade the steep frequency scaling of conventional axion haloscopes. In this follow-up work, we see that a generalized conic geometry enables robust frequency-tuning for these large-volume cm-wave cavities. The frequency-defining dimension of a conic shell-cavity changes symmetrically and uniformly during tuning, maintaining a high axion coupling efficiency (the "form factor") to an external solenoid field. It is further shown that such tunable geometry is not restricted to circular cones. A general prescription for arbitrary volume-filling conic shell-cavities is developed and direct solutions are obtained for the created numerical models. Additionally, the largest of the realized designs is a meandering "brain" cavity that is tunable over a frequency range of 20%. The scan rate of this cavity is three orders of magnitude larger than that of a scaled cylindrical cavity used in the current generation experiments. The prospects for such a large improvement in the scan rate should motivate R& D efforts in fabrication and other implementation techniques. If these engineering challenges can be met, cavity-based axion haloscopes can stay competitive at frequencies higher than a few GHz. We propose an experimental configuration at 20 GHz (~80 μeV) using an array of brain cavities and compare it with other proposals for similar frequencies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Asymmetric crossing of the attractive and repulsive magnetic potential by Abrikosov vortices

Here, we present studies of Abrikosov vortex motion across the magnetically charged edges of long thin ferromagnetic stripes placed above or under a thin superconducting film. The magnetic charges at the stripe edges form attractive or repulsive potential wells for vortices. Using a relatively small in-plane magnetic field to polarize the stripe edges and normal-to-plane magnetic field to induce up or down polarized vortices, we tune the attractive or repulsive stripe-vortex interactions. Imaging of the vortex dynamics reveals that the repulsive magnetic potential U m + acts as a robust vortex pinning barrier, while the attractive U m - has practically no effect on the vortex motion irrespective of the position of the stripes located above or underneath the superconducting film. We analyze the observed asymmetry using equations of the overdamped vortex motion. The formal analytical solution yields an asymmetry, but the numerical modeling with and without noise terms does not confirm it. Instead, we find that the asymmetry is caused by the creation of spontaneous vortices at the maxima of U m , which depends on the edge polarity and suppress U m - . Furthermore, our experiment and time-dependent Ginzburg-Landau simulations demonstrate that magnetic pinning dominates over vortex pinning due to corrugations of the superconducting layer deposited on top of the ferromagnetic stripes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Fast and robust all-electron density functional theory calculations in solids using orthogonalized enriched finite elements

Here, we present a computationally efficient approach to perform systematically convergent real-space all-electron Kohn-Sham density functional theory calculations for solids using an enriched finite element (FE) basis. The enriched FE basis is constructed by augmenting the classical FE basis with atom-centered numerical basis functions, comprising of atomic solutions to the Kohn-Sham problem. Notably, to improve the conditioning, we orthogonalize the enrichment functions with respect to the classical FE basis, without sacrificing the locality of the resultant basis. In addition to improved conditioning, this orthogonalization procedure also renders the overlap matrix block diagonal, greatly simplifying its inversion. Subsequently, we use a Chebyshev polynomial based filtering technique to efficiently compute the occupied eigenspace in each self-consistent field iteration. We demonstrate the accuracy and efficiency of the proposed approach on periodic unit cells and supercells. The benchmark studies show a staggering 130× speedup of the orthogonalized enriched FE basis over the classical FE basis. We also present a comparison of the orthogonalized enriched FE basis with the linearized augmented plane-wave + local orbitals basis, both in terms of accuracy and efficiency. Notably, we demonstrate that the orthogonalized enriched FE basis can handle large system sizes of ~10 000 electrons. Finally, we observe good parallel scalability of our implementation with 92% efficiency at 22× speedup for a system with 620 electrons

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

PMU-Based Decoupled State Estimation for Unsymmetrical Power Systems

Modal decomposition of measurement equations has already been shown to simplify the formulation and resulting computational complexity of three-phase state estimation of systems where all the transmission lines are three-phase and fully transposed. When there are non-transposed and/or mixed-phase lines, modal decomposition can no longer fully decouple the threephase measurement equations. Here, this paper addresses the above shortcoming by proposing a simple yet practical solution based on the commonly used numerical compensation techniques. Thus, it enables application of the powerful decoupling approach to any type of three-phase networks which may contain non-transposed or mixed-phase lines and are fully observable by PMUs. The proposed procedure modifies the measurement set by deriving additive terms that compensate for the neglected unsymmetrical effects. It will be shown that unbalanced systems including nontransposed and mixed-phase elements, can still be transformed into three decoupled subsystems and solved in parallel by the proposed approach. Performance of the proposed algorithm is validated against several IEEE test cases.

42 ENGINEERING↗
DOE OSTI2020

Analytical expressions of the surface shape of ‘diaboloid’ mirrors

Modern deterministic polishing processes allow fabrication of x-ray optics with almost any arbitrary aspherical surface shape. Among these optics, the so called "diaboloid"mirror is of special interest. The diaboloid mirror that converts a cylindrical wave to a spherical wave would improve focusing in x-ray beamlines implementing a diffraction element between a parabolic cylinder and a toroidal mirror. The replacement of the toroidal mirror in existing beamlines by the diaboloid mirror would mitigate aberrations. The shape of the diaboloid mirror is usually calculated numerically based on a truncated polynomial solution of the optical path problem. Here, we present an exact analytical solution for the shape of a diaboloid mirror as a function of the conjugate parameters of the mirror placed in a beamline. The derived analytical expressions for the diaboloid mirror in both the canonical and mirror-based coordinate systems are implemented in ray-tracing simulations to verify the beamline performances.

Yashchuk, Valeriy V.↗
DOE OSTI2022

High-dimensional discrete Fourier transform gates with a quantum frequency processor

The discrete Fourier transform (DFT) is of fundamental interest in photonic quantum information, yet the ability to scale it to high dimensions depends heavily on the physical encoding, with practical recipes lacking in emerging platforms such as frequency bins. In this article, we show that d -point frequency-bin DFTs can be realized with a fixed three-component quantum frequency processor (QFP), simply by adding to the electro-optic modulation signals one radio-frequency harmonic per each incremental increase in d . We verify gate fidelity F W > 0.9997 and success probability P W > 0.965 up to d = 10 in numerical simulations, and experimentally implement the solution for d = 3, utilizing measurements with parallel DFTs to quantify entanglement and perform tomography of multiple two-photon frequency-bin states. Our results furnish new opportunities for high-dimensional frequency-bin protocols in quantum communications and networking.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Pulse Detonation Engine for Advanced Oxy-Combustion of Coal-Based Fuel for Direct Power Extraction Applications

Pressure gain combustion and magnetohydrodynamic (MHD) systems have the potential to provide a step increase in the efficiency of combined-cycle power plants. Specific advantages include a net pressure increase to the system instead of a pressure drop, the high temperature of the detonation waves can increase the efficiency of power extraction compared to other processes, significant thermal energy can be released in a compact region, and the high velocities of the flow increase extraction of electrical power. In summary, a pressure gain combustor coupled with a MHD has the potential to be transformative. Despite the potential advantages, relatively little research has been conducted considering coupled pressure gain combustion systems with MHD systems. With this background and motivation, the overall goal of this effort was to advance the knowledge, technology, and computational tools associated with coupled detonation and MHD systems. A joint experimental and computational approach was used while seeking to accomplish the goals of this work. Specifically, two pulse detonation engines were developed and used for the experiments to produce detonations. Detonation speeds were measured for a variety of flow and fuel conditions (e.g., methane, propane, with coal particles). Preliminary electrical conductivity measurements were collected. An extensive amount of research was performed to identify sensitivities of detonation behaviors to the presence of combustion products. Computationally, a twofold approach was used in this work. First, a solver was developed for solving the governing equations for a reactive flow with coupled detonation and ionization chemistry. The solver was applied to study the impacts of seed material ionization on detonation. Second, a conservation element-solution element (CE-SE) based numerical solver for detonation studies with a reduced reaction mechanism for oxy-methane combustion was developed and verified on standard test cases. Key findings and contributions from this work are as follows. A system was developed for injecting powderized coal, or other seeding material, into a pulse-detonation system. The influence of a combustion product (i.e., CO 2 ) on detonation behavior was identified. Knowledge gained from this work is applicable to devices such as rotating detonation engines, where combustion products mix with fresh reactants. A system for measuring the electrical conductivity of the exhaust from a pulse-detonation engine was developed. The open-source solver, Clawpack, was extended to solve the reactive Euler equations for simulating detonations. A coupled combustion and ionization chemistry was developed in a single chemical kinetic model for methane oxidation. This model can be used to solve coupled MHD and detonation simulations. It was found that parasitic interactions from ionization chemistry with the magnetic field can reduce the detonation velocity by up to 8%, with a potential impact on power extraction of 15%. It is recommended that interactions between the detonation front and MHD field be considered. Higher gas temperatures and velocities were achieved owing to oxy-fuel detonations. Use of radical dissociation reactions in the reduced reaction mechanism, was found to be critical in predicting detonation temperature and velocity accurately.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Substation-Level Grid Topology Optimization Using Bus Splitting

Operations of substation circuit breakers are important for maintenance needs and topology reconfiguration in power systems. Bus splitting is one type of topology change where the two bus bars at a substation can become electrically disconnected under certain actions of circuit breakers. Because these events involve detailed substation modeling, they are typically not considered in routine power system operation and control. In this paper, an improved substation-level topology optimization framework is developed by expanding traditional line switching decisions by breaker-level bus splitting, which can further reduce grid congestion and generation costs. A tight McCormick relaxation is proposed to reformulate the bilinear terms in the resultant optimization problem to linear inequality constraints. Thus, a tractable mixed-integer linear program reformulation is attained that allows for efficient solutions in real-time operations. Numerical studies on the IEEE 14-bus and 118-bus systems demonstrate the computational performance and economic benefits of the proposed topology optimization approach.

bus split↗
DOE OSTI2025

Nonlinear Optimal Control of Electron Dynamics Within Hartree-Fock Theory

Consider the problem of determining the optimal applied electric field to drive a molecule from an initial state to a desired target state. For even moderately sized molecules, solving this problem directly using the exact equations of motion—the time-dependent Schrödinger equation (TDSE)—is numerically intractable. Here, we present a solution of this problem within time-dependent Hartree-Fock (TDHF) theory, a mean field approximation of the TDSE. Optimality is defined in terms of minimizing the total control effort while maximizing the overlap between desired and achieved target states. We frame this problem as an optimization problem constrained by the nonlinear TDHF equations; we solve it using trust region optimization with gradients computed via a custom-built adjoint state method. For three molecular systems, we show that with very small neural network parametrizations of the control, our method yields solutions that achieve desired targets within acceptable constraints and tolerances.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Adaptive Mesh Refinement Simulations for Turbulent Reacting Flow

With the increased availability of exascale computing hardware, detailed simulations of realistic devices can be performed at practically relevant time and length scales. Insights into the multiscale driving mechanisms in compressible reacting flow systems with complex geometry, such as combustors, can be used for design optimization and technology improvements. However, to effectively perform these simulations, advanced numerical algorithms must be used to maintain solution accuracy without incurring undue computational costs. PeleC, part of the Pele suite of codes, leverages block-structured adaptive mesh refinement (AMR) through the AMReX library to capture fine-scale flow features in compressible reacting flows. In this talk, we discuss recent improvements to the numerical algorithms, particularly in regard to describing flows at complex boundary structures, and PeleC's performance on exascale computing hardware. We will demonstrate that PeleC is well-suited for modern, extreme-scale, heterogenous compute platforms.

combustion↗
DOE OSTI2020

Facile and quantitative estimation of strain in nanobubbles with arbitrary symmetry in 2D semiconductors verified using hyperspectral nano-optical imaging

When layers of van der Waals materials are deposited via exfoliation or viscoelastic stamping, nanobubbles are sometimes created from aggregated trapped fluids. Though they can be considered a nuisance, nanobubbles have attracted scientific interest in their own right owing to their ability to generate large in-plane strain gradients that lead to rich optoelectronic phenomena, especially in the semiconducting transition metal dichalcogenides. Determination of the strain within the nanobubbles, which is crucial to understanding these effects, can be approximated using elasticity theory. However, the Föppl–von Kármán equations that describe strain in a distorted thin plate are highly nonlinear and often necessitate assuming circular symmetry to achieve an analytical solution. Here, we present an easily implemented numerical method to solve for strain tensors of nanobubbles with arbitrary symmetry in 2D crystals. The method only requires topographic information from atomic force microscopy and the Poisson ratio of the 2D material. We verify that this method reproduces the strain for circularly symmetric nanobubbles that have known analytical solutions. Finally, we use the method to reproduce the Grüneisen parameter of the E' mode for 1L-WS2 nanobubbles on template-stripped Au by comparing the derived strain with measured Raman shifts from tip-enhanced Raman spectroscopy, demonstrating the utility of our method for estimating localized strain in 2D crystals.

2D materials↗
DOE OSTI2026

Semi-analytic solutions to the Noh problem with a black box EoS

The objective of this paper is to derive a method of constructing semi-analytic solutions to the Noh problem when the equation of state is a black box. Such solutions can be used for verification tests of hydrodynamics codes. We present the underlying theory, the method for finding solutions, and several examples of derived semi-analytic solutions. We end by performing a classic verification convergence test comparing numerical results from a hydrodynamics code against a non-trivial semi-analytic solution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Numerically exact generalized Green's function cluster expansions for electron-phonon problems

We generalize the family of approximate momentum average methods to formulate a numerically exact, convergent hierarchy of equations whose solution provides an efficient algorithm to compute the Green's function of a particle dressed by bosons suitable in the entire parameter regime. We use this approach to extract ground-state properties and spectral functions. Our approximation-free framework, dubbed the generalized Green's function cluster expansion (GGCE), allows access to exact numerical results in the extreme adiabatic limit, where many standard methods struggle or completely fail. We showcase the performance of the method, specializing three important models of charge-boson coupling in solids and molecular complexes: the molecular Holstein model, which describes coupling between charge density and local distortions, the Peierls model, which describes modulation of charge hopping due to intersite distortions, and a more complex Holstein + Peierls system with couplings to two different phonon modes, paradigmatic of charge-lattice interactions in organic crystals. Furthermore, the GGCE serves as an efficient approach that can be systematically extended to different physical scenarios, thus providing a tool to model the frequency dependence of dressed particles in realistic settings.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Fabrication of Biofunctionalized PVA/PEI Nanofibers for Rare Earth Element Extraction from Coal Mine Tailings

The need for rare earth elements (REEs) for new technology is becoming increasingly stiff with demand at a record high. Biological ligands like lanmodulin and lanthanide binding tags (LBT) have gained noticeable attention for their selective extraction of valuable lanthanide metals in aqueous solution. Building on previous work using microbeads as the extraction media, we propose that nanofibers can improve extraction capabilities. Nanofibers could provide more efficient flow and higher extraction yields due to their high surface area to volume ratio. Relying on thiolmaleimide click chemistry as the biological ligand immobilization route, electrospinning as the nanofiber manufacturing method, and numerous characterization techniques, we produced electrospinning polymer solutions and nanofibers with a maximum neodymium adsorbent capacity of 12.1 μmol/g of extraction media immobilized with LBT. Although the prior work produced an extraction media with an equivalent capacity of ~40 μmol/g, the electrospinning process is much more economical and scalable. This work paves the way for an industrially relevant, environmentally-friendly REE extraction material. Further work is required to produce higher quality nanofibers, to investigate common ion interference with REE adsorbance, and to scale production up to a relevant scale.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Transport in Stochastic Media with Random Chord Length Distributions

Thermal radiation transport computations in binary Markovian random mixtures rely almost exclusively on the Levermore-Pomraning (LP) model which is obtained by applying a heuristic closure to the ensemble averaged random medium transport equation. The validity of this model has been extensively tested by comparing numerical results over a broad parameter range (material types and mixing parameters) against benchmark solutions in planar geometry. The conditions under which the LP-model provides useful results and when it breaks down are now well established, but work to date has been largely restricted to homogeneous mixing statistics, i.e., the mean chord lengths of both materials are taken to be spatially constant. In recent work, this limitation was relaxed by allowing the mean chord lengths and, in a consistent fashion, the volume fractions in the LP-model to vary spatially and in a follow-up investigation benchmark solutions were obtained by ensemble averaging results over material realizations sampled from a nonhomogeneous Poisson process (NHPP). Numerical experiments in rod geometry with specifically linear and quadratic spatial dependence of chord lengths showed that the material averaged radiation intensities vary nonmonotonically with depth into the medium, in stark contrast to solutions obtained assuming uniform chord lengths. Moreover, depending on the local optical thickness and strength of scattering, the LP-model results showed locally more nuanced deviations from the benchmark solutions than was the case with constant chord lengths. These limited numerical investigations highlight the nontrivial qualitative and quantitative consequences of nonhomogeneous mixing statistics, in particular that closure approximations may not be uniformly valid or invalid over the problem domain.

42 ENGINEERING↗
DOE OSTI2020

Newton versus the machine: solving the chaotic three-body problem using deep neural networks

ABSTRACT Since its formulation by Sir Isaac Newton, the problem of solving the equations of motion for three bodies under their own gravitational force has remained practically unsolved. Currently, the solution for a given initialization can only be found by performing laborious iterative calculations that have unpredictable and potentially infinite computational cost, due to the system’s chaotic nature. We show that an ensemble of converged solutions for the planar chaotic three-body problem obtained using an arbitrarily precise numerical integrator can be used to train a deep artificial neural network (ANN) that, over a bounded time interval, provides accurate solutions at a fixed computational cost and up to 100 million times faster than the numerical integrator. In addition, we demonstrate the importance of training an ANN using converged solutions from an arbitrary precise integrator, relative to solutions computed by a conventional fixed precision integrator, which can introduce errors in the training data, due to numerical round-off and time discretization, that are learned by the ANN. Our results provide evidence that, for computationally challenging regions of phase space, a trained ANN can replace existing numerical solvers, enabling fast and scalable simulations of many-body systems to shed light on outstanding phenomena such as the formation of black hole binary systems or the origin of the core collapse in dense star clusters.

Breen, Philip G.↗
DOE OSTI2020

Implicit fast sweeping method for hyperbolic systems of conservation laws

Implicit time-accurate methods are often used to integrate stiff problems where explicit schemes impose severe time step restrictions. This paper presents an efficient numerical framework based on the Fast Sweeping Method (FSM) for solving linear and nonlinear hyperbolic systems of conservation laws. The solution at each discrete location is computed by sweeping the numerical domain in several predetermined directions that follow the causality of the characteristic families. The use of a fractional step strategy eliminates the need for a solution selection criterion while one-sided stencils limit the number of sweeps to at most 2 d for d space dimensions. This work focuses on the first-order implicit upwind method since it constitutes the building block for high-order conservative schemes. For problems where the degree of stiffness evolves over time, implicit-explicit hybridization can be accomplished with the same algorithm by simply switching the stencil at each time level. As opposed to traditional implicit solvers, the sweeping method does not require a local time linearization of the fluxes thereby preserving the nonlinear stability properties of the original implicit scheme. It also avoids the large computational and memory requirements associated with solving large block-diagonal systems of equations. Here, a series of one- and two-dimensional test cases are presented for the inviscid Burgers' equation and the reactive Euler equations. The results indicate that the implicit FSM can allow a major reduction in the number of time steps even in the presence of discontinuous solution profiles.

74 ATOMIC AND MOLECULAR PHYSICS↗