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 1,027 records · Page 57

What Dominates the Coronal Emission Spectrum During the Cycle of Impulsive Heating and Cooling?

The smoking gun of small-scale, impulsive events heating the solar corona is expected to be the presence of a hot ( > 5 MK) plasma component. Evidence for this has been scarce, but has gradually begun to accumulate due to recent studies designed to constrain the high temperature part of the emission measure distribution. However, the detected hot component is often weaker than models predict and this is due in part to the common modeling assumption that the ionization balance remains in equilibrium. The launch of the latest generation of space-based observing instrumentation aboard Hinode and the Solar Dynamics Observatory (SDO) has brought the matter of the ionization state of the plasma firmly to the forefront. It is timely to consider exactly what emission current instruments would detect when observing a corona heated impulsively on small-scales by nanoflares. Only after we understand the full effects of nonequilibrium ionization can we draw meaningful conclusions about the plasma that is (or is not) present. We have therefore performed a series of hydrodynamic simulations for a variety of different nanoflare properties and initial conditions. Our study has led to several key conclusions. 1. Deviations from equilibrium are greatest for short-duration nanoflares at low initial coronal densities. 2. Hot emission lines are the most affected and are suppressed sometimes to the point of being invisible. 3. The emission detected in all of the SDO-AIA channels is generally dominated by warm, over-dense, cooling plasma. 4. It is difficult not to create coronal loops that emit strongly at 1.5 MK and in the range 2 to 5 MK, which are the most commonly observed kind, for a broad range of nanoflare scenarios. 5. The Fe XV (284.16 ) emission in most of our models is about 10 times brighter than the Ca XVII (192.82 ) emission, consistent with observations. Our overarching conclusion is that small-scale, impulsive heating inducing a nonequilibrium ionization state leads to predictions for observable quantities that are entirely consistent with what is actually observed.

Bradshaw, Stephen J.↗

Identifying and Characterizing Impact Melt Outcrops in the Nectaris Basin

The Nectaris Basin is an 820-km diameter, multi-ring impact basin located on the near side of the Moon. Nectaris is a defining stratigraphic horizon based on relationships between ejecta units, giving its name to the Nectarian epoch of lunar history. Lunar basin chronology based on higher resolution LRO imagery and topography, while assigning some important basins like Serenitatis to pre-Nectarian time, were generally consistent with those previously derived. Based on this stratigraphy, at least 11 large basins formed in the time between Nectaris and Imbrium. The absolute age of Nectaris, therefore, is a crucial marker in the lunar time-stratigraphic sequence for understanding the impact flux on the Moon, and by extension, the entire inner solar system. For several decades, workers have attempted to constrain the age of the Nectaris basin through radiometric dating of lunar samples. However, there is little agreement on which samples in our collection represent Nectaris, if any, and what the correct radiometric age of such samples is. The importance of the age of Nectaris goes far beyond assigning a stratigraphic marker to lunar chronology. Several dynamical models use Nectaris as their pin date, so that this date becomes crucial in understanding the time-correlated effects in the rest of the solar system. The importance of the Nectaris basin age, coupled with its nearside, mid-latitude location, make remnants of the impact-melt sheet an attractive target for a future mission, either for in-situ dating or for sample return. We have started exploring this possibility. We have begun a consortium data-analysis effort bringing multiple datasets and analysis methods to bear on these putative impact-melt deposits to characterize their extent, elemental composition and mineralogy, maturity and geologic setting, and to identify potential landing sites that meet both operational safety and science requirements.

Nectaris Basin↗

Consequences and benefits of utilizing continuous vibro-impact representations in constrained pipeline conveying fluid systems

