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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Seasonal-to-Interannual Variability in Antarctic Sea-Ice Dynamics, and Its Impact on Surface Fluxes and Water Mass Production

Strong seasonal and interannual signals in Antarctic bottom-water outflow remain unexplained yet are highly correlated with anomalies in net sea-ice growth in coastal polynyas. The mechanisms responsible for driving salination and replenishment and rejuvenation of the dense shelf "source" waters likely also generate pulses of bottom water outflow. The objective of this research is to investigate time-scales of variability in the dynamics of sea-ice in the Southern Ocean in order to determine the primary sites for production of dense shelf waters. We are using a merged satellite/buoy sea-ice motion data set for the period 1978-present day to compute the dynamics of opening and closing of coastal polynyas over the continental shelf. The Ocean Circulation and Climate Advanced Model (OCCAM) ocean general circulation model with coupled sea-ice dynamics is presently forced using National Center for Environmental Prediction (NCEP) data to simulate fluxes and the salination impact of the ocean shelf regions. This work is relevant in the context of measuring the influence of polar sea-ice dynamics upon polar ocean characteristics, and thereby upon global thermohaline ocean circulation. Interannual variability in simulated net freezing rate in the Southern Weddell Sea is shown for the period 1986-1993. There is a pronounced maximum of ice production in 1988 and minimum in 1991 in response to anomalies in equatorward meridional wind velocity. This follows a similar approximate 8-year interannual cycle in Sea Surface Temperature (SST) and satellite-derived ice-edge anomalies reported elsewhere as the "Antarctic Circumpolar Wave." The amplitude of interannual fluctuations in annual net ice production are about 40% of the mean value, implying significant interannual variance in brine rejection and upper ocean heat loss. Southward anomalies in wind stress induce negative anomalies in open water production, which are observed in passive microwave satellite images. Thus, cycles of enhanced poleward wind stress reduce ice growth by compacting the ice along the coastline and closing open water in leads and polynyas. Model simulations confirm that years of low ice production, such as 1991, coincide with years of lower than normal bottom water outflow. Future plans include the assimilation of satellite ice concentrations and ice drift dynamics to more accurately constrain boundary conditions in the model.

Drinkwater, Mark R.↗
NASA NTRS2013

Explicit Low-Thrust Guidance for Reference Orbit Targeting

The problem of a low-thrust spacecraft controlled to a reference orbit is addressed in this paper. A simple and explicit low-thrust guidance scheme with constrained thrust magnitude is developed by combining the fundamental equations of motion for constrained systems from analytical dynamics with a Lyapunov-based method. Examples are given for a spacecraft controlled to a reference trajectory in the circular restricted three body problem.

low-thrust guidance↗
DOE OSTI2025

Physics-Informed Graph Neural Networks for Collaborative Dynamic Reconfiguration and Voltage Regulation in Unbalanced Distribution Systems

Network reconfiguration has long been employed as a strategic approach to minimize power distribution system losses and effectively regulate voltage levels. Tap-changing voltage regulators are also critical for controlling bus voltages, especially in accommodating the increasing integration of distributed energy resources (DERs) with intermittent outputs. This paper introduces novel methodologies to address the challenges of dynamic reconfiguration and optimal tap setting in unbalanced three-phase distribution systems. We propose an approximated mixed-integer quadratically constrained program (MIQCP) to model dynamic reconfiguration, along with a pioneering formulation for voltage regulator (VR) tap-setting based on Special Ordered Set type 1 (SOS1). To mitigate computational complexity, we propose a physics-informed spatial-temporal graph convolutional network (STGCN) with an integrated link classifier. The proposed approach enables efficient solution generation by fixing specific variables in the MIQCP instance and solving the simplified sub-MIP using an MIP solver. Numerical studies demonstrate the superior prediction accuracy of our STGCN model compared to baseline neural network models, resulting in reduced DER curtailment and voltage deviation with shorter computation time.

dynamic reconfiguration↗
DOE OSTI2026

Real-Space Constrained Density Functional Theory Investigation of Site-Specific, Interfacial Charge Recombination Dynamics Across the Au Nanoparticle/TiO 2 Heterojunction

