Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Constrained control”

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 415 records · Page 23

Biophysical principles of choanoflagellate self-organization

Inspired by the patterns of multicellularity in choanoflagellates, the closest living relatives of animals, we quantify the biophysical processes underlying the morphogenesis of rosette colonies in the choanoflagellate Salpingoeca rosetta . We find that rosettes reproducibly transition from an early stage of 2-dimensional (2D) growth to a later stage of 3D growth, despite the underlying variability of the cell lineages. Our perturbative experiments demonstrate the fundamental importance of a basally secreted extracellular matrix (ECM) for rosette morphogenesis and show that the interaction of the ECM with cells in the colony physically constrains the packing of proliferating cells and, thus, controls colony shape. Simulations of a biophysically inspired model that accounts for the size and shape of the individual cells, the fraction of ECM, and its stiffness relative to that of the cells suffices to explain our observations and yields a morphospace consistent with observations across a range of multicellular choanoflagellate colonies. Overall, our biophysical perspective on rosette development complements previous genetic perspectives and, thus, helps illuminate the interplay between cell biology and physics in regulating morphogenesis.

59 BASIC BIOLOGICAL SCIENCES↗

CAFE: A NEW RELATIVISTIC MHD CODE

We introduce CAFE, a new independent code designed to solve the equations of relativistic ideal magnetohydrodynamics (RMHD) in three dimensions. We present the standard tests for an RMHD code and for the relativistic hydrodynamics regime because we have not reported them before. The tests include the one-dimensional Riemann problems related to blast waves, head-on collisions of streams, and states with transverse velocities, with and without magnetic field, which is aligned or transverse, constant or discontinuous across the initial discontinuity. Among the two-dimensional (2D) and 3D tests without magnetic field, we include the 2D Riemann problem, a one-dimensional shock tube along a diagonal, the high-speed Emery wind tunnel, the Kelvin–Helmholtz (KH) instability, a set of jets, and a 3D spherical blast wave, whereas in the presence of a magnetic field we show the magnetic rotor, the cylindrical explosion, a case of Kelvin–Helmholtz instability, and a 3D magnetic field advection loop. The code uses high-resolution shock-capturing methods, and we present the error analysis for a combination that uses the Harten, Lax, van Leer, and Einfeldt (HLLE) flux formula combined with a linear, piecewise parabolic method and fifth-order weighted essentially nonoscillatory reconstructors. We use the flux-constrained transport and the divergence cleaning methods to control the divergence-free magnetic field constraint.

79 ASTRONOMY AND ASTROPHYSICS↗

Benchmark study of the two-dimensional Hubbard model with auxiliary-field quantum Monte Carlo method

Ground-state properties of the Hubbard model on a two-dimensional square lattice are studied by the auxiliary-field quantum Monte Carlo method. Accurate results for energy, double occupancy, effective hopping, magnetization, and momentum distribution are calculated for interaction strengths of U/t from 2 to 8, for a range of densities including half-filling and n = 0.3,0.5,0.6, 0.75, and 0.875. At half-filling, the results are numerically exact. Away from half-filling, the constrained path Monte Carlo method is employed to control the sign problem. Our results are obtained with several advances in the computational algorithm, which are described in this paper in detail. We discuss the advantages of generalized Hartree-Fock trial wave functions and its connection to pairing wave functions, as well as the interplay with different forms of Hubbard-Stratonovich decompositions. We study the use of different twist angle sets when applying the twist averaged boundary conditions. We propose the use of quasirandom sequences, which improves the convergence to the thermodynamic limit over pseudorandom and other sequences. With it and a careful finite size scaling analysis, we are able to obtain accurate values of ground-state properties in the thermodynamic limit. Detailed results for finite-sized systems up to 16 × 16 are also provided for benchmark purposes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Dark Energy Survey Year 3 results: Constraints on extensions to Λ CDM with weak lensing and galaxy clustering

