Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “high order”

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

Conversion of CO 2 into Synthetic Graphite

Carbon Dioxide (CO 2 ) can be utilized as a source for producing synthetic graphite. Our method converts CO 2 into solid carbon. These CO 2 -sourced solid carbon have an anisotropic property that allow them to be graphitized into synthetic graphite. Our patented thermal catalytic process uses carbon dioxide and hydrogen as the feedstocks for the reaction, producing solid carbon material, with distilled water as the sole byproduct. The process is designed as a closed loop, where all feedstock gases entering the reactor are converted into durable bulk carbons and water, with no emissions. These carbons materials can be sintered into a solid structure and then thermally annealed, taking the carbon atoms from a disordered state to a highly ordered state of graphitic material – producing synthetic graphite structures. Phase 1 will prove that these bulk carbons can be sintered into a machinable solid carbon structure, which will then be annealed into a fully graphitic solid structure. Standard analytics will be done to verify the extent of graphitization (Raman, XRD). Proving this CO 2 -sourced carbon pathway to synthetic graphite structures will give the US a domestic supply route for this critical material and remove reliance on foreign sources. In addition, since the synthetic graphite structure is made from CO 2 -sourced carbon, this is also a method to mitigate CO 2 . Graphite structures are used in high temperature and high friction applications, for example, as anodes for steel refining, for brake pad linings, as neutron reflectors in nuclear reactors.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Highly Symmetric, Self-Assembling 3D DNA Crystals with Cubic and Trigonal Lattices

The rational design of nanoscopic DNA tiles has yielded highly ordered crystalline matter in 2D and 3D. The most well-studied 3D tile is the DNA tensegrity triangle, which is known to self-assemble into macroscopic crystals. However, contemporary rational design parameters for 3D DNA crystals nearly universally invoke integer numbers of DNA helical turns and Watson–Crick (WC) base pairs. Here, in this study, 24-bp edges are substituted into a previously 21-bp (two helical turns of DNA) tensegrity triangle motif to explore whether such unconventional motif can self-assemble into 3D crystals. The use of noncanonical base pairs in the sticky ends results in a cubic arrangement of tensegrity triangles with exceedingly high symmetry, assembling a lattice from winding helical axes and diamond-like tessellation patterns. Reverting this motif to sticky ends with Watson–Crick pairs results in a trigonal hexagonal arrangement, replicating this diamond arrangement in a hexagonal context. These results showcase that the authors can generate unexpected, highly complex, pathways for materials design by testing modifications to 3D tiles without prior knowledge of the ensuing symmetry. This study expands the rational design toolbox for DNA nanotechnology; and it further illustrates the existence of yet-unexplored arrangements of crystalline soft matter.

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

Quantum prediction of ultra-low thermal conductivity in lithium intercalation materials

Lithium-intercalated layered transition-metal oxides, Li x TMO 2, brought about a paradigm change in rechargeable batteries in recent decades and show promise for use in memristors, a type of device for future neural computing and on-chip storage. Thermal transport properties, although being a crucial element in limiting the charging/discharging rate, package density, energy efficiency, and safety of batteries as well as the controllability and energy consumption of memristors, are poorly managed or even understood yet. Here, for the first time, we employ quantum calculations including high-order lattice anharmonicity and find that the thermal conductivity κ of Li x TMO 2 materials is significantly lower than hitherto believed. More specifically, the theoretical upper limit of κ of LiCoO 2 is 6 W/m-K, 2–6 times lower than the prior theoretical predictions. Delithiation further reduces κ by 40–70% for LiCoO 2 and LiNbO 2 . Grain boundaries, strain, and porosity are yet additional causes of thermal-conductivity reduction, while Li-ion diffusion and electrical transport are found to have only a minor effect on phonon thermal transport. The results elucidate several long-standing issues regarding the thermal transport in lithium-intercalated materials and provide guidance toward designing high-energy-density batteries and controllable memristors.

25 ENERGY STORAGE↗
DOE OSTI2024

A Verification Suite of Test Cases for the Barotropic Solver of Ocean Models