Au nanoparticle (NP)/TiO 2 heterojunction is a representative system to study interfacial charge transfer in photocatalysis and photovoltaics, where suppressing recombination from TiO 2 to Au can enhance hot carrier extraction. We apply real-space constrained density functional theory (CDFT) with Marcus theory to quantify charge recombination time scales across Au/TiO 2 . This approach enables direct control and visualization of charge-separated states, aligning with site-specific probes like time-resolved X-ray photoelectron spectroscopy (trXPS). We find that the charge-separated state features a bipolaron, with recombination dominated by TiO 2 LUMO to Au HOMO transitions, primarily at interfacial Au sites. Marcus rate predictions are benchmarked with surface hopping methods, quantifying differences in time scales and computational efficiency. Lastly, we examine how the Au cluster size affects the free energy change (ΔG) and reorganization energy (λ), explaining trends in closed-shell systems and highlighting challenges for open-shell extrapolations. Overall, CDFT + Marcus theory provides efficient, mechanistically transparent interfacial charge transfer modeling, and we clearly defined its applicability and limitation.

Glenna, Drew M. [Univ. of Idaho, Idaho Falls, ID (↗
NASA NTRS1982

Constrained and unconstrained modes - Some modeling aspects of flexible spacecraft

Spacecraft that are partially rigid and partially flexible may be dynamically modeled in terms of either 'constrained' modes of vibration, for which the rigid part is held motionless, or the 'unconstrained' modes of the whole vehicle. A model 'completeness index' is defined for each expansion and used as a truncation criterion. Using as an example a large flexible platform with a rigid body at the center, it is shown that as the rigid portion of the vehicle becomes very small, a great many constrained modes are needed to achieve satisfactory model completeness. Otherwise, one finds an erroneous distribution of angular momentum among the modes. It is also shown that the importance of unconstrained modes is not necessarily ordered by frequency.

Hablani, H. B.↗
DOE OSTI2023

An inexact semismooth Newton method with application to adaptive randomized sketching for dynamic optimization

In many applications, one can only access the inexact gradients and inexact hessian times vector products. Thus it is essential to consider algorithms that can handle such inexact quantities with a guaranteed convergence to solution. An inexact adaptive and provably convergent semismooth Newton method is considered to solve constrained optimization problems. In particular, dynamic optimization problems, which are known to be highly expensive, are the focus. A memory efficient semismooth Newton algorithm is introduced for these problems. The source of efficiency and inexactness is the randomized matrix sketching. Further, applications to optimization problems constrained by partial differential equations are also considered.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2012

Direct Detection and Orbit Analysis of the Exoplanets HR 8799 bcd from Archival 2005 Keck/NIRC2 Data

We present previously unpublished July 2005 H-band coronagraphic data of the young, planet-hosting star HR 8799 from the newly-released Keck/NIRC2 archive. Despite poor observing conditions, we detect three of the planets (HR 8799 bcd), two of them (HR 8799 bc) without advanced image processing. Comparing these data with previously published 1998-2011 astrometry and that from re-reduced October 2010 Keck data constrains the orbits of the planets. Analyzing the planets' astrometry separately, HR 8799 d's orbit is likely inclined at least 25 deg from face-on and the others may be on in inclined orbits. For semimajor axis ratios consistent with a 4:2:1 mean-motion resonance our analysis yields precise values for HR 8799 bcd's orbital parameters and strictly constrains the planets' eccentricities to be less than 0.18-0.3. However, we find no acceptable orbital solutions with this resonance that place the planets in face-on orbits; HR 8799 d shows the largest deviation from such orbits. Moreover, few orbits make HR 8799 d coplanar with b and c, whereas dynamical stability analyses used to constrain the planets' masses typically assume coplanar and/or fare.on orbits. This paper illustrates the significant science gain enabled with the release of the NIRC2 archive.

Currie, Thayne↗
DOE OSTI2023

Flow Directions and Ages of Subsurface Water in a Salt Marsh System Constrained by Isotope Tracing

Salt marshes are dynamic hydrologic systems where terrestrial groundwater, terrestrial surface water, and seawater mix due to bi-directional flows and pressure gradients. Due to the counteracting terrestrial and marine forcings that control these environments, we do not comprehensively understand water fluxes in these complex coastal systems. To understand the water sources, flow directions, and velocities in salt marsh porewater, we employed a combination of geochemical tracers and analytical models across a hillslope-to-salt marsh continuum in a salt marsh experiencing daily inundation of estuarine surface water (SW) from tides and mixing of fresh seasonal groundwater. We used tritium ( 3 H) as a hydrologic tracer to assess porewater ages and stable water isotope ($\delta$ 2 H and $\delta$ 18 O) analyses to separate isotopically distinct estuarine and terrestrial groundwater across different depths and landscape positions in the study transect. We employed electrical conductivity to constrain the role of source mixing and evapotranspiration in salt marsh hydrology. Salinity and stable isotopes revealed that transpiration, rather than evaporation, increased subsurface water salinity to concentrations above estuarine SW during summer. Elevated salinity at depth indicated that salt marsh subsurface water is recharged during the dry growing season. Seasonal recharge patterns drive long-term deep subsurface water dynamics across the salt marsh, with 3 H ages of 3–7 years, and daily tidal cycles drive short-term shallow porewater dynamics with 3 H ages of 0 ± 3.6 years. Furthermore, our conceptual understanding of the spatiotemporal changes in SW-subsurface water interactions at the terrestrial-marine interface quantifies the hydrological constraints we are missing to improve our understanding of biogeochemical cycles within the salt marsh.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Role of Surface Features on the Initial Dissolution of CH 3 NH 3 PbI 3 Perovskite in Liquid Water: An Ab Initio Molecular Dynamics Study

In this study, the degradation of CH 3 NH 3 PbI 3 (MAPbI 3 ) hybrid organic inorganic perovskite (HOIP) by water has been the major issue hampering its use in commercial perovskites solar cells (PSCs) as MAPbI 3 HOIP have been known to easily degrade in the presence of water molecules. However, even though there have been numerous studies investigating this phenomenon, there is still no consensus on the mechanisms of initial stages of dissolution. Here, we attempt to consolidate differing mechanistic interpretations previously reported in the literature through the use of the first-principles constrained ab-initio molecular dynamics (AIMD) to study the mechanisms, kinetics, and thermodynamics that accompany the degradation of MAPbI 3 HOIP in liquid water. We consider not only the dissolution of the species found on the pristine MAPbI 3 HOIP surfaces, but also the dissolution from defect sites to imitate the successive degradation steps, and propose a sequence of events in the initial phase of MAPbI 3 HOIP dissolution. By comparing the dissolution free energy barrier between surface species of different surficial types, we find that the dominant dissolution mechanisms of surface species varies widely based on the specific surface features adjacent to the dissolving ion. The high sensitivity of dissolution mechanism to surface features has contributed to the many dissolution mechanisms proposed in the literature. In contrast, the dissolution free energy barriers are mainly determined by the dissolving species rather than the type of surfaces, and the type of surfaces the ions are dissolving from are inconsequential toward the dissolution free energy barrier. However, the presence of surface defects such as vacancy sites are found to significantly lower the dissolution free energy barriers. Based on the estimated dissolution free energy barriers from different types of surfaces that we investigated in this study, we proposed that the dissolution of MAPbI 3 HOIP in liquid water originates from surface defect sites that propagate laterally along the surface layer of the MAPbI 3 HOIP crystal.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2007

Influences of Local Sea-Surface Temperatures and Large-scale Dynamics on Monthly Precipitation Inferred from Two 10-year GCM-Simulations

Two parallel sets of 10-year long: January 1, 1982 to December 31, 1991, simulations were made with the finite volume General Circulation Model (fvGCM) in which the model integrations were forced with prescribed sea-surface temperature fields (SSTs) available as two separate SST-datasets. One dataset contained naturally varying monthly SSTs for the chosen period, and the oth& had the 12-monthly mean SSTs for the same period. Plots of evaporation, precipitation, and atmosphere-column moisture convergence, binned by l C SST intervals show that except for the tropics, the precipitation is more strongly constrained by large-scale dynamics as opposed to local SST. Binning data by SST naturally provided an ensemble average of data contributed from disparate locations with same SST; such averages could be expected to mitigate all location related influences. However, the plots revealed: i) evaporation, vertical velocity, and precipitation are very robust and remarkably similar for each of the two simulations and even for the data from 1987-ENSO-year simulation; ii) while the evaporation increased monotonically with SST up to about 27 C, the precipitation did not; iii) precipitation correlated much better with the column vertical velocity as opposed to SST suggesting that the influence of dynamical circulation including non-local SSTs is stronger than local-SSTs. The precipitation fields were doubly binned with respect to SST and boundary-layer mass and/or moisture convergence. The analysis discerned the rate of change of precipitation with local SST as a sum of partial derivative of precipitation with local SST plus partial derivative of precipitation with boundary layer moisture convergence multiplied by the rate of change of boundary-layer moisture convergence with SST (see Eqn. 3 of Section 4.5). This analysis is mathematically rigorous as well as provides a quantitative measure of the influence of local SST on the local precipitation. The results were recast to examine the dependence of local rainfall on local SSTs; it was discernible only in the tropics. Our methodology can be used for computing relationship between any forcing function and its effect(s) on a chosen field.