We constrain six possible extensions to the Λ cold dark matter (CDM) model using measurements from the Dark Energy Survey’s first three years of observations, alone and in combination with external cosmological probes. The DES data are the two-point correlation functions of weak gravitational lensing, galaxy clustering, and their cross-correlation. We use simulated data vectors and blind analyses of real data to validate the robustness of our results to astrophysical and modeling systematic errors. In many cases, constraining power is limited by the absence of theoretical predictions beyond the linear regime that are reliable at our required precision. The Λ CDM extensions are dark energy with a time-dependent equation of state, nonzero spatial curvature, additional relativistic degrees of freedom, sterile neutrinos with eV-scale mass, modifications of gravitational physics, and a binned σ 8 ( z ) model which serves as a phenomenological probe of structure growth. For the time-varying dark energy equation of state evaluated at the pivot redshift we find ( w p , w a ) = ( - 0.9 9 - 0.17 + 0.28 , - 0.9 ± 1.2 ) at 68% confidence with z p = 0.24 from the DES measurements alone, and ( w p , w a ) = ( - 1.0 3 - 0.03 + 0.04 , - 0. 4 - 0.3 + 0.4 ) with z p = 0.21 for the combination of all data considered. Curvature constraints of Ω k = 0.0009 ± 0.0017 and effective relativistic species N eff = 3.1 0 - 0.16 + 0.15 are dominated by external data, though adding DES information to external low-redshift probes tightens the Ω k constraints that can be made without cosmic microwave background observables by 20%. For massive sterile neutrinos, DES combined with external data improves the upper bound on the mass m eff by a factor of 3 compared to previous analyses, giving 95% limits of ( Δ N eff , m eff ) ≤ ( 0.28 , 0.20 eV ) when using priors matching a comparable Planck analysis. For modified gravity, we constrain changes to the lensing and Poisson equations controlled by functions Σ ( k , z ) = Σ 0 Ω Λ ( z ) / Ω Λ , 0 and μ ( k , z ) = μ 0 Ω Λ ( z ) / Ω Λ , 0 , respectively, to Σ 0 = 0. 6 - 0.5 + 0.4 from DES alone and ( Σ 0 , μ 0 ) = ( 0.04 ± 0.05 , 0.0 8 - 0.19 + 0.21 ) for the combination of all data, both at 68% confidence. Overall, we find no significant evidence for physics beyond Λ CDM .

79 ASTRONOMY AND ASTROPHYSICS↗

Stochastic Multi-Timescale Power System Operations With Variable Wind Generation

This paper describes a novel set of stochastic unit commitment and economic dispatch models that consider stochastic loads and variable generation at multiple operational timescales. The stochastic model includes four distinct stages: stochastic day-ahead security-constrained unit commitment (SCUC), stochastic real-time SCUC, stochastic real-time security-constrained economic dispatch (SCED), and deterministic automatic generation control (AGC). These sub-models are integrated together such that they are continually updated with decisions passed from one to another. The progressive hedging algorithm (PHA) is applied to solve the stochastic models to maintain the computational tractability of the proposed models. Comparative case studies with deterministic approaches are conducted in low wind and high wind penetration scenarios to highlight the advantages of the proposed methodology, one with perfect forecasts and the other with current state-of-the-art but imperfect deterministic forecasts. The effectiveness of the proposed method is evaluated with sensitivity tests using both economic and reliability metrics to provide a broader view of its impact.

17 WIND ENERGY↗

Search for Majorana neutrinos in same-sign WW scattering events from pp collisions at $\sqrt{s}=13$ TeV

A search for Majorana neutrinos in same-sign WW scattering events is presented. The analysis uses $\sqrt{s}=13$ TeV proton–proton collision data with an integrated luminosity of 140 fb -1 recorded during 2015–2018 by the ATLAS detector at the Large Hadron Collider. The analysis targets final states including exactly two same-sign muons and at least two hadronic jets well separated in rapidity. The modelling of the main backgrounds, from Standard Model same-sign WW scattering and WZ production, is constrained with data in dedicated signal-depleted control regions. The distribution of the transverse momentum of the second-hardest muon is used to search for signals originating from a heavy Majorana neutrino with a mass between 50 GeV and 20 TeV. No significant excess is observed over the background expectation. The results are interpreted in a benchmark scenario of the Phenomenological Type-I Seesaw model. In addition, the sensitivity to the Weinberg operator is investigated. Upper limits at the 95% confidence level are placed on the squared muon-neutrino–heavy-neutrino mass-mixing matrix element |V μN | 2 as a function of the heavy Majorana neutrino’s mass m N , and on the effective μμ Majorana neutrino mass |m μμ |.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Enabling Real-Time Communication in Multi-Agent Systems: A Graph Neural Network Based Approach

