Search NASA⌕ Search

SEARCH · Search NASA

Results for “Constrained parameter spaces”

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 361 records · Page 20

Dynamic Rupture Modeling of the M7.1, 2019 Ridgecrest, California, Earthquake

An important part of the simulation technique implemented in the strong motion prediction recipe developed in Japan, originally proposed by Irikura and Miyake (2011) is the earthquake rupture characterization that is based on established empirical scaling relationships between source parameters and seismic moment. These scaling relationships and rupture characterization are derived from recorded data and are tested by simulations. Recent analysis of results obtained with the Recipe for the 2016 M7.1 Kumamoto earthquake performed by several researchers, indicate that in order to reproduce near-fault long-period ground motions with permanent displacements, typical for earthquakes with surface rupturing, the recipe needs modifications. In particular, improvements are required in characterizing the rupture kinematics at shallow depths, above the seismogenic zone (Irikura et al., 2019). Modeling of rupture dynamics for recorded earthquakes, and numerical simulations for different rupture scenarios, provides physical constrains on slip evolution in time and space, which can be used in adjustments of empirical source scaling relationships for large crustal earthquakes. In this study we mainly focused on rupture dynamics modeling of the 2019 M7 Ridgecrest, California earthquake. The characterization of the kinematic rupture model derived from the dynamic rupture models of this earthquake can inform further refinements of rupture modeling proposed by the Recipe. Through different rupture scenarios we investigated potential relationships between characteristics of Long Period Ground Motion Generation Area (LMGA), used in the Recipe, and asperity depth and rupture initiation.

58 GEOSCIENCES↗

Exposure and Terrestrial Histories of New Lunar and Martian Meteorites

Cosmogenic nuclide studies of lunar and Martian meteorites have contributed significantly to our understanding of these objects. By measuring a combination of cosmogenic stable- and radionuclides, we can determine a number of important properties of those meteorites. Most lunar meteorites have complex cosmic ray exposure histories, having been exposed both at some depth on the lunar surface (2 irradiation) before their ejection and as small bodies in space (4 irradiation) during transport from the Moon to the Earth. On the other hand, we have not observed evidence of complex exposure history for any Martian meteorites, so far. These exposures were then followed by residence on Earth s surface, a time commonly referred to as the terrestrial age. In addition to their complement of galactic cosmic ray (GCR) produced nuclides some lunar and Martian meteorites contain nuclides produced by solar cosmic rays (SCR). Unraveling the complex history of these objects requires the measurement of at least four cosmogenic nuclides. The specific goals of these measurements are to constrain or set limits on the following shielding or exposure parameters: (1) the depth of the sample at the time of ejection from the Moon or Mars; (2) the transit time (4 exposure age) from ejection off the lunar or Martian surface to the time of capture by the Earth and (3) the terrestrial residence time. The sum of the transit time and residence time yield an ejection age. The ejection age, in conjunction with the sample depth on the Moon or Mars, can then be used to model impact and ejection mechanisms.

Nishiizumi, K.↗

A variational approach to dynamics of flexible multibody systems

This paper presents a variational formulation of constrained dynamics of flexible multibody systems, using a vector-variational calculus approach. Body reference frames are used to define global position and orientation of individual bodies in the system, located and oriented by position of its origin and Euler parameters, respectively. Small strain linear elastic deformation of individual components, relative to their body references frames, is defined by linear combinations of deformation modes that are induced by constraint reaction forces and normal modes of vibration. A library of kinematic couplings between flexible and/or rigid bodies is defined and analyzed. Variational equations of motion for multibody systems are obtained and reduced to mixed differential-algebraic equations of motion. A space structure that must deform during deployment is analyzed, to illustrate use of the methods developed.

Wu, Shih-Chin↗

The Importance of Earth as a Meteor Detector

