Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data 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 919 records · Page 51

TOPEX/POSEIDON tides estimated using a global inverse model

Altimetric data from the TOPEX/POSEIDON mission will be used for studies of global ocean circulation and marine geophysics. However, it is first necessary to remove the ocean tides, which are aliased in the raw data. The tides are constrained by the two distinct types of information: the hydrodynamic equations which the tidal fields of elevations and velocities must satisfy, and direct observational data from tide gauges and satellite altimetry. Here we develop and apply a generalized inverse method, which allows us to combine rationally all of this information into global tidal fields best fitting both the data and the dynamics, in a least squares sense. The resulting inverse solution is a sum of the direct solution to the astronomically forced Laplace tidal equations and a linear combination of the representers for the data functionals. The representer functions (one for each datum) are determined by the dynamical equations, and by our prior estimates of the statistics or errors in these equations. Our major task is a direct numerical calculation of these representers. This task is computationally intensive, but well suited to massively parallel processing. By calculating the representers we reduce the full (infinite dimensional) problem to a relatively low-dimensional problem at the outset, allowing full control over the conditioning and hence the stability of the inverse solution. With the representers calculated we can easily update our model as additional TOPEX/POSEIDON data become available. As an initial illustration we invert harmonic constants from a set of 80 open-ocean tide gauges. We then present a practical scheme for direct inversion of TOPEX/POSEIDON crossover data. We apply this method to 38 cycles of geophysical data records (GDR) data, computing preliminary global estimates of the four principal tidal constituents, M(sub 2), S(sub 2), K(sub 1) and O(sub 1). The inverse solution yields tidal fields which are simultaneously smoother, and in better agreement with altimetric and ground truth data, than previously proposed tidal models. Relative to the 'default' tidal corrections provided with the TOPEX/POSEIDON GDR, the inverse solution reduces crossover difference variances significantly (approximately 20-30%), even though only a small number of free parameters (approximately equal to 1000) are actually fit to the crossover data.

Egbert, Gary D.↗

The Effects of Assumed Source Depth and Shear-Wave Velocity on Moment Tensors Estimated for Small, Contained Chemical Explosions in Granite

ABSTRACT The Source Phenomenology Experiment (SPE-Arizona) included of a series of chemical explosions detonated within a copper mine in Arizona. This study focuses on ground motions from detonations in the copper mine, which are analyzed to assess the uniqueness of the resulting source representation when the source region propagation characteristics have a range of possible models. P-wave velocities are well constrained by refraction data with less constraint of the S-wave velocities. The effects of explosion source depth and VS are assessed with Green’s functions for a range of models in which VP is held constant. Propagation models with a Poisson’s value of 0.25 and a source depth 30–60 m most accurately replicate the data. The explosion was detonated at a centroid depth of 30 m, so trade-offs in depth are demonstrated. The compensated linear vector dipole and explosion components of the Green’s functions convolved with a Mueller–Murphy source function are compared. Both produce significant energy in the 2–12 Hz band, due to surface-wave contributions with no clear depth dependencies above 20 Hz. The range of propagation models is used with the observational data to invert for the frequency-domain moment tensor. Fits to the data from these inversions have cross-correlation values of 0.64, demonstrating effectiveness in replicating the observations with the assumed propagation path effects and resulting source function. Inversions produce horizontal dipoles (Mxx and Myy), roughly half the maximum amplitude of Mzz, consistent with a compensated linear vector dipole source, which is frequency dependent. Denny and Johnson, Mueller–Murphy, Walter and Ford, and the revised Mueller–Murphy source models, parameterized for granite, are compared to the moment tensors. Despite a nonisotropic moment tensor source, the revised Mueller–Murphy isotropic source model best replicates the long-period moments, overshoot, and corner frequency.

Geochemistry & Geophysics↗

Application of fiber bridging models to fatigue crack growth in unidirectional titanium matrix composites

