Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constrained regions”

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 379 records · Page 21

Climate Change Adaptation Science Activities at NASA Johnson Space Center

The Johnson Space Center (JSC), located in the southeast metropolitan region of Houston, TX is the prime NASA center for human spaceflight operations and astronaut training, but it also houses the unique collection of returned extraterrestrial samples, including lunar samples from the Apollo missions. The Center's location adjacent to Clear Lake and the Clear Creek watershed, an estuary of Galveston Bay, puts it at direct annual risk from hurricanes, but also from a number of other climate-related hazards including drought, floods, sea level rise, heat waves, and high wind events all assigned Threat Levels of 2 or 3 in the most recent NASA Center Disaster/Risk Matrix produced by the Climate Adaptation Science Investigator Working Group. Based on prior CASI workshops at other NASA centers, it is recognized that JSC is highly vulnerable to climate-change related hazards and has a need for adaptation strategies. We will present an overview of prior CASI-related work at JSC, including publication of a climate change and adaptation informational data brochure, and a Resilience and Adaptation to Climate Risks Workshop that was held at JSC in early March 2012. Major outcomes of that workshop that form a basis for work going forward are 1) a realization that JSC is embedded in a regional environmental and social context, and that potential climate change effects and adaptation strategies will not, and should not, be constrained by the Center fence line; 2) a desire to coordinate data collection and adaptation planning activities with interested stakeholders to form a regional climate change adaptation center that could facilitate interaction with CASI; 3) recognition that there is a wide array of basic data (remotely sensed, in situ, GIS/mapping, and historical) available through JSC and other stakeholders, but this data is not yet centrally accessible for planning purposes.

Stefanov, William L.↗

A Regional CO2 Observing System Simulation Experiment for the ASCENDS Satellite Mission

Top-down estimates of the spatiotemporal variations in emissions and uptake of CO2 will benefit from the increasing measurement density brought by recent and future additions to the suite of in situ and remote CO2 measurement platforms. In particular, the planned NASA Active Sensing of CO2 Emissions over Nights, Days, and Seasons (ASCENDS) satellite mission will provide greater coverage in cloudy regions, at high latitudes, and at night than passive satellite systems, as well as high precision and accuracy. In a novel approach to quantifying the ability of satellite column measurements to constrain CO2 fluxes, we use a portable library of footprints (surface influence functions) generated by the WRF-STILT Lagrangian transport model in a regional Bayesian synthesis inversion. The regional Lagrangian framework is well suited to make use of ASCENDS observations to constrain fluxes at high resolution, in this case at 1 degree latitude x 1 degree longitude and weekly for North America. We consider random measurement errors only, modeled as a function of mission and instrument design specifications along with realistic atmospheric and surface conditions. We find that the ASCENDS observations could potentially reduce flux uncertainties substantially at biome and finer scales. At the 1 degree x 1 degree, weekly scale, the largest uncertainty reductions, on the order of 50 percent, occur where and when there is good coverage by observations with low measurement errors and the a priori uncertainties are large. Uncertainty reductions are smaller for a 1.57 micron candidate wavelength than for a 2.05 micron wavelength, and are smaller for the higher of the two measurement error levels that we consider (1.0 ppm vs. 0.5 ppm clear-sky error at Railroad Valley, Nevada). Uncertainty reductions at the annual, biome scale range from 40 percent to 75 percent across our four instrument design cases, and from 65 percent to 85 percent for the continent as a whole. Our uncertainty reductions at various scales are substantially smaller than those from a global ASCENDS inversion on a coarser grid, demonstrating how quantitative results can depend on inversion methodology. The a posteriori flux uncertainties we obtain, ranging from 0.01 to 0.06 Pg C yr-1 across the biomes, would meet requirements for improved understanding of long-term carbon sinks suggested by a previous study.

remote sensing↗

Spatially-Explicit Monitoring of Crop Photosynthetic Capacity Through the Use of Space-Based Chlorophyll Fluorescence Data

