Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constrained dynamics”

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

Fast methods for multisite charge transfer. Processes II. Analytic nuclear gradients and nonadiabatic dynamics for cCASSCF(1,n) and cCASSCF(2n-1,n) wavefunctions

In this work we derive and implement analytic nuclear gradients and derivative couplings for a constrained complete active space self-consistent field with a small active space designed to model electron or hole transfer. Using a Lagrangian formalism, we are able to differentiate both the CASSCF energy and the constraint (which is required for smooth surfaces over a wide range of parameter space), and the resulting efficient algorithm can be immediately applied to nonadiabatic dynamics simulations of charge transfer processes. Here, we run initial surface-hopping simulations of a proton coupled electron transfer event for a phenoxyl–phenol system.

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

Evolution of radiative thermal instability in a confined medium

Thermally bistable fluid tends to self-organize into clouds of hot and cold material, which are internally uniform and separated by thin conduction fronts. The evolution of these clouds has been studied for isobaric systems, but when pressure is instead treated as a dynamical quantity and allowed to evolve self-consistently, fundamentally different dynamics appear. Such a treatment is necessary in some laboratory plasmas, whose volume is constrained but whose pressure can vary. Here, solutions are derived for the evolution of clouds, accounting for pressure variation and interactions between conduction fronts. Additional stable configurations and secondary instabilities are derived, which may be relevant to fusion plasmas and to the study of photoionized plasma in the laboratory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Charge Symmetry Violation Quark Distributions Via Precise Measurement of p+/p- Rations in Semi-Inclusive Deep Inelastic Scattering

Charge symmetry in the nucleon parton distributions assumes the distribution of quarks in the proton is related to those in the neutron. Indirect experimental evidence constrains Charge Symmetry Violation (CSV) to be less than 9%. In Quantum Chromo-Dynamics (QCD), charge symmetry is broken by the mass difference between the up and down quarks. CSV in the valence region can be extracted from precision measurements of the cross section ratio of charged pion production in semi-inclusive deep inelastic scattering (SIDIS) from deuterium. Such measurements were proposed and carried out in experiment E12-09-002 at Jefferson Lab. The experiment was conducted in Hall C from fall 2018 to spring 2019 using the upgraded 10.6 GeV incident electron beam. In this SIDIS experiment, charged pions are detected in coincidence with scattered electrons covering the four-momentum transfer of the virtual photon Q2 from 4 to 5.5 GeV2, the Bjorken variable x for 0.35-0.65, and the fraction of energy transfer carried by the outgoing pion z from 0.4-0.7.

Jia, Shuo↗
DOE OSTI2022

The 3D Kinematics of the Orion Nebula Cluster: NIRSPEC-AO Radial Velocities of the Core Population

The kinematics and dynamics of stellar and substellar populations within young, still-forming clusters provide valuable information for constraining theories of formation mechanisms. Using Keck II NIRSPEC+AO data, we have measured radial velocities for 56 low-mass sources within 4' of the core of the Orion Nebula Cluster (ONC). We also remeasure radial velocities for 172 sources observed with SDSS/APOGEE. These data are combined with proper motions measured using HST ACS/WFPC2/WFC3IR and Keck II NIRC2, creating a sample of 135 sources with all three velocity components. The velocities measured are consistent with a normal distribution in all three components. We measure intrinsic velocity dispersions of (${\sigma }_{{v}_{\alpha }}$, ${\sigma }_{{v}_{\delta }}$, ${\sigma }_{{v}_{r}}$) = (1.64 ± 0.12, 2.03 ± 0.13, ${2.56}_{-0.17}^{+0.16}$) km s –1 . Our computed intrinsic velocity dispersion profiles are consistent with the dynamical equilibrium models from Da Rio et al. (2014) in the tangential direction but not in the line-of-sight direction, possibly indicating that the core of the ONC is not yet virialized, and may require a nonspherical potential to explain the observed velocity dispersion profiles. We also observe a slight elongation along the north–south direction following the filament, which has been well studied in previous literature, and an elongation in the line-of-sight to tangential velocity direction. These 3D kinematics will help in the development of realistic models of the formation and early evolution of massive clusters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Pauli energy contribution to the nucleus-nucleus interaction

