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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Highly cascaded first-order fiber Bragg gratings in highly multimode optical fibers for distributed temperature sensing under harsh environment conditions

This study presents a pioneering technique for fabricating highly cascaded first-order fiber Bragg gratings (FBGs) using a femtosecond laser-assisted point-by-point inscription method in highly multimode optical fibers, specifically Sapphire crystalline fiber, and pure silica coreless fiber. Notably, it marks the first successful demonstration of a distributed array comprising 10 FBGs within highly multimode fibers. This achievement is facilitated by a high-power laser technique that yields larger reflectors characterized by a Gaussian intensity profile. These first-order FBGs offer various advantages, including enhanced reflectivity, reduced fabrication time, and simplified spectral characteristics, enhancing their accessibility for interpretation when contrasted with higher-order FBGs. In addition to that it encompasses a comprehensive analysis of the robustness and efficacy of these FBGs, with particular emphasis on their ability to endure extreme temperatures. These FBGs demonstrate an advantageous capability for localized multi-point temperature monitoring, reaching temperatures up to 1500°C with sapphire crystalline fiber and 1100°C with pure silica coreless fiber. This resilience makes them suitable for deployment in harsh environmental conditions. This innovative approach substantially broadens the potential applications of highly multimode optical fibers, particularly in the arena of sensing and communication, where challenges related to thermal gradients and harsh environments prevail. Furthermore, these groundbreaking first-order FBGs signify a substantial advancement in the realm of distributed temperature sensing, offering supreme capabilities for temperature monitoring and signal stability. As such, our work holds the promise of a substantial impact on industries and applications that demand unwavering reliability under extreme conditions.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

A Condensed Constrained Nonconforming Mortar-Based Approach for Preconditioning Finite Element Discretization Problems

This paper presents and studies an approach for constructing auxiliary space preconditioners for finite element problems using a constrained nonconforming reformulation that is based on a proposed modified version of the mortar method. The well-known mortar finite element discretization method is modified to admit a local structure, providing an element-by-element or subdomain-by-subdomain assembly property. This is achieved via the introduction of additional trace finite element spaces and degrees of freedom (unknowns) associated with the interfaces between adjacent elements or subdomains. The resulting nonconforming formulation and a reduced-via-static-condensation Schur complement form on the interfaces are used in the construction of auxiliary space preconditioners for a given conforming finite element discretization problem. Overall, the properties of these preconditioners are studied and their performance is illustrated on model second order scalar elliptic problems utilizing high order elements.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A note on higher-order and nonlinear limiting approaches for continuously bounds-preserving discontinuous Galerkin methods

In Dzanic (2024), a limiting approach for high-order discontinuous Galerkin schemes was introduced which allowed for imposing constraints on the solution continuously (i.e., everywhere within the element). While exact for linear constraint functionals, this approach only imposed a sufficient (but not the minimum necessary) amount of limiting for nonlinear constraint functionals. This short note shows how this limiting approach can be extended to allow exactness for general nonlinear quasiconcave constraint functionals through a nonlinear limiting procedure, reducing unnecessary numerical dissipation. Finally, some examples are shown for nonlinear pressure and entropy constraints in the compressible gas dynamics equations, where both analytic and iterative approaches are used.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

High harmonic generation in Haldane model quantum dots

We study theoretically the nonlinear electron dynamics of Haldane model quantum dots placed in the field of an ultrashort optical pulse. One of the tuning parameters of the Haldane model is the phase accumulated by an electron during its transfer between the next nearest neighbor sites. We study how this parameter affect the nonlinear electron dynamics and the generation of high optical harmonics. With increasing the phase from its zero value to 90°, the low-energy electron states in the conduction and valence bands become more localized near the edges of the quantum dot resulting in suppression of the band gap, enhancement of the dipole inter-band coupling, and strong suppression of the average low-energy electron density of the states. As a result, the nonlinear response of the electron system of a Haldane model quantum dot is the strongest at intermediate values of the phase, ≈ 30− 40°. At these values of the phase, the electron dynamics is irreversible and a few first high-order harmonics have the largest intensities. Here, when the phase approaches 90° value the generation of high harmonics is strongly suppressed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Tackling P3HT:Y‐Series Miscibility Through Advanced Processing for Tunable Aggregation