Sud, Y. C.↗
NASA NTRS2022

Electromagnetic Time and Wavelength Domain Predictions for Accreting Binary Black Holes

Accreting massive binary black holes in gaseous environments are keymulti-messenger sources for the LISA mission. We will present some ofthe latest predictions of the observational signatures of accreting binary black holes across the electromagnetic spectrum from the fieldand provide some of our own. In particular, we will provide a briefsummary of the progress made in our group to understand these systemstheoretically using high-performance simulations of generalrelativistic magnetohydrodynamics (MHD) and dynamic spacetimes.Simulating these systems realistically is challenging asradiation-coupled MHD must be considered over large dynamic ranges inspace and time. Constraining our view to the circumbinary disk regionfirst, we will show how the evolution duration, accretion disk size,mass ratio, and thermodynamics model affect the structure andvariability of the accretion flow. We particularly emphasize how theseparameters influence the over-density feature, or "lump", which orbitsthe binary near the edge of the cavity, since it is responsible formost of the electromagnetic emission's variability---a key signatureof a system being a binary. Extending to smaller length scales, wewill report on simulations following accretion all the down to theevent horizons so that we may begin to investigate how black hole spinaffects mini-disk dynamics, accretion rate, and jet power. We willalso report on our radiative transfer predictions for how black holespin alters the spectral energy distribution of these systems.Wrapping up, we will itemize the remaining key theoreticaltechnologies and directions needed to be accomplished to betterunderstand these sources, and assess strategies for searching for themin the coming decade.

