Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model fitting”

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,333 records · Page 74

THE DETECTION RATE OF EARLY UV EMISSION FROM SUPERNOVAE: A DEDICATED GALEX/PTF SURVEY AND CALIBRATED THEORETICAL ESTIMATES

The radius and surface composition of an exploding massive star, as well as the explosion energy per unit mass, can be measured using early UV observations of core-collapse supernovae (SNe). We present the first results from a simultaneous GALEX/PTF search for early ultraviolet (UV) emission from SNe. Six SNe II and one Type II superluminous SN (SLSN-II) are clearly detected in the GALEX near-UV (NUV) data. We compare our detection rate with theoretical estimates based on early, shock-cooling UV light curves calculated from models that fit existing Swift and GALEX observations well, combined with volumetric SN rates. We find that our observations are in good agreement with calculated rates assuming that red supergiants (RSGs) explode with fiducial radii of 500 R{sub ⊙}, explosion energies of 10{sup 51} erg, and ejecta masses of 10 M{sub ⊙}. Exploding blue supergiants and Wolf–Rayet stars are poorly constrained. We describe how such observations can be used to derive the progenitor radius, surface composition, and explosion energy per unit mass of such SN events, and we demonstrate why UV observations are critical for such measurements. We use the fiducial RSG parameters to estimate the detection rate of SNe during the shock-cooling phase (<1 day after explosion) for several ground-based surveys (PTF, ZTF, and LSST). We show that the proposed wide-field UV explorer ULTRASAT mission is expected to find >85 SNe per year (∼0.5 SN per deg{sup 2}), independent of host galaxy extinction, down to an NUV detection limit of 21.5 mag AB. Our pilot GALEX/PTF project thus convincingly demonstrates that a dedicated, systematic SN survey at the NUV band is a compelling method to study how massive stars end their life.

79 ASTRONOMY AND ASTROPHYSICS↗

Accretion flow properties of MAXI J1543-564 during 2011 outburst from the TCAF solution

We derive accretion flow properties of the transient black hole candidate (BHC) MAXI J1543–564 using the RXTE data. We use the two-component advective flow (TCAF) solution to fit the data of the initial rising phase of outburst (from 2011 May 10 to 15). The 2.5–25 keV spectra are fitted using the TCAF solution fits file as a local additive table model in XSPEC. We extract physical flow parameters such as the two-component (Keplerian disk and sub-Keplerian halo) accretion rates and size and the property of the Compton cloud (post-shock region close to a black hole). Similar to other classical transient BHCs, monotonic evolution of low-frequency quasi-periodic oscillations (QPOs) is observed during the rising phase of the outburst, which is fitted with the propagating oscillatory shock (POS) model, which describes how the Compton cloud properties change from day to day. From the nature of variations of TCAF model fitted physical flow parameters and QPOs, we only found hard-intermediate and soft-intermediate spectral states during this phase of the outburst under study. We also calculated the frequency of the dominating QPOs from the TCAF model fitted shock parameters and found that they roughly match with the observed and POS model fitted values. From our spectro-temporal study of the source with TCAF and POS models, the most probable mass of the BHC is found to be 12.6–14.0 M {sub ⊙}, or 13{sub −0.4}{sup +1.0} M{sub ⊙}.

79 ASTRONOMY AND ASTROPHYSICS↗

The Next Generation Virgo Cluster Survey. V. modeling the dynamics of M87 with the made-to-measure method

We study the dynamics of the giant elliptical galaxy M87 from the central to the outermost regions with the made-to-measure (M2M) method. We use a new catalog of 922 globular cluster line-of-sight velocities extending to a projected radius of 180 kpc (equivalent to 25 M87 effective radii), and SAURON integral field unit data within the central 2.4 kpc. There are 263 globular clusters, mainly located beyond 40 kpc, newly observed by the Next Generation Virgo Survey. For the M2M modeling, the gravitational potential is taken as a combination of a luminous matter potential with a constant stellar mass-to-light ratio and a dark matter potential modeled as a logarithmic potential. Our best-fit dynamical model returns a stellar mass-to-light ratio in the I band of M/L{sub I} = 6.0 ± 0.3 M{sub ⊙} L{sub ⊙}{sup −1} with a dark matter potential scale velocity of 591 ± 50 km s{sup –1} and scale radius of 42 ± 10 kpc. We determine the total mass of M87 within 180 kpc to be (1.5 ± 0.2) × 10{sup 13} M {sub ☉}. The mass within 40 kpc is smaller than previous estimates determined using globular cluster kinematics that did not extend beyond ∼45 kpc. With our new globular cluster velocities at much larger radii, we see that globular clusters around 40 kpc show an anomalously large velocity dispersion which affected previous results. The mass we derive is in good agreement with that inferred from ROSAT X-ray observation out to 180 kpc. Within 30 kpc our mass is also consistent with that inferred from Chandra and XMM-Newton X-ray observations, while within 120 kpc it is about 20% smaller. The model velocity dispersion anisotropy β parameter for the globular clusters in M87 is small, varying from –0.2 at the center to 0.2 at ∼40 kpc, and gradually decreasing to zero at ∼120 kpc.