Several fiber bridging models were reviewed and applied to study the matrix fatigue crack growth behavior in center notched (0)(sub 8) SCS-6/Ti-15-3 and (0)(sub 4) SCS-6/Ti-6Al-4V laminates. Observations revealed that fatigue damage consisted primarily of matrix cracks and fiber matrix interfacial failure in the (0)(sub 8) SCS-6/Ti-15-3 laminates. Fiber-matrix interface failure included fracture of the brittle reaction zone and cracking between the two carbon rich fiber coatings. Intact fibers in the wake of the matrix cracks reduce the stress intensity factor range. Thus, an applied stress intensity factor range is inappropriate to characterize matrix crack growth behavior. Fiber bridging models were used to determine the matrix stress intensity factor range in titanium metal matrix composites. In these models, the fibers in the wake of the crack are idealized as a closure pressure. An unknown constant frictional shear stress is assumed to act along the debond or slip length of the bridging fibers. The frictional shear stress was used as a curve fitting parameter to available data (crack growth data, crack opening displacement data, and debond length data). Large variations in the frictional shear stress required to fit the experimental data indicate that the fiber bridging models in their present form lack predictive capabilities. However, these models provide an efficient and relatively simple engineering method for conducting parametric studies of the matrix growth behavior based on constituent properties.

Bakuckas, J. G., Jr.↗

Exocentric judgements in real environments and stereoscopic displays

Spatial direction errors during interpretation of perspective images, such as 3D map displays, may originate from misjudgment of the orientation of the viewing direction used to make the display. One source of these errors could be perceptual evidence of the display surface. Two experiments are reported in which the same judgment exocentric task was presented, but the cues to the picture surface were reduced or eliminated by presenting the task as a stereoscopic, virtual image or by a geometrically matched physical model. A theory developed to model exocentric direction errors on perspective displays has been fitted to the data from these two experiments. The parameters estimated from the fit in both experiments indicate that the subjects may be more correctly estimating the viewing direction than in ordinary perspective displays. Consequently, in some real world or stereo viewing conditions, errors in estimating the viewing direction are not likely to dominate exocentric direction errors.

Ellis, Stephen R.↗

Gamma-Ray Burst Afterglows as Probes of Environment and Blastwave Physics II: The Distribution of p and Structure of the Circumburst Medium

We constrain blastwave parameters and the circumburst media of a subsample of BeppoSAX Gamma-Ray Bursts. For this sample we derive the values of the injected electron energy distribution index, p, and the density structure index of the circumburst medium, k, from simultaneous spectral fits to their X-ray, optical and nIR afterglow data. The spectral fits have been done in count space and include the effects of metallicity, and are compared with the previously reported optical and X-ray temporal behaviour. Assuming the fireball model, we can find a mean value of p for the sample as a whole of 2.035. A statistical analysis Of the distribution demonstrates that the p values in this sample are inconsistent with a single universal value for p at the 3sigma level or greater. This approach provides us with a measured distribution of circumburst density structures rather than considering only the cases of k = 0 (homogeneous) and k = 2 (wind-like). We find five GRBs for which k can be well constrained, and in four of these cases the circumburst medium is clearly wind-like. The fifth source has a value of 0 less than or equal to k less than or equal to 1, consistent with a homogeneous circumburst medium.

Starling, R. L. C.↗

Gamma-Ray Burst Afterglows as Probes of Environment and Blast Wave Physics. II. The Distribution of rho and Structure of the Circumburst Medium

We constrain blast wave parameters and the circumburst media ofa subsample of 10 BeppoSAX gamma-ray bursts (GRBs). For this sample we derive the values of the injected electron energy distribution index, p, and the density structure index of the circumburst medium, k, from simultaneous spectral fits to their X-ray, optical, and NIR afterglow data. The spectral fits have been done in count space and include the effects ofmetallicity, and are compared with the previously reported optical and X-ray temporal behavior. Using the blast wave model and some assumptions which include on-axis viewing and standard jet structure, constant blast wave energy, and no evolution of the microphysical parameters, we find a mean value ofp for the sample as a whole of 9.... oa -0.003.0" 2 a_ statistical analysis of the distribution demonstrates that the p-values in this sample are inconsistent with a single universal value forp at the 3 _ level or greater, which has significant implications for particle acceleration models. This approach provides us with a measured distribution ofcircumburst density structures rather than considering only the cases of k ----0 (homogeneous) and k - 2 (windlike). We find five GRBs for which k can be well constrained, and in four of these cases the circumburst medium is clearly windlike. The fifth source has a value of 0 < k < 1, consistent with a homogeneous circumburst medium.

Starling, R. L. C.↗

Episodes of Cross-Polar Transport in the Arctic Troposphere During July 2008 as Seen from Models, Satellite, and Aircraft Observations

