Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “discretization”

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 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Critical amino acid residues in the N-terminal domain of NADPH-dependent assimilatory sulfite reductase flavoprotein mediate octameric assembly

How large, flexible enzymes assemble into defined oligomeric architectures remains a central question in biology. NADPH-dependent assimilatory sulfite reductase (SiR) forms a heterododecamer built on an octameric flavoprotein (SiRFP) core, yet the molecular basis for this assembly has been unresolved because of its disordered N-terminus. Here, we use ion mobility mass spectrometry, small-angle neutron scattering, and mutagenesis to define the mechanism of SiRFP oligomerization. We show that SiRFP forms a discrete, stable octamer in solution. We also report that its N-terminal 52-residue segment is necessary and sufficient to mediate assembly, also mediating oligomerization when fused to a heterologous protein. Structure-guided mutagenesis identifies four residues (Gln22, Tyr39, Phe40, and Gln47) whose substitution disrupts the octamer, producing concentration-dependent lower-order species while retaining catalytic activity. These findings define the determinants of SiRFP assembly with broader implications for engineering homomeric protein complexes.

Nagy, Gergely [ORNL] (ORCID:0000000327420198)↗
DOE OSTI2021

The anomaly that was not meant IIB

Type IIB supergravity enjoys a discrete non-Abelian duality group, which has potential quantum anomalies. In this paper we explicitly compute these, and present the bordism group that controls them, modulo some physically motivated assumptions. Quite surprisingly, we find that they do not vanish, which naively would signal an inconsistency of F-theory. Remarkably, a subtle modification of the standard 10d Chern-Simons term cancels these anomalies, a fact which relies on the specific field content of type IIB supergravity. We also discover other ways to cancel this anomaly, via a topological analog of the Green-Schwarz mechanism. These alternative type IIB theories have the same low energy supergravity limit as ordinary type IIB, but a different spectrum of extended objects. They could either be part of the Swampland, or connect to the standard theory via domain walls.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Protein model quality assessment using rotation–equivariant transformations on point clouds

Machine learning research concerning protein structure has seen a surge in popularity over the last years with promising advances for basic science and drug discovery. Working with macromolecular structure in a machine learning context requires an adequate numerical representation, and researchers have extensively studied representations such as graphs, discretized 3D grids, and distance maps. As part of CASP14, we explored a new and conceptually simple representation in a blind experiment: atoms as points in 3D, each with associated features. These features—initially just the basic element type of each atom—are updated through a series of neural network layers featuring rotation-equivariant convolutions. Starting from all atoms, we further aggregate information at the level of alpha carbons before making a prediction at the level of the entire protein structure. We find that this approach yields competitive results in protein model quality assessment despite its simplicity and despite the fact that it incorporates minimal prior information and is trained on relatively little data. As a result, its performance and generality are particularly noteworthy in an era where highly complex, customized machine learning methods such as AlphaFold 2 have come to dominate protein structure prediction.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Lattice Dynamics of Sb 2 Se 3 from Inelastic Neutron and X-Ray Scattering

The lattice dynamics of orthorhombic Sb 2 Se 3 is studied by a combination of inelastic neutron and 121 Sb nuclear inelastic scattering giving access to the total and Sb partial density of phonon states (DPS). The Se partial DPS is determined from the difference between the total and Sb partial DPS. The total DPS is determined at 39, 150, and 300 K, and an analysis of the temperature-induced mode shifts in combination with low-temperature powder diffraction data is provided. Using an earlier reported theoretical approach, the corresponding total and partial DPS of Sb 2 Se 3 are calculated by first-principles calculations. Herein, a detailed analysis of the Grüneisen parameter, element-specific and bulk Debye temperatures, and the mean force constants as derived from the experimental data and discrete Fourier transform calculations is provided. In general, the calculations underestimate the strength of the covalent Sb-Se bonds.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

A two-fluid single-column model of turbulent shallow convection. Part II: Single-column model formulation and numerics

