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 1,207 records · Page 67

Passage of Bolides Through the Atmosphere

Different fragmentation models are applied to a number of events, including the entry of TC3 2008 asteroid in order to reproduce existing observational data. Keywords meteoroid entry fragmentation modeling 1 Introduction Fragmentation is a very important phenomenon which occurs during the meteoroid entry into the atmosphere and adds more drastic effects than mere deceleration and ablation. Modeling of bolide fragmentation (100 106 kg in mass) may be divided into several approaches. Detail fitting of observational data (deceleration and/or light curves) allows the determination of some meteoroid parameters (ablation and shape-density coefficients, fragmentation points, amount of mass loss) (Ceplecha et al. 1993; Ceplecha and ReVelle 2005). Observational data with high accuracy are needed for the gross-fragmentation model (Ceplecha et al. 1993), which is used for the analysis of European and Desert bolide networks data. Hydrodynamical models, which describe the entry of the meteoroid including evolution of its material, are applied mainly for large bodies (>106 kg) (Boslough et al. 1994; Svetsov et al. 1995; Shuvalov and Artemieva 2002, and others). Numerous papers were devoted to the application of standard equations for large meteoroid entry in the attempts to reproduce dynamics and/or radiation for different bolides and to predict meteorite falls. These modeling efforts are often supplemented by different fragmentation models (Baldwin and Sheaffer, 1971; Borovi.ka et al. 1998; Artemieva and Shuvalov, 2001; Bland and Artemieva, 2006, and others). The fragmentation may occur in different ways. For example, few large fragments are formed. These pieces initially interact through their shock waves and then continue their flight independently. The progressive fragmentation model suggests that meteoroids are disrupted into fragments, which continue their flight as independent bodies and may be disrupted further. Similar models were suggested in numerous papers, beginning with Levin (1956) and initial interaction of fragments started to be taken into account after the paper by Passey and Melosh (1980). The progressive fragmentation model with lateral spreading of formed fragments is widely used (Artemieva and Shuvalov, 1996; Nemtchinov and Popova, 1997; Borovi.ka et al. 1998; Bland and Artemieva, 2006).

Popova, O.↗

PolyFit: A C++ code for Polynomial Curve Fit with Calculation of Error Bars

In radiobiology, many dose-response results are modeled using the so-called linear-quadratic(LQ)model, which means that results are modeled as a function of dose Das 𝑅(𝐷)=𝛽0+𝛽1𝐷+𝛽2𝐷2. The coefficients 𝛽0, 𝛽1and 𝛽2are obtained from fitting a series of data points (𝑥𝑖,𝑦𝑖), which is usually done using a least-square method. The LQ and more generally the polynomial fit capability is implemented in many software that analyzes data. However, it is often convenient to do the fitting programmatically, especially when a large number of datasets should be analyzed. Furthermore, depending on the software used, some features may not be implemented. In this mini-review, I discuss the basis of polynomial fitting, including the calculation of errors on the coefficients and results, use of weighting and fixing the intercept value (the coefficient 𝛽0). A simple C++ code to perform the polynomial curve fitting is also provided. This code should be useful not only in radiobiology but in other fields of science as well.

Ianik Plante↗

Application of community data to surface complexation modeling framework development: Iron oxide protolysis

This study presents a comprehensive community data-driven surface complexation modeling framework for simulating potentiometric titration of mineral surfaces. Compiled community data for ferrihydrite, goethite, hematite, and magnetite are fit to produce representative protolysis constants that can reproduce potentiometric titration data collected from multiple literature sources. Using this framework, the impact of surface complexation model type and surface site density (SSD) on the fit quality and protolysis constants can be readily evaluated. For example, the non-electrostatic model yielded a poor data fit compared to diffuse double layer model and constant capacitance models due to the absence of known surface charge effects. Regardless of the choice of iron oxide mineral, pK a1 decreased with increasing SSD while the opposite tendency was observed for pK a2 . This newly developed framework demonstrates a method to reconcile community data-wide potentiometric titration data using Findable, Accessible, Interoperable, Reusable data principles to produce mineral protolysis constants that improve robustness of surface complexation models for applications in metal sorption and reactive transport modeling. The framework is readily expandable (as community data increase) and extensible (as the number of minerals increase). The framework provides a path forward for developing self-consistent, comprehensive, and updateable surface complexation databases for surface complexation and reactive transport modeling.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