Plant functional traits such as photosynthetic capacity are critical parameters for terrestrial biosphere models. However, their spatial and temporal characteristics are still poorly represented. In this study, we used satellite observations of sun-induced fluorescence (SIF) to estimate top-of-canopy photosynthetic capacity (maximum carboxylation rate, V(sub cmax) at a reference temperature of 25 deg C) for crops, which was in turn utilized to simulate regional gross primary production (GPP). We first estimate the key parameter, V(sub cmax), in the widely-used FvCB photosynthesis model using field measurements of CO2 and water fluxes during 2007-2012 at seven crop eddy covariance flux sites over the US Corn Belt. The results showed that satellite far-red SIF retrievals have a stronger link to V(sub cmax) at the seasonal scale (R(squared)=0.70 for C4 and R(squared)=0.63 for C3 crop) as compared with widely-used vegetation indices. We calibrate an empirical model linking V(sub cmax) with SIF that was used to estimate spatially and temporally varying crop V(sub cmax) for the US Corn Belt region. The resulting V(sub cmax) maps are used together with meteorological data from MERRA reanalysis data and vegetation structural parameters derived from the satellite- based spectral reflectance data to constrain the Soil-Canopy Observation of Photosynthesis and Energy (SCOPE) balance model in order to estimate regional crop GPP. Our results show a substantial improvement in the seasonal and spatial patterns of cropland GPP when compared with crop yield inventory data. The evaluation with tall tower atmospheric CO2 measurements further supports our estimation of spatiotemporal V(sub cmax) from space-borne SIF. Considering that SIF has a direct link to photosynthetic activity, our findings highlight the potential to infer regional V(sub cmax) using remotely sensed SIF data and to use this information for a better quantification of regional cropland carbon cycles.

Drought↗

Determining Stress Orientation in Rock Valley, Nevada, Using Ambient Seismic Noise

The stress field and the mechanical properties of rocks are important to consider for nuclear explosion monitoring due to their effect on seismic wave radiation from earthquakes and explosions. At the Rock Valley Direct Comparison site, the regional orientation of the maximum horizontal compressive stress (SH max ) is well constrained, but it is unknown whether there are local heterogeneities. Here, I show that stress-induced anisotropy in nonlinear elasticity can be used to estimate the orientation of SH max . Rocks have compliant internal contacts, such as fractures and mineral grain boundaries, that respond to applied strains more strongly than individual mineral crystals. This strain response is asymmetric between compressive and dilatational strains and is affected by anisotropy in the ambient stress field. Traditional seismic velocity measurements are of strain-averaged velocity, which is less sensitive than nonlinear elasticity to the behavior of compliant internal contacts and therefore the stress field and fracture behavior. My results show that the orientation of SH max measured using ambient seismic noise aligns with regional estimates, with some heterogeneity that can also be explained by the limitations of the seismic array. I demonstrate the potential of a passive technique for monitoring the stress field in places that currently lack local measurements.

58 GEOSCIENCES↗

Remote sensing-based vegetation and soil moisture constraints reduce irrigation estimation uncertainty

Understanding the human water footprint and its impact on the hydrological cycle is essential to inform water management under climate change. Despite efforts in estimating irrigation water withdrawals in earth system models, uncertainties and discrepancies exist within and across modeling systems conditioned by model structure, irrigation parameterization, and the choice of input datasets. Achieving model reliability could be much more challenging for data-sparse regions, given limited access to ground truth for parameterization and validation. Here, we demonstrate the potential of utilizing remotely sensed vegetation and soil moisture observations in constraining irrigation estimation in the Noah-MP land surface model. Results indicate that the two constraints together can effectively reduce model sensitivity to the choice of irrigation parameterization by 7%–43%. It also improves the characterization of the spatial patterns of irrigation and its impact on evapotranspiration and surface soil moisture by correcting for vegetation conditions and irrigation timing. This study highlights the importance of utilizing remotely sensed soil moisture and vegetation measurements in detecting irrigation signals and correcting for vegetation growth. Integrating the two remote sensing datasets into the model provides an effective and less feature engineered approach to constraining the uncertainty of irrigation modeling. Such strategies can be potentially transferred to other modeling systems and applied to regions across the globe.

Wanshu Nie↗

Constraint energy minimizing generalized multiscale finite element method for multi-continuum Richards equations

In fluid flow simulation, the multi-continuum model is a useful strategy. When the heterogeneity and contrast of coefficients are high, the system becomes multiscale, and some kinds of reduced order methods are demanded. Combining these techniques with nonlinearity, we will consider in this paper a dual-continuum model which is generalized as a multi-continuum model for a coupled system of nonlinear Richards equations as unsaturated flows, in complex heterogeneous fractured porous media; and we will solve it by a novel multiscale approach utilizing the constraint energy minimizing generalized multiscale finite element method (CEM-GMsFEM). In particular, such a nonlinear system will be discretized in time and then linearized by Picard iteration (whose global convergence is proved theoretically). Subsequently, we tackle the resulting linearized equations by the CEM-GMsFEM and obtain proper offline multiscale basis functions to span the multiscale space (which contains the pressure solution). More specifically, we first introduce two new sources of samples, and the GMsFEM is used over each coarse block to build local auxiliary multiscale basis functions via solving local spectral problems, that are crucial for detecting high-contrast channels. Second, per oversampled coarse region, local multiscale basis functions are created through the CEM as constrainedly minimizing an energy functional. Various numerical tests for our approach reveal that the error converges with the coarse-grid size and that only few oversampling layers as well as basis functions are needed.