The two-fluid single-column model of Thuburn et al. (Quart. J. R. Meteorol. Soc., 2019, 145, 1535–1550) is extended to include moisture and horizontal wind shear. Turbulent kinetic energy is introduced as a prognostic variable, dependence on a diagnosed boundary-layer height is removed, and subfilter fluxes are approximated using a two-fluid version of a Mellor–Yamada scheme. Three mechanisms for entrainment and detrainment processes are introduced, which represent entrainment of unstable air at the surface, forced detrainment of air at the top of the boundary/cloud layers, and turbulent mixing that relaxes the convective fluid to a reference profile. A semi-implicit Eulerian discretization replaces the semi-implicit semi-Lagrangian implementation of Thuburn et al. (Quart. J. R. Meteorol. Soc., 2019, 145, 1535–1550) to improve numerical stability and conservation. The equations for the implicit time step are solved using a quasi-Newton method, which is shown to perform well in numerical tests for conservation and convergence. The two-fluid single-column model presented in this article will be applied to simulations of shallow cumulus convection in Part III.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Learning compact physics‐aware delayed photocurrent models using dynamic mode decomposition

Abstract Radiation‐induced photocurrent in semiconductor devices can be simulated using complex physics‐based models, which are accurate, but computationally expensive. This presents a challenge for implementing device characteristics in high‐level circuit simulations where it is computationally infeasible to evaluate detailed models for multiple individual circuit elements. In this work we demonstrate a procedure for learning compact delayed photocurrent models that are efficient enough to implement in large‐scale circuit simulations, but remain faithful to the underlying physics. Our approach utilizes dynamic mode decomposition (DMD), a system identification technique for learning reduced‐order discrete‐time dynamical systems from time series data based on singular value decomposition. To obtain physics‐aware device models, we simulate the excess carrier density induced by radiation pulses by solving numerically the ambipolar diffusion equation, then use the simulated internal state as training data for the DMD algorithm. Our results show that the significantly reduced‐order delayed photocurrent models obtained via this method accurately approximate the dynamics of the internal excess carrier density—which can be used to calculate the induced current at the device boundaries—while remaining compact enough to incorporate into larger circuit simulations.

Hanson, Joshua↗
DOE OSTI2025

Hierarchically Structured Vitrimer Biocomposites for Sustainable Manufacturing

Polymers containing dynamic covalent bonds (DCBs) exhibit thermoplastic-like flow above their topology freezing temperature (T v ) while maintaining thermoset-like properties below it, making them promising for sustainable manufacturing. However, their large-scale adoption remains limited due to challenges in accurately determining T v and achieving efficient fiber-matrix bonding in composite applications. Here, hierarchically structured epoxy-anhydride-based polyester vitrimer composites reinforced with cellulosic filaments is demonstrated, where hydroxyl groups on fiber surfaces participate directly in transesterification with the matrix. Further, this dynamic interfacial bonding delivers exceptional mechanical properties, including ≈70 MPa shear strength and >10% strain-to-failure, while enabling thermal malleability. Using a combination of nuclear magnetic resonance and nano-infrared spectroscopies, direct evidence is provided that chemical bond exchange begins well below the conventionally measured T v , supporting the hypothesis that rheologically determined T v reflects a combination of chemical exchange and frictional dynamics rather than a discrete transition. The composites demonstrate excellent processability through vacuum-assisted resin transfer molding and maintain >90% of their mechanical properties after multiple thermal reforming cycles. These findings advance both the fundamental understanding of vitrimeric transitions and the practical development of sustainable, high-performance composite materials.

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

Impact of genotype‐calling methodologies on genome‐wide association and genomic prediction in polyploids

Abstract Discovery and analysis of genetic variants underlying agriculturally important traits are key to molecular breeding of crops. Reduced representation approaches have provided cost‐efficient genotyping using next‐generation sequencing. However, accurate genotype calling from next‐generation sequencing data is challenging, particularly in polyploid species due to their genome complexity. Recently developed Bayesian statistical methods implemented in available software packages, polyRAD, EBG, and updog, incorporate error rates and population parameters to accurately estimate allelic dosage across any ploidy. We used empirical and simulated data to evaluate the three Bayesian algorithms and demonstrated their impact on the power of genome‐wide association study (GWAS) analysis and the accuracy of genomic prediction. We further incorporated uncertainty in allelic dosage estimation by testing continuous genotype calls and comparing their performance to discrete genotypes in GWAS and genomic prediction. We tested the genotype‐calling methods using data from two autotetraploid species, Miscanthus sacchariflorus and Vaccinium corymbosum , and performed GWAS and genomic prediction. In the empirical study, the tested Bayesian genotype‐calling algorithms differed in their downstream effects on GWAS and genomic prediction, with some showing advantages over others. Through subsequent simulation studies, we observed that at low read depth, polyRAD was advantageous in its effect on GWAS power and limit of false positives. Additionally, we found that continuous genotypes increased the accuracy of genomic prediction, by reducing genotyping error, particularly at low sequencing depth. Our results indicate that by using the Bayesian algorithm implemented in polyRAD and continuous genotypes, we can accurately and cost‐efficiently implement GWAS and genomic prediction in polyploid crops.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Modeling the effect of wind speed and direction shear on utility‐scale wind turbine power production