High-pressure structural systematics of dysprosium metal compressed in a neon pressure medium to 182 GPa

Here we present an experimental and theoretical study of dysprosium metal compressed in the soft pressure transmitting medium Ne up to 182 and 300 GPa, respectively. Angle-dispersive x-ray powder diffraction data from each of the high-pressure polymorphs shows anisotropic compression behavior indicating changes to the electron density distribution throughout its polymorphic landscape. We compare the monoclinic (mC4) and orthorhombic (oF16) structures for the collapsed structure for Dy above 82 GPa and verify that the oF16 structure offers a better fit to our data than the previously reported mC4 structure. Further, we have found that the oF16 structure undergoes similar anisotropic compression of its lattice parameters, with a turning point above 160 GPa; suggesting a potential phase transition at pressures much higher than achieved in this study. Density functional theory calculations show the likely candidate for this new high-pressure phase is the isosymmetric oF8 structure, which is predicted to be lower in energy than the oF16 structure above 275 GPa.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Long-term, hard X-ray observations of Scorpius X-1 from HEAO 1

The results of observations of Sco X-1 from three separate scanning epochs during 1977 and 1978 by the UCSD/MIT instruments aboard HEAO 1 in the 15-200 keV range are reported. Intensity and spectral variations have been studied on three time scales: 33-minute scans, daily averages, and weekly averages. No relationship was found between intensity and temperature over the very narrow range of temperatures observed, but comparison of the variability of the intensity, emission measure, and temperature implies at least two different emission mechanisms are at work. It is found that the spectral data are best fitted by either a thin-thermal-bremsstrahlung model of about 5 keV or a Comptonized thermal model of temperature about 3 keV. At higher energies (50-200 keV), no flux is found in excess of the thermal (as claimed by previous observers), and a 2-sigma upper limit from scanning data of 10 to the -5th photons/sq cm sec keV is set. This is an order of magnitude lower than reported detections. Using the distance derived from optical and UV measurements, X-ray luminosities from Sco X-1 incompatible with predictions for spherical accretion onto a white dwarf and emission volumes inconsistent with polar caps of white dwarfs are found.

Soong, Y.↗

Limits on coronal reflection using high-frequency solar oscillations

Acoustic waves in the Sun with frequencies above about 5.3 mHz can propagate in the chromosphere. We examine imaged solar intensity data for evidence of reflection of these waves in the upper chromosphere, where the temperature increases by a large factor over a short distance. Our method is to compare the observed and theoretically derived frequency spacings between peaks in the power spectrum. We find that our theoretical frequencies provide the best fit to the data when the reflection in the upper atmosphere is eliminated. In particular, the model of Kumar (1993b), which includes the source depth, and radiative damping, in the calculation of power spectra but ignores chromospheric reflection, gives peak frequencies that are in good agreement with the observations. For acoustic waves of frequency greater than 6 mHz we put an upper limit to the reflectivity of chromosphere and corona, using our method, of about 10%. At a given spherical harmonic degree, the frequency spacing between peaks in the data generally decreases with increasing frequency, because the lower turning point of the waves is moving inward. However, between 5 and 5.5 mHz the frequency spacing increases slightly. This feature is probably associated with the acoustic cutoff frequency in the solar atmosphere, i.e., it indicates a transition from trapped waves to propagating waves. We are able to reproduce the observed behavior by a crude modeling of the solar atmosphere. Further study of these peaks should provide an independent way of exploring the mean structure of the solar atmosphere, particularly around the temperature minimum region.

Kumar, P.↗

Nonlinear Aerodynamic Modeling From Flight Data Using Advanced Piloted Maneuvers and Fuzzy Logic