97 MATHEMATICS AND COMPUTING↗

Vegetation responses to air dryness amplify future land surface warming

Temperature exerts a first-order control on vegetation photosynthesis and transpiration. Yet most studies investigating temperature impacts on plants rely on near-surface air temperature, rather than canopy temperature—the temperature plants actually experience. Because canopy temperature directly regulates ecosystem function, it provides a more accurate measure of vegetation–climate interactions. Combining Earth System Model (ESM) simulations and satellite observations in a dual emergent constraint, here we show that canopy temperature is projected to increase substantially more than air temperature (~0.11-degrees more or a 16% increase in their difference) over the 21st century. The ESM ensemble median fails to capture these stronger increases in the majority of vegetated regions. We find that the largest projected increases in the difference between canopy and air temperature are predicted to occur in regions where elevating moisture stress—particularly rising vapor pressure deficit—increasingly constrains vegetation growth and transpiration. This implies that future warming will impose stronger constraints on plant function than currently estimated. Relying on air temperature alone will therefore lead to systematic underestimation of temperature effects on photosynthesis, vegetation growth, and the land carbon sink. Accurate representation of canopy temperature in ESMs is thus essential to improve projections of ecosystem responses and feedbacks to climate change.

Green, Julia K↗

OzDES Reverberation Mapping Program: Hβ lags from the 6-yr survey

ABSTRACT Reverberation mapping measurements have been used to constrain the relationship between the size of the broad-line region and luminosity of active galactic nuclei (AGN). This R–L relation is used to estimate single-epoch virial black hole masses, and has been proposed to use to standardize AGN to determine cosmological distances. We present reverberation measurements made with Hβ from the 6-yr Australian Dark Energy Survey (OzDES) Reverberation Mapping Program. We successfully recover reverberation lags for eight AGN at 0.12 < z < 0.71, probing higher redshifts than the bulk of Hβ measurements made to date. Our fit to the R–L relation has a slope of α = 0.41 ± 0.03 and an intrinsic scatter of σ = 0.23 ± 0.02 dex. The results from our multi-object spectroscopic survey are consistent with previous measurements made by dedicated source-by-source campaigns, and with the observed dependence on accretion rate. Future surveys, including LSST, TiDES, and SDSS-V, which will be revisiting some of our observed fields, will be able to build on the results of our first-generation multi-object reverberation mapping survey.

79 ASTRONOMY AND ASTROPHYSICS↗

Measurement of photonuclear jet production in ultraperipheral Pb + Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV with the ATLAS detector

In ultrarelativistic heavy ion collisions at the LHC, each nucleus acts a sources of high-energy real photons that can scatter off the opposing nucleus in ultraperipheral photonuclear (𝛾 + 𝐴) collisions. Hard scattering processes initiated by the photons in such collisions provide a novel method for probing nuclear parton distributions in a kinematic region not easily accessible to other measurements. ATLAS has measured production of dijet and multijet final states in ultraperipheral Pb + Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV using a dataset recorded in 2018 with an integrated luminosity of 1.72 nb −1 . Photonuclear final states are selected by requiring a rapidity gap in the photon direction; this selects events where one of the outgoing nuclei remains intact. Jets are reconstructed using the anti-𝑘 t algorithm with radius parameter, 𝑅 = 0.4. Triple-differential cross sections, unfolded for detector response, are measured and presented using two sets of kinematic variables. The first set consists of the total transverse momentum (𝐻 T ), rapidity, and mass of the jet system. The second set uses 𝐻 T and particle-level nuclear and photon parton momentum fractions, 𝑥 A and 𝑧 𝛾 , respectively. The results are compared with leading-order perturbative QCD calculations of photonuclear jet production cross sections, where all leading order predictions using existing fits fall below the data in the shadowing region. More detailed theoretical comparisons will allow these results to strongly constrain nuclear parton distributions, and these data provide results from the LHC directly comparable to early physics results at the planned Electron-Ion Collider.

Parton distribution functions↗

Coulomb excitation of 80 Sr and the limits of the 𝑁 ≈ 𝑍 ≈ 40 island of deformation

