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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 433 records · Page 24

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Reinforcement Learning for Load-balanced Parallel Particle Tracing

We explore an online reinforcement learning (RL) paradigm to dynamically optimize parallel particle tracing performance in distributed-memory systems. Our method combines three novel components: (1) a work donation algorithm, (2) a high-order workload estimation model, and (3) a communication cost model. First, we design an RL-based work donation algorithm. Our algorithm monitors workloads of processes and creates RL agents to donate data blocks and particles from high-workload processes to low-workload processes to minimize program execution time. The agents learn the donation strategy on the fly based on reward and cost functions designed to consider processes' workload changes and data transfer costs of donation actions. Second, we propose a workload estimation model, helping RL agents estimate the workload distribution of processes in future computations. Third, we design a communication cost model that considers both block and particle data exchange costs, helping RL agents make effective decisions with minimized communication costs. We demonstrate that our algorithm adapts to different flow behaviors in large-scale fluid dynamics, ocean, and weather simulation data. Our algorithm improves parallel particle tracing performance in terms of parallel efficiency, load balance, and costs of I/O and communication for evaluations with up to 16,384 processors.

Distributed and parallel particle tracing↗

DOE OSTI2021

Situating the Vector Density Approach Among Contemporary Continuum Theories of Dislocation Dynamics

Abstract For the past century, dislocations have been understood to be the carriers of plastic deformation in crystalline solids. However, their collective behavior is still poorly understood. Progress in understanding the collective behavior of dislocations has primarily come in one of two modes: the simulation of systems of interacting discrete dislocations and the treatment of density measures of varying complexity that are considered as continuum fields. A summary of contemporary models of continuum dislocation dynamics is presented. These include, in order of complexity, the two-dimensional statistical theory of dislocations, the field dislocation mechanics treating the total Kröner–Nye tensor, vector density approaches that treat geometrically necessary dislocations on each slip system of a crystal, and high-order theories that examine the effect of dislocation curvature and distribution over orientation. Each of theories contain common themes, including statistical closure of the kinetic dislocation transport equations and treatment of dislocation reactions such as junction formation. An emphasis is placed on how these common themes rely on closure relations obtained by analysis of discrete dislocation dynamics experiments. The outlook of these various continuum theories of dislocation motion is then discussed.

Engineering↗

DOE OSTI2023

Turbulence Modeling with Nek5000/RS, SOD2D and Alya

We present Validation and Verification (V & V) study of two high-order spectral element method (SEM) based Computational Fluid Dyanimcs (CFD) codes that we will utilize for work on turbulence modeling: Nek5000 and SOD2D. While the former solves the incompressible form of Navier-Stokes equation, the latter works with compressible set of equations and uses an entropy-viscosity formulation to account for the discontinuities for high Mach number flows. We demonstrate the accuracy of these codes for two benchmark sub-sonic turbulent flows: periodic channel and pipe flow, by carrying out first and second order statistical analysis, grid convergence and turbulent structure analysis using wall-resolved large eddy simulations (WRLES). Later, we report on the implementation and testing of various wall-modeling strategies for large eddy simulation of turbulent flows in Nek5000. These include both classical log-law based and decision-tree based machine-learning models. Accuracy of these closure strategies are analyzed and necessary future work is outlined.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Nanostructure‐Derived Antireflectivity in Leafhopper Brochosomes

Understanding how insect‐derived biomaterials interact with light has led to new advances and interdisciplinary insights in entomology and physics. Leafhoppers are insects that coat themselves with highly ordered biological nanostructures known as brochosomes. Brochosomes are thought to provide a range of protective properties to leafhoppers, such as hydrophobicity and antireflectivity, which has inspired the development of synthetic brochosomes that mimic their structures. Despite recent progress, the high antireflective properties of brochosome structures are not fully understood. Herein, a combination of experiments and computational modeling is used to understand the structure‐, material‐, and polarization‐dependent optical properties of brochosomes modeled on the geometries found in three leafhopper species. The results qualitatively represent that light interference interaction with nanostructures naturally occurring in brochosomes is responsible for the spectral tuning and the asymmetric line shape of the reflectance spectra. Whereas prior work has focused on the computational modeling of idealized pitted particles, this work shows that light–matter interactions with brochosome structures can be tuned by varying the geometry of their cage‐like nanoscale features and by changing the arrangement of multiparticle assemblies. Broadly, this work establishes principles for the guided design of new optically active materials inspired by these unique insect nanostructures.