Results of the Aeronautics Research Mission Directorate Seedling Project Phase I research project entitled "Nonlinear Aerodynamics Modeling using Fuzzy Logic" are presented. Efficient and rapid flight test capabilities were developed for estimating highly nonlinear models of airplane aerodynamics over a large flight envelope. Results showed that the flight maneuvers developed, used in conjunction with the fuzzy-logic system identification algorithms, produced very good model fits of the data, with no model structure inputs required, for flight conditions ranging from cruise to departure and spin conditions.

Brandon, Jay M.↗

A Double-Wiebe Function for Reactivity Controlled Compression Ignition Combustion Using Reformate Diesel

Abstract Reactivity controlled compression ignition (RCCI) combustion has previously been proposed as a method to achieve high fuel conversion efficiency and reduce engine emissions. A single-fuel RCCI combustion strategy can have decreased fuel system complexity by using a reformate fuel for port fuel injection and the parent fuel (diesel) for direct injection. This paper presents a one-dimensional computational model of a compression ignition engine with single-fuel RCCI. A Wiebe function is used to predict the combustion process by representing the mass fraction burned (MFB) on a crank angle resolved basis. One single-Wiebe function (SWF) and two double-Wiebe functions(DWFs) were fitted to experimentally derive MFB data using the least-square method. The fitted results were compared with MFBs calculated from experimental data to verify the accuracy. The SWF did not fully capture the MFB curve with high fidelity while the detailed DWF captured the MFB curve within a root mean square error of 1.4%. The reduced double-Wiebe function (RDWF) also resulted in a predicted combustion profile with similar accuracy. Hence, the RDWF was used in a GT-power thermodynamic study to understand the effects of the low-temperature heat release (LTHR) fraction and combustion phasing on combustion characteristics. At optimum phasing of 5–10 crank angle degree after the top dead center, increasing the LTHR fraction from 20% to 60% resulted in the fuel conversion efficiency increasing from 39.5% to 41.1%, thus suggesting that the reformate fuel-based RCCI strategy is viable to unlock improved combustion performance.

Energy & Fuels↗

Methods—Analyzing Electrochemical Kinetic Parameters in Deep Eutectic Solvents Using an Extended Butler-Volmer Equation

Deep eutectic solvents (DESs) are promising electrolytes for electrochemical redox reactions, which can be used in redox flow batteries (RFBs). However, in some systems like the Fe 2+/3+ redox reaction in ethaline, traditional Tafel-based kinetic analysis generates unreasonable kinetic parameters (i.e., large anodic/cathodic charge transfer coefficients ( α a , α c ) along with low exchange current densities ( i 0 )). This hinders a comprehensive kinetic and kinetic mechanism study. Here, we perform a detailed investigation of the Tafel analysis using a series of synthetic rotating disk electrode (RDE) data. We find the Tafel analysis only works well when i 0 < 0.57 i lim (limiting current density) in our scenario and leads to abnormal kinetic values once i 0 exceeds this limiting value. Thus, we propose an extended Butler-Volmer (ex-BV) analysis based on modern non-linear fitting techniques to obtain the actual kinetic parameters for such systems. The results show that this method fits the RDE data closely and generates reliable α a , α c and i 0 values, demonstrating that it is a good replacement for traditional Tafel analysis for kinetic studies in high-viscosity electrolytes such as DES systems.

Electrochemistry↗

Effect of contrast on the active control of a moving line