The region of 𝑁 ≈ 𝑍 ≈ 40 has long been associated with strongly deformed nuclear configurations. The presence of this strong deformation was recently confirmed through lifetime measurements in 𝑁 ≈ 𝑍 Sr and Zr nuclei. Theoretically, however, these nuclei present a challenge due to the vast valence space required to incorporate all necessary orbitals. Recent state-of-the-art predictions indicate a near axial prolate deformation for 𝑁 = 𝑍 and 𝑁 = 𝑍+2 nuclei between 𝑁 = 𝑍 = 36 and 𝑁 = 𝑍 = 40. In this work we investigate the shores of this island of deformation through a sub-barrier Coulomb excitation study of the 𝑁 = 𝑍+4 nucleus, 80 Sr . Extracting a spectroscopic quadrupole moment of 𝑄 𝑠 ⁡($2^+_1$)= $0.5^{+0.8}_{−0.9}$𝑒⁢b, we find that 80 Sr is inconsistent with significant axial prolate deformation with a significance of 1.5⁢𝜎. This result, albeit with a large uncertainty, indicates that the predicted region of strong prolate deformation around 𝑁 = 𝑍 = 40 is tightly constrained to the quartet of nuclei: 76,78 Sr and 78,80 Zr .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nonlinear Optimal Control of Electron Dynamics Within Hartree-Fock Theory

Consider the problem of determining the optimal applied electric field to drive a molecule from an initial state to a desired target state. For even moderately sized molecules, solving this problem directly using the exact equations of motion—the time-dependent Schrödinger equation (TDSE)—is numerically intractable. Here, we present a solution of this problem within time-dependent Hartree-Fock (TDHF) theory, a mean field approximation of the TDSE. Optimality is defined in terms of minimizing the total control effort while maximizing the overlap between desired and achieved target states. We frame this problem as an optimization problem constrained by the nonlinear TDHF equations; we solve it using trust region optimization with gradients computed via a custom-built adjoint state method. For three molecular systems, we show that with very small neural network parametrizations of the control, our method yields solutions that achieve desired targets within acceptable constraints and tolerances.

97 MATHEMATICS AND COMPUTING↗

Curie isotherm surfaces inferred from high-altitude magnetic anomaly data

Two-dimensional Curie depth models along two sections across the western United States are derived using an equivalent layer magnetization model derived from Magsat data, and the results are used to constrain finite element thermal models developed by Mayhew and Majer (1980). Regional heat flow variations predicted by the models compare favorably with those inferred from direct measurements. The methodology is applied to upward-continued aeromagnetic data, and the results are found to be in good agreement with previous Curie depth estimates and average measured heat flow.

Mayhew, M. A.↗

On theories of rotating convection zones

It is shown that the time rate of change brought about by turbulent convective motions in the angular momentum of a thin spherical shell is such as to increase the angular velocity of the lower part of the solar convection zone (SCZ) and to decrease the angular velocity of the upper part. Arguments are presented in favor of the following very tentative model of rotation in the SCZ: (1) the lower region is in weaker differential rotation than the surface and not constrained by the Taylor-Proudman theorem, and (2) the observed solar differential rotation at the surface is generated as the SCZ relaxes from the state in the lower part to the state at the surface. In the upper and lower layers of the SCZ, angular-momentum conservation between the turbulent motions and viscous stresses leads to an angular velocity increasing upward.

Durney, B. R.↗

Thematic mapper study of Alaskan ophiolites

The combinations of Thematic Mapper (TM) bands that best distinguish basalts of the Brooks Range ophiolites were determined. Geochemical analyses, including major, trace, and rare earth elements (REE), are being done in order to study the significance of TM spectral variations that were observed within some of the sampled rock units. An image of the topography of the western Brooks Range and Colville Basin was constructed. Elevation data for the rest of Northern Alaska are being acquired to expand the area covered by the topography image. Two balanced cross sections (one along the eastern margin, the other along the western margin of the Brooks Range) are being constructed, using the techniques of fault-bend and fault-propagation folding. These are being used to obtain regional shortening estimates for the Brooks Range in an attempt to constrain tectonic models for the evolution of Northern Alaska. The TM data are being used to confirm reconnaissance maps and to obtain structural data where no maps exist. Along with the TM data, digital topography, seismic reflection profiles, and magnetic and gravity surveys are examined to better understand the evolution of the Colville Basin, north of the Brooks Range.

Bird, J. M.↗

Dynamo magnetic-field generation in turbulent accretion disks

