Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical stability”

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 OSTI2022

Glycan shield of the ebolavirus envelope glycoprotein GP

The envelope glycoprotein GP of the ebolaviruses is essential for host cell entry and the primary target of the host antibody response. GP is heavily glycosylated with up to 17 N-linked sites, numerous O-linked glycans in its disordered mucin-like domain (MLD), and three predicted C-linked mannosylation sites. Glycosylation is important for host cell attachment, GP stability and fusion activity, and shielding from neutralization by serum antibodies. Here, we use glycoproteomics to profile the site-specific glycosylation patterns of ebolavirus GP. We detect up to 16 unique O-linked glycosylation sites in the MLD, and two O-linked sites in the receptor-binding GP1 subunit. Multiple O-linked glycans are observed within N-linked glycosylation sequons, suggesting crosstalk between the two types of modifications. We confirmed C-mannosylation of W288 in full-length trimeric GP. We find complex glycosylation at the majority of N-linked sites, while the conserved sites N257 and especially N563 are enriched in unprocessed glycans, suggesting a role in host-cell attachment via DC-SIGN/L-SIGN. Our findings illustrate how N-, O-, and C-linked glycans together build the heterogeneous glycan shield of GP, guiding future immunological studies and functional interpretation of ebolavirus GP-antibody interactions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

A Game-Theoretic Quantum Algorithm for Solving Magic Squares

Variational quantum algorithms (VQAs) offer a promising near-term approach to finding optimal quantum strategies for playing non-local games. These games test quantum correlations beyond classical limits and enable entanglement verification. In this work, we present a variational framework for the Magic Square Game (MSG), a two-player non-local game with perfect quantum advantage. We construct a value Hamiltonian that encodes the game’s parity and consistency constraints, then optimize parameterize quantum circuits to minimize this cost. Our approach build on the stabilizer formalism, leverages commutation structure for circuit design, and is hardware-efficient. Compared to existing work, our contribution emphasizes algebraic structure an interpretability. We validate our method through numerical experiments and outline generalizations to larger games.

Chehade, Sarah [ORNL]↗
DOE OSTI2026

Hydrodynamic Modelling for Tidal Energy Kites: Cooperative Research and Development Final Report, CRADA Number CRD-24-30123

The primary goal of the work is to improve accuracy of hydrodynamic loads acting on the lifting surface of the kite turbine due to complex air-water-lifting body interaction in a low-order numerical model using high-fidelity numerical input. The result of this study will help the Participant better design control algorithms for the kite turbine so that the turbine can harvest more energy while maintaining stability of the platform as well as its reliability. The use of hydrodynamic coefficients/load components from computational fluid dynamic simulation will inform better inputs for a medium-fidelity dynamic model (specifically an OrcaFlex model) which is intended for use in rapid design iterations in the design process.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

Challenges in Continuous In-Field Critical Current Testing of High-Temperature Superconducting Tapes: Thermal and Mechanical Perspectives

High-temperature superconductors (HTS) are essential for ultra-high-field applications requiring exceptional current-carrying capacity under extreme conditions. However, systematic characterization of critical current in long-length conductors remains challenging due to complex thermal, electromag netic, and mechanical interactions during continuous testing. This study reports the development of a continuous in-field magnetization testing system for position-dependent critical current measurement in HTS tapes at 20 K under 7.5 T fields applied normal to the tape plane, enabling identification of performance-limiting regions that could compromise magnet stability. Here, the system addresses two fundamental challenges inherent to cryogenic reel to-reel testing. First, thermal management requires continuous cooling of a moving conductor to 20 K, achieved through liquid nitrogen precooling combined with a 100 W@20 K Gifford McMahon cryocooler. Second, screening currents in high fields generate Lorentz forces that induce twisting, bowing, and potential delamination. To mitigate these risks, we propose mechanical reinforcement and active current density suppression strategies. Numerical simulations using the stream function formulation reveal four primary failure modes: frictional heating at guide interfaces, unstable equilibria causing deformation, transverse current-induced stresses at guide transitions, and unsupported forces in vertical spans. Our mitigation strategies include PTFE coated guides to minimize friction, spring-loaded stabilization mechanisms to maintain tape alignment, controlled pre-heating using the liquid nitrogen thermal jacket to suppress critical current at stress points, and optimized guide positioning to minimize force accumulation. The experimental system is nearing completion, with testing planned to commence within two months. Preliminary validation at 65 K under 0.5 T demonstrates strong correlation between simulation-predicted mechanical instabilities and observed critical current variations during conductor tran sitions through the measurement region. These findings establish a robust foundation for quality assurance protocols essential to next-generation superconducting magnet applications.