In many passive visual tasks, human perceptual judgments are contrast dependent. To explore whether these contrast dependencies of visual perception also affect closed-loop manual control tasks, we examined visuomotor performance as humans actively controlled a moving luminance-defined line over a range of contrasts. Four subjects were asked to use a joystick to keep a horizontal line centered on a display as its vertical position was perturbed by a sum of sinusoids under two control regimes. The total root mean square (RMS) position error decreased quasi-linearly with increasing log contrast across the tested range (mean slope across subjects: -8.0 and -7.7% per log2 contrast, for the two control regimes, respectively). Frequency-response (Bode) plots showed a systematic increase in open-loop gain (mean slope: 1.44 and 1.30 dB per log2 contrast, respectively), and decrease in phase lag with increasing contrast, which can be accounted for by a decrease in response time delay (mean slope: 32 and 40 ms per log2 contrast, respectively). The performance data are well fit by a Crossover Model proposed by McRuer and Krendel, which allowed us to identify both visual position and motion cues driving performance. This analysis revealed that the position and motion cues used to support manual control under both control regimes appear equally sensitive to changes in stimulus contrast. In conclusion, our data show that active control of a moving visual stimulus is as dependent on contrast as passive perception and suggest that this effect is attributed to a shared contrast sensitivity early in the visual pathway, before any specialization for motion processing.

Contrast Sensitivity/physiology↗

Kids Making Noise: Presenting Noise Spectra of Aluminum MKIDS for Astronomy

Since the field has identified microwave kinetic inductance detectors (MKIDs) as a competitive candidate for future applications in infrared and sub-millimeter wave detection, it is increasingly important to establish the materials and techniques that produce highly sensitive, well-understood, and consistent detectors. We present new noise power spectral density (PSD) measurements for several thin film (25nm) aluminum coplanar waveguide (CPW) microwave resonators. We include noise PSD measurements for several bath temperatures and photon occupation numbers. We examine the trends exhibited by the fractional frequency noise spectra of our detectors at 1 kHz as a function of bath temperature and power. For our lowest temperature data (58mK) and a readout power of -84 dBm we observe fractional frequency noise in the range Sбf / f r 2 = a few x10 -17 . We observe a roughly Sбf / f r 2 (1 kHz) ∝ Pint -0.5 relationship between internal power and fractional frequency noise at 1 kHz. We attempted more complex fits to the PSDs, but concluded that more measurements are needed to converge on meaningful fits to the data.

MKIDs↗

Velocity and Temperature Measurements in High-speed Flows with Naturally Present Dust Particles Using Rayleigh and Mie Scattering

Dust particles and occasional moisture condensations are unavoidable reality of all wind tunnels. On the path to pursue a goal of velocity and temperature measurements in large transonic and supersonic wind tunnels we have created a tabletop, spectrally resolved, Rayleigh-Mie scattering setup around a small jet fed by ambient and lightly seeded air to determine the viability and accuracy of the technique. The other reality of a wind tunnel setup is the background scattering or the glare at laser frequency, which contaminates the Rayleigh-Mie scattered light. This is simulated by backgrounds with different reflectivity towards the collection optics. Light from a CW laser is delivered via an optical fiber and the scattered light is spectrally resolved using a stabilized Fabry-Perot interferometer, followed by imaging on an EMCCD camera. A model of the of the combined background glare, Mie scattering, and the Rayleigh spectrum was fitted to the camera image using maximum likelihood estimation. Since the background glare occurs at the known frequency of the incident light and the Mie scattering peak corresponds to the Doppler shift from the bulk velocity, both were easily identified, and provided a measure of flow velocity. Preliminary results show that the Rayleigh spectrum can also be resolved, which provides a measure of temperature. It is observed that a slight drift of the laser frequency during data collection affected fitting of the model function leading to larger error. A feedback loop-based stabilization system is on development to take advantage of slight tunability of the laser via a piezo-control. Preliminary results are presented in the abstract. More extensive data from a systematic survey will be presented in the final paper.

Rayleigh scattering↗

Real-Time Detection of Dust Devils from Pressure Readings