In order to understand the complex meteoroid environment that threatens spacecraft in near Earth space, ground based observational techniques provide us with an effective option for evaluating the important parameters of flux, speed, mass, density and directionality. While no observational method (in ]situ or ground based detectors) provides bias free measurements of direct quantities, there are several benefits to using the Earth as a detector for calibrating meteoroid models that extend to the inner Solar System and beyond Earth fs orbit. The Earth constantly interacts with meteoritic particles from comets and asteroids as they evolve from their parent objects through dynamical processes. Modeling these dynamical processes with a physics ]based approach can help reveal the causes of observational features but developing those models involves constraining sequences of poorly defined parameters. These parameters must be tuned to a particular detector that offers direct or indirect measurements of mass, density, flux and speed - the quantities used in evaluating penetration risk. The following topics will be discussed: a) Consistency between ground based and in.situ measurements of flux in Earth orbit. b) Consistency between radar and electro.optical speed determinations and agreement of these with dynamical studies. c) The fact that there are six sporadic sources which cannot be determined by integrated line of sight emissivity: Observations using IRAS or COBE cannot give you all important directionality needed for risk evaluations. d) The large collecting area of the Earth s atmosphere compared to in.situ detectors gives better number statistics in the threat regime. The NASA Meteoroid Environment Office utilizes ground based radar and electro.optical measurements to calibrate the Meteoroid Engineering Model (MEM). As MEM is an engineering design model our main emphasis is on threat size particle, speed, flux and directionality. However, MEM has acknowledged limitations particularly with regard to source strength as a function of mass and work is proceeding with improved distributions that will rectify many of these shortcomings.

Koehler, Heather↗

Domain Adaptive Graph Neural Networks for Constraining Cosmological Parameters Across Multiple Data Sets

Deep learning models have been shown to outperform methods that rely on summary statistics, like the power spectrum, in extracting information from complex cosmological data sets. However, due to differences in the subgrid physics implementation and numerical approximations across different simulation suites, models trained on data from one cosmological simulation show a drop in performance when tested on another. Similarly, models trained on any of the simulations would also likely experience a drop in performance when applied to observational data. Training on data from two different suites of the CAMELS hydrodynamic cosmological simulations, we examine the generalization capabilities of Domain Adaptive Graph Neural Networks (DA-GNNs). By utilizing GNNs, we capitalize on their capacity to capture structured scale-free cosmological information from galaxy distributions. Moreover, by including unsupervised domain adaptation via Maximum Mean Discrepancy (MMD), we enable our models to extract domain-invariant features. We demonstrate that DA-GNN achieves higher accuracy and robustness on cross-dataset tasks. Using data visualizations, we show the effects of domain adaptation on proper latent space data alignment. This shows that DA-GNNs are a promising method for extracting domain-independent cosmological information, a vital step toward robust deep learning for real cosmic survey data.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantifying the propagation of parametric uncertainty on flux balance analysis

Flux balance analysis (FBA) and associated techniques operating on stoichiometric genome-scale metabolic models play a central role in quantifying metabolic flows and constraining feasible phenotypes. At the heart of these methods lie two important assumptions: (i) the biomass precursors and energy requirements neither change in response to growth conditions nor environmental/genetic perturbations, and (ii) metabolite production and consumption rates are equal at all times (i.e., steady-state). Despite the stringency of these two assumptions, FBA has been shown to be surprisingly robust at predicting cellular phenotypes. In this paper, we formally assess the impact of these two assumptions on FBA results by quantifying how uncertainty in biomass reaction coefficients, and departures from steady-state due to temporal fluctuations could propagate to FBA results. In the first case, conditional sampling of parameter space is required to re-weigh the biomass reaction so as the molecular weight remains equal to 1 g mmol –1 , and in the second case, metabolite (and elemental) pool conservation must be imposed under temporally varying conditions. Results confirm the importance of enforcing the aforementioned constraints and explain the robustness of FBA biomass yield predictions.

59 BASIC BIOLOGICAL SCIENCES↗

Dark energy survey year 3 results: Cosmology with peaks using an emulator approach

