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 955 records · Page 53

A Spatiotemporal-Aware Weighting Scheme for Improving Climate Model Ensemble Predictions

Multimodel ensembling has been widely used to improve climate model predictions, and the improvement strongly depends on the ensembling scheme. In this work, we propose a Bayesian neural network (BNN) ensembling method, which combines climate models within a Bayesian model averaging framework, to improve the predictive capability of model ensembles. Our proposed BNN approach calculates spatiotemporally varying model weights and biases by leveraging individual models' simulation skill, calibrates the ensemble prediction against observations by considering observation data uncertainty, and quantifies epistemic uncertainty when extrapolating to new conditions. More importantly, the BNN method provides interpretability about which climate model contributes more to the ensemble prediction at which locations and times. Thus, beyond its predictive capability, the method also brings insights and understanding of the models to guide further model and data development. In this study, we design experiments using an ensemble of CMIP6 climate model simulations to illustrate the BNN ensembling method's capability with respect to prediction accuracy, interpretability, and uncertainty quantification (UQ). We demonstrate that BNN can correctly assign larger weights to the regions and seasons where the individual model fits the observation better. Moreover, its offered interpretability is consistent with our understanding of localized climate model performance. Additionally, BNN shows an increasing uncertainty when the prediction is farther away from the period with constrained data, which appropriately reflects our trustworthiness of the models in the changing climate.

54 ENVIRONMENTAL SCIENCES↗

The effects of missing data on global ozone estimates

The effects of missing data and model truncation on estimates of the global mean, zonal distribution, and global distribution of ozone are considered. It is shown that missing data can introduce biased estimates with errors that are not accounted for in the accuracy calculations of empirical modeling techniques. Data-fill techniques are introduced and used for evaluating error bounds and constraining the estimate in areas of sparse and missing data. It is found that the accuracy of the global mean estimate is more dependent on data distribution than model size. Zonal features can be accurately described by 7th order models over regions of adequate data distribution. Data variance accounted for by higher order models appears to represent climatological features of columnar ozone rather than pure error. Data-fill techniques can prevent artificial feature generation in regions of sparse or missing data without degrading high order estimates over dense data regions.

Drewry, J. W.↗

Pattern Search Algorithms for Bound Constrained Minimization

We present a convergence theory for pattern search methods for solving bound constrained nonlinear programs. The analysis relies on the abstract structure of pattern search methods and an understanding of how the pattern interacts with the bound constraints. This analysis makes it possible to develop pattern search methods for bound constrained problems while only slightly restricting the flexibility present in pattern search methods for unconstrained problems. We prove global convergence despite the fact that pattern search methods do not have explicit information concerning the gradient and its projection onto the feasible region and consequently are unable to enforce explicitly a notion of sufficient feasible decrease.

Lewis, Robert Michael↗

NASA Supercomputer Improves Prospects for Ocean Climate Research

Estimates of ocean circulation constrained by in situ and remotely sensed observations have become routinely available during the past five years, and they are being applied to myriad scientific and operational problems [Stammer et al.,2002]. Under the Global Ocean Data Assimilation Experiment (GODAE), several regional and global estimates have evolved for applications in climate research, seasonal forecasting, naval operations, marine safety, fisheries,the offshore oil industry, coastal management, and other areas. This article reports on recent progress by one effort, the consortium for Estimating the Circulation and Climate of the Ocean (ECCO), toward a next-generation synthesis of ocean and sea-ice data that is global, that covers the full ocean depth, and that permits eddies.

Estimating the Circulation and Climate of the Ocea↗

A Test of Sensitivity to Convective Transport in a Global Atmospheric CO2 Simulation