79 ASTRONOMY AND ASTROPHYSICS↗

Measurements of the time-dependent cosmic-ray Sun shadow with seven years of IceCube data: Comparison with the Solar cycle and magnetic field models

Observations of the time-dependent cosmic-ray Sun shadow have been proven as a valuable diagnostic for the assessment of solar magnetic field models. Here, seven years of IceCube data are compared to solar activity and solar magnetic field models. A quantitative comparison of solar magnetic field models with IceCube data on the event rate level is performed for the first time. Additionally, a first energy-dependent analysis is presented and compared to recent predictions. We use seven years of IceCube data for the moon and the Sun and compare them to simulations on data rate level. The simulations are performed for the geometrical shadow hypothesis for the moon and the Sun and for a cosmic-ray propagation model governed by the solar magnetic field for the case of the Sun. We find that a linearly decreasing relationship between Sun shadow strength and solar activity is preferred over a constant relationship at the 6.4σ level. We test two commonly used models of the coronal magnetic field, both combined with a Parker spiral, by modeling cosmic-ray propagation in the solar magnetic field. Both models predict a weakening of the shadow in times of high solar activity as it is also visible in the data. We find tensions with the data on the order of 3σ for both models, assuming only statistical uncertainties. The magnetic field model CSSS fits the data slightly better than the PFSS model. This is generally consistent with what is found previously by the Tibet AS-γ Experiment; a deviation of the data from the two models is, however, not significant at this point. Regarding the energy dependence of the Sun shadow, we find indications that the shadowing effect increases with energy during times of high solar activity, in agreement with theoretical predictions.

79 ASTRONOMY AND ASTROPHYSICS↗

SU-F-I-63: Relaxation Times of Lipid Resonances in NAFLD Animal Model Using Enhanced Curve Fitting

Purpose: The objective of this study is to evaluate the relaxation time of methylene resonance in comparison with other lipid resonances. Methods: The examinations were performed on a 3.0T MRI scanner using a four-channel animal coil. Eight more Sprague-Dawley rats in the same baseline weight range were housed with ad libitum access to water and a high-fat (HF) diet (60% fat, 20% protein, and 20% carbohydrate). In order to avoid large blood vessels, a voxel (0.8×0.8×0.8 cm{sup 3}) was placed in a homogeneous area of the liver parenchyma during free breathing. Lipid relaxations in NC and HF diet rats were estimated at a fixed repetition time (TR) of 6000 msec, and multi echo time (TEs) of 40–220 msec. All spectra for data measurement were processed using the Advanced Method for Accurate, Robust, and Efficient Spectral (AMARES) fitting algorithm of the Java-based Magnetic Resonance User Interface (jMRUI) package. Results: The mean T2 relaxation time of the methylene resonance in normal-chow diet was 37.1 msec (M{sub 0}, 2.9±0.5), with a standard deviation of 4.3 msec. The mean T2 relaxation time of the methylene resonance was 31.4 msec (M{sub 0}, 3.7±0.3), with a standard deviation of 1.8 msec. The T2 relaxation times of methylene protons were higher in normal-chow diet rats than in HF rats (p<0.05), and the extrapolated M{sub 0} values were higher in HF rats than in NC rats (p<0.005). The excellent linear fit with R{sup 2}>0.9971 and R{sup 2}>0.9987 indicates T2 relaxation decay curves with mono-exponential function. Conclusion: In in vivo, a sufficient spectral resolution and a sufficiently high signal-to-noise ratio (SNR) can be achieved, so that the data measured over short TE values can be extrapolated back to TE = 0 to produce better estimates of the relative weights of the spectral components. In the short term, treating the effective decay rate as exponential is an adequate approximation.