ABSTRACT We constrain the matter density Ωm and the amplitude of density fluctuations σ8 within the ΛCDM cosmological model with shear peak statistics and angular convergence power spectra using mass maps constructed from the first three years of data of the Dark Energy Survey (DES Y3). We use tomographic shear peak statistics, including cross-peaks: peak counts calculated on maps created by taking a harmonic space product of the convergence of two tomographic redshift bins. Our analysis follows a forward-modelling scheme to create a likelihood of these statistics using N-body simulations, using a Gaussian process emulator. We take into account the uncertainty from the remaining, largely unconstrained ΛCDM parameters (Ωb, ns, and h). We include the following lensing systematics: multiplicative shear bias, photometric redshift uncertainty, and galaxy intrinsic alignment. Stringent scale cuts are applied to avoid biases from unmodelled baryonic physics. We find that the additional non-Gaussian information leads to a tightening of the constraints on the structure growth parameter yielding $S_8~\equiv ~\sigma _8\sqrt{\Omega _{\mathrm{m}}/0.3}~=~0.797_{-0.013}^{+0.015}$ (68 per cent confidence limits), with a precision of 1.8 per cent, an improvement of 38 per cent compared to the angular power spectra only case. The results obtained with the angular power spectra and peak counts are found to be in agreement with each other and no significant difference in S8 is recorded. We find a mild tension of $1.5 \, \sigma$ between our study and the results from Planck 2018, with our analysis yielding a lower S8. Furthermore, we observe that the combination of angular power spectra and tomographic peak counts breaks the degeneracy between galaxy intrinsic alignment AIA and S8, improving cosmological constraints. We run a suite of tests concluding that our results are robust and consistent with the results from other studies using DES Y3 data.

79 ASTRONOMY AND ASTROPHYSICS↗

Planck constraints and gravitational wave forecasts for primordial black hole dark matter seeded by multifield inflation

Here, we perform a Markov Chain Monte Carlo analysis of a simple yet generic multifield inflation model characterized by two scalar fields coupled to each other and nonminimally coupled to gravity, fit to Planck 2018 cosmic microwave background (CMB) data. In particular, model parameters are constrained by data on the amplitude of the primordial power spectrum of scalar curvature perturbations on CMB scales A s , the spectral index n s , and the ratio of power in tensor to scalar modes r, with a prior that the primordial power spectrum should also lead to primordial black hole (PBH) production sufficient to account for the observed dark matter abundance. We find that n s in particular controls the constraints on our model. Whereas previous studies of PBH formation from an ultraslow-roll phase of inflation have highlighted the need for at least one model parameter to be highly fine-tuned, we identify a degeneracy direction in parameter space such that shifts by ~10% of one parameter can be compensated by comparable shifts in other parameters while preserving a close fit between model predictions and observations. Furthermore, we find this allowed parameter region produces observable gravitational wave signals in the frequency ranges to which upcoming experiments are projected to be sensitive, including Advanced LIGO and Virgo, the Einstein Telescope, Cosmic Explorer, DECIGO, and LISA.

79 ASTRONOMY AND ASTROPHYSICS↗

Using a Genetic Algorithm to Model Broadband Regional Waveforms for Crustal Structure in the Western United States

In this study, we analyze regional seismograms to obtain the crustal structure in the eastern Great Basin and western Colorado plateau. Adopting a for- ward-modeling approach, we develop a genetic algorithm (GA) based parameter search technique to constrain the one-dimensional crustal structure in these regions. The data are broadband three-component seismograms recorded at the 1994-95 IRIS PASSCAL Colorado Plateau to Great Basin experiment (CPGB) stations and supplemented by data from U.S. National Seismic Network (USNSN) stations in Utah and Nevada. We use the southwestern Wyoming mine collapse event (M(sub b) = 5.2) that occurred on 3 February 1995 as the seismic source. We model the regional seismograms using a four-layer crustal model with constant layer parameters. Timing of teleseismic receiver functions at CPGB stations are added as an additional constraint in the modeling. GA allows us to efficiently search the model space. A carefully chosen fitness function and a windowing scheme are added to the algorithm to prevent search stagnation. The technique is tested with synthetic data, both with and without random Gaussian noise added to it. Several separate model searches are carried out to estimate the variability of the model parameters. The average Colorado plateau crustal structure is characterized by a 40-km-thick crust with velocity increases at depths of about 10 and 25 km and a fast lower crust while the Great Basin has approximately 35- km-thick crust and a 2.9-km-thick sedimentary layer.