During the POLARCAT summer campaign in 2008, two episodes (2 5 July and 7 10 July 2008) occurred where low-pressure systems traveled from Siberia across the Arctic Ocean towards the North Pole. The two cyclones had extensive smoke plumes from Siberian forest fires and anthropogenic sources in East Asia embedded in their associated air masses, creating an excellent opportunity to use satellite and aircraft observations to validate the performance of atmospheric transport models in the Arctic, which is a challenging model domain due to numerical and other complications. Here we compare transport simulations of carbon monoxide (CO) from the Lagrangian transport model FLEXPART and the Eulerian chemical transport model TOMCAT with retrievals of total column CO from the IASI passive infrared sensor onboard the MetOp-A satellite. The main aspect of the comparison is how realistic horizontal and vertical structures are represented in the model simulations. Analysis of CALIPSO lidar curtains and in situ aircraft measurements provide further independent reference points to assess how reliable the model simulations are and what the main limitations are. The horizontal structure of mid-latitude pollution plumes agrees well between the IASI total column CO and the model simulations. However, finer-scale structures are too quickly diffused in the Eulerian model. Applying the IASI averaging kernels to the model data is essential for a meaningful comparison. Using aircraft data as a reference suggests that the satellite data are biased high, while TOMCAT is biased low. FLEXPART fits the aircraft data rather well, but due to added background concentrations the simulation is not independent from observations. The multi-data, multi-model approach allows separating the influences of meteorological fields, model realisation, and grid type on the plume structure. In addition to the very good agreement between simulated and observed total column CO fields, the results also highlight the difficulty to identify a data set that most realistically represents the actual pollution state of the Arctic atmosphere.

Sodemann, H.↗

Enhancements of Bayesian Blocks; Application to Large Light Curve Databases

Bayesian Blocks are optimal piecewise linear representations (step function fits) of light-curves. The simple algorithm implementing this idea, using dynamic programming, has been extended to include more data modes and fitness metrics, multivariate analysis, and data on the circle (Studies in Astronomical Time Series Analysis. VI. Bayesian Block Representations, Scargle, Norris, Jackson and Chiang 2013, ApJ, 764, 167), as well as new results on background subtraction and refinement of the procedure for precise timing of transient events in sparse data. Example demonstrations will include exploratory analysis of the Kepler light curve archive in a search for "star-tickling" signals from extraterrestrial civilizations. (The Cepheid Galactic Internet, Learned, Kudritzki, Pakvasa1, and Zee, 2008, arXiv: 0809.0339; Walkowicz et al., in progress).

Kepler light curve archive↗

Characteristics of DSD Bulk Parameters: Implication to Radar Rain Retrieval

With the use of 213,456 one-minute measured data of droplet-size distribution (DSD) of rain collected during several National Aeronautics and Space Administration (NASA)-sponsored field campaigns, the relationships between rainfall rate R, mass-weighted diameter D(m) and normalized intercept parameter N(w) of the gamma DSD are studied. It is found, based on the simulations of the gamma DSD model, that R, D(m) and N(w) are closely interrelated, and that the ratio of R to N(w) is solely a function of D(m), independent of the shape factor μ of the gamma distribution. Furthermore, the model-produced ratio agrees well with those from the DSD data. When a power-law equation is applied to fit the model data, we have: R=aN(w)D(sup b, sub m) , where a=1.588×10(exp −4) , b=4.706 . Analysis of two-parameter relationships such as R–D(m), N(w)–R and N(w)–D(m) reveals that R and D(m) are moderately correlated while N(w) and D(m) are negatively correlated. N(w) and R, however, are uncorrelated. The gamma DSD model also reveals that variation of R–D(m) relation is caused primarily by N(w). For the application of the Ku- and Ka-band dual-frequency radar for the retrieval of the DSD bulk parameters as well as the specific radar attenuations, the study is carried out to relate the dual-frequency radar reflectivity factors to the DSD and attenuation parameters.

Liang Liao↗

A Misperception in Reliability Growth Modelling

The Duane reliability growth model is n(t)/t = k t^-alpha (1) The reliability growth rate is alpha, the downward slope of n(t)/t versus t. It usually varies from 0.2 to 0.6. k is a constant. Crow used a 56-failure data set to illustrate reliability growth.1 A graphical Duane model fit to this data gives n(t)/t = 0.640 t^-0.283 (2) A problem in using the Duane-Crow reliability growth model is that it assumes that reliability growth continues and the failure rate decreases throughout the test period. It is more usual that reliability growth stops when the failure rated is low enough. Growth testing is often followed by testing with a low constant failure rate due to rare or uncorrectable failure modes. As more and more low constant rate acceptable failures accumulate after the period of reliability growth, the reliability growth time exponent alpha decreases toward zero. This occurs if constant rate failures are treated as occurring during the reliability growth period. It is more accurate to model a period of initial reliability growth followed by testing without repair to more accurately determine the final constant failure rate. This is done in the abcd model. n(t)/t = a t^-b + c from t = 0 to td (3) = c + d after td, where d = a td^-b (4) The term a t^-b describes the continuous reliability growth that continues out to time td and c is the constant uncorrected failure rate. The parameter d represents an additional constant failure rate due to correctable but uncorrected failure modes. After the reliability growth process is terminated, the failure rate n(t)/t = c + d.