Global connectivity enables effective coordination in Multi-Agent Systems (MAS). Solving these connection problems under hardware constraints is an NP-hard non-Euclidean Degree Constrained Minimum Spanning Tree (DCMST) problem. Prior MAS controllers coordinate team movement for task completion and collision avoidance; some considering Line-of-Sight (LOS) maintenance but prioritizing flexibility over guarantees. Evolutionary Algorithms (EA) have been shown to find good solutions for DCMST, but their performance degrades with larger populations required to support a large MAS. We present a method based on edge graph attention networks, trained offline to reduce online computation times. Empirical comparisons with greedy polynomial-time solvers and EA show that our method leverages latent graph information to consistently find constraint-satisfying solutions in less time.

connectivity maintenance↗

On the correlation between the stress exponent for creep determined by nanoindentation and the mechanism of action enabling stress relief in indium

Instrumented indentation performed at room temperature with a Berkovich and 10 μm radius sphere has been used to measure the stress exponent for creep before and after the strain burst observed in well-annealed, high-purity indium. Before the strain burst, the measured values are successfully rationalized using a new model based on stress directed diffusional flow along the interface between the indenter tip and test specimen. After the strain burst, the measured stress exponents are found to be representative of dislocation glide and climb assisted glide. Here these results are compared and contrasted to the previous experimental investigations and modeling efforts of Feng et al., Lucas et al., and Li et al. Collectively, the experimental observations and rationalization presented here provide significant new insight into the mechanisms of action that control the competition for stress relief in small, constrained volumes of crystalline metals subjected to high homologous temperatures.

36 MATERIALS SCIENCE↗

A Regional Phase Amplitude Model of 2-D Attenuation for North America

We analyzed seismic attenuation patterns across the North America using over 70,000,000 Lg wave amplitudes recorded by the various networks at frequencies of 0.1-32 Hz. Our inversion solved for laterally varying attenuation, site terms, moments, and apparent stress following Phillips et al., (2016). The inversion was anchored by independently constrained: corner frequencies (via coda spectral ratios) to control the attenuation-stress tradeoff and moment measurements of teleseismic (GCMT, USGS) and regional (St. Louis University and UC Berkeley) earthquakes to provided absolute scaling. The quality factor (Q) shows clear regional patterns: low values in coastal, volcanic, and tectonically active regions, and high values in stable areas like the Great Plains and major plateaus throughout North America. These 2-D Q models enable improved regional source characterization, magnitude estimation, and yield determination.

58 GEOSCIENCES↗

The Effects of the Uncertainty of Departures on Multi-Center Traffic Management Advisor (TMA) Scheduling

The Multi-center Traffic Management Advisor (McTMA) provides a platform for regional or national traffic flow management, by allowing long-range cooperative time-based metering to constrained resources, such as airports or air traffic control center boundaries. Part of the demand for resources is made up of proposed departures, whose actual departure time is difficult to predict. For this reason, McTMA does not schedule the departures in advance, but rather relies on traffic managers to input their requested departure time. Because this happens only a short while before the aircraft's actual departure, McTMA is unable to accurately predict the amount of delay airborne aircraft will need to take in order to accommodate the departures. The proportion of demand which is made up by such proposed departures increases as the horizon over which metering occurs gets larger. This study provides an initial analysis of the severity of this problem in a 400-500 nautical mile metering horizon and discusses potential solutions to accommodate these departures. The challenge is to smoothly incorporate departures with the airborne stream while not excessively delaying the departures.' In particular, three solutions are reviewed: (1) scheduling the departures at their proposed departure time; (2) not scheduling the departures in advance; and (3) scheduling the departures at some time in the future based on an estimated error in their proposed time. The first solution is to have McTMA to automatically schedule the departures at their proposed departure times. Since the proposed departure times are indicated in their flight times in advance, this method is the simplest, but studies have shown that these proposed times are often incorrect2 The second option is the current practice, which avoids these inaccuracies by only scheduling aircraft when a confirmed prediction of departure time is obtained from the tower of the departure airport. Lastly, McTMA can schedule the departures at a predicted departure time based on statistical data of past departure time performance. It has been found that departures usually have a wheels-up time after their indicated proposed departure time, as shown in Figure 1. Hence, the departures were scheduled at a time in the future based on the mean error in proposed departure times for their airport.