Background: The Pauli exclusion principle plays a crucial role as a building block of many-body quantal systems comprised of fermions. It also induces a “Pauli repulsion” in the interaction between di-nuclear systems. It has been shown in [Phys. Rev. C 95, 031601(R) (2017)] that the Pauli repulsion widens the nucleus-nucleus potential barrier, thus hindering sub-barrier fusion. Purpose: To investigate the proton and neutron contributions to the Pauli repulsion, both in the bare potential neglecting shape polarization and transfer between the reactants, as well as in the dynamical potential obtained by accounting for such dynamical rearrangements. Methods: As the basis of our study we utilize the Pauli kinetic energy (PKE) obtained by studying the nuclear localization function (NLF). Recently this approach has been generalized to incorporate all of the dynamical and time-odd terms present in the nuclear energy density functional. This approach is employed in the density constrained frozen Hartree-Fock (DCFHF) and in the density constrained time-dependent Hartree-Fock (DC- TDHF) microscopic methods. Results: The PKE spatial distribution shows that a repulsion occurs in the neck between the nuclei when they first touch. Inside the barrier, neutrons can contribute significantly more to the Pauli repulsion in neutron-rich systems. Dynamical effects tend to lower the Pauli repulsion near the barrier. Proton and neutron dynamical contributions to the PKE significantly differ inside the barrier for asymmetric collisions, which is interpreted as an effect of multinucleon transfer. Conclusions: The PKE is shown to make a significant contribution to nuclear interaction potentials. Protons and neutrons can play very different roles in both the bare potential and in the dynamical rearrangement. Further microscopic studies are required to better understand the role of transfer and to investigate the effect of pairing and deformation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Probabilistic learning on manifolds constrained by nonlinear partial differential equations for small datasets

A novel extension of the Probabilistic Learning on Manifolds (PLoM) is presented. It makes it possible to synthesize solutions to a wide range of nonlinear stochastic boundary value problems described by partial differential equations (PDEs) for which a stochastic computational model (SCM) is available and which depend on a vector-valued random control parameter. The cost of a single numerical evaluation of this SCM is assumed to be such that only a limited number of points can be computed for constructing the training dataset (small data). Each point of the training dataset is made up of realizations from a vector-valued stochastic process (the stochastic solution) and the associated random control parameter on which it depends. The presented PLoM constrained by PDE allows for generating a large number of learned realizations of the stochastic process and its corresponding random control parameter. These learned realizations are generated so as to minimize the vector-valued random residual of the PDE in the mean-square sense. Appropriate novel methods are developed to solve this challenging problem. Three applications are presented. The first one is a simple uncertain nonlinear dynamical system with a nonstationary stochastic excitation. The second one concerns the 2D nonlinear unsteady Navier–Stokes equations for incompressible flows in which the Reynolds number is the random control parameter. Here, the last one deals with the nonlinear dynamics of a 3D elastic structure with uncertainties. The results obtained make it possible to validate the PLoM constrained by stochastic PDE but also provide further validation of the PLoM without constraint.

Machine learning↗
DOE OSTI2025

Seismicity-constrained fault detection and characterization with a multitask machine learning model

Geological fault detection and characterization are crucial for understanding subsurface dynamics across scales. While methods for fault delineation based on either seismicity location analysis or seismic image reflector discontinuity are well-established, a systematic approach that integrates both data types remains absent. We develop a novel machine learning model that unifies seismic reflector images and seismicity location information to automatically identify geological faults and characterize their geometrical properties. The model encodes a seismic image and a seismicity location image separately, and fuses the encoded features with a spatial-channel attention fusion module to improve the learning of important features in both inputs. We design an automated strategy to generate high-quality synthetic training data and labels. To improve the realism of the seismicity location image, we include random seismicity noise and missing seismicity location associated with some of the faults. We validate the model’s efficacy and accuracy using synthetic data examples and two field data examples. Moreover, we show that fine-tuning the trained model with a small, domain-specific dataset enhances its fidelity for field data applications. The results demonstrate that integrating seismicity location and seismic images into a unified framework allows the end-to-end neural network to achieve higher fidelity and accuracy in delineating subsurface faults and their geometrical properties compared with image-only fault detection methods. Our approach offers an adaptive data-driven tool for geological fault characterization and seismic hazard mitigation, bridging the gap between seismicity location and image-based fault detection methods.

58 GEOSCIENCES↗
DOE OSTI2025

Beam optics ramping in underconstrained lattice design: Application to Electron-Ion Collider hadron storage ring cooling section

This paper presents the lattice design and optics ramping strategy for the cooling section of the hadron storage ring at the electron-ion collider. The main challenge is that available tuning knobs exceed beam-optics constraints. Independently optimized injection and top-energy optics often yield disconnected solutions, making interpolation impossible. To address this, we propose two new methods. The first is an intermediate-penalty scheme that ensures ramping path continuity by penalizing constraint violations at intermediate states. The second is a continuation-based approach that adapts high-energy optics to low energy, guided by an adaptive weighting scheme to balance injection and ramping constraints. The solutions meet all beam dynamics and hardware limits. The two methods offer a general strategy for ramping in systems where the solution space is under-constrained and the starting and target configurations are far apart.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