Banerjee, Progna↗

DOE OSTI2021

Dynamic Phasor Modeling of Multi-Converter Systems

Generalized form of Dynamic Phasor (DP)-based modeling of multi-converter systems containing high order harmonics has not been proposed in the literature yet due to the complexity and large number of variables that should be included in the model. In this work a generalized form of a system-level model comprising high number of single-phase converters is presented. The proposed method can be applied to any number of single-phase voltage source converters (VSIs). Test cases of this modeling method containing 5 to 100 parallel connected single-phase voltage source inverters are modeled and simulation times are recorded. The results of the proposed method are compared and validated with conventional average models as well as detailed switching models. A systematic approach for comparing the accuracy and timestep between dynamic phasor modeling method and detailed switching model is illustrated. Advantage of DP models over conventional average models for stability assessment are also discussed at the end of the paper.

Xue, Yaosuo↗

DOE OSTI2025

Revealing and Engineering Assembly Pathways of 3D DNA Origami Crystals

Recent developments in nanomaterial self-assembly demonstrate the capability to create tailored nanostructures by engineering both the binding coordination and specificity of interactions between material subunits. DNA origami frames allow for the design and fabrication of a broad variety of ordered 3D nanoscale architectures through self-assembly, facilitated by frame-to-frame bonds with designable strength and specificity. While the bond design is critical to lattice formation, the assembly process itself is often dependent on a thermal pathway. Highly ordered nanoscale frameworks, assembled from DNA frames, are predominantly crystallized through thermal annealing pathways that typically follow a “slow” cooling approach, with experiments on the time scale of days yielding DNA origami crystals in the range of 1−10 μm. This extended assembly time scale hinders the study of crystal formation pathways, necessitating a deeper understanding of factors governing successful annealing. Lack of insight into time scale also presents a practical limitation for material fabrication. Here, we investigate key factors affecting lattice assembly pathways and demonstrate that precise engineering of assembly conditions greatly reduces assembly times by up to nearly 2 orders of magnitude. We evaluate the nucleation and growth of crystals via optical and electron microscopy, and small-angle X-ray scattering techniques, mapping the time−temperature-transformation of superlattices from the melt through single-crystal optical tracking. The results show that origami frame assembly can be described by classical nucleation and growth theory, which can, in turn, be used to prescribe the growth of the crystals. Lastly, these findings are applied to demonstrate thermal pathway-dependent assembly, forming distinct assemblies based on different thermal annealing profiles.

36 MATERIALS SCIENCE↗

DOE OSTI2026

N -Terminal Octylated Peptoid Hydrogels as 3D-Printable Cell Scaffolds and Proteolytically Robust Cargo Depots

Supramolecular hydrogels that mimic the extracellular matrix (ECM) represent promising materials for tissue engineering and drug delivery. However, conventional hydrogels formed via the self-assembly of natural or synthetic building blocks often face a trade-off between biological functionality and biochemical stability, limiting their utility in long-term or protease-rich environments. Peptoids, a class of peptide-inspired, sequence-defined polymers, offer a compelling alternative due to their exceptional proteolytic resistance and bioactivity. Despite this potential, the development of supramolecular peptoid hydrogels has been hindered by the absence of backbone hydrogen bond donors, which limits long-range ordering necessary for efficient hydrogel formation. This work describes a short peptoid functionalized at the N-terminus with an octyl chain that readily self-assembles into hydrogels. Hydrophobic interactions among pendant octyl groups promote directional peptoid packing into highly ordered nanosheets, which interconnect to form a porous hydrogel network. These hydrogels exhibit tunable viscoelasticity, shear-thinning, and self-healing properties, enabling their use as inks for extrusion-based 3D printing. They support NIH-3T3 fibroblast adhesion, spreading, and proliferation, maintaining greater than 95% cell viability over 4 days. Moreover, the hydrogels retain their macroscopic integrity under protease-rich conditions, enabling sustained cargo release and uniform cellular uptake. Together, this study demonstrates a class of supramolecular peptoid hydrogelators that integrate biocompatibility, 3D printability, and proteolytic stability, providing a versatile platform for ECMmimetic scaffolds in regenerative medicine and long-term therapeutic delivery.

cargo delivery↗

DOE OSTI2024

Role of gas composition in weakened nonlinear standing wave excitation and improved plasma radial uniformity in very-high-frequency asymmetric capacitive Ar/CF 4 discharges