A method for real-time detection of dust devils at a given location is based on identifying the abrupt, temporary decreases in atmospheric pressure that are characteristic of dust devils as they travel through that location. The method was conceived for use in a study of dust devils on the Martian surface, where bandwidth limitations encourage the transmission of only those blocks of data that are most likely to contain information about features of interest, such as dust devils. The method, which is a form of intelligent data compression, could readily be adapted to use for the same purpose in scientific investigation of dust devils on Earth. In this method, the readings of an atmospheric- pressure sensor are repeatedly digitized, recorded, and processed by an algorithm that looks for extreme deviations from a continually updated model of the current pressure environment. The question in formulating the algorithm is how to model current normal observations and what minimum magnitude deviation can be considered sufficiently anomalous as to indicate the presence of a dust devil. There is no single, simple answer to this question: any answer necessarily entails a compromise between false detections and misses. For the original Mars application, the answer was sought through analysis of sliding time windows of digitized pressure readings. Windows of 5-, 10-, and 15-minute durations were considered. The windows were advanced in increments of 30 seconds. Increments of other sizes can also be used, but computational cost increases as the increment decreases and analysis is performed more frequently. Pressure models were defined using a polynomial fit to the data within the windows. For example, the figure depicts pressure readings from a 10-minute window wherein the model was defined by a third-degree polynomial fit to the readings and dust devils were identified as negative deviations larger than both 3 standard deviations (from the mean) and 0.05 mbar in magnitude. An algorithm embodying the detection scheme of this example was found to yield a miss rate of just 8 percent and a false-detection rate of 57 percent when evaluated on historical pressure-sensor data collected by the Mars Pathfinder lander. Since dust devils occur infrequently over the course of a mission, prioritizing observations that contain successful detections could greatly conserve bandwidth allocated to a given mission. This technique can be used on future Mars landers and rovers, such as Mars Phoenix and the Mars Science Laboratory.

Wagstaff, Kiri↗

Evaluation of models to predict insolation on tilted surfaces

An empirical study was performed to evaluate the validity of various insolation models which employ either an isotropic or an anisotropic distribution approximation for sky light when predicting insolation on tilted surfaces. Data sets of measured hourly insolation values were obtained over a 6-month period using pyranometers which received diffuse and total solar radiation on a horizontal plane and total radiation on surfaces tilted toward the equator at 37 degrees and 60 degrees angles above the horizon. Data on the horizontal surfaces were used in the insolation models to predict insolation on the tilted surface; comparisons of measured vs calculated insolation on the tilted surface were examined to test the validity of the sky light approximations. It was found that the Liu-Jordan isotropic distribution model provides a good fit to empirical data under overcast skies but underestimates the amount of solar radiation incident on tilted surfaces under clear and partly cloudy conditions.

Klucher, T. M.↗

Adsorption of terbium (III) on DGA and LN resins: Thermodynamics, isotherms, and kinetics

Two commercially available extraction chromatography (EXC) resins containing N,N,N’,N’-tetra-n-octyldiglycolamide (DGA Resin, Normal, 50 – 100 μm) and Bis(2-ethylhexyl) phosphate (LN Resin, 100 – 150 μm) were used as adsorbents to study fundamental adsorption properties such as thermodynamic values, equilibrium isotherms, and kinetic uptake models for terbium(III) adsorption. Weight distribution ratios (D w ) for terbium on DGA and LN resins were measured using a [ 160 Tb]Tb 3+ radiometric tracer in nitric acid as a function of acidity, temperature, initial analyte concentration, and equilibrium time. The D w values showed increasing binding affinity for DGA resin at high nitric acid concentrations and decreasing binding affinity for LN resins. Thermodynamic studies for DGA and LN resins revealed that the Gibbs free energy (ΔG) increased consistently with temperature. To model equilibrium data, increasingly higher parameter equilibrium isotherm models (Henry (1) < Langmuir, Freundlich (2) < Redlich-Peterson (3) < Fritz-Schluender (4)) were compared on their root mean squared errors (RMSE) and adjusted determination coefficients to determine the most applicable model. In all cases, the empirical four-parameter Fritz-Schluender isotherm demonstrated a superior fit. Similar comparisons for reaction-based kinetic models (Pseudo-first-order < Pseudo-second-order < Pseudo-n-order) revealed that the higher-order PNO model yielded a superior fit of kinetic data for both resins. Furthermore, in some cases, adsorption isotherms and kinetic models could also be modeled by a lower-order model with minimal change in error parameters. Weber-Morris plots revealed that two linear sections are observed for each resin, where the first linear segment is attributed to fast (film diffusion) adsorption of terbium, followed by slower intraparticle diffusion of terbium through the pores as the rate-limiting step. Based on the Weber-Morris plot, both film and intraparticle diffusion are involved in controlling the kinetic rate of adsorption for DGA and LN resins.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Magneto-structural studies of an unusual [Mn III Mn II Gd III (OR) 4 ] 4– partial cubane from 2,2'-bis- p- t Bu-calix[4]arene