Magnetic fields can play important roles in the dynamics and evolution of accretion disks. The presence of strong differential rotation and vertical density gradients in turbulent disks allows the alpha-omega dynamo mechanism to offset the turbulent dissipation and maintain strong magnetic fields. It is found that MHD dynamo magnetic-field normal modes in an accretion disk are highly localized to restricted regions of a disk. Implications for the character of real, dynamically constrained magnetic fields in accretion disks are discussed. The magnetic stress due to the mean magnetic field is found to be of the order of a viscous stress. The dominant stress, however, is likely to come from small-scale fluctuating magnetic fields. These fields may also give rise to energetic flares above the disk surface, providing a possible explanation for the highly variable hard X-ray emission from objects like Cyg X-l.

Stepinski, T. F.↗

Large scale nonlinear numerical optimal control for finite element models of flexible structures

This paper discusses the development of large scale numerical optimal control algorithms for nonlinear systems and their application to finite element models of structures. This work is based on our expansion of the optimal control algorithm (DDP) in the following steps: improvement of convergence for initial policies in non-convex regions, development of a numerically accurate penalty function method approach for constrained DDP problems, and parallel processing on supercomputers. The expanded constrained DDP algorithm was applied to the control of a four-bay, two dimensional truss with 12 soft members, which generates geometric nonlinearities. Using an explicit finite element model to describe the structural system requires 32 state variables and 10,000 time steps. Our numerical results indicate that for constrained or unconstrained structural problems with nonlinear dynamics, the results obtained by our expanded constrained DDP are significantly better than those obtained using linear-quadratic feedback control.

Shoemaker, Christine A.↗

Topography and tectonics of the central New Madrid seismic zone: Results of numerical experiements using a three-dimensional boundary element program

We present results of a series of numerical experiments designed to test hypothetical mechanisms that derive deformation in the New Madrid seismic zone. Experiments are constrained by subtle topography and the distribution of seismicity in the region. We use a new boundary element algorithm that permits calcuation of the three-dimensional deformation field. Surface displacement fields are calculated for the New Madrid zone under both far-field (plate tectonics scale) and locally derived driving strains. Results demonstrate that surface displacement fields cannot distinguish between either a far-field simple or pure shear strain field or one that involves a deep shear zone beneath the upper crustal faults. Thus, neither geomorphic nor geodetic studies alone are expected to reveal the ultimate driving mechanism behind the present-day deformation. We have also tested hypotheses about strain accommodation within the New Madrid contractional step-over by including linking faults, two southwest dipping and one vertical, recently inferred from microearthquake data. Only those models with step-over faults are able to predict the observed topography. Surface displacement fields for long-term, relaxed deformation predict the distribution of uplift and subsidence in the contractional step-over remarkably well. Generation of these displacement fields appear to require slip on both the two northeast trending vertical faults and the two dipping faults in the step-over region, with very minor displacements occurring during the interseismic period when the northeast trending vertical faults are locked. These models suggest that the gently dippling central step-over fault is a reverse fault and that the steeper fault, extending to the southeast of the step-over, acts as a normal fault over the long term.

Gomberg, Joan↗

The Influence of surface binding energy on sputtering in models of the sodium exosphere of Mercury

We have simulated the sodium (Na) exosphere of Mercury to show how the exosphere is affected by the assumed surface binding energy (SBE) of Na in the sputtered component. We constrained ion precipitation onto the surface using distributions for the cusp regions that are consistent with measurements by the MErcury Surface, Space ENvironment, GEochemistry, and Ranging Fast Imaging Plasma Spectrometer instrument. We have simulated sputtering with SBEs of 0.27, 2.6, 4.4, and 7.9 eV, with the lowest value commonly used in exosphere models and the highest from recent molecular dynamics calculations for the Na-bearing feldspar end-member, albite. A gradual change in the exosphere is seen as the yield decreases and the ejecta energy increases with increasing SBE. We describe the corresponding exosphere source functions for ion sputtering (IS), as well as for the previously studied processes of micrometeoroid impact vaporization and photon-stimulated desorption (PSD), along with their release energy distributions and spatial distributions. We have summed the contributions of the various source processes to explain how and when the different sources can be distinguished by observations. The modeled exosphere scale heights range from 72 km for PSD to over 1000 km for IS using a SBE of 7.9 eV. We find that the processes responsible for generating Mercury's Na exosphere are separable by measuring line-of-sight column densities tangent to the planet at various altitudes and positions around the planet. Our initial results are consistent with the Na being sputtered from a high-SBE material such as feldspar, which has been predicted to be abundant on the Mercury's surface.

Rosemary M. Killen↗