Abstract The higher harmonics generated by nonlinear sheath motion would enhance the standing wave effect, and thus lead to center-peaked plasma density profile in very-high-frequency (VHF) capacitive discharges. In this work, a nonlinear transmission line (NTL) model introduced in Zhouet al(2021Plasma Sources Sci. Technol.30125017) has been extended, with radial transport of various particles and nonlinear sheath motion into account, to investigate the effects of CF 4 fraction α on the nonlinear standing wave excitation and plasma radial uniformity in VHF (60 MHz) capacitively coupled Ar/CF 4 plasmas at 3 Pa. The results indicate that for pure Ar discharges (i.e. α = 0%), the nonlinearly excited harmonics with short wavelength significantly enhance the electron power absorption at the radial center, resulting in a pronounced central-high plasma density profile. As α increases, the high-order harmonics are gradually damped due to the increase of resistance, as well as the longer wavelength caused by thicker sheath thickness. Thus, the radial profile of the electron absorbed power density shifts from center-peak to edge-high. Besides, at the radial center, the electron density and Ar + ion density decrease with α , CF 3 + ion density shows an increasing trend, while F − ion density initially rises and then decreases. Moreover, the density profiles of all the species become more uniform at higher CF 4 fraction, due to the suppressed nonlinear standing wave excitation and the longer wavelength of the nonlinear harmonics.

Physics↗

DOE OSTI2021

Magnetic probe diagnostics of nonlinear standing waves and bulk ohmic electron power absorption in capacitive discharges

It is recognized that standing wave effects appearing in large-area, very-high-frequency capacitively coupled plasma (CCP) reactors cause center-high plasma non-uniformity. Using a high-frequency magnetic probe, we present a direct experimental diagnostic of the nonlinear standing waves and bulk ohmic electron power absorption dynamics in low pressure CCP discharges for different driving frequencies of 13.56, 30, and 60 MHz. The design, principle, calibration, and validation of the probe are described in detail. Spatial structures of the harmonics of the magnetic field, determined by the magnetic probe, were used to calculate the distributions of the harmonic current and the corresponding ohmic electron power deposition, providing insights into the behavior of nonlinear harmonics. At a low driving frequency, i.e. 13.56 MHz, no remarkable nonlinear standing waves were identified and the bulk ohmic electron power absorption was observed to be negligible. Here, the harmonic magnetic field/current was found to increase dramatically with the driving frequency, due to decreased sheath reactance and more remarkable nonlinear standing waves at a higher driving frequency, leading to the enhancements of the ohmic heating and the plasma density in the bulk, specifically at the electrode center. At a high driving frequency, i.e. 60 MHz, the high-order harmonic current density and the corresponding ohmic electron power absorption exhibited a similar node structure, with the main peak on axis, and one or more minor peaks between the electrode center and the edge, contributing to the center-high profile of the plasma density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Spectrally Stabilized Interface Capturing Formulation and Implementation in Nek5000/NekRS

This report documents the formulation of a novel level-set method for incompressible two-phase flows in the continuous Galerkin (CG) high order spectral element framework. The overall method hinges on a novel implementation of the spectral vanishing viscosity (SVV) operator for the stabilization of linear/non-linear hyperbolic problems. The multidimensional SVV convolution kernels, which in essence, have a similar effect as a high pass filter applied to the derivatives, are formulated by exploiting the tensor product form, analogous to the construction of the usual stiffness matrix system. The resulting kernels are directionally decoupled and ensure a linear, symmetric positive definite, elliptic matrix operator. The SVV formulation is demonstrated to provide a robust stabilizing mechanism through challenging linear and non-linear hyperbolic problems, including problems pertinent to the level-set formulation. The two-phase framework conceptualized herein is based on the conservative level-set (CLS) method which represents the interface between the fluids by the 0.5 iso-contour of the smoothed Heaviside function. The CLS method is augmented with a preconditioning procedure for interface normals using the signed distance function which precludes the manifestation of spurious oscillations in the vicinty of the interface. Further, the existing mixed explicit-implicit approach for the solution of Navier-Stokes equations in Nek5000, as described in Tomboulides et al, is augmented with a pressure coefficient splitting approach for the Poisson equation, which greatly accelerated the convergence of pressure solver for two-phase systems with large density ratio. The robustness and accuracy of the overall two-phase method is demonstrated through canonical challenging problems involving high density and viscosity ratios, with and without surface tension. The two-phase formulation is wholly implemented in Nek5000 and the SVV stabilization method is implemented in NekRS, which is the essential precursor to the two-phase framework, undergoing active development.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