Two approximations to convective transport have been implemented in an offline chemistry transport model (CTM) to explore the impact on calculated atmospheric CO2 distributions. GlobalCO2 in the year 2000 is simulated using theCTM driven by assimilated meteorological fields from the NASA s Goddard Earth Observation System Data Assimilation System, Version 4 (GEOS-4). The model simulates atmospheric CO2 by adopting the same CO2 emission inventory and dynamical modules as described in Kawa et al. (convective transport scheme denoted as Conv1). Conv1 approximates the convective transport by using the bulk convective mass fluxes to redistribute trace gases. The alternate approximation, Conv2, partitions fluxes into updraft and downdraft, as well as into entrainment and detrainment, and has potential to yield a more realistic simulation of vertical redistribution through deep convection. Replacing Conv1 by Conv2 results in an overestimate of CO2 over biospheric sink regions. The largest discrepancies result in a CO2 difference of about 7.8 ppm in the July NH boreal forest, which is about 30% of the CO2 seasonality for that area. These differences are compared to those produced by emission scenario variations constrained by the framework of Intergovernmental Panel on Climate Change (IPCC) to account for possible land use change and residual terrestrial CO2 sink. It is shown that the overestimated CO2 driven by Conv2 can be offset by introducing these supplemental emissions.

Bian, H.↗

Robust Path Planning and Feedback Design Under Stochastic Uncertainty

Autonomous vehicles require optimal path planning algorithms to achieve mission goals while avoiding obstacles and being robust to uncertainties. The uncertainties arise from exogenous disturbances, modeling errors, and sensor noise, which can be characterized via stochastic models. Previous work defined a notion of robustness in a stochastic setting by using the concept of chance constraints. This requires that mission constraint violation can occur with a probability less than a prescribed value.In this paper we describe a novel method for optimal chance constrained path planning with feedback design. The approach optimizes both the reference trajectory to be followed and the feedback controller used to reject uncertainty. Our method extends recent results in constrained control synthesis based on convex optimization to solve control problems with nonconvex constraints. This extension is essential for path planning problems, which inherently have nonconvex obstacle avoidance constraints. Unlike previous approaches to chance constrained path planning, the new approach optimizes the feedback gain as wellas the reference trajectory.The key idea is to couple a fast, nonconvex solver that does not take into account uncertainty, with existing robust approaches that apply only to convex feasible regions. By alternating between robust and nonrobust solutions, the new algorithm guarantees convergence to a global optimum. We apply the new method to an unmanned aircraft and show simulation results that demonstrate the efficacy of the approach.

autonomuys vehicles↗

Constraints on Vesta's Interior Structure Using Gravity and Shape Models from the Dawn Mission

We use the shape and gravity field of Vesta determined from observations of the Dawn spacecraft to place constraints on the asteroid's interior structure. We compute a three-layer interior structure model by minimizing the power of the residual gravity anomaly. The densities of the mantle and crust are based on constraints derived from the Howardite-Eucrite-Diogenite (HED) meteorites. Vesta's present-day shape is not in hydrostatic equilibrium. The Rheasilvia and Veneneia impact basins have a large effect on Vesta's shape and are the main source of deviation from hydrostatic shape. Constraining a pre-giant-impact rotation rate and orientation of the spin axis from an ellipsoidal fit to the parts of Vesta unaffected by the giant impacts, and using the theory of figure, we can constrain the shape of the core. Our solution for Vesta's crust-mantle interface reveals a belt of thick crust around Rheasilvia and Veneneia. The thinnest crust is in the floor of the two basins and in the Vestalia Terra region. Our solution does not reveal an uplift of the crust-mantle boundary to the surface in the largest basins. This, together with the lack of olivine detected by the Visible and Infrared Spectrometer (VIR) data in Rheasilvia and Veneneia, indicates that Vesta's presumed olivine mantle was either not brought to the surface by these large impacts or was covered by ejecta from subsequent impacts.

Interiors↗

Key Differences in Operating a Rover on the Moon vs. Mars

