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

On optimal control of hybrid dynamical systems using complementarity constraints

Optimal control for switch-based dynamical systems is a challenging problem in the process control literature. In this study, we model these systems as hybrid dynamical systems with finite number of unknown switching points and reformulate them using non-smooth and non-convex complementarity constraints as a mathematical program with complementarity constraints (MPCC). We utilize a moving finite element based strategy to discretize the differential equation system to accurately locate the unknown switching points at the finite element boundary and achieve high-order accuracy at intermediate non-collocation points. We propose a globalization approach to solve the discretized MPCC problem using a mixed NLP/MILP-based strategy to converge to a non-spurious first-order optimal solution. The method is tested on three dynamic optimization examples, including a gas–liquid tank model and an optimal control problem with a sliding mode solution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Collocation methods for nonlinear differential equations on low-rank manifolds

We introduce new methods for integrating nonlinear differential equations on low-rank manifolds. These methods rely on interpolatory projections onto the tangent space, enabling low-rank time integration of vector fields that can be evaluated entry-wise. A key advantage of our approach is that it does not require the vector field to exhibit low-rank structure, thereby overcoming significant limitations of traditional dynamical low-rank methods based on orthogonal projection. To construct the interpolatory projectors, we develop a sparse tensor sampling algorithm based on the discrete empirical interpolation method (DEIM) that parameterizes tensor train manifolds and their tangent spaces with cross interpolation. Using these projectors, we propose two time integration schemes on low-rank tensor train manifolds. The first scheme integrates the solution at selected interpolation indices and constructs the solution with cross interpolation. The second scheme generalizes the well-known orthogonal projector-splitting integrator to interpolatory projectors. We demonstrate the proposed methods with applications to several tensor differential equations arising from the discretization of partial differential equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

CLPNets: Coupled Lie–Poisson neural networks for multi-part Hamiltonian systems with symmetries

To accurately compute data-based prediction of Hamiltonian systems, it is essential to utilize methods that preserve the structure of the equations over time. We consider a particularly challenging case of systems with interacting parts that do not reduce to pure momentum evolution. Such systems are essential in scientific computations, such as discretization of a continuum elastic rod, which can be viewed as the group of rotations and translations $SE(3)$. The evolution involves not only the momenta but also the relative positions and orientations of the particles. The presence of Lie group-valued elements, such as relative positions and orientations, poses a problem for applying previously derived methods for data-based computing. We develop a novel method of data-based computation and complete phase space learning of such systems. We follow the original framework of SympNets (Jin et al., 2020) and LPNets (Eldred et al., 2024), building the neural network from phase space mappings that preserve the Lie–Poisson structure. We derive a novel system of mappings that are built into neural networks describing the evolution of such systems. We call such networks Coupled Lie–Poisson Neural Networks, or CLPNets. We consider increasingly complex examples for the applications of CLPNets, starting with the rotation of two rigid bodies about a common axis, progressing to the free rotation of two rigid bodies, and finally to the evolution of two connected and interacting $SE(3)$ components, describing the discretization of an elastic rod into two elements. Our method preserves all Casimir invariants to machine precision, preserves energy to high accuracy, and shows good resistance to the curse of dimensionality, requiring only a few thousand data points for all cases studied (three to eighteen dimensions). Additionally, the method is highly economical in memory requirements, requiring only about 200 parameters for the most complex case considered.

Data-based modeling↗
DOE OSTI2020

The CLAS12 Forward Time-of-Flight system

The Forward Time-of-Flight system for the large-acceptance CLAS12 spectrometer in Hall B at the Thomas Jefferson National Accelerator Facility is described. The system is positioned at distances in the range from 6.2 m to 7.2 m from the beam–target interaction point and spans laboratory polar angles from 5 deg → 45 deg and nearly the full azimuth. The system consists of 540 individual scintillation counters with double-ended readout that range in length from 17 cm to 426 cm of discrete widths of 6 cm, 15 cm, and 22 cm, and of discrete thicknesses of 5 cm and 6 cm. We see the effective counter time resolution for passing charged particles varies from 50 ps for the shortest counters at small angles to 200 ps for the longest counters at large angles. The detectors are part of the forward-angle particle identification system for CLAS12 during offline event reconstruction and are a component of the online data acquisition trigger to select final state event topologies with forward-going charged particles.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

A novel methodology for gamma-ray spectra dataset procurement over varying standoff distances and source activities