60 APPLIED LIFE SCIENCES↗

Modeling the effects of vorinostat in vivo reveals both transient and delayed HIV transcriptional activation and minimal killing of latently infected cells

Recent efforts to cure human immunodeficiency virus type-1 (HIV-1) infection have focused on developing latency reversing agents as a first step to eradicate the latent reservoir. The histone deacetylase inhibitor, vorinostat, has been shown to activate HIV RNA transcription in CD4+ T-cells and alter host cell gene transcription in HIV-infected individuals on antiretroviral therapy. In order to understand how latently infected cells respond dynamically to vorinostat treatment and determine the impact of vorinostat on reservoir size in vivo , we have constructed viral dynamic models of latency that incorporate vorinostat treatment. We fitted these models to data collected from a recent clinical trial in which vorinostat was administered daily for 14 days to HIV-infected individuals on suppressive ART. The results show that HIV transcription is increased transiently during the first few hours or days of treatment and that there is a delay before a sustained increase of HIV transcription, whose duration varies among study participants and may depend on the long term impact of vorinostat on host gene expression. Parameter estimation suggests that in latently infected cells, HIV transcription induced by vorinostat occurs at lower levels than in productively infected cells. Lastly, the estimated loss rate of transcriptionally induced cells remains close to baseline in most study participants, suggesting vorinostat treatment does not induce latently infected cell killing and thus reduce the latent reservoir in vivo .

59 BASIC BIOLOGICAL SCIENCES↗

Technical Report on Waveform Fit Metrics for Global Models

The new WAVEFORMS Initiative in the Ground-based Nuclear Detonation Detection (GNDD) program includes an increased emphasis on the development of Earth models and methods to predict entire seismic and acoustic waveforms more accurately. In general, this increased emphasis is predicated on the need to better characterize seismic events and provide improved model-based discrimination between event types including earthquakes and explosions. More specifically, while current moment tensor inversion methods tend to work well for larger events (M>~4) using tuned 1-D Earth models, the development of state-of-the-art 3-D models and methods is required for the prediction of shorter period waves over large areas for discrimination of smaller events. There is no standard metric for model-based waveform prediction accuracy used in the waveform modeling/inversion community. However, there are several popular waveform misfit definitions; and minimizing the corresponding objective functions is the goal of waveform inversion. Some example misfit definitions employed for adjoint waveform tomography include measures of simple travel time differences (e.g. Tape et al., 2010), cross-correlation travel time differences (e.g. Luo and Schuster, 1991), multi-taper frequency dependent methods (e.g. Lei et al., 2020), time-frequency phase misfit functions (e.g. Fichtner 2010; Rodgers et al., 2022), normalized cross-correlation methods (e.g. Tao et al., 2018), and others. In some cases, these misfit definitions also involve complicated weighting schemes and summations over multiple frequency bands making it difficult to duplicate the misfit measurement with alternative models and datasets. Although each of the misfit definitions mentioned above are useful for developing waveform models, the actual misfit values are not usually meaningful outside of a given project, model, and/or dataset. Therefore, it is difficult to understand and communicate model performance for predicting waveforms and comparing to other models and/or new model iterations with a different dataset. Therefore, there is a need for a generalized method for evaluating overall model performance that is independent from the specific misfit chosen to develop the waveform models that is also intuitive and meaningful. In this report, we describe a new metric we refer to as ‘Percent of Correlated Signal’. The following sections describe and demonstrate the metric with a case study event and a more rigorous test using a random selection of globally distributed events. While the focus here is on global tomography models, the metric is meant to applicable to regional ‘wiggle-for-wiggle’ waveform models/studies as well.

58 GEOSCIENCES↗

The unusually high halo concentration of the fossil group NGC 6482: evidence for weak adiabatic contraction