Abstract Wind speed and direction variations across the rotor affect power production. As utility‐scale turbines extend higher into the atmospheric boundary layer (ABL) with larger rotor diameters and hub heights, they increasingly encounter more complex wind speed and direction variations. We assess three models for power production that account for wind speed and direction shear. Two are based on actuator disc representations, and the third is a blade element representation. We also evaluate the predictions from a standard power curve model that has no knowledge of wind shear. The predictions from each model, driven by wind profile measurements from a profiling LiDAR, are compared to concurrent power measurements from an adjacent utility‐scale wind turbine. In the field measurements of the utility‐scale turbine, discrete combinations of speed and direction shear induce changes in power production of −19% to +34% relative to the turbine power curve for a given hub height wind speed. Positive speed shear generally corresponds to over‐performance and increasing magnitudes of direction shear to greater under‐performance, relative to the power curve. Overall, the blade element model produces both higher correlation and lower error relative to the other models, but its quantitative accuracy depends on induction and controller sub‐models. To further assess the influence of complex, non‐monotonic wind profiles, we also drive the models with best‐fit power law wind speed profiles and linear wind direction profiles. These idealized inputs produce qualitative and quantitative differences in power predictions from each model, demonstrating that time‐varying, non‐monotonic wind shear affects wind power production.

Energy & Fuels↗
DOE OSTI2026

Multipoint Aerostructural Optimization of Wind Turbine Rotors Using a Coupled Blade‐Resolved Aerostructural Solver

Physics‐based design optimization workflows thread the needle between computational cost limitations and simulation complexity, often compromising between modeling detail and the range of operating design conditions. Multipoint aerostructural optimization of wind turbine rotors has so far been confined to low‐fidelity analyses or to high‐fidelity studies with simplified structural models, leaving the most complex design trade‐offs unexplored. We close this gap by performing the first tightly coupled gradient‐based multipoint aerostructural rotor optimization using 3D aerodynamic and structural solvers with discrete coupled adjoints. The optimizer simultaneously varies blade planform, airfoil shapes, and structural thickness through more than 270 design variables, minimizing a weighted combination of rotor mass and power across multiple wind speeds. Applied to a modified DTU 10‐MW benchmark under conservative structural and aerodynamic constraints, our multipoint optimization reduces rotor mass by up to 36% and increases power by 12%–15% across the main operating conditions; biasing the objective toward power yields power gains up to 18% and a 17% mass reduction. For a nominal wind distribution, 3‐point rotor designs accounting for low RPM and high thrust conditions capture dominant trade‐offs and outperform single‐point designs. Adding two off‐design points changes individual‐condition power by less than 3% but leaves the weighted average within 0.5%, and the mass‐power bias has a stronger effect on the final design than the operating‐point weighting itself. Our framework extends naturally to richer load cases and site‐specific wind distributions, providing a basis for high‐fidelity multipoint design earlier in industrial workflows.

17 WIND ENERGY↗
DOE OSTI2024

A taxonomy of automatic differentiation pitfalls

Automatic differentiation is a popular technique for computing derivatives of computer programs. While automatic differentiation has been successfully used in countless engineering, science, and machine learning applications, it can sometimes nevertheless produce surprising results. In this paper, we categorize problematic usages of automatic differentiation, and illustrate each category with examples such as chaos, time-averages, discretizations, fixed-point loops, lookup tables, linear solvers, and probabilistic programs, in the hope that readers may more easily avoid or detect such pitfalls. We also review debugging techniques and their effectiveness in these situations.