The adoption of machine learning approaches for gamma-ray spectroscopy has received considerable attention in the literature. Many studies have investigated the deployment of various algorithm architectures to a specific task. However, little attention has been afforded to the development of the datasets leveraged to train the models. Such training datasets typically span a set of environmental or detector parameters to encompass a problem space of interest to a user. Variations in these measurement parameters will also induce fluctuations in the detector response, including expected pile-up and ground scatter effects. Fundamental to this work is the understanding that 1) the underlying spectral shape varies as the measurement parameters change and 2) the statistical uncertainties associated with two spectra impact their level of similarity. While previous studies attribute some arbitrary discretization to the measurement parameters for the generation of their synthetic training data, this work introduces a principled methodology for efficient spectral-based discretization of a problem space. A signal-to-noise ratio (SNR) respective spectral comparison measure and a Gaussian Process Regression (GPR) model are used to predict the spectral similarity across a range of measurement parameters. This innovative approach effectively showcased its capability by dividing a problem space, ranging from 5 cm to 100 cm standoff distances and 5 μCi–100 μCi of 137 Cs, into three unique combinations of measurement parameters. The findings from this work will aid in creating more robust datasets, which incorporate many possible measurement scenarios, reduce the number of required experimental test set measurements, and possibly enable experimental training data collection for gamma-ray spectroscopy.

data science↗
DOE OSTI2022

A whole-core steady-state thermal-hydraulic model for annular fuel type fluoride-salt-cooled reactors

A whole-core, steady-state thermal-hydraulic model is developed for the fluoride-salt-cooled small modular advanced high-temperature reactor (SmAHTR) that employs an annular fuel configuration. This pre-conceptual reactor design by Oak Ridge National Laboratory (ORNL) has the annular fuel and moderator pins arranged in a hexagonal layout. The FLiBe coolant flows from the bottom to the top of the core, parallel to the hexagonal bundle. The fuel and moderator pins in the core are discretized into finite volumes and the 3-D heat conduction equation is solved to obtain the temperature profile. Inter-fuel assembly conduction is also addressed. For this fuel assembly configuration, the coolant flows through two distinct regions – the hexagonal pin bundle and the annulus between the fuel pin and the tie rod. The fluid flow through the hexagonal bundles is modeled using the subchannel approach, in which the coolant region is discretized into corner, edge and interior subchannels and the resulting conservation equations are systematically solved. The 1-D mass, momentum and energy equations are solved for the annulus channels between the fuel pin and the tie rod. Pertinent closure models from the literature are employed to close the system of equations. We also performed a preliminary code-to-code comparison between the present model and a CFD model.. The resulting thermal-hydraulic model can provide temperature, flow rate and pressure drop profiles for the different solid and fluid regions throughout the entire core. Whole-core thermal-hydraulic results for a representative power profile are presented and discussed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

The divergence of nearby trajectories in soft-sphere DEM

The n-body instability is investigated with the soft-sphere discrete element method. The divergence of nearby trajectories is quantified by the dynamical memory time. Using the inverse proportionality between the dynamical memory time and the largest Lyapunov exponent, the soft-sphere discrete element method results are compared to previous hard-sphere molecular dynamics data for the first time. Good agreement is observed at low concentrations and the degree of instability is shown to increase asymptotically with increasing spring stiffness. At particle concentrations above 30%, the soft-sphere Lyapunov exponents increase faster than the corresponding hard-sphere data. Finally, this paper concludes with a demonstration of how this case study may be used in conjunction with regression testing and code verification activities.

42 ENGINEERING↗
DOE OSTI2023

Deep reinforcement learning assisted co-optimization of Volt-VAR grid service in distribution networks

With the increasing penetration of distributed energy resources in distribution networks, Volt-VAR control and optimization (VVC/VVO) have become very important to ensure an acceptable quality of service to all customers. System operators can rely on slow-responding utility devices, including capacitor banks and on-load tap changing transformers, along with fast-responding battery and photovoltaic (PV) inverters for the VVC/VVO implementation. Because of variations in response time of these two classes of devices, and different control actions (discrete versus continuous), coordinated and optimal scheduling and operation have become of utmost importance. Here, this paper develops a look-ahead deep reinforcement learning (DRL)-based multi-objective VVO technique to improve the voltage profile of active distribution networks, decrease network and inverter power loss, and save the operational cost of the grid. It proposes a deep deterministic policy gradient (DDPG)-based approach to schedule the optimal reactive and/or active power set-points of fast-responding inverters, and a deep Q-network (DQN)-based DRL agent to schedule the discrete decisions variables of slow-responding assets. The reactive power output of PV and battery smart inverters are scheduled at 30-minute intervals and the capacitors’ commitment status is scheduled with several hour intervals. The proposed framework is validated on the modified IEEE 34-bus and 123-bus test cases with embedded PV and PV-plus-storage. To validate the efficacy of the proposed VVO, it is compared with several scenarios, including the base case without VVO, localized droop control of DERs, DDPG-only, and twin delayed DDPG (TD3) agent-based DRL techniques. The results justify the superior performance of the proposed method to improve the voltage profile, reduce network power loss, and minimize the look-ahead grid operational cost while minimizing the undesirable power losses in inverters as a result of power factor adjustments.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Receiver willingness to participate in off-hour service programs