The command and control model for spacecraft operations, as well as the distribution of tasks between ground assets and in space assets, whether with a crew or solely robotic, is fundamentally constrained by the round trip light time between the space asset and the control facility (presumably on Earth, though not required). For an asset on Mars, the round trip light time varies, from roughly fourteen minutes to up to forty minutes. For a Lunar asset the round-trip light time is measured in only a few seconds, but current communications systems may more than double the latency with system overhead. For a Lunar Asset the total command latency may range from six seconds to more than forty, depending on communications overhead and data rates. Further, these variables are not always predictable, thus complicating operations. There are several differentiating factors for Lunar vs. Mars operations, Round trip light time/Atmosphere/Lighting and ShadowsTerrain type and knowledge/Round trip light time has implications for the distribution of tasks between ground and in space assets. Even at Lunar Distances, the combination of round trip light time plus communications systems overhead does not enable joy stick driving of a rover. The best that can be done, if driving from Earth, is near real time command and control. By 2030, driving from in space may be possible. Productivity on Mars requires either long operational sequences of commands, as is done for current rovers such as Curiosity, significant autonomous capability or, as may be possible by 2030, command and control support from space. Another implication of the long round trip light time from Earth to Mars, is that flight software functions must be resident on the in space asset. On the Moon, there is considerably more flexibility, enabling processing functions, to be resident on Earth or in space. This provides the opportunity to take advantage of the considerable processing power available on the ground, but may be constrained by data rates. On the Moon, for practical operational purposes, there is no atmosphere. Hence there is no scattering of light in the shadows. This has implications for image interpretation and driving near the poles. The Moon has permanently shadowed regions (PSR), unique terrain with unknown surface properties. With no scattering of light in shadows, driving on the Moon, particularly at the poles, where we have strong evidence of water, may prove to be hazardous and complex, requiring non-optical sensors, such as LIDAR.

Trimble, Jay↗

A Global Analysis of Wildfire Smoke Injection Heights Derived from Space-Based Multi-Angle Imaging

We present an analysis of over 23,000 globally distributed wildfire smoke plume injection heights derived from Multi-angle Imaging SpectroRadiometer (MISR) space-based, multi-angle stereo imaging. Both pixel-weighted and aerosol optical depth (AOD)-weighted results are given, stratified by region, biome, and month or season. This offers an observational resource for assessing first-principle plume-rise modelling, and can provide some constraints on smoke dispersion modelling for climate and air quality applications. The main limitation is that the satellite is in a sun-synchronous orbit, crossing the equator at about 10:30 a.m. local time on the day side. Overall, plumes occur preferentially during the northern mid-latitude burning season, and the vast majority inject smoke near-surface. However, the heavily forested regions of North and South America, and Africa produce the most frequent elevated plumes and the highest AOD values; some smoke is injected to altitudes well above 2 km in nearly all regions and biomes. Planetary boundary layer (PBL) versus free troposphere injection is a critical factor affecting smoke dispersion and environmental impact, and is affected by both the smoke injection height and the PBL height; an example assessment is made here, but constraining the PBL height for this application warrants further work.

Val Martin, Maria↗

Untangling causality in midlatitude aerosol–cloud adjustments

Aerosol–cloud interactions represent the leading uncertainty in our ability to infer climate sensitivity from the observational record. The forcing from changes in cloud albedo driven by increases in cloud droplet number (Nd) (the first indirect effect) is confidently negative and has narrowed its probable range in the last decade, but the sign and strength of forcing associated with changes in cloud macrophysics in response to aerosol (aerosol–cloud adjustments) remain uncertain. This uncertainty reflects our inability to accurately quantify variability not associated with a causal link flowing from the cloud microphysical state to the cloud macrophysical state. Once variability associated with meteorology has been removed, covariance between the liquid water path (LWP) averaged across cloudy and clear regions (here characterizing the macrophysical state) and Nd (characterizing the microphysical) is the sum of two causal pathways linking Nd to LWP: Nd altering LWP (adjustments) and precipitation scavenging aerosol and thus depleting Nd. Only the former term is relevant to constraining adjustments, but disentangling these terms in observations is challenging. We hypothesize that the diversity of constraints on aerosol–cloud adjustments in the literature may be partly due to not explicitly characterizing covariance flowing from cloud to aerosol and aerosol to cloud. Here, we restrict our analysis to the regime of extratropical clouds outside of low-pressure centers associated with cyclonic activity. Observations from MAC-LWP (Multisensor Advanced Climatology of Liquid Water Path) and MODIS are compared to simulations in the Met Office Unified Model (UM) GA7.1 (the atmosphere model of HadGEM3-GC3.1 and UKESM1). The meteorological predictors of LWP are found to be similar between the model and observations. There is also agreement with previous literature on cloud-controlling factors finding that increasing stability, moisture, and sensible heat flux enhance LWP, while increasing subsidence and sea surface temperature decrease it. A simulation where cloud microphysics are insensitive to changes in Nd is used to characterize covariance between Nd and LWP that is induced by factors other than aerosol–cloud adjustments. By removing variability associated with meteorology and scavenging, we infer the sensitivity of LWP to changes in Nd. Application of this technique to UM GA7.1 simulations reproduces the true model adjustment strength. Observational constraints developed using simulated covariability not induced by adjustments and observed covariability between Nd and LWP predict a 25 %–30 % overestimate by the UM GA7.1 in LWP change and a 30 %–35 % overestimate in associated radiative forcing.