Although fossil galaxy systems are thought to be very old, and thus should possess above-average halo concentrations, typically fossils have unexceptional concentrations of their masses. We revisit the massive isolated elliptical galaxy/fossil group NGC 6482 for which previous X-ray studies of a modest Chandra observation obtained a very uncertain, but also possibly very high, halo concentration. We present new measurements of the hot gas surface brightness, temperature, and iron abundance using the modest Chandra observation and a previously unpublished Suzaku observation, the latter of which allows the measurements of the gas properties to be extended out to ∼r{sub 2500}. By constructing hydrostatic equilibrium models of the gas with separate components for the gas, the stellar mass of the large central galaxy (BCG), and the dark matter (DM), we measure c{sub 200}=32.2±7.1 and M{sub 200}=(4.5±0.6)×10{sup 12} M{sub ⊙} using a Navarro-Frenk-White (NFW) DM profile. For a halo of this mass, the measured concentration c{sub 200} exceeds the mean value (7.1) expected for relaxed ΛCDM halos by 3.5σ in terms of the observational error, and by 6σ considering the intrinsic scatter in the ΛCDM c–M relation, which situates NGC 6482 as the most extreme outlier known for a fossil system. We explored several variants of adiabatic contraction (AC) models and, while the AC models provide fits of the same quality as the uncontracted models, they do have the following advantages: (1) lower c{sub 200} that is less of an outlier in the ΛCDM c–M relation, and (2) baryon fractions (f{sub b,200}) that agree better with the mean cosmic value. While the standard AC prescriptions yield a BCG stellar mass that is uncomfortably low compared to results from stellar population synthesis (SPS) models, a weaker AC variant that artificially shuts off cooling and star formation at z = 2 yields the same stellar mass as the uncontracted models. These are probably the reasons why our X-ray analysis prefers this weaker AC variant applied to either an NFW or Einasto DM halo. Finally, the stellar mass we measure for the BCG from the hydrostatic analysis strongly favors results from SPS models with a Chabrier or Kroupa initial mass function (IMF) over a Salpeter IMF.

79 ASTRONOMY AND ASTROPHYSICS↗

η Carinae's Dusty Homunculus Nebula from Near-infrared to Submillimeter Wavelengths: Mass, Composition, and Evidence for Fading Opacity

Infrared observations of the dusty, massive Homunculus Nebula around the luminous blue variable η Carinae are crucial to characterize the mass-loss history and help constrain the mechanisms leading to the great eruption. We present the 2.4–670 μm spectral energy distribution, constructed from legacy Infrared Space Observatory observations and new spectroscopy obtained with the Herschel Space Observatory. Using radiative transfer modeling, we find that the two best-fit dust models yield compositions that are consistent with CNO-processed material, with iron, pyroxene and other metal-rich silicates, corundum, and magnesium-iron sulfide in common. Spherical corundum grains are supported by the good match to a narrow 20.2 μm feature. Our preferred model contains nitrides AlN and Si{sub 3}N{sub 4} in low abundances. Dust masses range from 0.25 to 0.44 M{sub ⊙}, but M{sub tot}⩾45 M{sub ⊙} in both cases, due to an expected high Fe gas-to-dust ratio. The bulk of dust is within a 5{sup ′′}×7{sup ′′} central region. An additional compact feature is detected at 390 μm. We obtain L{sub IR} = 2.96 × 10{sup 6} L{sub ⊙}, a 25% decline from an average of mid-IR photometric levels observed in 1971–1977. This indicates a reduction in circumstellar extinction in conjunction with an increase in visual brightness, allowing 25%–40% of optical and UV radiation to escape from the central source. We also present an analysis of {sup 12}CO and {sup 13}CO J = 5–4 through 9–8 lines, showing that the abundances are consistent with expectations for CNO-processed material. The [{sup 12}C ii] line is detected in absorption, which we suspect originates in foreground material at very low excitation temperatures.

79 ASTRONOMY AND ASTROPHYSICS↗

A SIRTF interdiciplinary scientist proposal

Segal's chronometric cosmology provides an adequate fit to the radio source counts only for an unrealistic choice of spectral index. Since the typical observed spectral index of 0.75 gives a completely unacceptable X squared = 136 with 24 (or fewer) degrees of freedom, it is concluded that the actual Universe does not fit the chronometric model. Counts of ultraviolet excess quasistellar objects also show a steep N(S) curve that the chronometric cosmology cannot explain. Claims to the contrary by Segal, Loncaric, and Segal (1980) and Segal and Nicoll (1986) depend on a seemingly innocuous assumption that in fact destroys the power of the N(S) test. Even though the chronometric model gives a better fit that other non-evolving models it must be ruled out along with all non-evolving cosmologies.

Wright, E. L.↗

Are existing magnetospheric models excessively stretched?