Service trips, a frequently overlooked segment of urban traffic, represent a disproportionately large share of the negative impacts associated with commercial activity. The reason is that although service trips are less frequent than freight trips, they often take longer and thus occupy a significant share of commercial parking. Demand management programs seeking to move service trips off-hours can potentially reduce congestion and emissions but they have been understudied. Here, this paper describes research conducted to investigate receivers’ willingness to participate in off-hour programs for planned service activity. The research reported in this paper analyzes results from a survey of 189 business establishments in New York City and the Capital Region in New York State using descriptive analyses of the data, and by estimating a discrete choice model to gain insight into how business characteristics influence willingness to participate in off-hour services. To assess the impacts of a hypothetical off-hour services program, the authors applied the discrete choice model, together with service trip attraction models, to quantify the hours of daytime parking that could be eliminated by the proposed program in certain ZIP Codes in New York City and the Albany area. The research conducted led to the identification of numerous policy implications that will help policymakers understand and maximize the potential benefits of implementing off-hour service programs.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2025

Application of Lagrangian techniques for calculating the on-axis rotational transform

The Floquet exponents of periodic field lines are studied through the variations of the magnetic action on the magnetic axis, which is assumed to be elliptical. The near-axis formalism developed by Mercier, Solov'ev and Shafranov is combined with a Lagrangian approach. The on-axis Floquet exponent is shown to coincide with the on-axis rotational transform. A discrete solution suitable for numerical implementation is introduced, which gives the Floquet exponents as solutions to an eigenvalue problem. This discrete formalism expresses the exponents as the eigenvalues of a $6$ X $6$ matrix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Crystal structure of aminopentamide hydrogen sulfate, (C 19 H 25 N 2 O)(HSO 4 )