Scott C Noble↗
NASA NTRS2013

Initial Assessment of a Variable-Camber Continuous Trailing-Edge Flap System on a Rigid Wing for Drag Reduction in Subsonic Cruise

In this paper, we describe an initial optimization study of a Variable-Camber Continuous Trailing-Edge Flap (VCCTEF) system. The VCCTEF provides a light-weight control system for aircraft with long flexible wings, providing efficient high-lift capability for takeoff and landing, and greater efficiency with reduced drag at cruising flight by considering the effects of aeroelastic wing deformations in the control law. The VCCTEF system is comprised of a large number of distributed and individually-actuatable control surfaces that are constrained in movement relative to neighboring surfaces, and are non-trivially coupled through structural aeroelastic dynamics. Minimzation of drag results in a constrained, coupled, non-linear optimization over a high-dimension search space. In this paper, we describe the modeling, analysis, and optimization of the VCCTEF system control inputs for minimum drag in cruise. The purpose of this initial study is to quantify the expected benefits of the system concept. The scope of this analysis is limited to consideration of a rigid wing without structural flexibility in a steady-state cruise condition at various fuel weights. For analysis, we developed an optimization engine that couples geometric synthesis with vortex-lattice analysis to automate the optimization procedure. In this paper, we present and describe the VCCTEF system concept, optimization approach and tools, run-time performance, and results of the optimization at 20%, 50%, and 80% fuel load. This initial limited-scope study finds the VCCTEF system can potentially gain nearly 10% reduction in cruise drag, provides greater drag savings at lower operating weight, and efficiency is negatively impacted by the severity of relative constraints between control surfaces.