Reaction of 2,2'-bis- p - t Bu-calix[4]arene (H 8 L) with MnCl 2 ·4H 2 O, GdCl 3 ·6H 2 O and 2,6-pyridinedimethanol (H 2 pdm) affords [Mn III Mn II Gd III (H 3 L)(pdmH)(pdm)(MeOH) 2 (dmf)]·3MeCN·dmf ( 3 ·3MeCN·dmf) upon vapour diffusion of MeCN into the basic dmf/MeOH mother liquor. 3 crystallises in the tetragonal space group P 4 1 2 1 2 with the asymmetric unit comprising the entire cluster. The highly unusual core contains a triangular arrangement of Mn III Mn II Gd III ions housed within a [Mn III Mn II Gd III (OR) 4 ] 4– partial cubane. Magnetic susceptibility and magnetisation data reveal best fit parameters J Mn(II)–Mn(III) = +0.415 cm –1 , J Mn(III)–Gd(III) = +0.221 cm –1 , J Mn(II)–Gd(III) = –0.258 cm –1 and D Mn(III) = –4.139 cm –1 . Theoretically derived magnetic exchange interactions, anisotropy parameters, and magneto-structural correlations for 3 are in excellent agreement with the experimental data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Constraining the $\overline{K}$N coupled channel dynamics using femtoscopic correlations at the LHC

The interaction of K – with protons is characterised by the presence of several coupled channels, systems like $\overline{K}$ 0 n and πΣ with a similar mass and the same quantum numbers as the K – p state. The strengths of these couplings to the K – p system are of crucial importance for the understanding of the nature of the Λ(1405) resonance and of the attractive K – p strong interaction. In this article, we present measurements of the K – p correlation functions in relative momentum space obtained in pp collisions at $\sqrt{s}$=13 Te, in p–Pb collisions at $\sqrt{s_{NN}}$= 5.02 Te, and (semi)peripheral Pb–Pb collisions at $\sqrt{s_{NN}}$= 5.02 Te. The emitting source size, composed of a core radius anchored to the K + p correlation and of a resonance halo specific to each particle pair, varies between 1 and 2 fm in these collision systems. The strength and the effects of the $\overline{K}$ 0 n and πΣ inelastic channels on the measured K – p correlation function are investigated in the different colliding systems by comparing the data with state-of-the-art models of chiral potentials. A novel approach to determine the conversion weights ω, necessary to quantify the amount of produced inelastic channels in the correlation function, is presented. In this method, particle yields are estimated from thermal model predictions, and their kinematic distribution from blast-wave fits to measured data. The comparison of chiral potentials to the measured K – p interaction indicates that, while the πΣ –K – p dynamics is well reproduced by the model, the coupling to the $\overline{K}$ 0 n channel in the model is currently underestimated.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Laboratory measurements and modeling of microwave absorption by ammonia in gas mixtures applicable to giant planet atmospheres

Accurate knowledge of the microwave absorption behavior of ammonia is critical to the correct interpretation of radio astronomical and radio occultation data from the giant planets. New cavity resonator techniques developed at the Stanford Center for Radar Astronomy have allowed accurate laboratory measurements of the microwave absorptivity and refractivity spectra of gas mixtures containing trace amounts of ammonia. A parameterized version of the modified Ben-Reuven formalism of Berge and Bulkis was optimized to fit the new data. The new formalism produced by this method predicts ammonia absorptivity much more accurately than previous formalism over a significant range of conditions.

Spilker, T. R.↗