Autodiff↗
DOE OSTI2021

A Sparse-Grid Probabilistic Scheme for Approximation of the Runaway Probability of Electrons in Fusion Tokamak Simulation

Runaway electrons (RE) generated during magnetic disruptions present a major threat to the safe operation of plasma nuclear fusion reactors. A critical aspect of understanding RE dynamics is to calculate the runaway probability, i.e., the probability that an electron in the phase space will runaway on, or before, a prescribed time. Such probability can be obtained by solving the adjoint equation of the underlying Fokker-Planck equation that controls the electron dynamics. In this effort, we present a sparse-grid probabilistic scheme for computing the runaway probability. The key ingredient of our approach is to represent the solution of the adjoint equation as a conditional expectation, such that discretizing the differential operator reduces to the approximation of a set of integrals. Adaptive sparse grid interpolation is utilized to approximate the map from the phase space to the runaway probability. The main novelties of this effort are the integration of the sparse-grid method into the probabilistic numerical scheme for computing escape probability, and the application of the proposed method in computing RE probabilities. Two numerical examples are given to illustrate that the proposed method can achieve O(Δt) convergence, and that the local anisotropic adaptive refinement strategy (M. Stoyanov, Adaptive sparse grid construction in a context of local anisotropy and multiple hierarchical parents. In: Sparse Grids and Applications-Miami 2016, Springer, Berlin, 2018, pp. 175–199) can effectively handle the sharp transition layer between the runaway and non-runaway regions.

Yang, Minglei↗
DOE OSTI2023

High-Order Mesh hr-adaptivity for Surface Fitting to Implicit Geometries

We present an ℎ𝑟-adaptivity framework for morphing a given mesh to fit a target surface prescribed as the zero isocontour of a discrete function. In this framework, high-order meshing is posed as a variational minimization problem that depends on the mesh quality prescribed via the target matrix optimization paradigm (TMOP) and position of a subset of mesh nodes with respect to the target surface. The proposed formulation ensures that the variational problem is converged and mesh quality degradation near the surface is limited, even when the mesh topology is incompatible with the target surface. Additionally, a mesh subset-based approach and ℎ-refinement is introduced to efficiently increase fitting accuracy while reducing the computational cost of the mesh morphing problem. The ℎ𝑟-adaptivity technique extends to different element types in two- and three-dimensions, and can be used in existing finite element and spectral element frameworks to obtain high-order body-fitted meshes. Various numerical experiments demonstrate the robustness and accuracy of the fitting approach for problems of practical interest such as Lagrangian hydrodynamics and topology optimization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Scalar-induced neutrinoless double beta decay in SU(5)

We discuss the role of heavy scalar fields in mediating neutrinoless double beta decay (0νββ) within the SU(5) Grand Unified Theory framework, extended suitably to include neutrino mass. In such a minimal realistic SU(5) setup for fermion masses, the scalar contributions to 0νββ are extremely suppressed as a consequence of the proton decay bound. We circumvent this problem by imposing a discrete Z 3 symmetry. However, the scalar contributions to 0νββ remain suppressed in this SU(5) x Z 3 model due to the neutrino mass constraint. We find that the 0νββ contribution can be enhanced by extending the scalar sector with an additional 15-dimensional scalar representation with suitable Z 3 charge. Such an extension not only yields realistic fermion mass spectra but also leads to experimentally testable predictions in upcoming ton-scale 0νββ searches, which can be used as a sensitive probe of the new scalars across a broad range, from LHC-accessible scales up to ∼ 10 10 GeV.

Baryon/Lepton Number Violation↗
DOE OSTI2021

Heisenberg spin chain as a worldsheet coordinate for lightcone quantized string