Chen, Siwei [Princeton Plasma Physics Laboratory (↗
DOE OSTI2021

Root Influences on Mobilization and Export of Mineral-bound Soil Organic Matter

Biogeochemical cycles within mountainous watersheds are key regulators of ecosystem carbon storage and downstream nutrient loadings, and they have shown to be particularly vulnerable to climate change impacts. Increasing temperature and persistent droughts have already dramatically changed vegetation cover across the mountainous western US, with unknown consequences for soil carbon and nutrient cycles belowground. What remains elusive is to what extent associated changes in root-soil interactions may mobilize the vast pool of organic matter (OM) that has been stabilized by associations with minerals for centuries or millennia. Although plant root-driven OM mobilization from minerals may be a central control on carbon loss and nutrient export, such mechanisms are currently missing from conceptual and numerical models. The overall objective of this Exploratory Project is to identify the biogeochemical mechanisms by which roots destabilize mineral-associated organic matter (MAOM) and the cumulative impact on carbon and nutrient fate. To accomplish this goal, we integrated well-controlled laboratory experiments with in-field measurements and a scalable modeling approach. First, we conducted model system experiments to assess the vulnerability of MAOM to exudate-mediated mechanisms. Our results show that common root exudates effectively destabilize MAOM not only through direct, ligand-driven mobilization mechanisms, but also indirect, microbially-mediated mechanisms relying on secondary metabolites and enzymes. We further found that OM bound to poorly crystalline Fe and Al (hydr)oxides is more vulnerable to exudate-induced destabilization than OM bound to more crystalline phases, particularly in response to direct, ligand-promoted mechanisms. These findings demonstrate that the stability of MAOM is not just a function of their inherent properties, but also will depend in large parts on the ability of plant roots and microbes to produce exudates capable of triggering suitable mobilization mechanisms. We further employed a well-controlled rhizobox approach, combining advanced microsensor and mass spectrometry techniques, to resolve spatiotemporal variations in the composition and availability of exudates along single growing roots of grasses. Our results show that the composition of functionally relevant exudate compounds varies at extremely short time scales, seemingly shifting from ligands such as aromatic acids around root tips to less reactive metabolites such as amino acids around mature root segments. These results suggest a prevalence of direct MAOM mobilization mechanisms around the root tip, while indirect MAOM mobilization strategies may dominate around more mature root segments. Additionally, we utilized microsensor measurements in rhizoboxes to parameterize a rhizosphere (hydro)biogeochemistry reactive transport model (eSTOMP-ROOTS). The resulting model was used to assess how root exudations affects the stability of MAOM. Model simulations show that diel pulses of root exudation are strong enough to cause oscillations in biogeochemical conditions—particularly in pH, oxygen concentrations, and microbial activity— that repeatedly disrupt MAOM. Finally, we assessed field-scale impacts of root-driven MAOM destabilization over short and long time scales. A combination of in-field incubations and field-based characterizations showed that MAOM may not only be vulnerable to seasonal variations in root activity, but also to root-driven weathering over pedogenic time scales. In sum, results from our Exploratory Award highlight the strong control plant roots exert on MAOM mobilization and, thus, on the potential for carbon and nutrient export from watersheds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Chapter 10: Crystallography of Metabolic Enzymes

The metabolic enzymes like any enzymes generally display globular architecture where secondary structure elements and interactions between them preserve the spatial organization of the protein. A typical enzyme features a well-defined active site, designed for selective binding of the reaction substrate and facilitating a chemical reaction converting the substrate into a product. While many chemical reactions could be facilitated using only the functional groups that are found in proteins, the large percentage or intracellular reactions require use of cofactors, varying from single metal ions to relatively large molecules like numerous coenzymes, nucleotides and their derivatives, dinucleotides or hemes. Quite often these large cofactors become important not only for the catalytic function of the enzyme but also for the structural stability of it, as those are buried deep in the enzyme.

BASIC BIOLOGICAL SCIENCES,BIOMASS FUELS↗
DOE OSTI2022

A Package of Momentum and Heat Transfer Coefficients for the Stable Surface Layer Extended by New Coefficients over Sea Ice

The bulk transfer coefficients of momentum, heat, and humidity belong to the main ingredients of numerical weather prediction and climate models. They are needed for the calculation of turbulent fluxes in the surface layer and often rely on the Monin–Obukhov similarity theory requiring universal stability functions. The problem of a derivation of transfer coefficients based on different stability functions has been considered by many researchers over the years but it remains to this day. In this work, dedicated to the memory of S.S. Zilitinkevich, we also address this task, and obtain transfer coefficients from three pairs of theoretically derived stability functions suggested by Zilitinkevich and co-authors for stable conditions. Additionally, we construct non-iterative parametrizations of these transfer coefficients based on earlier work. Results are compared with state-of-the-art coefficients for land, ocean, and sea ice. The combined parametrizations form a package in a universal framework relying on a semi-analytical solution of the Monin-Obukhov similarity theory equations. A comparison with data of the Surface Heat Budget of the Arctic Ocean campaign (SHEBA) over sea ice reveals large differences between the coefficients for land conditions and the measurements over sea ice. However, two schemes of Zilitinkevich and co-authors show, after slight modification, good agreement with SHEBA although they had not been especially developed for sea ice. One pair of the modified transfer coefficients is superior and is compatible to earlier SHEBA-based parametrizations. Finally, an algorithm for practical use of all transfer coefficients in climate models is given.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Mathematical solutions in internal dose assessment: A comparison of Python-based differential equation solvers in biokinetic modeling

Abstract In biokinetic modeling systems employed for radiation protection, biological retention and excretion have been modeled as a series of discretized compartments representing the organs and tissues of the human body. Fractional retention and excretion in these organ and tissue systems have been mathematically governed by a series of coupled first-order ordinary differential equations (ODEs). The coupled ODE systems comprising the biokinetic models are usually stiff due to the severe difference between rapid and slow transfers between compartments. In this study, the capabilities of solving a complex coupled system of ODEs for biokinetic modeling were evaluated by comparing different Python programming language solvers and solving methods with the motivation of establishing a framework that enables multi-level analysis. The stability of the solvers was analyzed to select the best performers for solving the biokinetic problems. A Python-based linear algebraic method was also explored to examine how the numerical methods deviated from an analytical or semi-analytical method. Results demonstrated that customized implicit methods resulted in an enhanced stable solution for the inhaled 60 Co (Type M) and 131 I (Type F) exposure scenarios for the inhalation pathway of the International Commission on Radiological Protection (ICRP) Publication 130 Human Respiratory Tract Model (HRTM). The customized implementation of the Python-based implicit solvers resulted in approximately consistent solutions with the Python-based matrix exponential method ( expm ). The differences generally observed between the implicit solvers and expm are attributable to numerical precision and the order of numerical approximation of the numerical solvers. This study provides the first analysis of a list of Python ODE solvers and methods by comparing their usage for solving biokinetic models using the ICRP Publication 130 HRTM and provides a framework for the selection of the most appropriate ODE solvers and methods in Python language to implement for modeling the distribution of internal radioactivity.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2023

Full Technical Report: High-Rate Carbide Growth by Evaporation

Boron carbide (B 4 C) is an attractive inertial confined fusion (ICF) ablator and a unique ultra-hard material with numerous current and potential applications. What makes B 4 C attractive as an ablator is its relatively high density and the ability to form a stable glassy phase combined with excellent chemical and mechanical stability. Despite these advantages, the fabrication of B 4 C ablators has remained a challenge. Magnetron sputtering, our current frontrunner approach for B 4 C deposition, exhibits low deposition rates (< 2 μm/h) and undesirable nodular defects within the film structure that are believed to generate from the dusty plasma around the sputtering source. The goals of this project were to (1) establish alternative plasma-free, electron beam evaporation (EBE) capabilities for the deposition of thin films that are needed in mission-critical programs and; (2) demonstrate the feasibility to deposit B 4 C films at high rates by EBE in a bottom-up geometry and characterize their stoichiometry. In this feasibility study, we initially developed the capability to convert a general-use vacuum chamber into an electron beam evaporator as needed, or on demand. Then, we successfully used this setup to deposit films made of boron carbide and various other materials. By adjusting the deposition parameters and the evaporation source material, we were able to demonstrate the feasibility of reaching deposition rates as high as 8 μm/h, which is about 4 times the maximum rates obtained through sputtering methods. We predict that these values can still be further improved by implementing enhanced thermal control measures for the target source and optimizing electron beam parameters (ie. current, sweeping and frequency) since ceramic materials with similar densities to B 4 C can exhibit deposition rates of about twice this value. We anticipate that the skills and capabilities developed throughout this project will prove relevant not only to the HED and ICF campaigns but also to numerous existing and prospective applications focused on the development of ultra-hard carbide coatings. The fresh set of tools here developed will also facilitate the preparation and research on future new target materials of interest beyond B 4 C.

36 MATERIALS SCIENCE↗
DOE OSTI2021

De Sitter decays to infinity

Bubbles of nothing are a class of vacuum decay processes present in some theories with compactified extra dimensions. We investigate the existence and properties of bubbles of nothing in models where the scalar pseudomoduli controlling the size of the extra dimensions are stabilized at positive vacuum energy, which is a necessary feature of any realistic model. We map the construction of bubbles of nothing to a four-dimensional Coleman-De Luccia problem and establish necessary conditions on the asymptotic behavior of the scalar potential for the existence of suitable solutions. We perform detailed analyses in the context of five-dimensional theories with metastable dS 4 x S 1 vacua, using analytic approximations and numerical methods to calculate the decay rate. We find that bubbles of nothing sometimes exist in potentials with no ordinary Coleman-De Luccia decay process, and that in the examples we study, when both processes exist, the bubble of nothing decay rate is typically faster. Our methods can be generalized to other stabilizing potentials and internal manifolds.

classical theories of gravity↗
DOE OSTI2024

Using GANs to predict milling stability from limited data

Milling is a key manufacturing process that requires the selection of operating parameters that provide efficient performance. However, the presence of chatter, a self-excited vibration causing poor surface finish and potential damage to the machine and cutting tool, makes it challenging to select the appropriate parameters. To predict chatter, stability maps are commonly used, but their generation requires expensive data, making it difficult to employ these maps in industry. Therefore, there is a pressing need for an approach that can accurately predict stability maps using limited experimental data. This study introduces the new Encoder GAN (EGAN) approach based on Generative Adversarial Networks (GANs) that predicts stability maps using limited experimental data. The approach consists of the encoder, generator, and discriminator subnetworks and uses the trained encoder and generator to predict the target stability map. This versatile method can be applied to various tool setups and can accurately predict stability maps with limited experimental data (five to 10 cutting tests) even when there is little information available for unknown parameters. In conclusion, the study evaluates the proposed approach using both numerical data and experiments and demonstrates its superior performance compared to state-of-the-art benchmarks.

42 ENGINEERING↗
DOE OSTI2026

Multilevel Conditional Disturbance Rejection Control for Satellite Attitude Tracking

Recently, the conditional disturbance rejection controller (CDRC) was proposed to improve control performance by leveraging disturbances that have beneficial effects on the system. However, it only considers disturbances acting on the state variable directly influenced by them. Although fast convergence of this state can be achieved with the CDRC, it may unintentionally affect the convergence of the output (i.e., the primary state). Here, in this article, a multilevel CDRC is proposed to enhance satellite attitude control performance by accounting for the effect of disturbances on both attitude (output) and angular velocity. The extended state observer is employed to estimate the lumped disturbance, including the modeling errors and external disturbances. Then, a backstepping-based controller with the multilevel disturbance rejection law (ML-DRL) is designed for attitude tracking. The ML-DRL is developed to improve the control performance by using a disturbance with a damping effect on both attitude and velocity. Faster convergence of attitude and velocity can be achieved by conditionally compensating for the disturbance. The stability of the proposed control method is analyzed by demonstrating that the errors are bounded as time tends to infinity. The attitude control performance of the proposed method is evaluated through numerical examples conducted using the MATLAB/Simulink Multibody tool.

Active disturbance rejection control (ADRC)↗
DOE OSTI2023

In situ microscopy and spectroscopy characterization of microsized Sn anode for sodium-ion batteries

Microsized Sn is a promising anode material for sodium-ion batteries in terms of cost, specific capacity, and volumetric energy density, which however suffers from huge volume changes and rapid cell degradation upon cycling. Despite recent advances via nanostructured electrode design and interface engineering, the correlation between mechanical stability, solid-electrolyte interphase (SEI) and reaction kinetics/reversibility remains controversial and elusive. Here, in this work, by combining in situ scanning electron microcopy and X-ray absorption spectroscopy as well as X-ray photoelectron spectroscopy, we have investigated the underlying electro-chemo-mechanical behavior and their coupling effects during charge/discharge of microsized Sn anode. Our results revealed that microsized Sn is pulverized into nanoparticles with simultaneous formation of numerous voids and pores upon the 1st charge/discharge, while the electrolytes composition plays a critical role on the consequent parasitic reactions and eventually the sodiation/de-sodiation reversibility. In contrast to carbonate-based electrolytes, ether-based electrolytes enabled formation of inorganic species dominated SEI with improved mechanical strength, thus leading to higher specific capacity and improved cycling stability. The present findings are crucial for future development of microsized anode materials for rechargeable batteries with high volumetric energy density.

25 ENERGY STORAGE↗
DOE OSTI2022

The many-body electronic interactions of Fe(II)–porphyrin

Fe(II)–porphyrin complexes exhibit a diverse range of electronic interactions between the metal and macrocycle. Herein, the incremental full configuration interaction method is applied to the entire space of valence orbitals of a Fe(II)–porphyrin model using a modest basis set. A novel visualization framework is proposed to analyze individual many-body contributions to the correlation energy, providing detailed maps of this complex’s highly correlated electronic structure. Furthermore, this technique is used to parse the numerous interactions of two low-lying triplet states ( 3 A 2 g and 3 E g ) and to show that strong metal d–d and macrocycle π–π orbital interactions preferentially stabilize the 3 A 2 g state. d–π interactions, on the other hand, preferentially stabilize the 3 E g state and primarily appear when correlating six electrons at a time. Ultimately, the Fe(II)–porphyrin model’s full set of 88 valence electrons are correlated in 275 orbitals, showing the interactions up to the 4-body level, which covers the great majority of correlations in this system.

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

Error-correcting codes for fermionic quantum simulation

Utilizing the framework of \mathbb{Z}_2 ℤ 2 lattice gauge theories in the context of Pauli stabilizer codes, we present methodologies for simulating fermions via qubit systems on a two-dimensional square lattice. We investigate the symplectic automorphisms of the Pauli module over the Laurent polynomial ring. This enables us to systematically increase the code distances of stabilizer codes while fixing the rate between encoded logical fermions and physical qubits. We identify a family of stabilizer codes suitable for fermion simulation, achieving code distances of d=2,3,4,5,6,7, allowing correction of any \lfloor \frac{d-1}{2} \rfloor ⌊ d − 1 2 ⌋ -qubit error. In contrast to the traditional code concatenation approach, our method can increase the code distances without decreasing the (fermionic) code rate. In particular, we explicitly show all stabilizers and logical operators for codes with code distances of d=3,4,5. We provide syndromes for all Pauli errors and invent a syndrome-matching algorithm to compute code distances numerically.

Physics↗
DOE OSTI2025

Block-Structured Operator Inference for Coupled Multiphysics Model Reduction

This work presents a block-structured formulation of Operator Inference as a way to learn structured reduced-order models for multiphysics systems. The approach specifies the governing equation structure for each physics component and the structure of the coupling terms. Once the multiphysics structure is specified, the reduced-order model is learned from snapshot data following the nonintrusive Operator Inference methodology. In addition to preserving physical system structure, which in turn permits preservation of system properties such as stability and second-order structure, the block-structured approach has the advantages of reducing the overall dimensionality of the learning problem and admitting tailored regularization for each physics component. The numerical advantages of the block-structured formulation over a monolithic Operator Inference formulation are demonstrated for aeroelastic analysis, which couples aerodynamic and structural models. For the benchmark test case of the AGARD 445.6 wing, block-structured Operator Inference provides an average 20% online prediction speedup over monolithic Operator Inference across subsonic and supersonic flow conditions in both the stable and fluttering parameter regimes while preserving the accuracy achieved with monolithic Operator Inference.

42 ENGINEERING↗
DOE OSTI2024

Assessment of Cementitious Composites for High-Temperature Geothermal Wells

High-temperature (HT) geothermal wells can provide green power 24 hours a day, 7 days a week. Under harsh environmental and operational conditions, the long-term durability requirements of such wells require special cementitious composites for well construction. This paper reports a comprehensive assessment of geothermal cement composites in cyclic pressure function laboratory tests and field exposures in an HT geothermal well (300–350 °C), as well as a numerical model to complement the experimental results. Performances of calcium–aluminate cement (CAC)-based composites and calcium-free cement were compared against the reference ordinary Portland cement (OPC)/silica blend. The stability and degradation of the tested materials were characterized by crystalline composition, thermo-gravimetric and elemental analyses, morphological studies, water-fillable porosity, and mechanical property measurements. All CAC-based formulations outperformed the reference blend both in the function and exposure tests. The reference OPC/silica lost its mechanical properties during the 9-month well exposure through extensive HT carbonation, while the properties of the CAC-based blends improved over that period. The Modified Cam-Clay (MCC) plasticity parameters of several HT cement formulations were extracted from triaxial and Brazilian tests and verified against the experimental results of function cyclic tests. These parameters can be used in well integrity models to predict the field-scale behavior of the cement sheath under geothermal well conditions.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Virtual element approximations of the time-fractional nonlinear convection-diffusion equation on polygonal meshes

We extend the Virtual Element Method to a two-dimensional unsteady nonlinear convection-diffusion equation characterized by a fractional-order derivative with respect to the time variable. Our methodology is based on three fundamental technical components: a fractional version of the Grunwald-Letnikov approximation, discrete maximal regularity, and the regularity theory associated with non-linearity. We prove the method's well-posedness, i.e., the approximate solution's existence and uniqueness to the time-fractional convection-diffusion equation with a Lipschitz nonlinear source term. The fully discrete scheme inherently maintains stability and consistency by leveraging the discrete maximal regularity and the energy projection operator. The convergence in the L 2 -norm and H 1 -norm to various mesh configurations is validated by numerical results, underlining the practical effectiveness of the proposed method.

97 MATHEMATICS AND COMPUTING↗