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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS2022

Direct-imaging Discovery and Dynamical Mass of a Substellar Companion Orbiting an Accelerating Hyades Sun-like Star with SCExAO/CHARIS

We present the direct-imaging discovery of a substellar companion in orbit around a Sun-like star member of the Hyades open cluster. So far, no other substellar companions have been unambiguously confirmed via direct imaging around main-sequence stars in Hyades. The star HIP 21152 is an accelerating star as identified by the astrometry from the Gaia and Hipparcos satellites. We detected the companion, HIP 21152 B, in multiple epochs using the high-contrast imaging from SCExAO/CHARIS and Keck/NIRC2. We also obtained the stellar radialvelocity data from the Okayama 188 cm telescope. The CHARIS spectroscopy reveals that HIP 21152 B’s spectrum is consistent with the L/T transition, best fit by an early T dwarf. Our orbit modeling determines the semimajor axis and the dynamical mass of HIP 21152 B to be 17 +7.2 −3.8 and 27.8 +8.4 −5.4 M Jub , respectively. The mass ratio of HIP 21152 B relative to its host is ≈2%, near the planet/brown dwarf boundary suggested by recent surveys. Mass estimates inferred from luminosity-evolution models are slightly higher (33–42 MJup). With a dynamical mass and a well-constrained age due to the system’s Hyades membership, HIP 21152 B will become a critical benchmark in understanding the formation, evolution, and atmosphere of a substellar object as a function of mass and age. Our discovery is yet another key proof of concept for using precision astrometry to select direct-imaging targets.

Masayuki Kuzuhara↗
DOE OSTI2023

Scaling Wind Power Innovation Assessment for Rapid Energy Transition with Artificial Intelligence

Planning for energy system decarbonization requires new insights into the potential of renewable technologies, deployed at unprecedented scale, to meet urgent sustainability goals. However, limited scalability of current wind energy research tools restricts characterization of innovation impacts to isolated reference sites, challenging investment and decision making under rapid growth. We demonstrate the transformative potential of artificial intelligence (AI) to inform future technology advancement and energy systems design by leveraging a state-of-the-art surrogate model to conduct a series of fleet-wide wind plant layout optimizations for greater than 6,800 projected U.S. onshore buildout locations. We show how innovative wake steering technology can address an array of barriers to large-scale deployment and integration of wind power. Specifically, wake steering reduces required plant area by an average of 18% and could preserve upwards of 13,000 km2 for future greenfield deployment, potentially easing siting challenges associated with wind energy infrastructure. Further, by enabling reduced turbine spacing and increased energy production, flexible operations of wake steering improve levelized cost of energy, particularly for large plants and in land-constrained settings. Finally, optimizations that consider dynamic energy prices can deliver increased power production and revenue capture during high-value (often low-wind) periods, further bolstering plant economics. Our computationally efficient approach offers a pathway to accelerate nationwide geographic evaluation of innovative technologies.

graph neural networks↗
DOE OSTI2024

Electric Vehicles Charging Time Constrained Deliverable Provision of Secondary Frequency Regulation

Aggregation of electric vehicles (EVs) is a promising technique for providing secondary frequency regulation (SFR) in highly renewable energy-penetrated power systems. Equipped with energy storage devices, EV aggregation can provide reliable SFR. However, the main challenge is to guarantee reliable intra-interval SFR capacities and inter-interval delivery following the automatic generation control (AGC) signal. Furthermore, aggregated EV SFR provision will be further complicated by the EV charging time anxiety because SFR provision might extend EV's charging time. This paper proposes a deliverable EV SFR provision with a charging-time-constrained control strategy. First, a charging-time-constrained EV aggregation is proposed to address the uncertainty of EV capacity based on the state-space model considering the charging-time restriction of EV owners. Second, a real-time economic dispatch and time domain simulation (RTED-TDS) cosimulation framework is proposed to verify financial results and the dynamic performance of the EV SFR provision. Last, the proposed charging time-constrained EV aggregation is validated on the IEEE 39-bus system. In conclusion, the results demonstrate that with charging time-constrained EV aggregation, the dynamic performance of the system can be improved with a marginal increase in total cost. More importantly, the charging time constraint can be respected in the proposed SFR provision of the EV aggregation.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Transmitted wave measurements in cold sprayed materials under dynamic compression

