Search NASA⌕ Search

SEARCH · Search NASA

Results for “Regularization”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A Block-Structured Adaptive Mesh Framework to Solve Radiation Transfer Equation in Irregular Embedded Geometries

Radiation transport arises in various scientific, industrial, and medical fields, and understanding its effect in applications is needed to make accurate predictions, safety assessments and performance optimizations. Solving the Radiation Transport Equation (RTE) is challenging due to its integro-differential nature, which involves both differential and integral terms. The differential term describes the change in radiation intensity due to absorption and emission, while the integral term accounts for scattering. The accurate modeling of radiation is further complicated in many applications due to the complex, irregular geometries. Various methods exist for solving the RTE, including the zonal, Monte Carlo, spherical harmonics, discrete ordinates, and finite volume methods. Traditional mesh-based approaches, which rely on structured or unstructured meshes, struggle with irregular geometries due to: a) the difficulty of conforming structured grids to irregular domains, b) challenges in enforcing boundary conditions correctly, and c) the additional computational cost of unstructured mesh methods. This work presents a second-order accurate method for solving the RTE in irregular geometries. The radiation intensity is discretized using the finite-volume method in both spatial and angular directions on regular Cartesian grid blocks. Leveraging the block-structured adaptive mesh refinement (AMR) framework provided by AMReX, our method refines the grid locally to reduce spatial discretization error, ensuring a converged numerical solution while minimizing computational costs elsewhere. A two-stage deferred correction approach is employed: First, a first-order discretization on grid blocks is solved using an algebraic multigrid method in HYPRE. Second, a correction term is applied explicitly to achieve second-order accuracy. The correction term is calculated by approximating the radiation flux on cell faces using a Total Variation Diminishing (TVD) scheme. This approach ensures quick convergence of the multigrid method while preserving higher-order accuracy of the numerical solution. Irregular geometries are resolved as embedded boundaries (EB), resulting in both cut cells and regular cells. In cut cells, we modify the fluxes using face fractions and incorporate additional contributions from EB boundary conditions. To ensure higher-order convergence near the EB interface, the correction term is modified by interpolating the radiation intensity to fictitious ghost points. The implementation takes advantage of modern supercomputers by leveraging AMReX’sMPI/X parallelization strategy where X can be MPI or a GPU accelerator including CUDA, HIP and DPC++. We validate our solver using classical test cases, both with and without EB, demonstrating accuracy and efficiency. Additionally, we analyze the impact of adaptive mesh refinement on solution accuracy and computational cost, highlighting the advantages of our approach for high-resolution radiation transport simulations.

computational fluid dynamics (CFD)↗
DOE OSTI2025

Next-Level Energy Management in Manufacturing: Facility-Level Energy Digital Twin Framework Based on Machine Learning and Automated Data Collection

This research introduces an energy prediction framework at the facility level supported by automated data collection and machine learning models. It investigates whether reducing the prediction time scale allows for applying more complex machine learning techniques and if those techniques improve the prediction accuracy. The primary advantages of this framework lie in its automation of the energy prediction process and its provision of real-time energy data suitable for use in energy dashboards or digital twins. A sitewide dataset was created by combining 15 min energy and daily production data of five shops—assembly, battery, body (electric), body (gas), and paint—from a globally recognized electric vehicle manufacturer. Various machine learning models were evaluated on daily, weekly, and monthly datasets, including, in increasingly complex order: naïve, simple linear regression, net regularized generalized linear regression, principal component regression, k-nearest neighbor, random forest, and Bayesian regularized neural network. Compared to the current state-of-the-art energy consumption prediction for the industrial facility level, this research investigates more complex models and smaller time intervals for higher accuracy. The findings revealed that the more complex monthly models require a minimum of a year and a half of data to operate, while weekly models demand a year of data to achieve improved accuracy. Daily models can operate with only six months of data but exhibit poor performance due to reduced prediction accuracy of production. Key challenges identified include access to reliable, high-quality energy and production data and the initial demand for human labor.