Daniel T. McCoy↗

Using Mineralogy to Interpret Martian Geologic History

The surface mineralogy of Mars is key to interpreting the planet’s geologic history and constraining when and where Mars may have been habitable to microbial life. Orbital mineralogical data from infrared spectrometers collected over the last two decades have demonstrated that Mars is more than a basaltic planet. Thermal infrared (TIR) emission spectrometers have identified regional variations in igneous minerals, iron oxides, and chloride salts, showing evidence for magmatic evolution and aqueous alteration. Orbital visible/shortwave infrared reflectance spectra provide detailed information about the aqueous history of Mars through the identification of clay minerals, sulfates, carbonates, zeolites, and amorphous materials in ancient ~3-4-billion-year-old terrains. Mineralogical data from rovers allow us to better characterize mineral formation mechanisms by identifying mineral assemblages and placing them in geologic context using outcrop- to grain-scale images. TIR and Mössbauer data from the Mars Exploration Rovers helped identify an ancient dune-interdune playa environment with multiple episodes of groundwater at Meridiani Planum and volcanism-induced hydrothermalism at Gusev crater. The CheMin X-ray diffractometer on the Mars Science Laboratory Curiosity rover in Gale crater is the first fully quantitative mineralogical instrument sent to another planetary surface. CheMin data allow the quantification of minerals and X-ray amorphous materials with a mineral detection limit of <1 wt.% and provide crystal chemistry of major phases from refined unit-cell parameters. CheMin has revealed local mineralogical changes in ancient sedimentary rocks not seen from orbit and helped identify habitable environments. Throughout the 600+ m of vertical stratigraphy studied so far, changes in clay mineralogy, Fe-oxides/oxyhydroxides, sulfates, and carbonates demonstrate a long history of surface and groundwater at Gale crater with variable salinity, pH, Eh, and temperature. Recent CheMin data document a decrease in clay minerals and an increase in sulfate minerals in younger strata that corresponds with a change from fluvio-lacustrine to eolian depositional environments, potentially documenting a change to a drier climate across the planet.

E. B. Rampe↗

Structure of the Human BK Ion Channel in Lipid Environment

Voltage-gated and ligand-modulated ion channels play critical roles in excitable cells. To understand the interplay among voltage sensing, ligand binding, and channel opening, the structures of ion channels in various functional states and in lipid membrane environments need to be determined. Here, the random spherically constrained (RSC) single-particle cryo-EM method was employed to study human large conductance voltage- and calcium-activated potassium (hBK or hSlo1) channels reconstituted into liposomes. The hBK structure was determined at 3.5 Å resolution in the absence of Ca 2+ . Instead of the common fourfold symmetry observed in ligand-modulated ion channels, a twofold symmetry was observed in hBK in liposomes. Compared with the structure of isolated hSlo1 Ca 2+ sensing gating rings, two opposing subunits in hBK unfurled, resulting in a wider opening towards the transmembrane region of hBK. In the pore gate domain, two opposing subunits also moved downwards relative to the two other subunits.

59 BASIC BIOLOGICAL SCIENCES↗

Variability in Ice Nucleating Particles Across Greater Houston Texas