Spray-formed materials have complex microstructures which pose challenges for microscale and mesoscale modeling. To constrain these models, experimental measurements of wave profiles when subjecting the material to dynamic compression are necessary. The use of a gas gun to launch a shock into a material is a traditional method to understand wave propagation and provide information of time-dependent stress variations due to complex microstructures. This data contains information on wave reverberations within a material and provides a boundary condition for simulation. Here we present measurements of the wavespeed and wave profile at the rear surface of tantalum, niobium, and a tantalum/niobium blend subjected to plate impact. Measurements of the Hugoniot elastic limit are compared to previous work and wavespeeds are compared to longitudinal sound velocity measurements to examine wave damping due to the porous microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Shifting groundwater fluxes in bedrock fractures: Evidence from stream water radon and water isotopes

Geologic features (e.g., fractures and alluvial fans) can play an important role in the locations and volumes of groundwater discharge and degree of groundwater-surface water (GW-SW) interactions. However, the role of these features in controlling GW-SW dynamics and streamflow generation processes are not well constrained. GW-SW interactions and streamflow generation processes are further complicated by variability in precipitation inputs from summer and fall monsoon rains, as well as declines in snowpack and changing melt dynamics driven by warming temperatures. Using high spatial and temporal resolution radon and water stable isotope sampling and a 1D groundwater flux model, we evaluated how groundwater contributions and GW-SW interactions varied along a stream reach impacted by fractures (fractured-zone) and below the fractured hillslope (non-fractured zone) in Coal Creek, a Colorado River headwater stream affected by summer monsoons. During early summer, groundwater contributions from the fractured zone dominated, but declined throughout the summer. Groundwater contributions from the non-fractured zone were constant throughout the summer and became proportionally more important later in the summer. We hypothesize that groundwater in the non-fractured zone is dominantly sourced from a high-storage alluvial fan at the base of a tributary that is connected to Coal Creek throughout the summer and provides consistent groundwater influx. Water isotope data revealed that Coal Creek responds quickly to incoming precipitation early in the summer, and summer precipitation becomes more important for streamflow generation later in the summer. We quantified the change in catchment dynamic storage and found it negatively related to stream water isotope values, and positively related to modeled groundwater discharge and the ratio of fractured zone to non-fractured zone groundwater. Importantly, we interpret these relationships as declining hydrologic connectivity throughout the summer leading to late summer streamflow supported predominantly by shallow flow paths, with variable response to drying from geologic features based on their storage. As groundwater becomes more important for sustaining summer flows, quantifying local geologic controls on groundwater inputs and their response to variable moisture conditions may become critical for accurate predictions of streamflow.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Introducing Kynema, an Open-Source Performance-Portable Flexible-Multibody-Dynamics Solver

In this talk we introduce Kynema, an open-source general flexible-multibody-dynamics solver that is well suited for simulating wind turbine structural dynamics. Kynema uses a Lie-group time integrator for constrained systems and runs on both CPUs and GPUs. Timing results for simulations are presented for the IEA 15-MW turbine with and without aerodynamic forces.

17 WIND ENERGY↗
NASA NTRS1994

Numerical methods for multibody systems

This article gives a brief summary of some results obtained by Nasser on modeling and simulation of inequality problems in multibody dynamics. In particular, the augmented Lagrangian method discussed here is applied to a constrained motion problem with impulsive inequality constraints. A fundamental characteristic of the multibody dynamics problem is the lack of global convexity of its Lagrangian. The problem is transformed into a convex analysis problem by localization (piecewise linearization), where the augmented Lagrangian has been successfully used. A model test problem is considered and a set of numerical experiments is presented.

Glowinski, Roland↗
NASA NTRS2006

MEASURING THE MASS OF 4UO900-40 DYNAMICALLY