DUET - Distribution of Understory using Elliptical Transport: A mechanistic model of leaf litter and herbaceous spatial distribution based on tree canopy structure

Heterogeneity in surface fuels produced by overstory trees and understory vegetation is a major driver of fire behavior and ecosystem dynamics. Previous attempts at predicting tree leaf and needle litter accumulation over time have been constrained in scope to probabilistic models that consider a limited number of key factors influencing tree litter dispersal patterns and decomposition processes. Here we present a mechanistic model for estimating variation in surface fuels called the Distribution of Understory using Elliptical Transport (DUET). DUET uses a pre-generated voxelated canopy array and distributes the leaf litter from the trees in an elliptical shape based on tree species characteristics, wind data, and location-specific features. DUET then calculates grass growth between trees and around the litter patterns and sets moisture levels based on relative humidity. In this study, we use DUET to simulate spatial distributions of leaf litter and illustrate model sensitivity to several parameters, providing inputs for testing fuel spatial heterogeneity impacts on fire behavior using process-based models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Measurements of nonequilibrium interatomic forces using time-domain x-ray scattering

Here, we demonstrate an experimental approach to determining the excited-state interatomic forces using femtosecond x-ray pulses from an x-ray free-electron laser. We determine experimentally the excited-state interatomic forces that connect photoexcited carriers to the nonequilibrium lattice dynamics in the prototypical Peierls-distorted material, bismuth. The forces are obtained by a constrained least-squares fit of a pairwise interatomic force model to the excited-state phonon dispersion relation as measured by the time- and momentum-resolved x-ray diffuse scattering. We find that photoexcited carriers weaken predominantly the nearest-neighbor forces, which drives the measured softening of the transverse acoustic modes throughout the Brillouin zone as well as the zone-center A 1g optical mode. This demonstrates a bond-selective approach to measuring electron-phonon coupling relevant to a broad range of photoinduced phase transitions and transient light-driven states in quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Energy-Energy Correlators for Precision QCD

In this contribution to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021) we review recent progress in the evaluation and application of the Energy-Energy Correlator (EEC) event shape observable in $e^+e^⁻$ annihilation, hadronic collisions, and deep inelastic scattering. The importance of EEC as a precision probe of the perturbative and non perturbative aspects of QCD dynamics is emphasized. It can be used to extract the strong coupling constant and to constrain TMD distribution functions. Closely related energy-correlation shape variables have also been used to tag boosted objects produced in high energy collisions. The opportunities to study EEC at the future Electron-Ion Collider are also highlighted.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Coordination and divergence in community assembly processes across co-occurring microbial groups separated by cell size

Setting the pace of life and constraining the role of members in food webs, body size can affect the structure and dynamics of communities across multiple scales of biological organization (e.g., from the individual to the ecosystem). However, its effects on shaping microbial communities, as well as underlying assembly processes, remain poorly known. Here, we analyzed microbial diversity in the largest urban lake in China and disentangled the ecological processes governing microbial eukaryotes and prokaryotes using 16S and 18S amplicon sequencing. We found that pico/nano-eukaryotes (0.22–20 μm) and micro-eukaryotes (20–200 μm) showed significant differences in terms of both community composition and assembly processes even though they were characterized by similar phylotype diversity. We also found scale dependencies whereby micro-eukaryotes were strongly governed by environmental selection at the local scale and dispersal limitation at the regional scale. Interestingly, it was the micro-eukaryotes, rather than the pico/nano-eukaryotes, that shared similar distribution and community assembly patterns with the prokaryotes. This indicated that assembly processes of eukaryotes may be coupled or decoupled from prokaryotes’ assembly processes based on eukaryote cell size. While the results support the important influence of cell size, there may be other factors leading to different levels of assembly process coupling across size classes. Additional studies are needed to quantitatively parse the influence of cell size versus other factors as drivers of coordinated and divergent community assembly processes across microbial groups. Regardless of the governing mechanisms, our results show that there are clear patterns in how assembly processes are coupled across sub-communities defined by cell size. These size-structured patterns could be used to help predict shifts in microbial food webs in response to future disturbance.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

The Effect of Acceleration on the Separation Force in Constrained-Surface Stereolithography

Constrained-surface-based stereolithography has recently attracted much attention from both academic and industrial communities. Despite numerous experimental, numerical and theoretical efforts, the fundamental need to reduce the separation force between the newly cured part and constrained surface has not yet been completely solved. In this paper, we develop a fluid dynamics approach, proposed in our previous work, to theoretically model the separation force in 3D printing of a cylindrical part for flat and patterned windows. We demonstrate the possibility of separation force reduction with an accelerated movement of the printing platform. In particular, we investigate behaviors of transient parameter, its reduction rate, and separation force reduction with respect to elevation speed and time. The proposed approach involves deceleration and acceleration stages and allows to achieve the force reduction for the entire printing process. Finally, we provide implicit analytical solutions for time moments when switching between the stages can be done without noticeable increase of separation force and explicit expression for separation force in case of patterned window.