Bhattacharyya, Joydeep↗

The Juno mission as a probe of long-range new physics

Orbits of celestial objects, especially the geocentric and heliocentric ones, have been well explored to constrain new long-range forces beyond the Standard Model (SM), often referred to as fifth forces. In this paper, for the first time, we apply the motion of a spacecraft around Jupiter to probe fifth forces that don’t violate the equivalence principle. The spacecraft is the Juno orbiter, and ten of its early orbits already allow a precise determination of the Jovian gravitational field. We use the shift in the precession angle as a proxy to test non-gravitational interactions between Juno and Jupiter. Requiring that the contribution from the fifth force does not exceed the uncertainty of the precession shift inferred from data, we find that a new parameter space with the mass of the fifth-force mediator around 10 −14 eV is excluded at 95% C.L.

new light particles↗

Renormalized classical theory of quantum magnets

Here, we derive a renormalized classical spin (RCS) theory for 𝑆 >1/2 quantum magnets by constraining a generalized classical theory that includes all multipolar fluctuations to a reduced CP 1 phase space of dipolar SU(2) coherent states. When the spin Hamiltonian $\hat{ℋ}$(𝑆) is linear in the spin operators $\hat{𝑺}$ 𝑗 for each lattice site 𝑗, the RCS Hamiltonian $\tilde{ℋ}$ cl coincides with the usual classical model ℋ cl = lim 𝑆→∞⁡ $\hat{ℋ}$(𝑆). In the presence of nonlinear terms, however, the RCS theory is more accurate than ℋ cl . For the many materials modeled by spin Hamiltonians with (nonlinear) single-ion anisotropy terms, the use of the RCS theory is essential to accurately model phase diagrams and to extract the correct Hamiltonian parameters from neutron-scattering data.

magnetic anisotropy↗

Origin and propagation of galactic cosmic rays

The study of systematic trends in elemental abundances is important for unfolding the nuclear and/or atomic effects that should govern the shaping of source abundances and in constraining the parameters of cosmic ray acceleration models. In principle, much can be learned about the large-scale distributions of cosmic rays in the galaxy from all-sky gamma ray surveys such as COS-B and SAS-2. Because of the uncertainties in the matter distribution which come from the inability to measure the abundance of molecular hydrogen, the results are somewhat controversial. The leaky-box model accounts for a surprising amount of the data on heavy nuclei. However, a growing body of data indicates that the simple picture may have to be abandoned in favor of more complex models which contain additional parameters. Future experiments on the Spacelab and space station will hopefully be made of the spectra of individual nuclei at high energy. Antiprotons must be studied in the background free environment above the atmosphere with much higher reliability and presion to obtain spectral information.

Cesarsky, Catherine J.↗

Constraints on Cosmic Rays, Magnetic Fields, and Dark Matter from Gamma-ray Observations of the Coma Cluster of Galaxies with VERITAS and FERMI