Accurate measurements of neutron star masses are needed to constrain the equation of state of neutron star matter - of importance to both particle physics and the astrophysics of neutron stars - and to identify the evolutionary track of the progenitor stars that form neutron stars. The best measured values of the mass of 4UO900-40 (= Vela XR-l), 1.86 +/- 0.16 Msun (Barziv et al. 2001) and 1.93 +/- 0.20 Msun (Abubekerov et al. 2004), make it a leading candidate for the most massive neutron star known. The direct relationship between the maximum mass of neutron stars and the equation of state of ultra-dense matter makes 4UO900-40 an important neutron star mass to determine accurately. The confidence interval on previous mass estimates, obtained from observations that include parameters determined by non-dynamical methods, are not small enough to significantly restrict possible equations of state. We describe here a purely dynamical method for determining the mass of 4UO900-40, an X-ray pulsar, using the reprocessed UV pulses emitted by its BO.5Ib companion. One can derive the instantaneous radial velocity of each component by simultaneous X-ray and UV observations at the two quadratures of the system. The Doppler shift caused by the primary's rotational velocity and the illumination pattern of the X-rays on the primary, two of the three principal contributors to the uncertainty on the derived mass of the neutron star, almost exactly cancel by symmetry in this method. A heuristic measurement of the mass of 4UO900-40 using observations obtained previously with the High Speed Photometer on HST is given in Appendix A.

Dolan, J. F.↗
NASA NTRS1993

Goddard problem in presence of a dynamic pressure limit

The Goddard problem is that of maximizing the final altitude for a vertically ascending, rocket-powered vehicle under the influence of an inverse square gravitational field and atmospheric drag. The present paper deals with the effects of two additional constraints, namely, a dynamic pressure limit and specified final time. Nine different switching structures involving zero-thrust arcs, full-thrust arcs, singular-thrust arcs, and state-constrained arcs are obtained when the value of the dynamic pressure limit is varied between zero and infinity and the final time is specified between the minimum possible time within which all of the fuel can be burned and the natural final time that emerges for the problem with final time unspecified. For all points in the aforementioned domain of dynamic pressure limit and prescribed final time, the associated optimal switching structure is clearly identified. Finally, a simple intuitive feedback law is presented for the free time problem. For all values of prescribed dynamic pressure limit, this strategy yields a loss in final altitude of less than 3 percent with respect to the associated optimal solution.

Seywald, Hans↗
NASA NTRS1990

Application of adjoint operators to neural learning

A technique for the efficient analytical computation of such parameters of the neural architecture as synaptic weights and neural gain is presented as a single solution of a set of adjoint equations. The learning model discussed concentrates on the adiabatic approximation only. A problem of interest is represented by a system of N coupled equations, and then adjoint operators are introduced. A neural network is formalized as an adaptive dynamical system whose temporal evolution is governed by a set of coupled nonlinear differential equations. An approach based on the minimization of a constrained neuromorphic energylike function is applied, and the complete learning dynamics are obtained as a result of the calculations.

Barhen, J.↗
DOE OSTI2020

An error-in-constitutive equations strategy for topology optimization for frequency-domain dynamics

This paper presents a topology optimization formulation for frequency-domain dynamics to reduce solution dependence upon initial guess and considered loading conditions. Due to resonance phenomena in undamped steady-state dynamics, objectives measuring dynamic response possess many local minima that may represent poor solutions to a design problem, an issue exacerbated for design with respect to multiple frequencies. In this work, we propose an extension of the modified error-in-constitutive-equations (MECE) method, used previously in material identification inverse problems, as a new approach for frequency-domain dynamics topology optimization to mitigate these issues. The main idea of the proposed framework is to incorporate an additional penalty-like term in the objective function that measures the discrepancy in the constitutive relations between stresses and strains and between inertial forces and displacements. Then, the design problem is cast within a PDE-constrained optimization formulation in which we seek displacements, stresses, inertial forces, and a density-field solution that minimize our new objective subject to conservation of linear momentum plus some additional constraints. We show that this approach yields superior designs to conventional gradient-based optimization approaches that solely use a functional of displacements as the objective, while strictly enforcing the constitutive equations. The MECE strategy integrates into a density-based topology optimization scheme for void–solid or two-phase material structural design. We highlight the merits of our approach in a variety of scenarios for direct frequency response design, considering multiple frequency load cases and structural objectives.