digital twin↗
DOE OSTI2025

Investigations on the Thermal Stability and Kinetics of Biolubricants Synthesized from Different Types of Vegetable Oils

Petroleum-based lubricants raise environmental concerns due to their non-biodegradability and toxicity, whereas biobased lubricants underperform owing to low thermal stability. This study examined and compared three vegetable oils, along with their chemically modified versions, to better understand their suitability as biolubricants. High oleic soybean oil (HOSOY), regular soybean oil (RSOY), and waste cooking oil (WCO) were subjected to chemical modification, where isopropyl groups were attached to the fatty acid chains of the oils to produce branched oils, i.e., b-HOSOY, b-RSOY, and b-WCO. The detailed kinetic study of each regular and modified sample was investigated using thermogravimetric analysis. The kinetic parameters, such as the activation energies, reaction rate, and pre-exponential factor, were generated via Friedman methods. The differential thermal gravimetric (DTG) analysis showed low volatilization at the onset temperature in each modified oil as compared with the unmodified samples under an oxidative environment. Furthermore, the comparative kinetic studies demonstrated the enhanced thermoxidative stability of the modified products relative to their unaltered counterparts. Among the tested oils, the b-RSOY showed an average activation energy of 325 kJ/mol, followed by the b-WCO: 300 kJ/mol and the b-HOSOY: 251 kJ/mol, indicating the most stable modified product under an oxidative environment. For all the samples, the pre-exponential factors were in good agreement with the activation energies, which validates that finding the pre-exponential components is crucial to the kinetic analysis.

Sarker, Majher I. (ORCID:0000000299509274)↗
DOE OSTI2025

Ginsparg-Wilson Hamiltonians with Improved Chiral Symmetry

We construct a family of Ginsparg-Wilson Hamiltonians with improved chiral properties, starting from a construction of Creutz-Horvath-Neuberger that provides a doubler-free Hamiltonian lattice regularization for Dirac fermions in even spacetime dimensions. We use a higher-order generalization of the Ginsparg-Wilson relation due to Fujikawa, which yields an order-$k$ Hamiltonian overlap operator for each integer $k \geq 0$, with an exactly conserved but nonquantized chiral charge that becomes quantized as $k \to \infty$. Our construction provides physical insight into how Fujikawa's higher-order Ginsparg-Wilson relation improves chiral symmetry while reproducing the anomaly, highlighting the trade-offs inherent in any Hamiltonian lattice realization of an anomalous chiral symmetry. This class of Hamiltonian lattice regularizations, with their tunable chiral symmetry properties, offers potential advantages for quantum and tensor-network simulations.

Singh, Hersh [Fermilab]↗
DOE OSTI2026

Multibody for Everybody (M4E) - A Linearization Approach to Enable Frequency Domain Analysis, Time Integration and Control Co-Design

1.1 Background/Objectives: Marine energy represents a promising yet underexploited source of power. To increase the harvested power, significant efforts have been made to improve wave energy converter (WEC) modeling capabilities and optimize power take-off (PTO) performance; however, these efforts have often treated WEC dynamics, PTO design, and controller development sequentially. In contrast, control co-design (CCD) is emerging as a promising strategy to address these issues directly, creating a growing need for fast analysis tools suitable for repeated simulation and parametric studies [1]. To support this need, this work presents the Multibody for Everybody (M4E) [2] linearization module, which employs a symbolic toolbox to provide deeper insight of WEC design parameters. The objective is to demonstrate that a minimal-coordinate linearization of articulated WEC dynamics can provide accurate wave response predictions and substantial computational savings relative to nonlinear time-domain simulation, while preserving compatibility with broader wave-energy analysis workflows, enabling CCD. 1.2 Approach/Activities: The proposed approach linearizes the equations of motion, generated by M4E, in minimal coordinates about a selected operating point and combines the resulting system with frequencydomain hydrodynamic terms to incorporate the reduced mass, damping, stiffness, and forcing operators. The linearized model is used for both impedance-based response amplitude operator (RAO) prediction and rapid regular-wave time integration. The methodology is demonstrated on a single-flap device and a FOSWEC configuration, with linearized M4E responses compared against the corresponding nonlinear M4E simulations and WEC-Sim results. Regular-wave time histories, RAO trends, and runtime differences are assessed. The framework is also compatible with broader wave-energy workflows, including coupling to WecOptTool, although that capability is not the focus of this work [3]. 1.3 Results/Lessons: The linearized M4E model reproduces key regularwave response characteristics such as integration and Response Amplitude over multiple frequencies. This module matches nonlinear M4E and WEC-Sim results while substantially reducing integration cost. Thus, the proposed framework can serve as a rapid analysis layer for articulated WEC design, parameter studies, and controls-oriented workflows. The analysis is most appropriate in the near-equilibrium regime, about the linearization point.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