High-Mobility Hole Transport in Single-Grain PbSe Quantum Dot Superlattice Transistors

Epitaxially-fused superlattices of colloidal quantum dots (QD epi-SLs) may exhibit electronic minibands and high-mobility charge transport, but electrical measurements of epi-SLs have been limited to large-area, polycrystalline samples in which superlattice grain boundaries and intragrain defects suppress/obscure miniband effects. Systematic measurements of charge transport in individual, highly-ordered epi-SL grains would facilitate the study of minibands in QD films. Here, we demonstrate the air-free fabrication of microscale field-effect transistors (μ-FETs) with channels consisting of single PbSe QD epi-SL grains (2–7 μm channel dimensions) and analyze charge transport in these single-grain devices. The eight devices studied show p-channel or ambipolar transport with a hole mobility as high as 3.5 cm 2 V –1 s –1 at 290 K and 6.5 cm 2 V –1 s –1 at 170–220 K, one order of magnitude larger than that of previous QD solids. The mobility peaks at 150–220 K, but device hysteresis at higher temperatures makes the true mobility–temperature curve uncertain and evidence for miniband transport inconclusive.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Crossbreeding Effect of Chalcogenation and Iodination on Benzene Additives Enables Optimized Morphology and 19.68% Efficiency of Organic Solar Cells

Volatile solid additives have attracted increasing attention in optimizing the morphology and improving the performance of currently dominated non-fullerene acceptor-based organic solar cells (OSCs). However, the underlying principles governing the rational design of volatile solid additives remain elusive. Herein, a series of efficient volatile solid additives are successfully developed by the crossbreeding effect of chalcogenation and iodination for optimizing the morphology and improving the photovoltaic performances of OSCs. Five benzene derivatives of 1,4-dimethoxybenzene (DOB), 1-iodo-4-methoxybenzene (OIB), 1-iodo-4-methylthiobenzene (SIB), 1,4-dimethylthiobenzene (DSB) and 1,4-diiodobenzene (DIB) are systematically studied, where the widely used DIB is used as the reference. The effect of chalcogenation and iodination on the overall property is comprehensively investigated, which indicates that the versatile functional groups provided various types of noncovalent interactions with the host materials for modulating the morphology. Among them, SIB with the combination of sulphuration and iodination enabled more appropriate interactions with the host blend, giving rise to a highly ordered molecular packing and more favorable morphology. As a result, the binary OSCs based on PM6:L8-BO and PBTz-F:L8-BO as well as the ternary OSCs based on PBTz-F:PM6:L8-BO achieved impressive high PCEs of 18.87%, 18.81% and 19.68%, respectively, which are among the highest values for OSCs.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Consistent Second Moment Methods with Scalable Linear Solvers for Radiation Transport

Second moment methods (SMMs) are developed that are consistent with the discontinuous Galerkin spatial discretization of the discrete ordinates (or S$_N$) transport equations. The low-order (LO) diffusion system of equations is discretized with fully consistent P$_1$, local discontinuous Galerkin (LDG), and interior penalty (IP) methods. A discrete residual approach is used to derive SMM correction terms that make each of the LO systems consistent with the high-order discretization. We show that the consistent methods are more accurate and have better solution quality than independently discretized LO systems, that they preserve the diffusion limit, and that the LDG and IP consistent SMMs can be scalably solved in parallel on a challenging, multimaterial benchmark problem.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

An adaptive scalable fully implicit algorithm based on stabilized finite element for reduced visco-resistive MHD