Gritsenko, Dmitry (ORCID:0000000248181512)↗
DOE OSTI2025

GR-Athena++: General-relativistic Magnetohydrodynamics Simulations of Neutron Star Spacetimes

We present the extension of GR-Athena++ to general-relativistic magnetohydrodynamics (GRMHD) for applications to neutron star spacetimes. The new solver couples the constrained transport implementation of Athena++ to the Z4c formulation of the Einstein equations to simulate dynamical spacetimes with GRMHD using oct-tree adaptive mesh refinement. We consider benchmark problems for isolated and binary neutron star spacetimes demonstrating stable and convergent results at relatively low resolutions and without grid symmetries imposed. The code correctly captures magnetic field instabilities in nonrotating stars with total relative violation of the divergence-free constraint of 10 –16 . It handles evolutions with a microphysical equation of state and black hole formation in the gravitational collapse of a rapidly rotating star. For binaries, we demonstrate correctness of the evolution under the gravitational radiation reaction and show convergence of gravitational waveforms. We showcase the use of adaptive mesh refinement to resolve the Kelvin–Helmholtz instability at the collisional interface in a merger of magnetised binary neutron stars. GR-Athena++ shows strong scaling efficiencies above 80% in excess of 10 5 CPU cores and excellent weak scaling is shown up to ~5 × 10 5 CPU cores in a realistic production setup. GR-Athena++ allows for the robust simulation of GRMHD flows in strong and dynamical gravity with exa-scale computers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Applications of Multibody for Everybody (M4E) in Marine Energy

Multibody for Everybody (M4E) is an open-source symbolic dynamics modeling framework designed to automate the derivation of equations of motion and simulation of constrained multibody systems using the joint coordinate formulation.

16 TIDAL AND WAVE POWER↗
DOE OSTI2020

Formation of helices with controllable chirality in gel-fiber composites

While microscale helices with specified chirality are pervasive in biological systems, driving purely synthetic materials to self-organize into such useful structures remains a significant challenge. Through theory and simulations, we model the behavior of initially flat ribbons of gel-fiber composites and show that these materials self-organize into three-dimensional helices with the application of heat. We specifically focus on thermo-responsive gels, which shrink when heated to a point near the lower critical solubility temperature (LCST). Within the composite, regions of the LCST gel that lie far from the fibers can readily shrink in response to the increased temperature; however, near the fibers the network is inhibited from undergoing the heat-induced collapse. This competition between the constrained and the unconstrained regions of the heated sample plays a vital role in the structural evolution of the material. Additionally, we find that this dynamic process can be controlled through the careful placement of the fibers on the outer surface of the gel. We specifically identify fiber arrangements that lead to the formation of helical structures with controlled chirality. Moreover, we isolate scenarios where constraining the movement of one or both the free ends of the sample promotes and stabilizes the formation of helical structures with specified handedness. Such structurally tailored composites are useful for creating mechanically robust actuators, and “appendages” for soft robots that can controllably curl and uncurl with variations in temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Connecting the Light Curves of Type IIP Supernovae to the Properties of Their Progenitors

Abstract Observations of core-collapse supernovae (CCSNe) reveal a wealth of information about the dynamics of the supernova ejecta and its composition but very little direct information about the progenitor. Constraining properties of the progenitor and the explosion requires coupling the observations with a theoretical model of the explosion. Here we begin with the CCSN simulations of Couch et al., which use a nonparametric treatment of the neutrino transport while also accounting for turbulence and convection. In this work we use the SuperNova Explosion Code to evolve the CCSN hydrodynamics to later times and compute bolometric light curves. Focusing on Type IIP SNe (SNe IIP), we then (1) directly compare the theoretical STIR explosions to observations and (2) assess how properties of the progenitor’s core can be estimated from optical photometry in the plateau phase alone. First, the distribution of plateau luminosities ( L 50 ) and ejecta velocities achieved by our simulations is similar to the observed distributions. Second, we fit our models to the light curves and velocity evolution of some well-observed SNe. Third, we recover well-known correlations, as well as the difficulty of connecting any one SN property to zero-age main-sequence mass. Finally, we show that there is a usable, linear correlation between iron core mass and L 50 such that optical photometry alone of SNe IIP can give us insights into the cores of massive stars. Illustrating this by application to a few SNe, we find iron core masses of 1.3–1.5 M ⊙ with typical errors of 0.05 M ⊙ . Data are publicly available online on Zenodo: doi: 10.5281/zenodo.6631964 .

79 ASTRONOMY AND ASTROPHYSICS↗