The concentration and cloud-forming potential of a region's ice nucleating particle (INP) population have uncertain impacts on deep convective clouds. Specifically, ice nucleating particles (INPs) may affect various cloud properties related to the formation, lifetime, and precipitation of deep convective clouds. As part of the U.S. Department of Energy's TRacking Aerosol and Convection interaction ExpeRiment (TRACER) campaign, researchers from Texas A&M University deployed three Davis Rotating-drum Universal-size-cut Monitoring (DRUM) samplers throughout Greater Houston, Texas from June through September 2022. Ambient particles, collected at the surface with the DRUM samplers in four aerodynamic diameter size ranges (>3, 3–1.2, 1.2–0.34, and 0.34–0.15 μm), were analyzed in offline cold-stage ice nucleation experiments. The INP population in Greater Houston is complex, varying by site and day, but can be generalized by a weak to moderately efficient mode of INPs at −24°C and an efficient mode at −15°C. Analysis reveals that supermicron particles are largely responsible for ice nucleation warmer than −20°C across the region while submicron particles dominate at temperatures colder than −20°C. Additionally, significant spatial diversity in the INP population was observed, with differences in mean nucleation temperature between sites for nearly every size cut. Although INP concentrations were typically ∼0.08 L −1 at −20°C throughout the campaign, a notable region-wide increase in INP concentration for particles freezing at temperatures warmer than −20°C occurred from mid-August to mid-September. This comprehensive characterization of Greater Houston's INP population, including spatial, temporal, and particle size variations, can help constrain ice microphysics parameterizations in weather and climate models.