Thipphavong, Jane↗

Absolute Stability Analysis of a Phase Plane Controlled Spacecraft

Many aerospace attitude control systems utilize phase plane control schemes that include nonlinear elements such as dead zone and ideal relay. To evaluate phase plane control robustness, stability margin prediction methods must be developed. Absolute stability is extended to predict stability margins and to define an abort condition. A constrained optimization approach is also used to design flex filters for roll control. The design goal is to optimize vehicle tracking performance while maintaining adequate stability margins. Absolute stability is shown to provide satisfactory stability constraints for the optimization.

Jang, Jiann-Woei↗

Monitoring Ocean CO2 Fluxes from Space: GOSAT and OCO-2

The ocean is a major component of the global carbon cycle, emitting over 330 billion tons of carbon dioxide (CO2) into the atmosphere each year, or about 10 times that emitted fossil fuel combustion and all other human activities [1, 2]. The ocean reabsorbs a comparable amount of CO2 each year, along with ~25% of the CO2 emitted by these human activities. The nature and geographic distribution of the processes controlling these ocean CO2 fluxes are still poorly constrained by observations. A better understanding of these processes is essential to predict how this important CO2 sink may evolve as the climate changes.While in situ measurements of ocean CO2 fluxes can be very precise, the sampling density is far too sparse to quantify ocean CO2 sources and sinks over much of the globe. One way to improve the spatial resolution, coverage, and sampling frequency is to make observations of the column averaged CO2 dry air mole fraction, XCO2, from space [4, 5, 6]. Such measurements could provide global coverage at high resolution (< 100 km) on monthly time scales. High precision (< 1 part per million, ppm) is essential to resolve the small, near-surface CO2 variations associated with ocean fluxes and to better constrain the CO2 transport over the ocean. The Japanese Greenhouse gases Observing Satellite (GOSAT) and the NASA Orbiting Carbon Observatory (OCO) were first two space based sensors designed specifically for this task. GOSAT was successfully launched on January 23, 2009, and has been returning measurements of XCO2 since April 2009. The OCO mission was lost in February 2009, when its launch vehicle malfunctioned and failed to reach orbit. In early 2010, NASA authorized a re-flight of OCO, called OCO-2, which is currently under development.

GOSAT↗

Evaluation of Intercomparisons of Four Different Types of Model Simulating TWP-ICE

Four model intercomparisons were run and evaluated using the TWP-ICE field campaign, each involving different types of atmospheric model. Here we highlight what can be learnt from having single-column model (SCM), cloud-resolving model (CRM), global atmosphere model (GAM) and limited-area model (LAM) intercomparisons all based around the same field campaign. We also make recommendations for anyone planning further large multi-model intercomparisons to ensure they are of maximum value to the model development community. CRMs tended to match observations better than other model types, although there were exceptions such as outgoing long-wave radiation. All SCMs grew large temperature and moisture biases and performed worse than other model types for many diagnostics. The GAMs produced a delayed and significantly reduced peak in domain-average rain rate when compared to the observations. While it was shown that this was in part due to the analysis used to drive these models, the LAMs were also driven by this analysis and did not have the problem to the same extent. Based on differences between the models with parametrized convection (SCMs and GAMs) and those without (CRMs and LAMs), we speculate that that having explicit convection helps to constrain liquid water whereas the ice contents are controlled more by the representation of the microphysics.