Here, the effectiveness of continuous vibro-impact forcing representations for the cantilevered pipe that conveys fluid is explored and analyzed. The previously accepted forcing model utilizing a smoothened trilinear spring is estimated using three continuous forcing representations, namely, polynomial, rational polynomial, and hyperbolic tangent. The accuracy of the estimated forcing functions is investigated and analyzed by calculating the root mean square error, and bifurcation diagrams are generated and compared to the nominal system. Additionally, the dynamic response of the system is further characterized using Poincare maps, power spectra, and basins of attraction. Once all continuous forcing representations are analyzed and compared to the nominal system, the computational cost of each method is examined, and further limitations of the hyperbolic tangent method are discovered. It is proved that the hyperbolic tangent forcing representation most accurately captures the dynamic response of the pipeline, and the least accurate representation is the rational polynomial representation. Additionally, considerable computational cost is saved when employing the hyperbolic tangent representation compared to the discontinuous representation.

42 ENGINEERING↗

An algorithm for constrained maximization of the trace of a matrix

The 'trace' of a rectangular matrix is defined as the trace of a square matrix obtained by appending null rows (or columns) at the bottom (or right) end. The problem of maximizing the trace of a matrix, by permutations and mergers of rows and columns with constraints on the resulting size of the matrix, is of interest in comparison of maps and image-change detection. This correspondence presents an algorithm based on dynamic programming for efficient maximization of trace.

Ramapriyan, H. K.↗

Modal identities for elastic bodies, with application to vehicle dynamics and control

It is a standard procedure to analyze a flexible vehicle in terms of its vibration frequencies and mode shapes. However, the entire mode shape is not needed per se, but two integrals of the mode shape, pi and hi, which correspond to the momentum and angular momentum in Mode i. Together with the natural frequencies omega-i, these modal parameters satisfy several important identities, 25 of which are derived in this paper. Expansions in terms of both constrained and unconstrained modes are considered. A simple illustrative example is included. The paper concludes with some remarks on the theoretical and practical utility of these results, and several potential extensions to the theory are suggested.

Hughes, P. C.↗

Response of the Current Climate to Land‐Ocean Contrasts in Parameterized Cumulus Entrainment

Abstract Cumulus entrainment substantially regulates the earth's climate but remains poorly constrained in global climate models. Recent studies have shown that cumulus bulk entrainment (or dilution) is particularly sensitive to continentality, with the entrainment rate in simulated maritime cumuli nearly double that of continental cumuli. This study examines the impacts of such land–ocean entrainment contrasts on the current climate using 21‐year simulations with the Geophysical Fluid Dynamics Laboratory's High‐Resolution Atmospheric Model (HIRAM). In response to a 25% reduction in the HIRAM entrainment parameter c 0 over land, precipitation over tropical land regions increases by up to 40%. Along with directly facilitating enhanced convective precipitation, this c 0 reduction induces an increase in soil moisture, which may contribute to a further enhancement of convective precipitation over land. A 25% c 0 reduction over the oceans leads to more widespread modifications of convection patterns, with the strongest signal in the tropical Pacific. Deep convection shifts upstream (eastward) there, inducing enhanced large‐scale ascent over the central Pacific with compensating subsidence and reduced humidity and precipitation over the western Pacific (WP). Land–ocean variations in c 0 project onto the Pacific Walker circulation, with the 25% land reduction strengthening it by 4% and the 25% ocean reduction weakening it by 14%. These changes are driven by variations in convective and large‐scale stratiform heating over the Pacific. While reduced c 0 over land enhances diabatic heating in the Maritime Continent to strengthen the Walker circulation, reduced c 0 over the oceans decreases diabatic heating in the WP to weaken the Walker circulation.

54 ENVIRONMENTAL SCIENCES↗

Viterbi decoding of CRES signals in Project 8

Abstract Cyclotron radiation emission spectroscopy (CRES) is a modern approach for determining charged particle energies via high-precision frequency measurements of the emitted cyclotron radiation. For CRES experiments with gas within the fiducial volume, signal and noise dynamics can be modelled by a hidden Markov model. We introduce a novel application of the Viterbi algorithm in order to derive informational limits on the optimal detection of cyclotron radiation signals in this class of gas-filled CRES experiments, thereby providing concrete limits from which future reconstruction algorithms, as well as detector designs, can be constrained. The validity of the resultant decision rules is confirmed using both Monte Carlo and Project 8 data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The Role of Seasonality and the ENSO Mode in Central and East Pacific ENSO Growth and Evolution