Abstract The development of any atmosphere or ocean model warrants a suite of test cases (TCs) to verify its spatial and temporal discretizations, order of accuracy, stability, reproducibility, portability, scalability, etc. In this paper, we present a suite of shallow water TCs designed to verify the barotropic solver of atmosphere and ocean models. These include the non‐dispersive coastal Kelvin wave; the dispersive inertia‐gravity wave; the dispersive planetary and topographic Rossby waves; the barotropic tide; and a non‐linear manufactured solution. These TCs check the implementation of the linear pressure gradient term; the linear constant or variable‐coefficient Coriolis and bathymetry terms; and the non‐linear advection terms. Simulation results are presented for a variety of time‐stepping methods as well as two spatial discretizations: a mimetic finite volume method based on the TRiSK scheme, and a high‐order discontinuous Galerkin spectral element method. The experimental procedure for conducting these numerical experiments is detailed. It underscores several key considerations that vary depending on the chosen spatial discretization method. Finally, convergence studies of every TC are conducted with refinement in both space and time, only in space, and only in time. The convergence slopes match the expected theoretical predictions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Retarded Charge-Carrier Recombination in Photoelectrochemical Cells from Plasmon-Induced Resonance Energy Transfer

N-type metal oxides such as hematite (α-Fe 2 O 3 ) and bismuth vanadate (BiVO 4 ) are promising candidate materials for efficient photoelectrochemical water splitting; however, their short minority carrier diffusion length and restricted carrier lifetime result in undesired rapid charge recombination. In this paper, a 2D arranged globular Au nanosphere (NS) monolayer array with a highly ordered hexagonal hole pattern (hereafter, Au array) is introduced onto the surface of photoanodes comprised of metal oxide films via a facile drying and transfer-printing process. Through plasmon-induced resonance energy transfer, the Au array provides a strong electromagnetic field in the near-surface area of the metal oxide film. The near-field coupling interaction and amplification of the electromagnetic field suppress the charge recombination with long-lived photogenerated holes and simultaneously enhance the light harvesting and charge transfer efficiencies. Consequently, an over 3.3-fold higher photocurrent density at 1.23 V versus reversible hydrogen electrode (RHE) is achieved for the Au array/α-Fe 2 O 3 . Furthermore, the high versatility of this transfer printing of Au arrays is demonstrated by introducing it on the molybdenum-doped BiVO 4 film, resulting in 1.5-fold higher photocurrent density at 1.23 V versus RHE. The tailored metal film design can provide a potential strategy for the versatile application in various light-mediated energy conversion and optoelectronic devices.

2D pattern array↗
DOE OSTI2024

Optimal local truncation error method for 3-D elasticity interface problems

The paper deals with a new effective numerical technique on unfitted Cartesian meshes for simulations of heterogeneous elastic materials. Here, we develop the optimal local truncation error method (OLTEM) with 27- point stencils (similar to those for linear finite elements) for the 3-D time-independent elasticity equations with irregular interfaces. Only displacement unknowns at each internal Cartesian grid point are used. The interface conditions are added to the expression for the local truncation error and do not change the width of the stencils. The unknown stencil coefficients are calculated by the minimization of the local truncation error of the stencil equations and yield the optimal second order of accuracy for OLTEM with the 27-point stencils on unfitted Cartesian meshes. A new post-processing procedure for accurate stress calculations has been developed. Similar to basic computations it uses OLTEM with the 27-point stencils and the elasticity equations. The post-processing procedure can be easily extended to unstructured meshes and can be independently used with existing numerical techniques (e.g., with finite elements). Numerical experiments show that at an accuracy of 0.1% for stresses, OLTEM with the new post-processing procedure significantly (by 10 5 -10 9 times) reduces the number of degrees of freedom compared to linear finite elements. OLTEM with the 27-point stencils yields even more accurate results than high-order finite elements with wider stencils.

42 ENGINEERING↗
DOE OSTI2025

Anodized Aluminum Oxide Membrane Ionic Memristors