Harry W Jones↗

Active region coronal evolution

Scaling relations between coronal base pressure and longitudinal photospheric magnetic field strength are tested for the case of a single active region observed for five solar rotations from Skylab. The evolution of measureable quantities, such as coronal thermal energy content, total longitudinal photospheric magnetic flux, region scale size, and peak energy density, is traced throughout the five rotations observed. The theoretically derived scaling law of Golub et al. (1980) is found to provide an acceptable fit to the data throughout the entire evolutionary history of the region from an age of about 3 days to the fully evolved state in which the mature active region merges into the general large-scale structure of the quiet corona. An alternative scaling law obtained by including the results of Galeev et al. (1981), however, is found to provide a somewhat better fit to the data. The study is seen as providing additional justification for the belief that magnetic field-related heating is the operative mechanism in the solar corona.

Golub, L.↗

ROSAT PSPC observations of two X-ray-faint early-type galaxies: NGC 4365 and NGC 4382

We present the results of ROSAT Positive Sensitive Proportional Counter (PSPC) observations of the two early-type galaxies NGC 4365 and NGC 4382. These galaxies are among those observed with Einstein to have the lowest X-ray to optical flux ratios of early-type galaxies. The PSCP data show that for radii r greater than 50 arcsec the radial distributions of the X-ray surface brightness are consistent with the optical distributions of King (1978). We also find that these galaxies have X-ray spectra significantly different from those observed in X-ray-bright ellipticals, with a relative excess of counts detected in the softest spectral channels. This confirms earlier Einstein results. The characteristics of the ROSAT PSPC do not allow us to discriminate between possible spectral models. If we adopt a two-component thermal model on the grounds of physical plausibility, we find that the spectral data can be fitted with a very soft optically thin component, with kT approximately 0.2 keV, and a hard component with kT greater than (1.0-1.5) keV. The hard component has a luminosity consistent with that expected from the integrated emission of a population of low mass-X-ray binaries in these galaxies; the nature of the very soft component is more speculative. Candidates include the coronal emission of late-type stars, supersoft X-ray sources, RS CVn, and perhaps a hot Interstellar Medium (ISM). Alternatively, the spectal data may be fitted with a 0.6-1 keV bremsstrahlung spectrum (expontential plus Gaunt), and may suggest the presence of a totally new population of X-ray sources.

Fabbiano, G.↗

Analysis of spin frustration in an Fe III 7 cluster using a combination of computational, experimental, and magnetostructural correlation methods

The synthesis, structure, and magnetic properties are reported for [Fe 7 O 3 (O 2 C t Bu) 9 (mda) 3 (H 2 O) 3 ] ( 1 ), where mdaH 2 is N -methyldiethanolamine. 1 was prepared from the reaction of [Fe 3 O(O 2 C t Bu) 6 (H 2 O) 3 ](NO 3 ) with mdaH 2 in a 1:~3 ratio in MeCN. The core of 1 consists of a central octahedral Fe III ion held within a non-planar Fe 6 loop by three μ 3 -O 2- and three μ 2 -RO - arms from the three mda 2- chelates. Variable-temperature dc and ac magnetic susceptibility studies revealed dominant antiferromagnetic coupling, leading to a ground state spin of S = 5 / 2 . The ground state was confirmed by a fit of magnetization data collected in the 0.1–7.0 T and 1.8–10.0 K ranges. The four Fe 2 pairwise exchange parameters ( J 1 - J 4 ) were estimated by independent methods: theoretical calculations using either broken symmetry energy differences (-46.3, -16.2, -3.9, and - 28.1 cm -1 , respectively) or Green’s function approximation methods (-41.4, -14.8, -13.2, and - 24.7 cm -1 ), and a magnetostructural correlation (MSC) previously developed for high nuclearity Fe III /O complexes (-39.5, -13.8, -6.7, and - 23.5 cm -1 ). Additionally, the J 1 - J 4 obtained from the MSC and theoretical methods were used with the program PHI to both simulate χ M T vs T as well as to serve as reasonable input values to fit the experimental data (-41.0, -11.4, -5.0, and - 27.3 cm -1 ). Analysis of the J ij led to identification of the spin frustration effects operative and the resultant spin vector alignments at each Fe III ion, thus allowing for the rationalization of the experimental ground state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fast and flexible analysis of direct dark matter search data with machine learning