Although the energy spectrum of the Heisenberg spin chain on a circle defined by H = 1 4 Σ k = 1 M σ k x σ k + 1 x + σ k y σ k + 1 y + Δ σ k z σ k + 1 z is well known for any fixed M, the boundary conditions vary according to whether Mϵ 4N + r, where r = -1,0,1,2, and also according to the parity of the number of overturned spins in the state, In string theory all these cases must be allowed because interactions involve a string with M spins breaking into strings with M 1 < M and M - M 1 spins (or vice versa). We organize the energy spectrum and degeneracies of H in the case Δ = 0 where the system is equivalent to a system of free fermions. In spite of the multiplicity of special cases, in the limit M → ∞ the spectrum is that of a free compactified worldsheet field. Such a field can be interpreted as a compact transverse string coordinate x(σ) ≡ x(σ) + R 0 . We construct the bosonization formulas explicitly in all separate cases, and for each sector give the Virasoro conformal generators in both fermionic and bosonic formulations. Furthermore from calculations in the literature for selected classes of excited states, there is strong evidence that the only change for Δ ≠ 0 is a change in the compactification radius R 0 → R Δ . As Δ→ -1 this radius goes to infinity, giving a concrete example of noncompact space emerging from a discrete dynamical system. Finally we apply our work to construct the three string vertex implied by a string whose bosonic coordinates emerge from this mechanism.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Vector boson dark matter from trinification

We show how trinification models based on the gauge group SU(3) C x SU(3) L x SU(3) R realized near the TeV scale can provide naturally a variety of dark matter (DM) candidates. These models contain a discrete T parity which may remain unbroken even after spontaneous symmetry breaking. The lightest T-odd particle, which could be a fermion, a scalar, or a gauge boson, can constitute the dark matter of the universe. This framework naturally admits a doublet-singlet fermionic DM, a singlet scalar DM, or a vector boson DM. Here we develop the vector boson DM scenario wherein the DM couples off-diagonally with the usual fermions and vector-like fermions present in the theory. We show consistency of this framework with dark matter relic abundance and direct detection limits as well as LHC constraints. We derive upper limits of 900 GeV on the vector gauge boson DM mass and 4.5 TeV on the vector-like quark masses. We also show the consistency of spontaneous gauge symmetry breaking down to Standard Model times an extra U(1) while preserving the T-parity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Solving the strong CP problem with massless grand-color quarks

We propose a solution to the strong CP problem that specifically relies on massless quarks and has no light axion. The QCD color group SU(3) c is embedded into a larger, simple gauge group (grand-color) where one of the massless, colored fermions enjoys an anomalous chiral symmetry, rendering the strong CP phase unphysical. The grand-color gauge group G GC is Higgsed down to SU(3) c × ${G}_{c^{\prime }}$, after which ${G}_{c^{\prime }}$ eventually confines at a lower scale, spontaneously breaking the chiral symmetry and generating a real, positive mass to the massless, colored fermion. Since the chiral symmetry has a ${G}_{c^{\prime }}$ anomaly, there is no corresponding light Nambu-Goldstone boson. The anomalous chiral symmetry can be an accidental symmetry that arises from an exact discrete symmetry without introducing a domain wall problem. Potential experimental signals of our mechanism include vector-like quarks near the TeV scale, pseudo Nambu-Goldstone bosons below the 10 GeV scale, light dark matter decay, and primordial gravitational waves from the new strong dynamics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Modified instanton sum in QCD and higher-groups

We consider the SU(N) Yang-Mills theory, whose topological sectors are restricted to the instanton number with integer multiples of p. We can formulate such a quantum field theory maintaining locality and unitarity, and the model contains both 2π-periodic scalar and 3-form gauge fields. This can be interpreted as coupling a topological theory to Yang-Mills theory, so the local dynamics becomes identical with that of pure Yang-Mills theory. The theory has not only $\mathbb{Z}$ N 1-form symmetry but also $\mathbb{Z}$ p 3-form symmetry, and we study the global nature of this theory from the recent ’t Hooft anomaly matching. The computation of ’t Hooft anomaly incorporates an intriguing higher-group structure. We also carefully examine that how such kinematical constraint is realized in the dynamics by using the large-N and also the reliable semiclassics on $\mathbb{R}$ 3 × S 1 , and we find that the topological susceptibility plays a role of the order parameter for the $\mathbb{Z}$ p 3-form symmetry. Introducing a fermion in the fundamental or adjoint representation, we find that the chiral symmetry becomes larger than the usual case by $\mathbb{Z}$ p , and it leads to the extra p vacua by discrete chiral symmetry breaking. No dynamical domain wall can interpolate those extra vacua since such objects must be charged under the 3-form symmetry in order to match the ’t Hooft anomaly.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