The crystal structure of aminopentamide hydrogen sulfate has been solved and refined using synchrotron X-ray powder diffraction data, and optimized using density functional techniques. Aminopentamide hydrogen sulfate crystallizes in space group P2 1 /c (#14) with a = 17.62255(14), b = 6.35534(4), c = 17.82499(10) Å, β = 96.4005(6)°, V = 1983.906(14) Å 3 , and Z = 4. The structure consists of layers parallel to the bc-plane with hydrogen sulfate anions at the core and aminopentamide cations on the outside. There is a strong charge-assisted O49–H53∙∙∙O52 hydrogen bond between the hydrogen sulfate anions. This hydrogen bond links the anions in a chain parallel to the b-axis. The cation forms a discrete N–H∙∙∙O hydrogen bond to the anion. The amide group also forms two weaker discrete hydrogen bonds to the anion. The three N–H∙∙∙O hydrogen bonds link the cations and anions into columns parallel to the b-axis. This commercial material from USP contained an unidentified impurity, the powder pattern of which could be indexed on a monoclinic unit cell. The powder pattern has been submitted to ICDD for inclusion in the Powder Diffraction File™ (PDF®).

36 MATERIALS SCIENCE↗
DOE OSTI2026

A Chemocatalytic Route to Stereoregular Poly(3-hydroxyhexanoate) and Its Statistical and Tri-Block Copolymers

Poly(3-hydroxyhexanoate) (P3HHx)-based poly(3- hydroxyalkanoate)s (PHAs) are biologically produced, commercially implemented PHAs, but little is known in the open literature about the structure and property characterizations of discrete, authentic homopolymer P3HHx, and its copolymers with poly(3- hydroxybutyrate) (P3HB). Here, we introduce a chemocatalytic route to stereoregular (both isotactic and syndiotactic) P3HHx and PHA copolymers with P3HB, including statistical copolymer P3HBHx with various levels of 3HHx incorporation and discrete hard−soft−hard triblock copolymer P3HB-b-P3HHx-b-P3HB, and present extensive characterizations of their structures, thermal properties, and mechanical performance. The synthesis is efficiently achieved by one-pot polymerization of eight-membered di(n-propyl) and dimethyl (for copolymerization) diolides catalyzed by chiral molecular catalysts. Notably, P3HBHx can be rapidly produced in quantitative yield and exhibits high molar mass (M n up to 551 kg/mol) as well as both high modulus (up to 1.39 GPa) and ductility (up to 445%), while P3HB-b-P3HHx-b-P3HB further extends application temperature windows by possessing a unique combination of low T g (−18 °C) and high T m (154 °C) values. These findings highlight the stereomicrostructural and architectural versatility of chemocatalytic routes to PHAs, which, in turn, can be utilized to largely tune the PHA thermal properties and mechanical performance.

Biopolymers↗
DOE OSTI2026

Optimal Design and Techno-Economic Analysis of 3D-Printed, Intensified Packings for Absorbers and Strippers in Solvent-Based CO 2 Capture

A potential technology for the CO 2 absorption process is utilizing intensified structured packing with embedded cooling/heating channels for continuous heat exchange, which can overcome limitations of discrete methods, such as discrete intercooling and centralized reboilers, to aid in reducing energy consumption and decreasing costs. This work investigates the modeling of intensified packing (IP) for the stripper tower, extending on previous work for the absorber, which distributes heat internally within the column, improving the thermodynamics for the solvent regeneration process. The model includes submodels for steam turbine extraction to produce steam at various qualities as well as a surrogate model for calculating steam enthalpy. A cost model for a plant-scale absorption capture process was developed, allowing for the design of the plant to be optimized, subject to minimizing capture cost using two different power plant flue gas sources. In this optimization, the placement of IP in both towers is optimized to balance the trade-off between enhanced heat transfer and reduced mass transfer volume. For natural gas combined cycle flue gas, the standard process configuration had a minimum cost of $\$$65.40/tonne CO 2 , and considering IP, the minimum capture cost is reduced to $\$$62.73/tonne, with utilization in the stripper column, which reduces yearly costs by up to $\$$2.67 MM/yr. Cooling the absorber through IP, or intercoolers, was only found to be beneficial at higher capture rates, with IP in both towers having a cost of capture of $\$$68.08/tonne at 99.9% capture, a reduction of $\$$12.64/tonne when using only intercoolers at the same capture rate. When capturing from pulverized-coal power plants, the minimum cost of capture when using IP in both towers is $\$$44.18/tonne (at 97% capture), while the standard configuration with and without intercoolers was $\$$45.69 and $\$$47.22 per tonne, respectively. This results in a reduction in yearly costs of $\$$16.98 MM/yr from the base-case configuration. At this higher CO 2 concentration, cooling in the absorber from the IP becomes extremely beneficial, reducing energy consumption by up to 6%.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2021

Adsorbate Partition Functions via Phase Space Integration: Quantifying the Effect of Translational Anharmonicity on Thermodynamic Properties

A new method for computing anharmonic thermophysical properties for adsorbates on metal surfaces is presented. Classical Monte Carlo phase space integration is performed to calculate the partition function for the motion of a hydrogen atom on Cu(111). Here, a minima-preserving neural network potential energy surface is used within the integration routine. Two different sampling schema for generating the training data are presented, and two different density functionals are used. The results are benchmarked against direct state counting results by using discrete variable representation. The phase space integration results are in excellent quantitative agreement with the benchmark results. Additionally, both the discrete variable representation and the phase space integration results confirm that the motion of H on Cu(111) is highly anharmonic. The results were applied to calculate the free energy of dissociative adsorption of H 2 and the resulting Langmuir isotherms at 400, 800, and 1200 K in a partial pressure range of 0–1 bar. It shows that the anharmonic effects lead to significantly higher predicted surface site fractions of hydrogen.

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

Probing the In Situ Pseudocapacitive Charge Storage in Ti 3 C 2 MXene Thin Films with X-ray Reflectivity

MXenes are a large family of two-dimensional materials that are attractive for energy storage due to their high-rate charging capabilities as well as for electrochemical actuators, water purification, and many other technologies. Ion intercalation during electrochemically driven charge and discharge processes is the fundamental process associated with MXene functionality, which we have characterized here using in situ and operando X-ray reflectivity (XRR). Experiments performed at the Advanced Photon Source at Argonne National Laboratory monitored the changes in the structure of a Ti 3 C 2 MXene film on a platinum current collector as a function of static applied potential between 0.3 and –0.7 V vs Ag/AgCl in an aqueous 0.1 M Li 2 SO 4 electrolyte. Negative potential sweeps lead to a contraction of 1.2 Å in the interlayer spacing and a loss of electron density between the layers, likely due to Li + ion insertion and water removal. The change in lattice spacing includes a continuous variation vs potential as well as an additional discrete contraction that occurs near –0.35 V that has the characteristics of a first-order transition. The continuous change in the MXene interlayer spacing is associated with the capacitive charge, while the discrete change in structure correlated to the weak feature in the cyclic voltammogram at –0.35 V can be interpreted as either a pseudocapacitive charging process or a potential-dependent change in capacity.

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

Direct Demonstration of Topological Stability of Magnetic Skyrmions via Topology Manipulation

Topological protection precludes a continuous deformation between topologically inequivalent configurations in a continuum. Motivated by this concept, magnetic skyrmions, topologically nontrivial spin textures, are expected to exhibit topological stability, thereby offering a prospect as a nanometer-scale nonvolatile information carrier. In real materials, however, atomic spins are configured as not continuous but discrete distributions, which raises a fundamental question if the topological stability is indeed preserved for real magnetic skyrmions. Answering this question necessitates a direct comparison between topologically nontrivial and trivial spin textures, but the direct comparison in one sample under the same magnetic fields has been challenging. Here we report how to selectively achieve either a skyrmion state or a topologically trivial bubble state in a single specimen and thereby experimentally show how robust the skyrmion structure is in comparison with the bubbles. We demonstrate that topologically nontrivial magnetic skyrmions show longer lifetimes than trivial bubble structures, evidencing the topological stability in a real discrete system. Our work corroborates the physical importance of the topology in the magnetic materials, which has hitherto been suggested by mathematical arguments, providing an important step toward ever-dense and more-stable magnetic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Friends or Foes: Fundamental Principles of Th-Organic Scaffold Chemistry Using Zr-Analogs as a Guide

The fundamental interest in actinide chemistry, particularly for the development of thorium-based materials, is experiencing a renaissance owing to the recent and rapidly growing attention to fuel cycle reactors, radiological daughters for nuclear medicine, and efficient nuclear stockpile development. Herein, we uncover fundamental principles of thorium chemistry on the example of Th-based extended structures such as metal–organic frameworks in comparison with the discrete systems and zirconium extended analogs, demonstrating remarkable over two-and-half-year chemical stability of Th-based frameworks as a function of metal node connectivity, amount of defects, and conformational linker rigidity through comprehensive spectroscopic and crystallographic analysis as well as theoretical modeling. Despite exceptional chemical stability, we report the first example of studies focusing on the reactivity of the most chemically stable Th-based frameworks in comparison with the discrete Th-based systems such as metal–organic complexes and a cage, contrasting multicycle recyclability and selectivity (>97%) of the extended structures in comparison with the molecular compounds. Altogether, the presented work not only establishes the conceptual foundation for evaluating the capabilities of Th-based materials but also represents a milestone for their multifaceted future and foreshadows their potential to shape the next era of actinide chemistry.

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

Advanced Newton methods for geodynamical models of Stokes flow with viscoplastic rheologies

Strain localization and resulting plasticity and failure play an important role in the evolution of the lithosphere. These phenomena are commonly modeled by Stokes flows with viscoplastic rheologies. The nonlinearities of these rheologies make the numerical solution of the resulting systems challenging, and iterative methods often converge slowly or not at all. Yet accurate solutions are critical for representing the physics. Moreover, for some rheology laws, aspects of solvability are still unknown. We study a basic but representative viscoplastic rheology law. The law involves a yield stress that is independent of the dynamic pressure, referred to as von Mises yield criterion. Two commonly used variants, perfect/ideal and composite viscoplasticity, are compared. We derive both variants from energy minimization principles, and we use this perspective to argue when solutions are unique. We propose a new stress-velocity Newton solution algorithm that treats the stress as an independent variable during the Newton linearization but requires solution only of Stokes systems that are of the usual velocity-pressure form. To study different solution algorithms, we implement 2-D and 3-D finite element discretizations, and we generate Stokes problems with up to 7 orders of magnitude viscosity contrasts, in which compression or tension results in significant nonlinear localization effects. Comparing the performance of the proposed Newton method with the standard Newton method and the Picard fixed-point method, we observe a significant reduction in the number of iterations and improved stability with respect to problem nonlinearity, mesh refinement, and the polynomial order of the discretization.

Rudi, Johann↗