trailing edge flaps↗
DOE OSTI2026

Dynamic Evolution of Mass and Physical Properties of Atmospheric Organic Aerosol under Solar Irradiance

Organic aerosol (OA) particles constitute a substantial fraction of sub-micron particulate mass in the atmosphere and play a critical role in climate system. OA undergoes dynamic aging processes in the atmosphere, with photolytic aging induced by ultraviolet solar irradiance being an important yet poorly characterized mechanism. Knowledge gaps persist regarding the role of volatility transformations during photolytic aging on the OA mass decay kinetics and the evolution of climate-relevant properties, such as hygroscopicity, hindering the model evaluation of OA spatiotemporal distributions and atmospheric budgets. In this study, we conduct isothermal photolytic aging experiments on both laboratory-generated secondary organic aerosols and ambient-collected particles from urban Atlanta, utilizing a high-sensitivity Quartz Crystal Microbalance. Our results reveal that photolytic aging reduces 40–66% of the low-volatility OA mass with lifetimes ranging from 8 to 200 hours under solar irradiance, and 44–92% of the photolytic mass loss is through slow evaporation of semi- or intermediate-volatile products, kinetically limited by their volatility. We observe up to ±50% changes in OA hygroscopicity with the transformation of fresh OA to photo-recalcitrant low-volatility products, associated with changes in oxygen-to-carbon ratio and molecular weight. A kinetic model incorporating photolytic volatility transformation provides a cohesive explanation for the observed photolysis-induced changes in mass, volatility, and hygroscopicity. Our results can help constrain model representation of the dynamic evolutions of mass and climate-relevant properties during photolytic aging processes of the ambient OA, improving our understanding of OA atmospheric behavior and climate impact.

Bai, Bin↗
NASA NTRS1988

Passive damping for space truss structures

Theoretical and experimental studies of passive damping techniques in truss-type structures are presented, with emphasis on the use of viscoelastic damping in the parallel load path. The constraining member length is shown to be a convenient design variable for enhancing damping performance. Results are presented for integral damping members made of thin-wall aluminum tubes, concentric constraining members, and viscoelastic materials in a six-bay truss structure at low frequency and low dynamic strain conditions. Integral members with graphite/epoxy constraining members exhibited relatively low damping values due to the possible polymer interaction during the cocure stage.

Chen, Gun-Shing↗
NASA NTRS2017

Implicit Formulation of Muscle Dynamics in OpenSim

Astronauts lose bone and muscle mass during spaceflight. Exercise countermeasure is the primary method for counteracting bone and muscle mass loss in space. New spacecraft exercise device concepts are currently being developed for the NASAs new crew exploration vehicle. The NASA Digital Astronaut Project (DAP) uses computational modeling to help determine if the new exercise devices will be effective as countermeasures. The NASA Digital Astronaut Project is developing the ability to utilize predictive simulation to provide insight into the change in kinematics and kinetics with a change in device and gravitational environment (1-g versus 0-g). For example, in space exercise the subject's body weight is applied in addition to the loads prescribed for musculoskeletal maintenance. How and where these loads are applied obviously directly impacts bone and tissue loads. Additionally, due to space vehicle structural requirements, exercise devices are often placed on vibration isolation systems. This changes the apparent impedance or stiffness of the device as seen by the user. Data collection under these conditions is often impractical and limited. Predictive modeling provides a means to have a virtual subject to test hypotheses. Predictive simulation provides a virtual subject for which we are able to perform studies such as sensitivity to device loading and vibration isolation without the need for laboratory kinematic or kinetic test data.Direct Collocation optimization provides an efficient means to perform task based optimization and predictive modeling. It is relatively straight forward to structure a physical exercise task in a Direct Collocation mathematical formulation: perform a motion such that you start at an initial pose, achieve a given amount of deflection i.e a squat, return to the initial pose, and minimize muscle activation cost. Direct Collocation is advantageous in that it does not require numerical integration to evaluate the objective function. Instead, the system dynamics are transformed to discrete time and the optimizer is constrained such that the solution is not considered to be a valid unless the dynamic equations are satisfied at all time points. The simulation and optimization are effectively done simultaneously. Due to the implicit integration, time steps can be more coarse than in a differential equation solver. In a gait scenario this means that that the model constraints and cost function are evaluated at 100 nodes in the gait cycle versus 10,000 integration steps in a variable-step forward dynamic simulation. Furthermore, no time is wasted on accurate simulations of movements that are far from the optimum. Constrained optimization algorithms require a Jacobian matrix that contains the partial derivatives of each of the dynamic constraints with respect to of each of the state and control variables at all time points. This is a large but sparse matrix. An implicit dynamics formulation requires computation of the dynamic residuals f as a function of the states x and their derivatives, and controls u:f(x, dxdt, u) 0If the dynamics of musculoskeletal system are formulated implicitly, the Jacobian elements are often available analytically, eliminating the need for numerical differentiation; this is obviously computationally advantageous. Additionally, implicit formulation of musculoskeletal dynamics do not suffer from singularities from low mass bodies, zero muscle activation, or other stiff system or