The magnetohydrodynamics (MHD) equations are continuum models used in the study of a wide range of plasma physics systems, including the evolution of complex plasma dynamics in tokamak disruptions. However, efficient numerical solution methods for MHD are extremely challenging due to disparate time and length scales, strong hyperbolic phenomena, and nonlinearity. Additionally, therefore the development of scalable, implicit MHD algorithms and high-resolution adaptive mesh refinement strategies is of considerable importance. In this work, we develop a high-order stabilized finite-element algorithm for the reduced visco-resistive MHD equations based on the MFEM finite element library (mfem.org). The scheme is fully implicit, solved with the Jacobian-free Newton-Krylov (JFNK) method with a physics-based preconditioning strategy. Our preconditioning strategy is a generalization of the physics-based preconditioning methods in Chacón et al. (2002) to adaptive, stabilized finite elements. Algebraic multigrid methods are used to invert sub-block operators to achieve scalability. A parallel adaptive mesh refinement scheme with dynamic load-balancing is implemented to efficiently resolve the multi-scale spatial features of the system. Our implementation uses the MFEM framework, which provides arbitrary-order polynomials and flexible adaptive conforming and non-conforming meshes capabilities. Results demonstrate the accuracy, efficiency, and scalability of the implicit scheme in the presence of large scale disparity. The potential of the AMR approach is demonstrated on an island coalescence problem in the high Lundquist-number regime (≥ 10 7 ) with the successful resolution of plasmoid instabilities and thin current sheets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Lipid vesicles induced ordered nanoassemblies of Janus nanoparticles

Since many advanced applications require specific assemblies of nanoparticles (NPs), considerable efforts have been made to fabricate nanoassemblies with specific geometries. Although nanoassemblies can be fabricated through top-down approaches, recent advances show that intricate nanoassemblies can also be obtained through self-assembly, mediated for example by DNA strands. Here, we show, through extensive molecular dynamics simulations, that highly ordered self-assemblies of NPs can be mediated by their adhesion to lipid vesicles (LVs). Specifically, Janus NPs are considered so that the amount by which they are wrapped by the LV is controlled. The specific geometry of the nanoassembly is the result of effective curvature-mediated repulsion between the NPs and the number of NPs adhering to the LV. Further, the NPs are arranged on the LV into polyhedra which satisfy the upper limit of Euler's polyhedral formula, including several deltahedra and three Platonic solids, corresponding to the tetrahedron, octahedron, and icosahedron.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

A multi-terawatt two-color beam for high-power field-controlled nonlinear optics

Two-color laser beams are instrumental in light-field control and enhancement of high-order harmonic, spectral supercontinuum, and terahertz radiation generated in gases, plasmas, and solids. Here we demonstrate a multi-terawatt two-color beam produced using a relativistic plasma mirror, with 110 mJ at 800 nm and 30 mJ at 400 nm. Both color components have high spatial quality and can be simultaneously focused, provided that the plasma mirror lies within a Rayleigh range of the driving fundamental beam. Favorable scaling of second-harmonic generation by plasma mirrors at relativistic intensities suggests them as an excellent tool for multi-color waveform synthesis beyond the petawatt level.

42 ENGINEERING↗

DOE OSTI2023

Spectral phase singularity and topological behavior in perfect absorption

Perfect absorbers, which can achieve total absorption of all incoming energy, have been extensively studied in the last decades for various important technologies in general wave systems. Here, we show that perfect absorption (PA) is generically associated with topological spectral phase singularity (SPS), carrying conserved and tunable quantized topological invariants in spectral space. The order of topological invariant $\tilde{v}$ depends on the number of degenerate outgoing channels. Two commonly studied absorbers, mirror-backed and all-dielectric structures, are reexamined from a topological perspective to reveal the generation, evolution, and annihilation of SPSs with $\tilde{v}$ = ±1 or even high orders (i.e., $\tilde{v}$ = ±2,±4). A strategy based on charge conservation of SPSs to design dual-band perfect absorbers has been established. Our findings establish the topological origin of the robust existence of PA. More broadly, in this letter, we highlight topology as the fundamental property of conventional scattering behaviors. Finally, this insight could lead to opportunities in applications such as biosensing, topological metasurfaces, and micro/nano thermal radiation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Progress with the 5D full-F continuum gyrokinetic code COGENT

COGENT is an Eulerian gyrokinetic code being developed for edge plasma modelling. The code is distinguished by the use of a high-order finite-volume (conservative) discretization combined with mapped multi-block grid technology. Our recent work is focused on the development of a 5D full-F COGENT version.Anumerical algorithm utilizing locally a field-aligned multi-block coordinate system is implemented to facilitate simulations of highly anisotropic microturbulence in the presence of a strong magnetic shear. In this approach, the toroidal direction is divided into blocks such that,within each block, the cells are field-aligned and a non-matching (non-conformal) grid interface is allowed at the block boundaries. In this paper we report on details of the numerical implementation and present preliminary results of verification studies performed for the case of the ion temperature gradient (ITG) instability in a sheared toroidal annulus geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