Thompson, Seth A. [Texas A & M Univ., College Stat↗

Where did the tumor start? An inverse solver with sparse localization for tumor growth models

In this work, we present a numerical scheme for solving an inverse problem for parameter estimation in tumor growth models for glioblastomas, a form of aggressive primary brain tumor. The growth model is a reaction–diffusion partial differential equation (PDE) for the tumor concentration. We use a PDE-constrained optimization formulation for the inverse problem. The unknown parameters are the reaction coefficient (proliferation), the diffusion coefficient (infiltration), and the initial condition field for the tumor PDE. Segmentation of magnetic resonance imaging (MRI) scans drive the inverse problem where segmented tumor regions serve as partial observations of the tumor concentration. Like most cases in clinical practice, we use data from a single time snapshot. Moreover, the precise time relative to the initiation of the tumor is unknown, which poses an additional difficulty for inversion. We perform a frozen-coefficient spectral analysis and show that the inverse problem is severely ill-posed. We introduce a biophysically motivated regularization on the structure and magnitude of the tumor initial condition. In particular, we assume that the tumor starts at a few locations (enforced with a sparsity constraint on the initial condition of the tumor) and that the initial condition magnitude in the maximum norm is equal to one. We solve the resulting optimization problem using an inexact quasi-Newton method combined with a compressive sampling algorithm for the sparsity constraint. Our implementation uses PETSc and AccFFT libraries. We conduct numerical experiments on synthetic and clinical images to highlight the improved performance of our solver over a previously existing solver that uses standard two-norm regularization for the calibration parameters. The existing solver is unable to localize the initial condition. Our new solver can localize the initial condition and recover infiltration and proliferation. In clinical datasets (for which the ground truth is unknown), our solver results in qualitatively different solutions compared to the two-norm regularized solver.

97 MATHEMATICS AND COMPUTING↗

Shell-model studies of the astrophysical r p -process reactions S 34 ( p , γ ) Cl 35 and Cl 34 g , m ( p , γ ) Ar 35

Background: Dust grains condensed in the outflows of presolar classical novae should have been present in the protosolar nebula. Candidates for such presolar nova grains have been found in primitive meteorites and can in principle be identified by their isotopic ratios, but the ratios predicted by state-of-the-art one-dimensional hydrodynamic models are uncertain due to nuclear-physics uncertainties. Purpose: To theoretically calculate the thermonuclear rates and uncertainties of the S 34 ( p , γ ) Cl 35 and Cl 34 g , m ( p , γ ) Ar 35 reactions and investigate their impacts on the predicted 34 S/ 32 S isotopic ratio for presolar nova grains. Method: A shell-model approach in a (0 +1 ) ℏ ω model space was used to calculate the properties of resonances in the 34 S(p ,γ ) 35 Cl and 34 g,m Cl (p ,γ ) 35 Ar reactions and their thermonuclear rates. Uncertainties were estimated using a Monte Carlo method. The implications of these rates and their uncertainties on sulfur isotopic nova yields were investigated using a postprocessing nucleosynthesis code. The rates for transitions from the ground state of 34 Cl as well as from the isomeric first excited state of 34 Cl were explicitly calculated. Results: At energies in the resonance region near the proton-emission threshold, many negative-parity states appear. Energies, spectroscopic factors, and proton-decay widths are reported. The resulting thermonuclear rates are compared with previous determinations. Conclusions: The shell-model calculations alone are sufficient to constrain the variation of the 34 S/ 32 S ratios to within about 30%. Uncertainties associated with other reactions must also be considered, but in general we find that the 34 S/ 32 S ratios are not a robust diagnostic to clearly identify presolar grains made from nova ejecta.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Search for a dark leptophilic scalar produced in association with τ + τ − pair in e + e − annihilation at center-of-mass energies near 10.58 GeV

A dark leptophilic scalar ( ϕ L ) is a hypothetical particle that couples only to leptons rather than quarks. We report on a search for ϕ L in the e + e − → τ + τ − ϕ L , ϕ L → ℓ + ℓ − ( ℓ = e , μ ) process using 626 fb − 1 of data collected by the Belle experiment near the ϒ ( 4 S ) resonance. We validate the backgrounds with multiple control regions in data, using a novel multiclass multivariate event classifier. In absence of a signal, we quote upper limits at the 90% confidence level on the coupling between ϕ L and leptons. Our bounds, obtained in a blinded approach, are 19% more constraining than the previous limits, averaged over the mass range 0.04 ≤ m ϕ L ≤ 6.5 GeV . We exclude the parameter space below 4 GeV favored by measurement of the anomalous magnetic moment of the muon. Published by the American Physical Society 2024

Biswas, D. (ORCID:0000000275433471)↗

Assessing the Behavior of Moderately Volatile Elements on the Moon in Order to Constrain Processes of Magmatic Evolution and Planetary Accretion

Heterogeneous isotopic compositions of Zn and Ga in lunar basalts and highlands rocks indicate that secondary processing affected the stable isotopic composition of major lunar reservoirs. It follows that these fractionating processes must be identified and characterized in order to make robust estimates for the isotopic composition of the mare basalt source regions or the bulk silicate Moon (BSM). Contrary to previous assertions (e.g., Kato & Moynier, 2017), MVEs such as Zn and Ga do not behave similarly on the Moon, in line with differences in their chemical behavior and volatility. There is clear evidence that impact reworking fractionates zinc isotopes but, in addition, zinc isotopes are also highly sensitive to the formation of Ti-bearing oxides such as ilmenite. In contrast, gallium is far less susceptible to volatile reworking, but its isotopic composition is fractionated during incorporation into plagioclase (Wimpenny et al., 2022; Render et al., 2023). Based on these findings we conclude that current isotopic compositions of Zn and Ga in lunar rocks, such as those measured for the mare basalt suite, are unlikely to be representative of the bulk Moon when it accreted. As such, estimates for the MVE isotopic composition of the BSM derived from mare basalts may not be accurate, which may have implications for the incorporation of these data into dynamical models of the Moon’s formation. In the following, we discuss the main findings from this project in more detail, focusing on the behaviour of zinc and gallium isotopic systematics in lunar rocks and the processes that control their isotopic fractionation on the Moon.

58 GEOSCIENCES↗

Flow distributions and discharge coefficient effects for jet array impingement with initial crossflow

To model the impingement cooled mid-chord region of gas turbine airfoils in cases where an initial crossflow is present, the paper presents experimentally determined flow distributions for jet arrays with ten spanwise rows of holes, which range from uniform to highly nonuniform. The jet flow after impingement is constrained to exit in a single direction along the channel formed by the jet orifice plate and the impingement surface. The streamwise distributions and crossflow velocities are presented for ratios of the initial crossflow rate to the total jet flow rate ranging from zero to unity. For crossflow to jet velocity ratios greater than a value somewhat less than unity, jet orifice discharge coefficients do not remain constant, but decrease significantly, showing a secondary dependence on z/d, where z is the channel height and d is the jet hole diameter.

Florschuetz, L. W.↗