We present the results from combining machine learning with the profile likelihood fit procedure, using data from the Large Underground Xenon (LUX) dark matter experiment. This approach demonstrates reduction in computation time by a factor of 30 when compared with the previous approach, without loss of performance on real data. We establish its flexibility to capture non-linear correlations between variables (such as smearing in light and charge signals due to position variation) by achieving equal performance using pulse areas with and without position-corrections applied. Its efficiency and scalability furthermore enables searching for dark matter using additional variables without significant computational burden. We demonstrate this by including a light signal pulse shape variable alongside more traditional inputs such as light and charge signal strengths. Furthermore, this technique can be exploited by future dark matter experiments to make use of additional information, reduce computational resources needed for signal searches and simulations, and make inclusion of physical nuisance parameters in fits tractable.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Transverse momentum dependent feed-down fractions for bottomonium production

We extract transverse-momentum-dependent feed-down fractions for bottomonium production using a data-driven approach. We use data published by the ATLAS, CMS, and LHCb Collaborations for √𝑠 =7 TeV proton-proton collisions. Based on this collected data, we produce fits to the differential cross sections for the production of both 𝑆- and 𝑃-wave bottomonium states. Combining these fits with branching ratios for excited state decays from the Particle Data Group, we compute the feed-down fractions for both the ϒ⁡(1⁢𝑆) and ϒ⁡(2⁢𝑆) as a function of transverse momentum. Our results indicate a strong dependence on transverse momentum, which is consistent with prior extractions of the feed-down fractions. When evaluated at the average momentum of the states, we find that approximately 75% of ϒ⁡(1⁢𝑆) and ϒ⁡(2⁢𝑆) states are produced directly. Our results for the transverse-momentum-dependent feed-down fractions are provided in tabulated form so that they can be used by other research groups.

Astronomy & Astrophysics↗

Finite width effects in nonleptonic D -meson decays

Many analyses of two-body nonleptonic decays of D -mesons rely on flavor S U ( 3 ) symmetry relations and fits of experimental data of decays rates to extract the universal transition amplitudes. Such fits assume that the final state mesons are well-defined asymptotic states of QCD. We develop a technique to take into account the finite width effects of the final state mesons and study their effects on the extracted values of transition amplitudes. Published by the American Physical Society 2025

Kumar, Girish (ORCID:0000000160512495)↗

Observations of noise bands associated with the upper hybrid resonance by the Imp 6 radio astronomy experiment.

The intense noise bands occurring near the upper hybrid resonance frequency have been observed with the Imp 6 GSFC radio astronomy experiment in the plasmasphere. The identification of the upper hybrid resonance provides an accurate measure of the local electron density and allows the observed noise data to be fit to the scale of characteristic frequencies in the plasma. The data are consistent with earlier theoretical interpretations in which noise is generated between the upper hybrid and plasma frequencies and propagates to region 4 of the CMA diagram, where it is reflected at the L = 0 plasma cutoff.

Mosier, S. R.↗

The correlation of Skylab L-band brightness temperatures with antecedent precipitation

The S194 L-band radiometer flown on the Skylab mission measured terrestrial radiation at the microwave wavelength of 21.4 cm. The terrain emissivity at this wavelength is strongly dependent on the soil moisture content, which can be inferred from antecedent precipitation. For the Skylab data acquisition pass from the Oklahoma panhandle to southeastern Texas on 11 June 1973, the S194 brightness temperatures are highly correlated with antecedent precipitation from the preceding eleven day period, but very little correlation was apparent for the preceding five day period. The correlation coefficient between the averaged antecedent precipitation index values and the corresponding S194 brightness temperatures between 230 K and 270 K, the region of apparent response to soil moisture in the data, was -0.97. The equation of the linear least squares line fitted to the data was: API (cm) = 31.99 -0.114 TB (K), where API is the antecedent precipitation index and TB is the S194 brightness temperature.

Mcfarland, M. J.↗