The distribution of the average north-south magnetic field component Bz in the vicinity of the neutral sheet has been investigated. This component is crucial for mappings between the nightside polar ionosphere and the equatorial magnetosphere. Data sets consisting of about 0.5 R(E) averages of magnetic field observations by the IMP/HEOS and ISEE spacecraft have been compared to the field predicted by the Tsyganenko (1987, 1989) models T87 and T89. It was found that both T87 and T89 underestimate Bz in the near tail region by as much as a factor of 2. Modified versions of the T87 model, incorporating plasma sheet warping, were obtained by fitting the model parameters via nonlinear least squares to the ISEE data set and yielded Bz values in agreement with the ISEE data. The study reveals an enormous scatter among the observed baseline values of Bz (on a time scale of 10-20 min), as well as intrinsic biases imposed by the mathematical structure of tail models, and these two factors (especially the first one) greatly limit the accuracy of model predictions of tail Bz.

Peredo, Mauricio↗

Parameterized Linear Longitudinal Airship Model

A parameterized linear mathematical model of the longitudinal dynamics of an airship is undergoing development. This model is intended to be used in designing control systems for future airships that would operate in the atmospheres of Earth and remote planets. Heretofore, the development of linearized models of the longitudinal dynamics of airships has been costly in that it has been necessary to perform extensive flight testing and to use system-identification techniques to construct models that fit the flight-test data. The present model is a generic one that can be relatively easily specialized to approximate the dynamics of specific airships at specific operating points, without need for further system identification, and with significantly less flight testing. The approach taken in the present development is to merge the linearized dynamical equations of an airship with techniques for estimation of aircraft stability derivatives, and to thereby make it possible to construct a linearized dynamical model of the longitudinal dynamics of a specific airship from geometric and aerodynamic data pertaining to that airship. (It is also planned to develop a model of the lateral dynamics by use of the same methods.) All of the aerodynamic data needed to construct the model of a specific airship can be obtained from wind-tunnel testing and computational fluid dynamics

Kulczycki, Eric↗

The Asymmetric Inner Disk of the Herbig Ae Star HD 163296 in the Eyes of VLTI/MATISSE: Evidence for a Vortex?

Context.A complex environment exists in the inner few astronomical units of planet-forming disks. High-angular-resolution observa-tions play a key role in our understanding of the disk structure and the dynamical processes at work.Aims.In this study we aim to characterize the mid-infrared brightness distribution of the inner disk of the young intermediate-massstar HD 163296 from early VLTI/MATISSE observations taken in theL- andN-bands. We put special emphasis on the detection ofpotential disk asymmetries.Methods.We use simple geometric models to fit the interferometric visibilities and closure phases. Our models include a smoothedring, a flat disk with an inner cavity, and a 2D Gaussian. The models can account for disk inclination and for azimuthal asymmetriesas well. We also perform numerical hydrodynamical simulations of the inner edge of the disk.Results.Our modeling reveals a significant brightness asymmetry in theL-band disk emission. The brightness maximum of the asym-metry is located at the NW part of the disk image, nearly at the position angle of the semimajor axis. The surface brightness ratio inthe azimuthal variation is3.5±0.2. Comparing our result on the location of the asymmetry with other interferometric measurements,we confirm that the morphology of ther<0.3au disk region is time-variable. We propose that this asymmetric structure, located in ornear the inner rim of the dusty disk, orbits the star. To find the physical origin of the asymmetry, we tested a hypothesis where a vortexis created by Rossby wave instability, and we find that a unique large-scale vortex may be compatible with our data. The half-lightradius of theL-band-emitting region is0.33±0.01au, the inclination is52◦+5◦−7◦, and the position angle is143◦±3◦. Our models predictthat a non-negligible fraction of theL-band disk emission originates inside the dust sublimation radius forμm-sized grains. Refractorygrains or large (&10μm-sized) grains could be the origin of this emission.N-band observations may also support a lack of smallsilicate grains in the innermost disk (r.0.6au), in agreement with our findings fromL-band data.

J Varga↗

On-board state of health monitoring of batteries using incremental capacity analysis

A method is provided for determining capacity of a battery. The method includes: defining a model for a battery, where the model relates terminal voltage of the battery to charged capacity of the battery; taking a plurality of voltage measures of the battery, where the voltage measures are taken through a range of states of charge and the range excludes the battery being fully charged and fully discharged; determining the parameters of the model by fitting the plurality of voltage measures to the model; determining an incremental capacity curve for the battery by taking derivative of the model; and quantifying a peak of the incremental capacity curve to thereby determine a capacity of the battery.

Sun, Jing↗