Abstract Polymer and small molecule blend thin films are of strong interest for organic electronics and particularly organic solar cells. The high miscibility in blends of ordinary P3HT and state‐of‐the‐art Y‐series non‐fullerene acceptors (NFAs) suppresses phase separation and aggregation challenging successful charge separation and transport. In a recent work, current‐induced doping (CID) is introduced, a method to precisely control the aggregation of Poly(3‐hexylthiophene) (P3HT) in solution. The highly ordered pre‐aggregation in solution is used here to control the P3HT aggregation in neat films and blends with Y12 (BTP‐4F‐12). This results in a 25‐fold increase in hole mobility in P3HT organic field‐effect transistor (OFET) devices and tunability of the P3HT aggregate quality in the presence of Y12 over large ranges. At the same time, particularly the Y12 long‐range ordering is heavily suppressed by increasing P3HT aggregation. However, solvent vapor annealing (SVA) leads to an extraordinarily high Y12 ordering, changes in the crystal orientation of Y12, and a further improvement of P3HT aggregation. A broad range of different degrees of aggregation of both materials can therefore be obtained in the final thin films solely by changing processing parameters without changing the composition of the material system.

Chemistry↗
DOE OSTI2025

Toward Higher-order Accuracy in Self-gravitating Hydrodynamics

High-order algorithms have emerged in numerical astrophysics as a promising avenue to reduce truncation error (proportional to a power of the linear resolution Δ x ) with only a moderate increase to computational expense. Significant effort has been placed in the development of finite-volume algorithms for (magneto)hydrodynamics; however, state-of-the-art astrophysical simulations tightly couple a plenitude of physics, additionally including gravity, photon transport, cosmic-ray transport, chemistry, and/or diffusion, to name a few. Algorithms frequently operator-split this additional physics (often a first-order error in time) and/or adopt a model wherein their evaluation is limited to second-order accuracy in space. In this work, we present a fourth-order-accurate finite-volume scheme for self-gravitating hydrodynamics on a uniform Cartesian grid. The method supplies source terms for the gravitational acceleration ( ρ g ) and gravitational energy release ( ρ v · g ) associated with fourth-order-accurate solutions to the Poisson equation. Our scheme (1) guarantees the conservation of total linear momentum while (2) decreasing (in proportion to Δ x 4 ) the effects of spurious heating and/or cooling associated with truncation error in the gravity. We demonstrate expected convergence rates for the algorithm by measuring errors in test problems evolving self-gravity modified linear waves and 3D polytropic equilibria. We test robustness of the algorithm by integrating an induced “inside-out” adiabatic collapse. We also discuss a method to smoothly downgrade the solution to second-order spatial accuracy to avoid spurious overshoots near steep density and/or pressure gradients.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A parallel p ‐adaptive discontinuous Galerkin method for the Euler equations with dynamic load‐balancing on tetrahedral grids

Abstract A novel p ‐adaptive discontinuous Galerkin (DG) method has been developed to solve the Euler equations on three‐dimensional tetrahedral grids. Hierarchical orthogonal basis functions are adopted for the DG spatial discretization while a third order TVD Runge‐Kutta method is used for the time integration. A vertex‐based limiter is applied to the numerical solution in order to eliminate oscillations in the high order method. An error indicator constructed from the solution of order and is used to adapt degrees of freedom in each computational element, which remarkably reduces the computational cost while still maintaining an accurate solution. The developed method is implemented with under the Charm++ parallel computing framework. Charm++ is a parallel computing framework that includes various load‐balancing strategies. Implementing the numerical solver under Charm++ system provides us with access to a suite of dynamic load balancing strategies. This can be efficiently used to alleviate the load imbalances created by p ‐adaptation. A number of numerical experiments are performed to demonstrate both the numerical accuracy and parallel performance of the developed p ‐adaptive DG method. It is observed that the unbalanced load distribution caused by the parallel p ‐adaptive DG method can be alleviated by the dynamic load balancing from Charm++ system. Due to this, high performance gain can be achieved. For the testcases studied in the current work, the parallel performance gain ranged from 1.5× to 3.7×. Therefore, the developed p ‐adaptive DG method can significantly reduce the total simulation time in comparison to the standard DG method without p ‐adaptation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Effect of Molecular Structure of Quinones and Carbon Electrode Surfaces on the Interfacial Electron Transfer Process