Dynamical Low-Rank Compression of Neural Networks with Robustness under Adversarial Attacks

Deployment of neural networks on resource-constrained devices demands models that are both compact and robust to adversarial inputs. However, compression and adversarial robustness often conflict. In this work, we introduce a dynamical low-rank training scheme enhanced with a novel spectral regularizer that controls the condition number of the low-rank core in each layer. This approach mitigates the sensitivity of compressed models to adversarial perturbations without sacrificing clean accuracy. The method is model- and data-agnostic, computationally efficient, and supports rank adaptivity to automatically compress the network at hand. Extensive experiments across standard architectures, datasets, and adversarial attacks show the regularized networks can achieve over 94 compression while recovering or improving adversarial accuracy relative to uncompressed baselines.

Schotthoefer, Steffen [ORNL] (ORCID:00000002156965↗
DOE OSTI2025

Nondestructive Evaluation of Concrete: Elastic Property Imaging Through Full-Waveform Inversion

Concrete is a vital material in construction—especially in the nuclear industry, where it is used in critical structures such as containment vessels. Over time, concrete can degrade due to harsh operational and environmental conditions, necessitating that its elastic properties be accurately evaluated to ensure structural integrity and safety. Traditional nondestructive evaluation methods such as ultrasound-based techniques often rely on simplifying assumptions that may not hold true for concrete. This paper presents an advanced ultrasound-based method that uses elastic full-waveform inversion (EFWI) to create detailed images of concrete’s mechanical properties. By accurately modeling wave behaviors such as scattering and reflection, we aim to overcome the limitations of conventional ultrasonic-based methods. In this work, the imaging problem involved reconstructing the various elastic properties of a heterogenous concrete block with three steel rebars embedded in it. The ultrasonic measurements were synthetically generated from multiple sources and receivers, and the reconstruction process was performed using a gradient-based optimization algorithm. Our approach leveraged EFWI to reconstruct high-resolution images of the pressure wave speed, shear wave speed, and density. Multiple misfit functions—including L2-norm, cross-correlation (CC), and L1-norm—combined with total variation (TV) regularization and parameter constraints using a Sigmoid function—were explored for the reconstruction. The results demonstrated that using the L1-norm misfit function in conjunction with TV regularization and Sigmoid constraints significantly improved the reconstruction quality in comparison to traditional methods. This approach provided clearer images with fewer artifacts and better captured background heterogeneity. Our findings highlight that, when properly designed, EFWI carries great potential for providing comprehensive, more accurate, and more reliable assessments of concrete conditions, as is crucial for the maintenance and safety of nuclear power plant structures.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Unveiling horizons in quantum critical collapse

Critical gravitational collapse offers a unique window into regimes of arbitrarily high curvature, culminating in a naked singularity arising from smooth initial data — thus providing a dynamical counterexample to weak cosmic censorship. Near the critical regime, quantum effects from the collapsing matter are expected to intervene before full quantum gravity resolves the singularity. Despite its fundamental significance, a self-consistent treatment has so far remained elusive. In this work, we perform a one-loop semiclassical analysis using the robust anomaly-based method in the canonical setup of Einstein gravity minimally coupled to a free, massless scalar field. Focusing on explicitly solvable near-critical solutions in both 2 + 1 and 3 + 1 dimensions, we analytically solve the semiclassical Einstein equations and obtain controlled, quantitative results for several long-standing questions within the dominant s-wave sector. We find that regularity uniquely selects a Boulware-like quantum state, encoding genuine vacuum polarization effects from the collapsing matter. Remarkably, the resulting quantum corrections manifest as a growing mode. Horizon-tracing analyses, incorporating both classical and quantum modes, reveal the emergence of a finite mass gap, signaling a phase transition from classical Type II to quantum-modified Type I behavior, thereby providing a quantum enforcement of the weak cosmic censorship. The most non-trivial aspect of our analysis involves dealing with non-conformal matter fields in explicitly time-dependent critical spacetimes. Along the way, we uncover intriguing and previously underexplored features of quantum field theory in curved spacetime.

2D Gravity↗
DOE OSTI2025

Renormalization-group equations of the LEFT at two loops: dimension-six baryon-number-violating operators

We present the second part of a systematic calculation of the two-loop anomalous dimensions for the low-energy effective field theory below the electroweak scale (LEFT): the baryon-number-violating sector at dimension six in the power counting. We obtain the results in two different schemes: in the algebraically consistent ’t Hooft-Veltman scheme for γ 5 , corrected for evanescent as well as chiral-symmetry-breaking effects through finite renormalizations; and in naive dimensional regularization, which in the considered sector of the theory does not lead to any ill-defined γ 5 -odd traces. Our results are of interest for a reanalysis of the constraints on physics beyond the Standard Model from proton-decay searches within an EFT framework at next-to-leading-logarithmic accuracy.

Baryon/Lepton Number Violation↗
DOE OSTI2025

Stereochemical Control of Water Transport Properties in Thiol‐yne Polymers

Barrier polymers underpin almost every commercial sector, yet the needs of several emerging areas remain unmet by commercially‐available materials, including temporary orthopedic implants, transient health monitors, neural implants, and other long‐term implants. The ability to tune polymer composition independently of polymer structure positions thiol‐yne click chemistry as a promising platform to serve these emerging technologies. Here, this work describes the differences in the hierarchical structure of stoichiometrically identical materials which differ only in the proportion of the cis versus trans backbone alkene stereochemistry. Varying the isomer content in this way directs different temperature and rate dependent crystallization behavior, which affords control over micron‐scale structure. This investigation focuses on how these stereochemical features affect the water vapor permeation process by several methods and develops an understanding of how this unique structural regularity improves barrier performance relative to a commercially available water barrier polymer, poly(ethylene terephthalate).

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Advancing Insights into Electrochemical Pre‐Treatments of Supported Nanoparticle Electrocatalysts by Combining a Design of Experiments Strategy with In Situ Characterization

Activation, break-in, and/or pre-treatment protocols are generally applied to energy conversion devices before regular operation to reach stable performance. There remains much to understand about the relationships among physical properties, performance, and electrochemical pre-treatments. Here, a design-of-experiments (DoE) strategy is employed to address this gap by demonstrating the influence of five pre-treatment parameters for carbon-supported Pt-nanoparticle catalysts on the electrocatalytic oxygen reduction reaction (ORR). A subset of pre-treatments, developed using a central composite design, are tested in a flow cell combined with an inductively-coupled plasma mass spectrometer (on-line ICP-MS). The DoE-based approach facilitates comprehensive insights from two orders of magnitude fewer experiments than a conventional grid search. The coupled on-line ICP-MS setup enables effective catalysis and real-time catalyst dissolution data. Leveraging insights from DoE for on-line ICP-MS and additional characterization, a model is built between the degradation of a multi-dimensional supported Pt surface, its performance, and applied electrochemical parameters. These investigations identify surface modifications, such as oxidation, and subsequent restructuring of Pt during pre-treatment as a primary cause of performance deterioration during ORR. By combining DoE with advanced characterization techniques, a powerful approach is demonstrated to gain a mechanistic understanding of pre-treatment protocols that can be broadly adapted to various reaction chemistries.

Platinum↗
DOE OSTI2024

Arbitrary Order Virtual Element Methods for High‐Order Phase‐Field Modeling of Dynamic Fracture

ABSTRACT Accurate modeling of fracture nucleation and propagation in brittle and ductile materials subjected to dynamic loading is important in predicting material damage and failure under extreme conditions. Phase‐field fracture models have garnered a lot of attention in recent years due to their success in representing damage and fracture processes in a wide class of materials and under a variety of loading conditions. Second‐order phase‐field fracture models are by far the most popular among researchers (and increasingly, among practitioners), but fourth‐order models have started to gain broader acceptance since their more recent introduction. The exact solution corresponding to these high‐order phase‐field fracture models has higher regularity. Thus, numerical solutions of the model equations can achieve improved accuracy and higher spatial convergence rates. In this work, we develop a virtual element framework for the high‐order phase‐field model of dynamic fracture. The virtual element method (VEM) can be regarded as a generalization of the classical finite element method. In addition to many other desirable characteristics, the VEM allows computing on polytopal meshes. Here, we use ‐conforming virtual elements and the generalized‐ time integration method for the momentum balance equation, and adopt ‐conforming virtual elements for the high‐order phase‐field equation. We verify our virtual element framework using classical quasi‐static benchmark problems and demonstrate its capabilities with the aid of numerical simulations of dynamic fracture in brittle materials.

42 ENGINEERING↗
DOE OSTI2024

Implementation and Validation of a Generalized Actuator Disk Parameterization for Wind Turbine Simulations Within the FastEddy Model

ABSTRACT Fast and accurate large‐eddy simulation (LES) of the atmospheric boundary layer plays a crucial role in advancing wind energy research. Long‐duration wind farm studies at turbine‐resolving scales have become increasingly important to understand the intricate interactions between large wind farms and the atmospheric boundary layer. However, the prohibitive computational cost of these turbulence‐ and turbine‐resolving simulations has precluded such modeling to be exercised on a regular basis. To that end, we implement and validate the generalized actuator disk (GAD) model in the computationally efficient, graphics processing unit (GPU)–resident, LES model FastEddy. We perform single‐turbine simulations under three atmospheric stabilities (neutral, unstable, and stable) and compare them against observations from the Scaled Wind Farm Technology (SWiFT) facility and other LES codes from the recent Wakebench turbine wake model benchmark. Our idealized LES results agree well with observed wake velocity deficit and downstream recovery across stability regimes. Turbine response in terms of rotational speed, generated power, torque, and thrust coefficient are well predicted across stability regimes and are consistent with the LES results from the benchmark. The FastEddy simulations are found to be at least two orders of magnitude more efficient than the traditional CPU‐based LES models, opening the door for realistic LES simulations of full wind plants as a viable standard practice.

17 WIND ENERGY↗
DOE OSTI2026

Tracking discontinuities in parameter space

We develop a geometric framework in Feynman-parameter space to determine constraints on the sequential discontinuities of Feynman integrals. Our method is based on tracking the deformation of the integration contour as external kinematics are analytically continued. This procedure imposes powerful constraints on the analytic structure of Feynman integrals, providing crucial inputs for their bootstrap. We demonstrate the usefulness of this framework by applying it to integrals in dimensional regularization, with higher propagator powers, and to examples with non-uniform transcendental weight. The method is illustrated with several one- and two-loop calculations.

Differential and Algebraic Geometry↗
DOE OSTI2026

Gaussian processes for inferring parton distributions

The extraction of parton distribution functions (PDFs) from experimental or lattice QCD data is an ill-posed inverse problem, where regularization strongly impacts both systematic uncertainties and the reliability of the results. We study a framework based on Gaussian Process Regression (GPR) to reconstruct PDFs from lattice QCD matrix elements. Within a Bayesian framework, Gaussian processes serve as flexible priors that encode uncertainties, correlations, and constraints without imposing rigid functional forms. We investigate a wide range of kernel choices, mean functions, and hyperparameter treatments. We quantify information gained from the data using the Kullback-Leibler divergence. Synthetic data tests demonstrate the consistency and robustness of the method. Our study establishes GPR as a systematic and non-parametric approach to PDF reconstruction, offering controlled uncertainty estimates and reduced model bias in lattice QCD analyses.

hadronic spectroscopy↗
DOE OSTI2024

Investigating two-dimensional adjoint QCD on the lattice

We present our investigations of SU(N) adjoint QCD in two dimensions with one Majorana fermion on the lattice. We determine the relevant parameter range for the simulations with Wilson fermions and present results for Polyakov loop, chiral condensate, and string tension. In the theory with massive fermions, all observables we checked show qualitative agreement between numerical lattice data and theory, while the massless limit is more subtle since chiral and non-invertible symmetry of the continuum theory are explicitly broken by lattice regularization. In thermal compactification, we observe N perturbative vacua for the holonomy potential at high-T with instanton events connecting them, and a unique vacuum at low-T. At finite-N, this is a cross-over and it turns to a phase transition at large-N thermodynamic limit. In circle compactification with periodic boundary conditions, we observe a unique center-symmetric minimum at any radius. In continuum, the instantons in the thermal case carry zero modes (for even N) and indeed, in the lattice simulations, we observe that chiral condensate is dominated by instanton centers, where zero modes are localized. We present lattice results on the issue of confinement vs. screening in the theory and comment on the roles of chiral symmetry and non-invertible symmetry.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Lattice realizations of topological defects in the critical (1+1)-d three-state Potts model

Topological/perfectly-transmissive defects play a fundamental role in the analysis of the symmetries of two dimensional conformal field theories (CFTs). In the present work, spin chain regularizations for these defects are proposed and analyzed in the case of the three-state Potts CFT. In particular, lattice versions for all the primitive defects are presented, with the remaining defects obtained from the fusion of the primitive ones. The defects are obtained by introducing modified interactions around two given sites of an otherwise homogeneous spin chain with periodic boundary condition. The various primitive defects are topological on the lattice except for one, which is topological only in the scaling limit. The lattice models are analyzed using a combination of exact diagonalization and density matrix renormalization group techniques. Low-lying energy spectra for different defect Hamiltonians as well as entanglement entropy of blocks located symmetrically around the defects are computed. The latter provides a convenient way to compute the g-function which characterizes various defects. Finally, the eigenvalues of the line operators in the “crossed channel” and fusion of different defect lines are also analyzed. The results are all in agreement with expectations from conformal field theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

On gauge amplitudes first appearing at two loops

We study scattering amplitudes in massless non-abelian gauge theory where all outgoing gluons have positive helicity. It has been argued recently by Costello that for a particular fermion representation (8 fundamentals plus one antisymmetric-tensor representation in SU(N)) the one-loop amplitudes vanish identically. We show that this vanishing leads to previously-observed identities among one-loop color-ordered partial amplitudes. We then turn to two loops, where Costello has computed the all-plus amplitudes for this theory, as rational functions of the kinematics for any number of gluons using the celestial chiral algebra (CCA) bootstrap. We show that in dimensional regularization, these two-loop amplitudes are not rational, and they are not even finite as ϵ → 0. However, the finite remainder for four gluons agrees with the formula by Costello. In addition, we provide a mass regulator for the infrared-divergent loop integrals; with this regulator, the CCA bootstrap formula is recovered exactly. Finally, we use the CCA bootstrap to compute the double-trace terms in the theory at two loops for an arbitrary number of gluons.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