Abstract A cyclostationary linear inverse model (CSLIM) is used to investigate the seasonal growth of tropical Pacific Ocean El Niño–Southern Oscillation (ENSO) events with canonical, central Pacific (CP), or eastern Pacific (EP) sea surface temperature (SST) characteristics. Analysis shows that all types of ENSO events experience maximum growth toward final states occurring in November and December. ENSO events with EP characteristics also experience growth into May and June, but CP events do not. A single dominant “ENSO mode,” growing from an equatorial heat content anomaly into a characteristic ENSO-type SST pattern in about 9 months (consistent with the delayed/recharge oscillator model of ENSO), is essential for the predictable development of all ENSO events. Notably, its seasonality is responsible for the late-calendar-year maximum in ENSO amplification. However, this ENSO mode alone does not capture the observed growth and evolution of diverse ENSO events, which additionally involve the seasonal evolution of other nonorthogonal Floquet modes. EP event growth occurs when the ENSO mode is initially “covered up” in combination with other Floquet modes. The ENSO mode’s slow seasonal evolution allows it to emerge while the other modes rapidly evolve and/or decay, leading to strongly amplifying and more predictable EP events. CP events develop when the initial state has a substantial contribution from Floquet modes with meridional mode–like SST structures. Thus, while nearly all ENSO events involve the seasonally varying ENSO-mode dynamics, the diversity and predictability of ENSO events cannot be understood without identifying contributions from the remaining Floquet modes. Significance Statement The purpose of this study is to identify structures that lead to seasonal growth of diverse types of El Niño–Southern Oscillation (ENSO) events. An important contribution from this study is that it uses an observationally constrained, empirically derived seasonal model. We find that processes affecting the evolution of diverse ENSO events are strongly seasonally dependent. ENSO events with eastern equatorial Pacific sea surface temperature (SST) characteristics are closely related to a single “ENSO mode” that resembles theoretical models of ENSO variability. ENSO events that have central equatorial Pacific SST characteristics include contributions from additional “meridional mode” structures that evolve via different physical processes. These findings are an important step in evaluating the seasonal predictability of ENSO diversity.

Vimont, Daniel J.↗

Photometry of Outer Solar System Objects from the Dark Energy Survey. II. A Joint Analysis of Trans-Neptunian Absolute Magnitudes, Colors, Light Curves and Dynamics

For the 696 trans-Neptunian objects (TNOs) with absolute magnitudes 5.5 < H r < 8.2 detected in the Dark Energy Survey, we characterize the relationships between their dynamical state and physical properties—namely H r , indicating size; colors, indicating surface composition; and flux variation semiamplitude A, indicating asphericity and surface inhomogeneity. We seek “birth” physical distributions that can recreate these parameters in every dynamical class. We show that the observed colors of these TNOs are consistent with two Gaussian distributions in griz space, “near-infrared bright” (NIRB) and “near-infrared faint” (NIRF), presumably an inner and outer birth population, respectively. We find a model in which both the NIRB and NIRF H r and A distributions are independent of current dynamical states, supporting their assignment as birth populations. All objects are consistent with a common rolling p(H r ), but NIRF objects are significantly more variable. Cold classicals (CCs) are purely NIRF, while hot classical (HC), scattered, and detached TNOs are consistent with ≈ 70% NIRB and the resonance NIRB fractions show significant variation. The NIRB components of the HCs and of some resonances have broader inclination distributions than the NIRFs, i.e. their current dynamics retains information about birth location. We find evidence for radial stratification within the birth NIRB population, in that HC NIRBs are on average redder than detached or scattered NIRBs; a similar effect distinguishes CCs from other NIRFs. We estimate total object counts and masses of each class within our H r range. These results will strongly constrain models of the outer solar system.

79 ASTRONOMY AND ASTROPHYSICS↗

Assimilation of multiple datasets results in large differences in regional- to global-scale NEE and GPP budgets simulated by a terrestrial biosphere model