Observations of radio halos and relics in galaxy clusters indicate efficient electron acceleration. Protons should likewise be accelerated and, on account of weak energy losses, can accumulate, suggesting that clusters may also be sources of very high energy (VHE; E greater than100 GeV) gamma-ray emission. We report here on VHE gamma-ray observations of the Coma galaxy cluster with the VERITAS array of imaging Cerenkov telescopes, with complementing Fermi Large Area Telescope observations at GeV energies. No significant gamma-ray emission from the Coma Cluster was detected. Integral flux upper limits at the 99 confidence level were measured to be on the order of (2-5) x 10(sup -8) photons m(sup -2) s(sup -1) (VERITAS,greater than 220 GeV) and approximately 2 x 10(sup -6) photons m(sup -2) s(sup -1) (Fermi, 1-3 GeV), respectively. We use the gamma-ray upper limits to constrain cosmic rays (CRs) and magnetic fields in Coma. Using an analytical approach, the CR-to-thermal pressure ratio is constrained to be less than 16% from VERITAS data and less than 1.7% from Fermi data (averaged within the virial radius). These upper limits are starting to constrain the CR physics in self-consistent cosmological cluster simulations and cap the maximum CR acceleration efficiency at structure formation shocks to be 50. Alternatively, this may argue for non-negligible CR transport processes such as CR streaming and diffusion into the outer cluster regions. Assuming that the radio-emitting electrons of the Coma halo result from hadronic CR interactions, the observations imply a lower limit on the central magnetic field in Coma of approximately (2-5.5)microG, depending on the radial magnetic field profile and on the gamma-ray spectral index. Since these values are below those inferred by Faraday rotation measurements in Coma (for most of the parameter space), this renders the hadronic model a very plausible explanation of the Coma radio halo. Finally, since galaxy clusters are dark matter (DM) dominated, the VERITAS upper limits have been used to place constraints on the thermally averaged product of the total self-annihilation cross section and the relative velocity of the DM particles, (sigma upsilon)

cosmic↗

Clarity through the neutrino fog: constraining new forces in dark matter detectors

The PandaX-4T and XENONnT experiments present indications of Coherent Elastic Neutrino Nucleus Scattering (CEνNS) from 8 B solar neutrinos at 2.6σ and 2.7σ, respectively. This constitutes the first observation of the neutrino “floor” or “fog”, an irreducible background that future dark matter searches in terrestrial detectors will have to contend with. Here, we first discuss the contributions from neutrino–electron scattering and from the Migdal effect in the region of interest of these experiments, and we argue that they are non-negligible. Second, we make use of the recent PandaX-4T and XENONnT data to derive novel constraints on light scalar and vector mediators coupling to neutrinos and quarks. We demonstrate that these experiments already provide world-leading laboratory constraints on new light mediators in some regions of parameter space.

neutrino interactions↗

Uncertainty quantification for equations of state: copper as an example

Equations of state are essential for providing a fundamental description of materials properties in thermodynamic equilibrium and are used to provide closure relations for hydrodynamics simulations. Generally, equations of state rely on simple physics-based parameterized materials models to inform on the free energy of a material through out a given thermodynamic state space. Historically the parameters of these models have been tuned by hand to fit various experimental data. However, modern optimization and uncertainty quantification techniques allow us to quickly test thousands of parameter combinations and obtain meaningful uncertainty estimates on the parameters, opening opportunities for assessing systematic uncertainties in experiments, assessing model adequacy, and more. In this report, we use Bayesian inference to fit the solid (fcc) equation of state of copper. We focus on fitting five different experimental datasets, including the isobaric density, isobaric heat capacity, room temperature isotherm, principal isentrope, and principal Hugoniot. We fit all five data types simultaneously, and then explore the extent to which combinations of 2 subsets of the 5 datasets can constrain the EOS parameters, as compared to the fit to all 5. This information is useful for investigating the extent to which different datasets can con strain EOS models and thereby help guide experimental investigations in order to best constrain the EOS. We also discuss ways that the methodologies can be used to investigate systematic discrepancies between experiments, as well as how the methods can be used to assess model uncertainty. The framework we develop is general, in that it can be used with a variety of optimization or uncertainty quantification techniques and with a variety of data sources, including both experimental and ab-inito data.

97 MATHEMATICS AND COMPUTING↗

ACCESS: Confirmation of No Potassium in the Atmosphere of WASP-31b