Quinones can undergo thermodynamically reversible proton-coupled electron transfer reactions and are being applied as electroactive compounds in aqueous organic batteries. However, the electrochemical reversibility of these compounds is affected not only by their molecular structure but also by the properties of a carbon-based electrode surface. This study combines experimental and theoretical approaches to understand this dependence. Specifically, we study the electron transfer kinetics of two synthesized quinone derivatives and two commercially available ones with a glassy carbon, a highly ordered pyrolytic graphite, and a high-edge-density graphite electrode (HEDGE). The electrochemical reversibility is notably improved on the HEDGE, which shows a higher density of defects and presents oxygenated functional groups at its surface. The electron transfer kinetics are controlled by adsorbed species onto the HEDGE. Molecular dynamics simulation and quantum mechanics calculations suggest defects with oxygen-containing functional groups, such as C-O and C$=$O, on HEDGE surfaces drive the interaction with the functional groups of the molecules, during physisorption from van der Waals forces. The presence of sulfonic acid side groups and a greater number of aromatic rings in the molecular structure may contribute to a higher stabilization of quinone derivatives on HEDGEs. We propose that high-performance carbon-based electrodes can be obtained without catalysts for organic batteries, by the engineering of carbon-based surfaces with edge-like defects and oxygenated functional groups.

25 ENERGY STORAGE↗
DOE OSTI2020

Auxiliary Space Preconditioning of Finite Element Equations Using a Nonconforming Interior Penalty Reformulation and Static Condensation

We modify the well-known interior penalty finite element discretization method so that it allows for element-by-element assembly. This is possible due to the introduction of additional unknowns associated with the interfaces between neighboring elements. The resulting bilinear form, and a Schur complement (reduced) version of it, are utilized in a number of auxiliary space preconditioners for the original conforming finite element discretization problem. Furthermore, these preconditioners are analyzed on the fine scale, and their performance is illustrated on model second order scalar elliptic problems discretized with high order elements.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Continuously bounds-preserving discontinuous Galerkin methods for hyperbolic conservation laws

For finite element approximations of transport phenomena, it is often necessary to apply a form of limiting to ensure that the discrete solution remains well-behaved and satisfies physical constraints. However, these limiting procedures are typically performed at discrete nodal locations, which is not sufficient to ensure the robustness of the scheme when the solution must be evaluated at arbitrary locations (e.g., for adaptive mesh refinement, remapping in arbitrary Lagrangian–Eulerian solvers, overset meshes, etc.). In this work, a novel limiting approach for discontinuous Galerkin methods is presented which ensures that the solution is continuously bounds-preserving (i.e., across the entire solution polynomial) for any arbitrary choice of basis, approximation order, and mesh element type. Through a modified formulation for the constraint functionals, the proposed approach requires only the solution of a single spatial scalar minimization problem per element for which a highly efficient numerical optimization procedure is presented. Here, the efficacy of this approach is shown in numerical experiments by enforcing continuous constraints in high-order unstructured discontinuous Galerkin discretizations of hyperbolic conservation laws, ranging from scalar transport with maximum principle preserving constraints to compressible gas dynamics with positivity-preserving constraints.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A short note on the accuracy of the discontinuous Galerkin method with reentrant faces

In this work, we study the convergence of the discontinuous Galerkin (DG) method applied to the advection–reaction equation on meshes with reentrant faces. On such meshes, the upwind numerical flux is not smooth, and so the numerical integration of the resulting face terms can only be expected to be first-order accurate. Despite this inexact integration, we prove that the DG method converges with order $\mathscr{O}$(h p+1/2 ), which is the same rate as in the case of exact integration. Consequently, specialized quadrature rules that accurately integrate the non-smooth numerical fluxes are not required for high-order accuracy. These results are numerically corroborated on examples of linear advection and discrete ordinates transport equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Real-fluid behavior in rapid compression machines: Does it matter?