42 ENGINEERING↗
DOE OSTI2022

SPARC: Structural properties associated with residue constraints

SPARC facilitates the generation of plausible hypotheses regarding underlying biochemical mechanisms by structurally characterizing protein sequence constraints. Such constraints appear as residues co-conserved in functionally related subgroups, as subtle pairwise correlations (i.e., direct couplings), and as correlations among these sequence features or with structural features. SPARC performs three types of analyses. First, based on pairwise sequence correlations, it estimates the biological relevance of alternative conformations and of homomeric contacts, as illustrated here for death domains. Second, it estimates the statistical significance of the correspondence between directly coupled residue pairs and interactions at heterodimeric interfaces. Third, given molecular dynamics simulated structures, it characterizes interactions among constrained residues or between such residues and ligands that: (a) are stably maintained during the simulation; (b) undergo correlated formation and/or disruption of interactions with other constrained residues; or (c) switch between alternative interactions. We illustrate this for two homohexameric complexes: the bacterial enhancer binding protein (bEBP) NtrC1, which activates transcription by remodeling RNA polymerase (RNAP) containing σ 54 , and for DnaB helicase, which opens DNA at the bacterial replication fork. Based on the NtrC1 analysis, we hypothesize possible mechanisms for inhibiting ATP hydrolysis until ADP is released from an adjacent subunit and for coupling ATP hydrolysis to restructuring of σ 54 binding loops. Based on the DnaB analysis, we hypothesize that DnaB ‘grabs’ ssDNA by flipping every fourth base and inserting it into cavities between subunits and that flipping of a DnaB-specific glutamine residue triggers ATP hydrolysis.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Faster Form of Electron Magnetic Reconnection with a Finite Length X-Line

Observations in Earth’s turbulent magnetosheath downstream of a quasiparallel bow shock reveal a prevalence of electron-scale current sheets favorable for electron-only reconnection where ions are not coupled to the reconnecting magnetic fields. In small-scale turbulence, magnetic structures associated with intense current sheets are limited in all dimensions. And since the coupling of ions are constrained by a minimum length scale, the dynamics of electron reconnection is likely to be 3D. Here, both 2D and 3D kinetic particle-in-cell simulations are used to investigate electron-only reconnection, focusing on the reconnection rate and associated electron flows. A new form of 3D electron-only reconnection spontaneously develops where the magnetic X-line is localized in the out-of-plane (z) direction. The consequence is an enhancement of the reconnection rate compared with two dimensions, which results from differential mass flux out of the diffusion region along z, enabling a faster inflow velocity and thus a larger reconnection rate. This outflow along z is due to the magnetic tension force in z just as the conventional exhaust tension force, allowing particles to leave the diffusion region efficiently along z unlike the 2D configuration.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Formation of Magnetically Truncated Accretion Disks in 3D Radiation-transport Two-temperature GRMHD Simulations

Multiwavelength observations suggest that the accretion disk in the hard and intermediate states of X-ray binaries (XRBs) and active galactic nucleus transitions from a cold, thin disk at large distances into a hot, thick flow close to the black hole (BH). However, the formation, structure, and dynamics of such truncated disks are poorly constrained due to the complexity of the thermodynamic, magnetic, and radiative processes involved. We present the first radiation-transport two-temperature general relativistic magnetohydrodynamic (GRMHD) simulations of truncated disks radiating at ~35% of the Eddington luminosity with and without large-scale poloidal magnetic flux. We demonstrate that when a geometrically thin accretion disk is threaded by large-scale net poloidal magnetic flux, it self-consistently transitions at small radii into a two-phase medium of cold gas clumps floating through a hot, magnetically dominated corona. This transition occurs at a well-defined truncation radius determined by the distance out to which the disk is saturated with magnetic flux. The average ion and electron temperatures in the semiopaque corona reach, respectively, $T_i$ ≳ 10 10 K and $T_e$ ≳ 5 × 10 8 K. The system produces radiation, powerful collimated jets, and broader winds at the total energy efficiency exceeding 90%, the highest ever energy extraction efficiency from a spinning BH by a radiatively efficient flow in a GRMHD simulation. This is consistent with jetted ejections observed during XRB outbursts. The two-phase medium may naturally lead to broadened iron line emission observed in the hard state.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1990