We present a new optical (400–950 nm) transmission spectrum of the hot Jupiter WASP-31b (M = 0.48 M {sub J}; R = 1.54 R {sub J}; P = 3.41 days), obtained by combining four transit observations. These transits were observed with IMACS on the Magellan Baade Telescope at Las Campanas Observatory as part of the ACCESS project. We investigate the presence of clouds/hazes in the upper atmosphere of this planet, as well as the contribution of stellar activity on the observed features. In addition, we search for absorption features of the alkali elements Na i and K i, with particular focus on K i, for which there have been two previously published disagreeing results. Observations with Hubble Space Telescope (HST)/STIS detected K i, whereas ground-based low- and high-resolution observations did not. We use equilibrium and nonequilibrium chemistry retrievals to explore the planetary and stellar parameter space of the system with our optical data combined with existing near-IR observations. Our best-fit model is that with a scattering slope consistent with a Rayleigh slope (α=5.3{sub −3.1}{sup +2.9}), high-altitude clouds at a log cloud top pressure of −3.6{sub −2.1}{sup +2.7} bars, and possible muted H{sub 2}O features. We find that our observations support other ground-based claims of no K i. Clouds are likely why signals like H{sub 2}O are extremely muted and Na or K cannot be detected. We then juxtapose our Magellan/IMACS transmission spectrum with existing VLT/FORS2, HST/WFC3, HST/STIS, and Spitzer observations to further constrain the optical-to-infrared atmospheric features of the planet. We find that a steeper scattering slope (α = 8.3 ± 1.5) is anchored by STIS wavelengths blueward of 400 nm and only the original STIS observations show significant potassium signal.

79 ASTRONOMY AND ASTROPHYSICS↗

On the compensation of geoid anomalies due to subducting slabs

Candidate models of the forces which oppose the sinking of slabs are all constrained to produce results consistent with the following observation: relative geoid highs, which one assumes are due to slabs, characteristically occur over subduction zones. Certain models which are otherwise plausible, such as those based on a Newtonian half-space mantle, yield geoid lows instead of highs. This study has extended a published model of viscous corner flow in subduction zones in order to demonstrate that it can, in certain cases, produce the requisite geoid highs. Specifically the relative geoid highs are produced if mantle flow is distinctly non-Newtonian (stress exponent n 2). Results in the form of deflection on vertical (or geoid slope) profiles are computed for typical values of the slab parameters; they are compared with a representative profile of geoid slopes derived from Seasat altimeter data in order to show qualitative similarities. It is concluded that the effect of non-Newtonian flow as opposed to Newtonian, is to spread out the induced surface deformation, thereby stretching out the regional compensation to wavelengths, (transverse to the trench) of several thousand kilometers.

Mcadoo, D. C.↗

Pickup Ion Phase Space Distributions at Titan in a Three Dimensional Exosphere

The composition and structure of neutral exospheres imbedded in moving plasmas can be determined by measurements of the velocity distributions of their pickup ion progeny. In turn, the velocity distributions are dependent on the spatial structure of the neutral source gases. Since Titan's neutral exosphere extends into the Saturn's magnetosphere (or solar wind) and well above its ionopause, it serves as a good place to analyze such characteristics. They are analyzed using pickup ion measurements made by the Cassini Plasma Spectrometer (CAPS) at Titan [e.g., Hartle et al., 2006] and an ion kinetic model. An early version of the model [Hartle and Sittler, 2007] is an expression describing the phase space density of pickup ions, which is derived from the Vlasov equation with an ion source that explicitly accounts for the velocity and spatial variation of the exosphere source gases. The current version used here includes exosphere source gases in three dimensions. A fundamental parameter of the phase space densities is the ratio of the gyroradius to the neutral scale height alpha, = r(sub g)/H. Titan's exosphere structure yields pickup ions whose phase space distributions are beam-like when alpha >> 1 and fluid-like when alpha << 1. Downstream from the source peak, the light pickup ions, with alpha << 1, are easily observed because their phase space densities are almost uniform over the orbit phases. On the other hand, the phase space distributions of the heavier ions, with alpha >> 1, peak over narrow velocity and spatial ranges. This beam-like nature makes it considerably more difficult to observe heavy ions because their downstream positions and viewing directions are narrowly constrained. Examples of these extremes will be discussed.

Hartle, Richard↗