Memory effect in ion transport (IT) at the solid–solution interface is uniquely attractive in that the conductance depends on or “memorizes” the previous states. Hysteretic and rectified transport properties offer exciting potential to developing advanced iontronics and neuromorphic functions, improving the efficiency of energy conversion and electrochemical processes, and overcoming the selectivity-throughput bottleneck in the enrichment of low abundant species for environment- and energy-friendly separations, among others. Herein, memory effects are discovered in the rectified electrokinetic IT through anodized aluminum oxide (AAO) membranes containing densely packed highly ordered nanochannels (10 10 per cm 2 ). Characteristic memristor responses of pinched current–potential loops are resolved in voltammetric experiments and successfully reproduced through finite element simulation. Excitatory and inhibitory conductance states are shown to arise from the enrichment and depletion of mobile charge carriers. Structurewise, the transport symmetry is broken by the barrier oxide layer (BOL) on the one end of the cylindrical nanochannels across the AAO membranes. Charge selectivity is attributed to the gradient(s) of the space charge density across the BOL characterized by depth profiling via X-ray photoelectron spectroscopy analysis. The space charge gradient(s) overcomes the fundamental limitation of widely exploited surface charge effects to enable intense rectification and hysteresis prevailing at very high ionic concentrations up to 1–2 M. A new strategy is developed for controlling the preferential IT direction and selectivity via counterion intercalation and extraction/exchange. Mechanistic understanding is further confirmed through parameter variations such as potential scan rate and ionic strength, which also demonstrates convenient controls of the related functions.

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

Two-dimensional phase cartography for high-harmonic spectroscopy

Extremely nonlinear spectroscopy based on high-order-harmonic generation has become a powerful investigation method for attosecond dynamics in gas and solid targets. In particular, the phase of harmonic emission was shown to carry profound insight into atomic and molecular structure and dynamics. However, current techniques offer phase measurements only along specific directions, thus providing partial characterization. Here we report on a new approach combining optical and quantum interferometers measuring along two dimensions the intensity and phase of harmonic emission from aligned molecules in the exact same experimental conditions. This two-dimensional cartography technique measures the phase with no arbitrary offset and no uncertainty on its sign. Measurements along different dimensions can be combined in two ways: either a single mapping or a redundant mapping allowing high-precision phase recovery using a Shack–Hartmann-like algorithm. We demonstrate both methods in a nitrogen test case, which allows disentangling structural and dynamical effects. Two-dimensional phase cartography paves the way to high-resolution high-harmonic spectroscopy for applications such as quantum orbital tomography and attosecond charge migration in molecules.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Slimmer Geminals For Accurate F12 Electronic Structure Models

The Slater-type F12 geminal length scales originally tuned for the second-order Mo̷ller-Plesset F12 method are too large for higher-order F12 methods formulated using the SP (diagonal fixed-coefficient spin-adapted) F12 ansatz. The new geminal parameters reported herein reduce the basis set incompleteness errors (BSIEs) of absolute coupled-cluster singles and doubles F12 correlation energies by a significant─and increase with the cardinal number of the basis─margin. The effect of geminal reoptimization is especially pronounced for the cc-pVXZ-F12 basis sets (specifically designed for use with F12 methods) relative to their conventional aug-cc-pVXZ counterparts. The BSIEs of relative energies are less affected, but substantial reductions can be obtained, especially for atomization energies and ionization potentials with the cc-pVXZ-F12 basis sets. The new geminal parameters are therefore recommended for all applications of high-order F12 methods, such as coupled-cluster F12 methods and transcorrelated F12 methods.