convection↗

Revision 1: Characterizing Basalt-Atmosphere Interactions on Venus: A Review of Thermodynamic and Experimental Results

The surface of Venus is in contact with a hot (~470° C), high pressure (92 bars), and caustic (CO 2 with S, but little H 2 O) atmosphere, which should cause progressive alteration of the crust in the form of sulfate and iron-oxide coatings; however, the exact rate of alteration and mineral species are not well constrained. Different experimental approaches, each with its own limitations, are currently being used to constrain mineralogy and alteration rate. One note is that no experimental approach has been able to fully replicate the necessary conditions and sustain them for a significant length of time. Further, geochemical modeling studies can also constrain surface alteration mineralogy, again with different assumptions and limitations. Here we review recent geochemical modeling and experimental studies to constrain the state of the art for alteration mineralogy, rate of alteration, open questions about the surface mineralogy of Venus, and what can be constrained before the fleet of missions arrive later this decade. Combining the new results confirm that basalt on the surface of Venus should react quickly and form coatings of sulfates and iron-oxides; however, the mineralogy and rate of alteration is dependent on physical properties of the protolith (including bulk composition, mineralogy, and crystallinity), as well as atmospheric composition, and surface temperature. Importantly, the geochemically modelling results show that the mineralogy is largely controlled by atmospheric oxygen fugacity, which is not well constrained for the near surface environment on Venus. Therefore, alteration experiments run over a range of oxygen and sulfur fugacities are needed across a wide range of Venus analog materials with varying mineralogy and crystallinity.

Justin Filiberto↗

Constrained optimum trajectories with specified range

The characteristics of optimum fixed-range trajectories whose structure is constrained to climb, steady cruise, and descent segments are derived by application of optimal control theory. The performance function consists of the sum of fuel and time costs, referred to as direct operating costs (DOC). The state variable is range-to-go and the independent variable is energy. In this formulation a cruise segment always occurs at the optimum cruise energy for sufficiently large range. At short ranges (500 n. mi. and less) a cruise segment may also occur below the optimum cruise energy. The existence of such a cruise segment depends primarily on the fuel flow vs thrust characteristics and on thrust constraints. If thrust is a free control variable along with airspeed, it is shown that such cruise segments will not generally occur. If thrust is constrained to some maximum value in climb and to some minimum in descent, such cruise segments generally will occur. The performance difference between free thrust and constrained thrust trajectories has been determined in computer calculations for an example transport aircraft.

Erzberger, H.↗

Optimum climb and descent trajectories for airline missions

The characteristics of optimum fixed-range trajectories whose structure is constrained to climb, steady cruise, and descent segments are derived by application of optimal control theory. The performance function consists of the sum of fuel and time costs, referred to as direct operating cost (DOC). The state variable is range to go and the independent variable is energy. In this formulation a cruise segment always occurs at the optimum cruise energy for sufficiently large range. At short ranges (400 n. mi. and less), a cruise segment may also occur below the optimum cruise energy. The existence of such a cruise segment depends primarily on the fuel flow vs thrust characteristics and on thrust constraints. If thrust is a free control variable along with airspeed, it is shown that such cruise segments will not generally occur. If thrust is constrained to some maximum value in climb and to some minimum in descent, such cruise segments generally will occur.

Erzberger, H.↗

Foundations of Molecular 'Isotomics'