In spite of the importance of land ecosystems in offsetting carbon dioxide emissions released by anthropogenic activities into the atmosphere, the spatiotemporal dynamics of terrestrial carbon fluxes remain largely uncertain at regional to global scales. Over the past decade, data assimilation (DA) techniques have grown in importance for improving these fluxes simulated by terrestrial biosphere models (TBMs), by optimizing model parameter values while also pinpointing possible parameterization deficiencies. Although the joint assimilation of multiple data streams is expected to constrain a wider range of model processes, their actual benefits in terms of reduction in model uncertainty are still under-researched, also given the technical challenges. In this study, we investigated with a consistent DA framework and the ORCHIDEE-LMDz TBM–atmosphere model how the assimilation of different combinations of data streams may result in different regional to global carbon budgets. To do so, we performed comprehensive DA experiments where three datasets (in situ measurements of net carbon exchange and latent heat fluxes, spaceborne estimates of the normalized difference vegetation index, and atmospheric CO 2 concentration data measured at stations) were assimilated alone or simultaneously. We thus evaluated their complementarity and usefulness to constrain net and gross C land fluxes. We found that a major challenge in improving the spatial distribution of the land C sinks and sources with atmospheric CO 2 data relates to the correction of the soil carbon imbalance.

54 ENVIRONMENTAL SCIENCES↗

Edge-Flames in Von Karman Swirling Flows

Classical understanding of diffusion flames dictates that they, unlike the premixed flames, do not possess a characteristic propagation velocity and are constrained by stoichiometric requirements at the flame surface. However, it has been commonly observed that when local extinction occurs within a diffusion flame sheet, the edges that are formed propagate with distinct speeds. In general, the propagation speed of these edges depend on their geometrical shape (concave, convex, or straight) among other factors. Recently, Buckmaster investigated the dynamics of straight diffusion flame edges separating burning and quenched regions using simplified one-dimensional models. He showed that these flame edges can have positive, negative, or zero velocity depending on the Damkoehler number of the equilibrium diffusion flame that support them. It was also shown that this unsteady flame-edge behavior is intrinsically linked to S-curve behavior of the diffusion flame with varying Damkoehler number. When the system Damkoehler number lies between the extinction and ignition limits, flame edges can propagate as an "ignition wave" or as a "failure wave," and for a critical Damkoehler number remain as a stationary flame-edge. We have extend Buckmaster's 1-d model to more general edge-flame configurations where the edges appear as "flame holes" or as "flame disks". These two configurations along with the straight-edge case cover the entire range of possible edge-flame geometry observable in planar diffusion-flame sheets. A generalized map of edge-flame propagation velocities as a function of the system Damkoehler number and the edge-flame radius is presented. Experimentally we show that edge flames can be created using diffusion flames embedded in von Karman boundary layers. In a von Karman boundary layer, the flow is generated by spinning a solid (fuel) disk in a quiescent ambient gas. Under normal gravity we were able to produce "flame disks" over a range of fuel-disk rotational velocities varying from 0 to 20 revolutions per second, by orienting the burning surface of the fuel disk facing downward.

Nayagam, Vedha↗

Burst Oscillation Probes of Neutron Stars and Nuclear Burning with LOFT

X-ray brightness oscillations during thermonuclear X-ray bursts--burst oscillations--have provided a new probe of neutron star spins as well as of the dependent nuclear burning processes. The frequency drift and amplitude evolution of the oscillations observed during bursts can in principle place constraints on the physics of thermonuclear flame spreading and the dynamics of the burning atmosphere. I use simulations appropriate to LOFT to explore the precision with which the time dependence of the oscillation frequency can be inferred. This can test, for example, different models for the frequency drift, such as up-lift versus geostrophic drift. I also explore the precision with which asymptotic frequencies can be constrained in order to estimate the capability for LOFT to detect the Doppler shifts induced by orbital motion of the neutron star from a sample of bursts at different orbital phases.

Strohmayer, Tod↗

Mesosphere and Lower Thermosphere Temperatures Simulated by WACCM‐X With NAVGEM‐HA Meteorological Analyses and Compared to SABER Observations