Powell, Samuel R. [Virginia Polytechnic Inst. and ↗
DOE OSTI2021

On the possibility of footprint compression with one lens in nonlinear accelerator lattice

Electromagnetic interaction of colliding beams along with other nonlinear fields often limits the beams' lifetimes and luminosities. Nonlinearities result in the spread of betatron frequencies (footprint) and, thus, may enhance dynamic diffusion of particles due to high order resonances. One of the possible ways to eliminate nonlinearities and overcome the corresponding difficulties is compensation of nonlinear forces, but, in practice, it is hardly possible to obtain exact linearity of the system. The compensation with a single nonlinear lens cannot cope with distributed nonlinearities, nonlinearities due to parasitic crossings, etc. Here, we present a method to compute parameters of nonlinear element (lens) that eliminates both the footprint and resonance strength without achieving full compensation.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Single-Crystal Polycationic Polymers Obtained by Single-Crystal-to-Single-Crystal Photopolymerization

The efficient preparation of single-crystalline ionic polymers and fundamental understanding of their structure–property relationships at the molecular level remains a challenge in chemistry and materials science. Herein, we observe the single-crystal structure of a highly ordered polycationic polymer (polyelectrolyte) and its proton conductivity. The polyelectrolyte single crystals can be prepared on a gram-scale in quantitative yield, by taking advantage of an ultraviolet/sunlight-induced topochemical polymerization, from a tricationic monomer—a self-complementary building block possessing a preorganized conformation. A single-crystal-to-single-crystal photopolymerization was revealed unambiguously by in situ single-crystal X-ray diffraction analysis, which was also employed to follow the progression of molecular structure from the monomer, to a partially polymerized intermediate, and, finally, to the polymer itself. Collinear polymer chains are held together tightly by multiple Coulombic interactions involving counterions to form two-dimensional lamellar sheets (1 nm in height) with sub-nanometer pores (5 Å). The polymer is extremely stable under 254 nm light irradiation and high temperature (above 500 K). The extraordinary mechanical strength and environmental stability—in combination with its impressive proton conductivity (~3 × 10 –4 S cm –1 )—endow the polymer with potential applications as a robust proton-conducting material. By marrying supramolecular chemistry with macromolecular science, the outcome represents a major step toward the controlled synthesis of single-crystalline polyelectrolyte materials with perfect tacticity.

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

Architected mesoporous crystalline magnesium silicates with ordered pore structures

Novel approaches to harness earth abundant silicates as building blocks for carbon dioxide removal, capture, utilization, and storage are gaining increasing attention in the context of sustainable and low carbon energy and resource recovery. Advancing a calibrated understanding of these fluid-silicate interactions is essential for developing scalable processes. One of the challenges in developing predictive controls over these interactions is the compositional and morphological heterogeneity of naturally occurring, heterogeneous magnesium silicate minerals. To address this challenge, the synthesis of architected mesoporous crystalline magnesium silicate (Mg 2 SiO 4 ) is proposed. While synthesis routes for producing amorphous mesoporous magnesium silicates have been developed via sol-gel methods, approaches to synthesize crystalline magnesium silicates with well-controlled pore size distributions have not been explored. The conventional approaches of converting matter that is amorphous to crystalline states at elevated temperatures results in a heterogeneous pore size distribution. To develop controls on pore size distribution, amorphous mesoporous magnesium silicates are coated with carbon. This approach preserves the pore size distributions during the amorphous to crystalline transformations of Mg-silicates at elevated temperatures. The carbon coating is removed on heating. Magnesium silicate particles produced using this approach have highly ordered pores around 2.58 nm and a specific surface area of 124.25 m 2 /g. In this study, we report the chemical compositions, morphologies, phase transitions, and pore structures of the intermediate materials produced during the synthesis of crystalline mesoporous magnesium silicates. The transitions, and pore structures of the intermediate materials produced during the synthesis of crystalline mesoporous magnesium silicates. Furthermore, the synthesis routes discussed in this study can be applied translationally to produce metal silicates with ordered mesoporous structures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Direct, simple and efficient computation of all components of the virtual-casing magnetic field in axisymmetric geometries with Kapur–Rokhlin quadrature

In a recent publication (Toler et al., J. Plasma Phys., vol. 89, issue 2, 2023, p. 905890210), we demonstrated that for axisymmetric geometries, the Kapur–Rokhlin quadrature rule provided an efficient and high-order accurate method for computing the normal component, on the plasma surface, of the magnetic field due to the toroidal current flowing in the plasma, via the virtual-casing principle. The calculation was indirect, as it required the prior computation of the magnetic vector potential from the virtual-casing principle, followed by the computation of its tangential derivative by Fourier differentiation, to obtain the normal component of the magnetic field. Our approach did not provide the other components of the virtual-casing magnetic field. In this letter, we show that a more direct and more general approach is available for the computation of the virtual-casing magnetic field. The Kapur–Rokhlin quadrature rule accurately calculates the principal value integrals in the expression for all the components of the magnetic field on the plasma boundary, and the numerical error converges at a rate nearly as high as the indirect method we presented previously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Continuum electromagnetic gyrokinetic simulations of turbulence in the tokamak scrape-off layer and laboratory devices

Here, we present algorithms and results from Gkeyll, a full-f continuum, electromagnetic gyrokinetic code, designed to study turbulence in the edge region of fusion devices. The edge is computationally very challenging, requiring robust algorithms that can handle large-amplitude fluctuations and stable interactions with plasma sheaths. We present an energy-conserving high-order discontinuous Galerkin scheme that solves gyrokinetic equations in Hamiltonian form. Efficiency is improved by a careful choice of basis functions and automatically generated computation kernels. Previous verification tests were performed in the straight-field-line large plasma device [Shi et al., J. Plasma Phys. 83, 905830304 (2017)] and the Texas Helimak, a simple magnetized torus [Bernard et al., Phys. Plasmas 26, 042301 (2019)], including the effect of end-plate biasing on turbulence. Results for the scrape-off layer for NSTX parameters with a model helical magnetic geometry with bad curvature have been obtained [Shi et al., Phys. Plasmas 26, 012307 (2019)]. In this paper, we present algorithms for the two formulations of electromagnetic gyrokinetics: the Hamiltonian and the symplectic. We describe each formulation and show results of benchmark tests. Although our scheme works for the Hamiltonian formulation, the presence of spurious numerical modes for high-β and large $k^2_⊥ρ^2_s$ regimes shows that the symplectic formulation is more robust. We then review our recent algorithm for the symplectic formulation [Mandell et al., J. Plasma Phys. 86, 905860109 (2020)], along with example application of this new capability. Maintaining positivity of the distribution function can be challenging, and we describe a new and novel exponential recovery based algorithm to address this.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Large-Eddy Simulation of a wind turbine using a Filtered Actuator Line Model

When dealing with multirotor devices such as quadcopters or wind farms, the cost of blade-resolved large-eddy simulation (LES) becomes prohibitive. Combining LES with a family of lower-fidelity models, called actuator line models (ALMs), has grown in popularity in the past decade. ALM replaces full blade resolution with an array of actuator points or lines parameterized by aerodynamic lift/drag polar plots along the blades. Body forces computed based on these actuator points are then projected onto the LES flow mesh, mimicking the effect of rotating blades on the flow. However, the optimal projection radius and the associated LES grid size is often too restrictive for multirotor simulations. Recently, a new tip-correction-based filtered ALM (F-ALM) was proposed by Martinez-Tossas and Meneveau (2019), which allows coarser-than-optimal grids by avoiding the associated overprediction of thrust. In this work, F-ALM is implemented into a high-order, in-house LES code to simulate National Renewable Energy Laboratory Phase VI wind turbine. It is then followed by a comparison between the baseline ALM and the newly implemented F-ALM in terms of instantaneous and time-averaged flow fields and blade loads, revealing the advantages of F-ALM in preventing the overprediction of power on coarse grids. Finally, this encourages accurate and affordable simulations of multirotor devices in the future.

17 WIND ENERGY↗
DOE OSTI2020

A chemically stabilized sulfur cathode for lean electrolyte lithium sulfur batteries

Lithium sulfur batteries (LSBs) are promising next-generation rechargeable batteries due to the high gravimetric energy, low cost, abundance, nontoxicity, and high sustainability of sulfur. However, the dissolution of high-order polysulfide in electrolytes and low Coulombic efficiency of Li anode require excess electrolytes and Li metal, which significantly reduce the energy density of LSBs. Quasi-solid-state LSBs, where sulfur is encapsulated in the micropores of carbon matrix and sealed by solid electrolyte interphase, can operate under lean electrolyte conditions, but a low sulfur loading in carbon matrix (<40 wt %) and low sulfur unitization (<70%) still limit the energy density in a cell level. Here, we significantly increase the sulfur loading in carbon to 60 wt % and sulfur utilization to ∼87% by dispersing sulfur in an oxygen-rich dense carbon host at a molecular level through strong chemical interactions of C–S and O–S. In an all-fluorinated organic lean electrolyte, the C/S cathode experiences a solid-state lithiation/delithiation reaction after the formation of solid electrolyte interphase in the first deep lithiation, completely avoiding the shuttle reaction. The chemically stabilized C/S composite retains a high reversible capacity of 541 mAh⋅g −1 (based on the total weight of the C/S composite) for 200 cycles under lean electrolyte conditions, corresponding to a high energy density of 974 Wh⋅kg −1 . The superior electrochemical performance of the chemical bonding-stabilized C/S composite renders it a promising cathode material for high-energy and long-cycle-life LSBs.

25 ENERGY STORAGE↗
DOE OSTI2023

Real-space solution to the electronic structure problem for nearly a million electrons

We report a Kohn–Sham density functional theory calculation of a system with more than 200 000 atoms and 800 000 electrons using a real-space high-order finite-difference method to investigate the electronic structure of large spherical silicon nanoclusters. Our system of choice was a 20 nm large spherical nanocluster with 202 617 silicon atoms and 13 836 hydrogen atoms used to passivate the dangling surface bonds. To speed up the convergence of the eigenspace, we utilized Chebyshev-filtered subspace iteration, and for sparse matrix–vector multiplications, we used blockwise Hilbert space-filling curves, implemented in the PARSEC code. For this calculation, we also replaced our orthonormalization + Rayleigh–Ritz step with a generalized eigenvalue problem step. We utilized all of the 8192 nodes (458 752 processors) on the Frontera machine at the Texas Advanced Computing Center. We achieved two Chebyshev-filtered subspace iterations, yielding a good approximation of the electronic density of states. Our work pushes the limits on the capabilities of the current electronic structure solvers to nearly 106 electrons and demonstrates the potential of the real-space approach to efficiently parallelize large calculations on modern high-performance computing platforms.

Chemistry↗
DOE OSTI2020

Ultrashort XUV pulse absorption spectroscopy of partially oxidized cobalt nanoparticles

High-order harmonic generation (HHG) based transient extreme ultraviolet (XUV) absorption spectroscopy is an emerging technique to trace photoinduced charge carrier dynamics in condensed phase materials with femtosecond and even attosecond temporal resolution and elemental specificity. However, its application to nanoparticulate samples that are relevant, for example, for novel photocatalytic light harvesting concepts, has been limited. This is in part due to the challenge to produce residual-free samples on ultrathin, XUV-transparent substrates as well as a widespread understanding that sparsely distributed nanoparticles do not provide sufficient contrast for XUV absorption measurements. Here, we present static XUV absorption spectra of partially oxidized Co nanowire-structures with diameters of approximately 4.5 nm and lengths between 10 and 40 nm, recorded with an ultrashort pulse HHG light source. Nanoparticles are synthesized by the agglomeration of Co atoms inside superfluid helium droplets, followed by surface deposition and oxidation in ambient air. The method is uniquely suited for residual-free synthesis of transition metal nanowires and their deposition on ultrathin substrates. Analysis by high-resolution transmission electron microscopy reveals the formation of CoO nanowires with regions of unoxidized Co in their interior. The nanoparticle samples are investigated in an HHG-driven ultrafast XUV absorption setup. Despite the low surface coverage of only 23%, the recorded spectrum exhibits a distinct absorption feature at the Co M 2,3 (2p) edge near 60 eV with a peak height of about 40 mOD. Finally, the results support the feasibility of table-top ultrafast transient XUV absorption studies of photoinduced dynamics in transition metal oxide nanoparticles with sub-monolayer surface coverage.

74 ATOMIC AND MOLECULAR PHYSICS↗