The naturally occurring rare isotopes are versions of common elements, such as hydrogen, carbon and oxygen, that contain a larger than usual number of neutrons in their atomic nuclei and therefore are higher in mass than the common atoms of that element. Isotopes exist for most elements and are found in most natural and synthetic materials, but are uneven in their distribution because chemical and physical processes are isotope-selective (e.g., a chemical reaction may proceed more rapidly for one isotope than for another). For this reason, abundances of isotopes in a material of interest can provide a record, or ‘signature’ of various features of that material’s origin and history. These signatures have been used in the geo, life, chemical and physical sciences in a wide variety of ways over close to 8 decades. However, many such applications struggle to reach unique interpretations of isotopic data because multiple factors combine to control a given sample’s overall isotopic content. That is, the factors controlling isotopic content are too numerous and complex to fully constrain from a simple measurement of a material’s isotope abundances. However, the distribution of isotopes within materials, at molecular scales potentially provides a vastly larger number and diversity of constraints on the chemical and physical processes that comprise a material’s history. The rare isotopes may be concentrated into one atomic position in a molecule relative to another, some proportion of molecules in a sample may contain two or more rare isotopes, and those multiply-isotope-substituted forms of molecules may also have uneven distributions of those isotopes across individual atomic sites. For these reasons, even small, seemingly simple molecules, such as sugars, amino acids or drug compounds, actually exist in a vast number of isotopically unique forms (often millions or more), and each one of those forms is in some sense an independent ‘vote’ on that sample’s history. This project has focused on opening this rich archive of information by enabling the creation of routinely and widely applicable ways of measuring and interpreting isotopic structures of molecules. This work has included the development of core technologies and analytical methods, advancing fundamental understanding of the physical and chemical properties of isotopic versions of molecules, and conducting proof of concept studies of illustrative geochemical, cosmochemical and forensic problems in order to show how these technologies, methods and principles come together to solve problems in new ways. A key to the success of this project was the adaptation of ‘Fourier transform mass spectrometry’ (FTMS) to the task of precisely measuring proportions of the rare, naturally occurring isotopic forms of molecules. FTMS is a highly specialized form of mass spectrometry that traps ions within magnetic or electrostatic cavities and, effectively, ‘listens’ (through registering of subtle electrical signals) to the harmonic signals they make while rapidly orbiting within those cavities. These signals have periods that are a function of their mass and strength (or ‘loudness’) that is proportional to their abundances. Thus, these signals constrain relative amounts of molecules that differ in their mass due to various isotopic substitutions. This technology has been essential to the identification of organic molecules in the life, chemical and environmental sciences for over 4 decades, but generally has lacked the control, stability and precision to meaningfully measure rare isotope forms of molecules. This project’s most fundamental contribution has been to modify FTMS, both in terms of hardware and methods, to enable such measurements. The raw data of molecular isotopic structure is tremendously voluminous and complex, so another important activity of this project has been developing the theoretical and data-science tools needed to interpret the data generated by this new form of isotopic measurement. A particularly challenging part of this task has been predicting molecular isotopic structure, as only through the comparison of measurements with predictions can we make progress on hypothesis driven research questions. We have attacked this this prediction task through a combination of first-principles chemical-physics models of the effects of isotope substitution on molecule properties and data-science models that permit us to generalize that chemical physics to cases that have not yet been studied by detailed chemical physics theory. The proof of concept applications we have pursued over the course of this study include biological reactions of amino acids and other biomolecules, non-biological synthesis of organic molecules in extra-terrestrial settings such as meteorites, petroleum geoscience questions concerning the origin and evolution of natural gas, oil and kerogen compounds, and forensic questions such as the sourcing of chemical weapons. The successes of these applications have laid the groundwork for the next phase of this field’s development, which will include larger scale and more ambitious studies of molecular isotopic structure as a means of diagnosing human diseases, such as cancer, and reconstructing detailed interpretations of the origin and evolution of organic molecules in modern and geological environments.

Cesar, Jaime↗

Nonlinear feedback control of multiple robot arms

Multiple coordinated robot arms are modeled by considering the arms: (1) as closed kinematic chains, and (2) as a force constrained mechanical system working on the same object simultaneously. In both formulations a new dynamic control method is discussed. It is based on a feedback linearization and simultaneous output decoupling technique. Applying a nonlinear feedback and a nonlinear coordinate transformation, the complicated model of the multiple robot arms in either formulation is converted into a linear and output decoupled system. The linear system control theory and optimal control theory are used to design robust controllers in the task space. The first formulation has the advantage of automatically handling the coordination and load distribution among the robot arms. In the second formulation, by choosing a general output equation, researchers can superimpose the position and velocity error feedback with the force-torque error feedback in the task space simultaneously.

Tarn, T. J.↗