Realistic modeling of the dynamics and variability in the upper mesosphere and lower thermosphere (UMLT) is critical to understand the coupling between different layers of the whole atmosphere system. Here we present simulations of the UMLT temperatures at ∼100 km altitude for one year during 2014 by the Whole Atmosphere Community Climate Model with thermosphere-ionosphere extension (WACCM-X) constrained below ∼90 km using meteorological analysis products of the high-altitude version of Navy Global Environmental Model (NAVGEM-HA). The model results are sampled at the same times and latitudes and longitudes as the satellite observations from Thermosphere Ionosphere and Mesosphere Electric Dynamics/Sounding of the Atmosphere using Broadband Emission Radiometry (TIMED/SABER). Comparisons show that the observed and modeled daily zonal mean temperatures are correlated (r ∼0.5–0.7) at most latitudes between ±50°. Both the observations and simulations show an annual variation at mid-latitudes in two hemispheres with the temperature maximum in summer and the minimum in winter, and at lower latitudes the semiannual variation becomes stronger having the temperature maximums at equinoxes and the minimums during solstices. However, the temperatures observed are on average ∼5–10 K (3%–5%) smaller than the model and the observations show a larger variability. Moreover, migrating tidal amplitudes are mostly overestimated by the model. Though differences are noticed, the WACCM-X simulations with NAVGEM-HA meteorological analyses are overall consistent with the SABER observations. These results support that whole atmosphere models informed by high altitude observations would help to simulate the UMLT variability and the atmosphere and ionosphere coupling.

Guiping Liu↗

Data-driven linear time advance operators for the acceleration of plasma physics simulation

In this study, we demonstrate the application of data-driven linear operator construction for time advance with a goal of accelerating plasma physics simulation. We apply dynamic mode decomposition (DMD) to data produced by the nonlinear SOLPS-ITER (Scrape-off Layer Plasma Simulator - International Thermonuclear Experimental Reactor) plasma boundary code suite in order to estimate a series of linear operators and monitor their predictive accuracy via online error analysis. We find that this approach defines when these dynamics can be represented by a sequence of approximate linear operators and is essential for providing consistent projections when compared to an unconstrained application. For linear diffusion and advection–diffusion fluid test problems, we construct and apply operators within explicit and implicit time advance schemes, demonstrating that stability can be robustly guaranteed in each case. We further investigate the use of the linear time advance operators within several integration methods including forward Euler, backward Euler, and the matrix exponential. The application of this method to simulation data from SOLPS-ITER, with varying levels of Markov chain Monte Carlo numerical noise, shows that constrained DMD operators yield a capability to identify, extract, and integrate a (slow) subset of the present timescales. Example applications show that for projected speedup factors of [Formula: see text], and [Formula: see text], a mean relative error of 3%, 5%, and 8% and maximum relative error less than 20% are achievable, which appears acceptable for typical SOLPS-ITER steady-state simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Crystalline Peptoid Nanofibers with a Single-Unit Cell Cross Section

Ultranarrow crystalline one-dimensional nanostructures formed from soft materials facilitate precise structural control in nanomaterial design, which is essential for biomedicine and nanotechnology applications. Systematic control of their hierarchical structure is challenging due to the complexities of simultaneously manipulating multiple noncovalent interactions at such small scales. We employed a polypeptoid crystal motif as a supramolecular synthon to engineer ultranarrow crystalline nanofibers constrained to a single lattice axis by incorporating a single ionizable side chain into the hydrophobic core of a nanosheet-forming peptoid. Cryogenic transmission electron microscopy of the nanofibers revealed detailed molecular arrangements of a unit-cell cross-section and the presence of distinct pH-dependent lattice isoforms that resulted in morphological transformations. Molecular dynamics simulations demonstrated that the ionizable side chain plays a critical role in changing the local conformation of the unit cell, which further impacts the dimensionality of hierarchical structures. Moreover, these fibers were readily functionalized with biological ligands to afford one-dimensional (1D) protein arrays. This approach for the high-precision bottom-up assembly of ultranarrow 1D nanostructures offers significant potential for developing novel biomimetic nanostructures.

Lee, Yen Jea↗

Drought impacts on microbial trait distribution and feedback to soil carbon cycling