Rapid compression machines (RCMs) have been extensively used to quantify fuel autoignition chemistry and validate chemical kinetic models at high-pressure conditions. Historically, the analyses of experimental and modeling RCM autoignition data have been conducted based on the adiabatic core hypothesis with ideal gas assumption, where real-fluid behavior has been completely overlooked, though this might be significant at common RCM test conditions. Here, this work presents a first-of-its-kind study that addresses two significant but overlooked questions for autoignition studies within RCMs in the fundamental combustion community: (i) experiment-wise, can unaccounted-for real-fluid behavior in RCMs affect the interpretation and analysis of RCM experimental data? and (ii) simulation-wise, can unaccounted-for real-fluid behavior in RCMs affect RCM autoignition modeling and the validation of chemical kinetic models? To this end, theories for real-fluid isentropic change are newly proposed and derived based on high-order Virial EoS, and are further incorporated into an effective-volume real-fluid autoignition modeling framework newly developed for RCMs. With detailed analyses, the strong real-fluid behavior in representative RCM tests is confirmed, which can greatly influence the interpretation of RCM autoignition experiments, particularly the determination of end-of-compression temperature and evolution of the adiabatic core in the reaction chamber. Furthermore, real-fluid RCM modeling results reveal that considerable error can be introduced into simulating RCM autoignition experiments when following the community-wide accepted effective-volume approach by assuming ideal-gas behavior, which can be as high as 64% in the simulated ignition delay time at compressed pressure of 125 bar and lead to contradictory validation results of chemical kinetic models. Therefore, we recommend the community to adopt frameworks with real-fluid behavior fully accounted for (e.g., the one developed in this study) to analyze and simulate past and future RCM experiments, so as to avoid misinterpretation of RCM autoignition experiments and eliminate the potential errors that can be introduced into the simulation results with the existing RCM modeling frameworks.

High-order Virial equation of state↗
DOE OSTI2023

Advanced Turbulence Models for Large-Scale Atmospheric Boundary Layer Flows

We present high-fidelity large-eddy-simulation (LES) modeling approaches for the turbulent atmospheric boundary layer (ABL) flows. Wind energy is a prime example of an application driven by ABL. Generation of electrical energy from farms of wind turbines at night in the stable ABL is a particularly interesting situation. In this report, we consider the well-known GEWEX (Global Energy and Water Cycle Experiment) Atmospheric Boundary Layer Study (GABLS) stably stratified benchmark LES case. We use a high-order spectral element code Nek5000/RS, which is supported under the DOE's Exascale Computing Project (ECP) Center for Efficient Exascale Discretizations (CEED) project, targeting application simulations on various acceleration-device based exascale computing platforms. In our earlier ANL report, we demonstrated our newly developed subgrid-scale (SGS) models based on high-pass filter (HPF), mean-field eddy viscosity (MFEV), and Smagorinsky (SMG) with no-slip and traction boundary conditions, provided with low-order statistics, convergence and turbulent structure analysis. In this report, we extend the range of our SGS modeling approaches in the context of the mean-field eddy viscosity (MFEV), to include the solution of an SGS turbulent kinetic energy equation (TKE). We demonstrate the model fidelity of Nek5000/RS in comparison to that of AMR-Wind, a block-structured second-order finite-volume code with adaptive-mesh-refinement capabilities, with which we studied scaling performance for both codes in comparison on DOE's leadership computing platforms.

17 WIND ENERGY↗
DOE OSTI2021

Anomalous ellipticity dependence of the generation of near-threshold harmonics in noble gases

We measured the ellipticity dependence of the harmonic yield of He, Ne, Ar, Kr, and Xe. It is found that the harmonic yield in general decreases monotonically with increasing laser ellipticity but anomalous dependence does occur for each atom for some harmonic orders. Additionally, we found that the degree of anomaly is weakest in helium but is stronger for the heavier atoms. The anomaly also depends on laser intensity. To explain these features, we developed a quantum trajectory Monte Carlo (QTMC) model based on the quantum path integral theory to study high-order harmonic generation (HHG). The model includes the effect of Coulomb potential from the target ion, thus it is capable of calculating near-threshold harmonics quantitatively. This model reveals that the presence of Coulomb potential would generate orbiting trajectories which are responsible for the anomalous ellipticity dependence near the threshold.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

A Field Enhancement Integration Design Featuring Misalignment Tolerance for Wireless EV Charging Using LCL Topology