Temperature conditions for chondrule formation

An attempt is made to constrain the chondrule-forming process from dynamic crystallization experiments on chondrule analogs and from correlations between texture and bulk composition in natural chondrules. Liquidus temperatures for chondrules from unequilibrated chondrites are calculated using Herzberg's (1979) method. The paper then infers whether each chondrule formed above, at, or below its liquidus based on textures produced in experiments. A range of temperatures to which chondrules may have been initially heated is then derived from the texture-composition-temperature relationships. Finally, the role of dust seeding, or external heterogeneous nucleation, in producing chondrule textures is examined.

Hewins, Roger H.↗
NASA NTRS1991

CVs and millisecond pulsar progenitors in globular clusters

The recent discovery of a large population of millisecond pulsars in globular clusters, together with earlier studies of both low luminosity X-ray sources and LMXBs in globulars, suggest there should be significant numbers of CVs in globulars. Although they have been searched for without success in selected cluster X-ray source fields, systematic surveys are lacking and would constrain binary production and both stellar and dynamical evolution in globular clusters. We describe the beginnings of such a search, using narrow band H-alpha imaging, and the sensitivities it might achieve.

Grindlay, J. E.↗
NASA NTRS2004

Meridional Variations of C2H2 and C2H6 in Jupiter's Atmosphere from Cassini CIRS Infrared Spectra

Hydrocarbons such as acetylene (C2H2) and ethane (C2H6) are important tracers in Jupiter's atmosphere, constraining our models of the chemical and dynamical processes. However, our knowledge of the vertical and meridional variations of their abundances has remained sparse. During the flyby of the Cassini spacecraft in December 2000, the Composite Infrared Spectrometer (CIRS) instrument was used to map the spatial variation of emissions from 10-1400 cm(sup -1) (1000-7 microns). In this paper we analyze a zonally-averaged set of CIRS spectra taken at the highest (0.5 cm(sup -1)) resolution, to infer atmospheric temperatures in the stratosphere at 0.5-20 mbar via the v4 band of CH4, and in the troposphere at 150-400 mbar, via the H2 absorption at 600-800 cm(sup -1). Simultaneously, we retrieve the abundances of C2H2 and C2H6 via the v5 and vg bands respectively. Tropospheric absorption and stratospheric emission are highly anti-correlated at the CIRS resolution, introducing a non-uniqueness into the retrievals, such that vertical gradient and column abundance cannot both be found without additional constraints. Assuming profile gradients from photochemical calculations, we show that the column abundance of C2H2 decreases sharply towards the poles by a factor approximately 4, while C2H6 is unchanged in the north and increasing in the south, by a factor approximately 1.8. An explanation for the meridional trends is proposed in terms of a combination of photochemistry and dynamics. Poleward, the decreasing UV flux is predicted to decrease the abundances of C2H2 and C2H6 by factors 2.7 and 3.5 respectively at a latitude 70 deg. However, the lifetime of C2H6 in the stratosphere (5 x 10(exp 9)) is much longer than the dynamical timescale for meridional motions inferred from SL-9 debris (5 x 10(exp 8 s)), and therefore the constant or rising abundance towards high latitudes likely indicates that meridional mixing dominates over photochemical effects. For C2H2, the opposite occurs, with the relatively short photochemical lifetime (3 x 10(exp 7 s)), compared to meridional mixing times, ensuring that the expected photochemical trends are visible.

Nixon, C. A.↗