Abstract Quantifying the impact of drought on microbial processes and its consequences for soil carbon cycling is hindered by the lack of underlying mechanistic understanding. Therefore, there is a need to scale up the physiological response to changing water status from individual soil microbes to collective communities across different ecosystems. Here we propose the use of a framework that incorporates trait‐based ecology to link drought‐impacted microbial processes to rates of soil carbon decomposition and stabilisation. We briefly synthesise existing knowledge on the effects of drought on microbial physiology at the individual to community scale, before integrating this understanding within a framework incorporating life‐history strategy, ecological strategy and biochemistry. This framework highlights a dynamic allocation to high yield (Y), resource acquisition (A) and stress tolerance (S) pathways as environmental conditions change. Y‐A‐S strategies represent sets of traits that tend to correlate due to physiological or evolutionary trade‐offs. This framework enables assessment of microbial processes along two key environmental gradients of water and resource availability, both of which are constrained by drought. The variable chemistry of biomass and necromass produced under different physiological strategies in response to drying–rewetting impacts organic matter decomposition and stabilisation in soils, and should also be considered when quantifying soil carbon balance. We highlight that diversion of resources away from microbial growth can alter soil organic matter chemistry and its persistence depending on the kind of microbial compounds produced. To advance such a framework, we highlight avenues of research that would enable the further identification and quantification of traits linked to Y‐A‐S strategies and the physiological outcomes at the community level under drought and rewetting, and conclude by hypothesising how ecosystem‐level changes might feedback on to the soil carbon cycle. A scalable understanding of microbial drought‐response mechanisms affecting soil carbon cycling will transform the way microbial physiology is represented in ecosystem studies. Read the free Plain Language Summary for this article on the Journal blog.

59 BASIC BIOLOGICAL SCIENCES↗

A Framework for a Supervisory Expert System for Robotic Manipulators with Joint-Position Limits and Joint-Rate Limits

This report addresses the problem of path planning and control of robotic manipulators which have joint-position limits and joint-rate limits. The manipulators move autonomously and carry out variable tasks in a dynamic, unstructured and cluttered environment. The issue considered is whether the robotic manipulator can achieve all its tasks, and if it cannot, the objective is to identify the closest achievable goal. This problem is formalized and systematically solved for generic manipulators by using inverse kinematics and forward kinematics. Inverse kinematics are employed to define the subspace, workspace and constrained workspace, which are then used to identify when a task is not achievable. The closest achievable goal is obtained by determining weights for an optimal control redistribution scheme. These weights are quantified by using forward kinematics. Conditions leading to joint rate limits are identified, in particular it is established that all generic manipulators have singularities at the boundary of their workspace, while some have loci of singularities inside their workspace. Once the manipulator singularity is identified the command redistribution scheme is used to compute the closest achievable Cartesian velocities. Two examples are used to illustrate the use of the algorithm: A three link planar manipulator and the Unimation Puma 560. Implementation of the derived algorithm is effected by using a supervisory expert system to check whether the desired goal lies in the constrained workspace and if not, to evoke the redistribution scheme which determines the constraint relaxation between end effector position and orientation, and then computes optimal gains.

Mutambara, Arthur G. O.↗

Multisystem feedbacks from a changing climate: Do altered hydrological dynamics control vadose zone carbon nutrient cycling and storage in shallow aquifer systems?

The vadose zone, which extends from upper soils to the subsurface water table, consists of many distinct habitats (including the critical zone), each with its own physical characteristics. Upper soils are typically richer in organic carbon chemical diversity and concentration, while deeper portions near the water table have less labile carbon and a greater percentage of humic acids and other long-lived organics. The availability of carbon and oxygen constrain the habitability of these zones. Typically, microorganisms (bacteria, archaea, and fungi) extend throughout the vadose zone and potentially deeper into the bedrock, while higher eukaryotes (i.e., arthropods and plants) are limited to the surficial soils. An exception to this is deep taproots of some tree species that can extend tens of meters into the subsurface. In subsurface systems, microbial metabolisms are constrained by the availability of carbon (organic and inorganic) and electron acceptors.

54 ENVIRONMENTAL SCIENCES↗