This article proposes a magnetic integration design for EV wireless power transfer (WPT) systems, where the compensation and transmitting coils overlap one on top of the other to share the ferrite layer without any decoupling consideration. The magnetic field generated by both the compensation and transmitting coils are exploited to transfer power. To this end, a compensation method is proposed to enable magnetic field enhancement without any reactive power flow between the compensation and the transmitting coils, and to achieve input zero phase angle and constant current output. In addition, an efficient finite element analysis-based coil optimization algorithm is proposed to improve the coil misalignment tolerance on the horizontal plane, in which a reversely connected inner coil is used to stabilize the system output under misalignment conditions. Analytical and simulation results confirm transmission flux density enhancement and reduced leakage field characteristics of the proposed coil design. Finally, a scaled-down WPT prototype is built and tested to verify the performance and effectiveness of the design. The proposed design achieves 91.17% efficiency and misalignment tolerance up to 200 mm in any XY-direction while maintaining the power transfer and its efficiency. The results of this work provide insights into magnetic integration design of high-order compensation topologies featuring higher compactness, less ferrite usage, magnetic field enhancement, and misalignment tolerance for WPT systems.

42 ENGINEERING↗
DOE OSTI2024

Insight of Synthesis of Single Crystal Ni-Rich LiNi 1-x-y Co x Mn y O 2 Cathodes

Single-crystal Ni-rich LiNi 1-x-y Co x Mn y O 2 (NCM) cathodes have garnered widespread attention in the lithium-ion battery community due to their unique advantages in mechanical performance and their ability to minimize interfacial electrochemical side reactions. The synthesis of single-crystal materials with monodisperse and appropriate size, minimal lattice defects, and highly ordered structures is the key for high-performance batteries. However, achieving this goal poses challenges due to the lack of in-depth understanding regarding specific experimental parameters and the solid reaction mechanism during the synthesis process. Further, in this review, the aim is to provide an in-depth analysis of the critical process parameters involved in the synthesis and their impact on crystal morphology, structure, and electrochemical performance. Consequently, the first section focuses on the effect of the precursor morphology, lithium salt, atmosphere, and sintering procedure. In the second section, the study delves into an in-depth discussion of the solid reaction and crystal growth mechanism. Lastly, it is concluded by highlighting the prospects and challenges associated with the synthesis and application of single-crystal Ni-rich NCM cathodes.

25 ENERGY STORAGE↗
DOE OSTI2021

Importance of charge self-consistency in first-principles description of strongly correlated systems

First-principles approaches have been successful in solving many-body Hamiltonians for real materials to an extent when correlations are weak or moderate. As the electronic correlations become stronger often embedding methods based on first-principles approaches are used to better treat the correlations by solving a suitably chosen many-body Hamiltonian with a higher level theory. The success of such embedding theories, often referred to as second-principles, is commonly measured by the quality of self-energy Σ which is either a function of energy or momentum or both. However, Σ should, in principle, also modify the electronic eigenfunctions and thus change the real space charge distribution. While such practices are not prevalent, some works that use embedding techniques do take into account these effects. In such cases, choice of partitioning, of the parameters defining the correlated Hamiltonian, of double-counting corrections, and the adequacy of low-level Hamiltonian hosting the correlated subspace hinder a systematic and unambiguous understanding of such effects. Further, for a large variety of correlated systems, strong correlations are largely confined to the charge sector. Then an adequate nonlocal low-order theory is important, and the high-order local correlations embedding contributes become redundant. Here we study the impact of charge self-consistency within two example cases, TiSe 2 and CrBr 3 , and show how real space charge re-distribution due to correlation effects taken into account within a first-principles Green’s function-based many-body perturbative approach is key in driving qualitative changes to the final electronic structure of these materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Explicit Runge–Kutta Methods that Alleviate Order Reduction

Explicit Runge–Kutta (RK) methods are susceptible to a reduction in the observed order of convergence when applied to an initial boundary value problem with time-dependent boundary conditions. We study conditions on explicit RK methods that guarantee high order convergence for linear problems; we refer to these conditions as weak stage order conditions. We prove a general relationship between the method’s order, weak stage order, and number of stages. Furthermore, we derive explicit RK methods with high weak stage order and demonstrate, through numerical tests, that they avoid the order reduction phenomenon up to any order for linear problems and up to order three for nonlinear problems.

explicit Runge–Kutta↗