physical exercise↗
NASA NTRS2011

Constraining the Absolute Orientation of eta Carinae's Binary Orbit: A 3-D Dynamical Model for the Broad [Fe III] Emission

We present a three-dimensional (3-D) dynamical model for the broad [Fe III] emission observed in Eta Carinae using the Hubble Space Telescope/Space Telescope Imaging Spectrograph (HST/STIS). This model is based on full 3-D Smoothed Particle Hydrodynamics (SPH) simulations of Eta Car's binary colliding winds. Radiative transfer codes are used to generate synthetic spectro-images of [Fe III] emission line structures at various observed orbital phases and STIS slit position angles (PAs). Through a parameter study that varies the orbital inclination i, the PA(theta) that the orbital plane projection of the line-of-sight makes with the apastron side of the semi-major axis, and the PA on the sky of the orbital axis, we are able, for the first time, to tightly constrain the absolute 3-D orientation of the binary orbit. To simultaneously reproduce the blue-shifted emission arcs observed at orbital phase 0.976, STIS slit PA = +38deg, and the temporal variations in emission seen at negative slit PAs, the binary needs to have an i approx. = 130deg to 145deg, Theta approx. = -15deg to +30deg, and an orbital axis projected on the sky at a P A approx. = 302deg to 327deg east of north. This represents a system with an orbital axis that is closely aligned with the inferred polar axis of the Homunculus nebula, in 3-D. The companion star, Eta(sub B), thus orbits clockwise on the sky and is on the observer's side of the system at apastron. This orientation has important implications for theories for the formation of the Homunculus and helps lay the groundwork for orbital modeling to determine the stellar masses.

Madura, T. I.↗
DOE OSTI2021

Measuring H 0 using X-ray and SZ effect observations of dynamically relaxed galaxy clusters

ABSTRACT We use a sample of 14 massive, dynamically relaxed galaxy clusters to constrain the Hubble constant, H0, by combining X-ray and Sunyaev–Zel’dovich (SZ) effect signals measured with Chandra, Planck, and Bolocam. This is the first such analysis to marginalize over an empirical, data-driven prior on the overall accuracy of X-ray temperature measurements, while our restriction to the most relaxed, massive clusters also minimizes astrophysical systematics. For a cosmological-constant model with Ωm = 0.3 and ΩΛ = 0.7, we find $H_0 = 67.3^{+21.3}_{-13.3}\, \mathrm{km}\, \mathrm{s}^{-1}\, \mathrm{Mpc}^{-1}$, limited by the temperature calibration uncertainty (compared to the statistically limited constraint of $H_0 = 72.3^{+7.6}_{-7.6}\, \mathrm{km}\, \mathrm{s}^{-1}\, \mathrm{Mpc}^{-1}$). The intrinsic scatter in the X-ray/SZ pressure ratio is found to be 13 ± 4 per cent (10 ± 3 per cent when two clusters with significant galactic dust emission are removed from the sample), consistent with being primarily due to triaxiality and projection. We discuss the prospects for reducing the dominant systematic limitation to this analysis, with improved X-ray calibration and/or precise measurements of the relativistic SZ effect providing a plausible route to per cent level constraints on H0.

79 ASTRONOMY AND ASTROPHYSICS↗