Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Curve 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

Mechanistic definition and prediction of the mass exchange coefficient between rivers and hyporheic zones: The $α$ of two $Ω$s

Solute transport in interconnected rivers and hyporheic zones is typically modeled through dual-domain models where first-order solute mass transfer between the two domains, Ω R and Ω HZ , is represented by a coefficient α. The transient storage model (TSM) is an example of such an approach. In practice, α is determined by fitting the tails of solute tracer breakthrough curves using a TSM. This approach has led to ambiguity regarding α’s physical meaning and transferability. Here, in this work, we investigated the physical basis for α and tested it with virtual experiments through the fully coupled multiphysics model hyporheicFoam for the Ω R – Ω HZ system. hyporheicFoam explicitly simulated coupled flow and solute transport over a kilometer with centimeter-scale resolution. Model results were analyzed to calculate α following its theoretical definition directly. Using the determined α within a TSM enables accurate reproduction of solute transport, underscoring α’s physical relevance and precision.

54 ENVIRONMENTAL SCIENCES↗

The Identity and Chemistry of C 7 H 7 Radicals Observed during Soot Formation

Here we used aerosol mass spectrometry coupled with tunable synchrotron photoionization to measure radical and closed-shell species associated with particle formation in premixed flames and during pyrolysis of butane, ethylene, and methane. We analyzed photoionization (PI) spectra for the C 7 H 7 radical to identify the isomers present during particle formation. For the combustion and pyrolysis of all three fuels, the PI spectra can be fit reasonably well with contributions from four radical isomers: benzyl, tropyl, vinylcyclopentadienyl, and o-tolyl. Although there are significant experimental uncertainties in the isomeric speciation of C7H 7 , the results clearly demonstrate that the isomeric composition of C 7 H 7 strongly depends on the combustion or pyrolysis conditions and the fuel or precursors. Fits to the PI spectra using reference curves for these isomers suggest that all of these isomers may contribute to m/z 91 in butane and methane flames, but only benzyl and vinylcyclopentadienyl contribute to the C 7 H 7 isomer signal in the ethylene flame. Only tropyl and benzyl appear to play a role during pyrolytic particle formation from ethylene, and only tropyl, vinylcyclopentadienyl, and o-tolyl appear to participate during particle formation from butane pyrolysis. There also seems to be a contribution from an isomer with an ionization energy below 7.5 eV for the flames but not for the pyrolysis conditions. Kinetic models with updated and new reactions and rate coefficients for the C 7 H 7 reaction network predict benzyl, tropyl, vinylcyclopentadienyl, and o-tolyl to be the primary C 7 H 7 isomers and predict negligible contributions from other C 7 H 7 isomers. These updated models provide better agreement with the measurements than the original versions of the models but, nonetheless, underpredict the relative concentrations of tropyl, vinylcyclopentadienyl, and o-tolyl in both flames and pyrolysis and overpredict benzyl in pyrolysis. Our results suggest that there are additional important formation pathways for the vinylcyclopentadienyl, tropyl, and o-tolyl radicals and/or loss pathways for the benzyl radical that are currently unaccounted for in the present models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multi-energy calibration of a PILATUS3 CdTe detector for hard x-ray measurements of magnetically confined fusion plasmas

A multi-energy hard x-ray pin-hole camera based on the PILATUS3 X 100K-M CdTe detector has been developed at the Princeton Plasma Physics Laboratory for installation on the Tungsten Environment in Steady State Tokamak. This camera will be employed to study thermal plasma features such as electron temperature as well as non-thermal effects such as fast electron tails produced by a lower hybrid radiofrequency current drive and the birth of runaway electrons. The innovative aspect of the system lies in the possibility of setting the threshold energy independently for each of the ~100k pixels of the detector. This feature allows for the measurement of the x-ray emission in multiple energy ranges with adequate space and time resolution (~1 cm, 2 ms) and coarse energy resolution. In this work, the energy dependence of each pixel was calibrated within the range 15 keV–100 keV using a tungsten x-ray tube and emission from a variety of fluorescence targets (from yttrium to uranium). The data corresponding to pairs of K α emission lines are fit to the characteristic responsivity (“S-curve”), which describes the detector sensitivity across the 64 possible energy threshold values for each pixel; this novel capability is explored by fine-tuning the voltage of a six-bit digital-analog converter after the charge-sensitive amplifier for each of the ~100k pixels. Here, this work presents the results of the calibration including a statistical analysis. It was found that the achievable energy resolution is mainly limited by the width of the S-curve to 3 keV–10 keV for threshold energies up to 50 keV, and to ≥20 keV for energies above 60 keV.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A low-luminosity core-collapse supernova very similar to SN 2005cs

ABSTRACT We present observations and analysis of PSN J17292918+7542390, a low-luminosity Type II-P supernova (LL SN IIP). The observed sample of such events is still low, and their nature is still under debate. Such SNe are similar to SN 2005cs, a well-observed LL Type II-P event, having low expansion velocities, and small ejected 56Ni mass. We have developed a robust and relatively fast Monte Carlo code that fits semi-analytic models to light curves of core-collapse SNe. This allows the estimation of the most important physical parameters, like the radius of the progenitor star, the mass of the ejected envelope, the mass of the radioactive nickel synthesized during the explosion, among others. PSN J17292918+7542390 has $R_0 = 91_{-70}^{+119} \times 10^{11} \, \text{cm}$, $M_\text{ej} = 9.89_{-1.00}^{+2.10} \, \mathrm{ M}_{\odot }$, $E_{\mbox{kin}} = 0.65_{-0.18}^{+0.19} \, \text{foe}$, and $v_{\mbox{exp}} = 3332_{-347}^{+216}$ km s−1, for its progenitor radius, ejecta mass, kinetic energy, and expansion velocity, respectively. The initial nickel mass of the PSN J17292918+7542390 turned out to be $1.55_{-0.70}^{+0.75} \times 10^{-3} \, \mathrm{M}_{\odot }$. The measured photospheric velocity at the earliest observed phase is 7000 km s−1. As far as we can tell based on the small population of observed LL SNe IIP, the determined values are typical for these events.

Jäger, Jr, Zoltán↗

Core-collapse supernovae in the Dark Energy Survey: luminosity functions and host galaxy demographics

ABSTRACT We present the luminosity functions and host galaxy properties of the Dark Energy Survey (DES) core-collapse supernova (CCSN) sample, consisting of 69 Type II and 50 Type Ibc spectroscopically and photometrically confirmed supernovae over a redshift range 0.045 < z < 0.25. We fit the observed DES griz CCSN light curves and K-correct to produce rest-frame R-band light curves. We compare the sample with lower redshift CCSN samples from Zwicky Transient Facility (ZTF) and Lick Observatory Supernova Search (LOSS). Comparing luminosity functions, the DES and ZTF samples of SNe II are brighter than that of LOSS with significances of 3.0σ and 2.5σ, respectively. While this difference could be caused by redshift evolution in the luminosity function, simpler explanations such as differing levels of host extinction remain a possibility. We find that the host galaxies of SNe II in DES are on average bluer than in ZTF, despite having consistent stellar mass distributions. We consider a number of possibilities to explain this – including galaxy evolution with redshift, selection biases in either the DES or ZTF samples, and systematic differences due to the different photometric bands available – but find that none can easily reconcile the differences in host colour between the two samples and thus its cause remains uncertain.

79 ASTRONOMY AND ASTROPHYSICS↗

SALT2 versus SALT3: updated model surfaces and their impacts on type Ia supernova cosmology

ABSTRACT For the past decade, SALT2 has been the most common model used to fit Type Ia supernova (SN Ia) light curves for dark energy analyses. Recently, the SALT3 model was released, which upgraded a number of model features but has not yet been used for measurements of dark energy. Here, we evaluate the impact of switching from SALT2 to SALT3 for a SN cosmology analysis. We train SALT2 and SALT3 on an identical training sample of 1083 well-calibrated Type Ia supernovae, ensuring that any differences found come from the underlying model framework. We publicly release the results of this training (the SALT ‘surfaces’). We then run a cosmology analysis on the public Dark Energy Survey 3-Yr Supernova data sample (DES-SN3YR), and on realistic simulations of those data. We provide the first estimate of the SN + CMB systematic uncertainty arising from the choice of SALT model framework (i.e. SALT2 versus SALT3), Δw = + 0.001 ± 0.005 – a negligible effect at the current level of dark energy analyses. We also find that the updated surfaces are less sensitive to photometric calibration uncertainties than previous SALT2 surfaces, with the average spectral energy density dispersion reduced by a factor of two over optical wavelengths. This offers an opportunity to reduce the contribution of calibration errors to SN cosmology uncertainty budgets.

79 ASTRONOMY AND ASTROPHYSICS↗

Leaf gas exchange and fitted parameters, two sites in Panama, 2022

Photosynthetic CO2 response curves (ACi curves), light response curves (AQ curves), dark adapted dark respiration, conductance curves, and survey measurements for leaves measured in the Parque Natural Metropolitano (PNM) and Gamboa, Panama, from January to April 2022 are presented. Measurements were made on leaves from 46 different tree, shrub and liana species, from sunlit canopy and understory locations. Leaf area index from 8 vertical profiles at PNM are also included. The aim of this measurement campaign was three-fold: (i) to improve our understanding of the vertical variation in leaf-level water use efficiency; (ii) to develop an understanding of the sensitivity of stomata to changes in environmental conditions, especially light, humidity, and temperature; (iii) to improve models which can predict leaf traits from leaf contact spectral measurements. All the gas exchange data and metadata are presented in .csv files and complete instrument output are included in .zip folders. Data and metadata meet the ESS-DIVE leaf-level gas exchange reporting format requirements. The protocol details are provided as pdf documents. In addition to gas exchange data reported here these samples were also used for measurement of leaf optical properties, carbon and nitrogen content, and leaf mass per unit leaf area (LMA). These data can be cross linked using the unique sample ID and are provided in separate data packages.

54 ENVIRONMENTAL SCIENCES↗

Yttrium solid phase equation of state with uncertainty quantification

We discuss the development of an equation of state (EOS) for solid yttrium up to roughly 80 GPa. The EOS makes use of different experimental datasets, including measurements of the isobaric density, isobaric specific heat, room temperature isotherm, and principal shock Hugoniot. The fitting procedure is carried out using Markov Chain Monte Carlo (MCMC), where we fit model parameters for both the cold curve and ion thermal models simultaneously. The results show close agreement with experimental data and provide meaningful uncertainty estimates on the model parameters. This work serves as a first step towards a multiphase EOS for yttrium, which will include higher pressure solid phases (> 80 GPa), as well as modeling of the liquid phase.

36 MATERIALS SCIENCE↗

Optimization and Multimachine Learning Algorithms to Predict Nanometal Surface Area Transfer Parameters for Gold and Silver Nanoparticles

Interactions between gold metallic nanoparticles and molecular dyes have been well described by the nanometal surface energy transfer (NSET) mechanism. However, the expansion and testing of this model for nanoparticles of different metal composition is needed to develop a greater variety of nanosensors for medical and commercial applications. In this study, the NSET formula was slightly modified in the size-dependent dampening constant and skin depth terms to allow for modeling of different metals as well as testing the quenching effects created by variously sized gold, silver, copper, and platinum nanoparticles. Overall, the metal nanoparticles followed more closely the NSET prediction than for Förster resonance energy transfer, though scattering effects began to occur at 20 nm in the nanoparticle diameter. To further improve the NSET theoretical equation, an attempt was made to set a best-fit line of the NSET theoretical equation curve onto the Au and Ag data points. An exhaustive grid search optimizer was applied in the ranges for two variables, 0.1≤C≤2.0 and 0≤α≤4, representing the metal dampening constant and the orientation of donor to the metal surface, respectively. Three different grid searches, starting from coarse (entire range) to finer (narrower range), resulted in more than one million total calculations with values C=2.0 and α=0.0736. The results improved the calculation, but further analysis needed to be conducted in order to find any additional missing physics. With that motivation, two artificial intelligence/machine learning (AI/ML) algorithms, multilayer perception and least absolute shrinkage and selection operator regression, gave a correlation coefficient, R2, greater than 0.97, indicating that the small dataset was not overfitting and was method-independent. This analysis indicates that an investigation is warranted to focus on deeper physics informed machine learning for the NSET equations.

Demers, Steven M. E. (ORCID:0000000192213246)↗

An analytical study of reduced-gravity flow dynamics

Addition of surface tension forces to a marker-and-cell code and the performance of four incompressible fluid simulations in reduced gravity, were studied. This marker-and-cell code has a variable grid capability with arbitrary curved boundaries and time dependent acceleration fields. The surface tension logic includes a spline fit of surface marker particles as well as contact angle logic for straight and curved wall boundaries. Three types of flow motion were simulated with the improved code: impulsive settling in a model Centaur LH2 tank, continuous settling in a model and full scale Centaur LO2 tank and mixing in a Centaur LH2 tank. The impulsive settling case confirmed a drop tower analysis which indicated more orderly fluid collection flow patterns with this method providing a potential savings in settling propellants. In the LO2 tank, fluid collection and flow simulation into the thrust barrel were achieved. The mixing simulation produced good results indicating both the development of the flow field and fluid interface behavior.

Bradshaw, R. D.↗

Extinction-sedimentation inversion technique for measuring size distribution of artificial fogs

In measuring the size distribution of artificial fog particles, it is important that the natural state of the particles not be disturbed by the measuring device, such as occurs when samples are drawn through tubes. This paper describes a method for carrying out such a measurement by allowing the fog particles to settle in quiet air inside an enclosure through which traverses a parallel beam of light for measuring the optical depth as a function of time. An analytic function fit to the optical depth time decay curve can be directly inverted to yield the size distribution. Results of one such experiment performed on artificial fogs are shown as an example. The forwardscattering corrections to the measured extinction coefficient are also discussed with the aim of optimizing the experimental design so that the error due to forwardscattering is minimized.

Deepak, A.↗

Brightness, density, and color profiles of three globular clusters - NGC 6440, 6541, and 7099

Using short-exposure SEC television camera observations, the radial distribution profiles of surface brightness, B-V colors and star densities are determined in three globular clusters, NGC 6440, 6541 and 7099, from about 2-200 arc sec from the cluster centers. Good agreement is found between the surface-brightness profiles and the star-density distributions; both can be fitted well with a King-type curve, and no significant evidence for a presence of cusp at r greater then 3 arc sec is observed. Also determined for the three clusters were central surface brightnesses, core radii, central mass densities, total cluster masses, relaxation times and central escape velocities. In addition, an upper limit of 5000 solar masses is placed on the mass of a possible central black hole in these clusters.

Williams, T. B.↗

Documentation of program AFTBDY to generate coordinate system for 3D after body using body fitted curvilinear coordinates, part 1

The computer program AFTBDY generates a body fitted curvilinear coordinate system for a wedge curved after body. This wedge curved after body is being used in an experimental program. The coordinate system generated by AFTBDY is used to solve 3D compressible N.S. equations. The coordinate system in the physical plane is a cartesian x,y,z system, whereas, in the transformed plane a rectangular xi, eta, zeta system is used. The coordinate system generated is such that in the transformed plane coordinate spacing in the xi, eta, zeta direction is constant and equal to unity. The physical plane coordinate lines in the different regions are clustered heavily or sparsely depending on the regions where physical quantities to be solved for by the N.S. equations have high or low gradients. The coordinate distribution in the physical plane is such that x stays constant in eta and zeta direction, whereas, z stays constant in xi and eta direction. The desired distribution in x and z is input to the program. Consequently, only the y-coordinate is solved for by the program AFTBDY.

Kumar, D.↗

Classification of the Allan Hills A77307 meteorite

Thermoluminescence (TL) measurements on Allan Hills A77307 (AH), a carbonaceous chondrite found in Antarctica, are compared with those on other chondrite and applied to its classification. Two lithologically different 250-mg samples were ground, freed of magnetic material, and ground again to pass a 100-micron sieve. Aliquots of 4 mg were heated to 500 C, exposed to beta radiation from a Sr-90 source, and heated at a rate of 7.3 C/sec in N2. TL was measured with a photomultiplier tube fitted with thermal and blue filters. Glow curves for AH and for seven other, established CO-type chondrites area presented, all exhibiting two major peaks of TL sensitivity. The peaks for the seven CO-type chondrites are found at 91 + or - 7 C and at 203 + or - 11 C; those for AH at 170 + or - 17 C and at approximately 250 C. This difference is considered significant and not due to random fluctuation or a typical sampling. From this comparison and from consideration of the weathering, preterrestrial-alteration, petrological and compositional evidence on AH, it is concluded that AH is a unique chondrite, possessing both similarities to and differences from the CO class.

Sears, D. W. G.↗

The energy spectra of solar energetic particles

A survey of recent results on the shapes and relative slopes of the spectra of various solar energetic particle populations is presented, with emphasis on the more extensive results currently available for protons, alphas and electrons. From previous work, it is found that proton spectra 0.8 to more than 400 MeV and alpha spectra 1.4 to 80 MeV/nucleon are best characterized, on average, by a functional form involving a Bessel function in momentum/nucleon. However, proton and alpha spectral slopes using this form are not equal, and there is significant variation from event to event. From other studies, electrons 0.02 to 20 MeV are also found to have curved spectra, but seem to be better fit with a double power law in energy. The spectral properties in both cases correlate with other measures of solar particle acceleration; e.g. gamma-ray line production, hard X-ray burst spectra and microwave fluxes.

Mcguire, R. E.↗

Interstellar absorption lines in the spectrum of sigma Sco using Copernicus observations

Since the launch of Copernicus in 1972, studies have been made of the depletion of gas-phase elements onto dust grains. A few stars have been studied in detail, resulting in a standard depletion pattern which has since been used for comparison. Recent developments, however, have suggested that this standard pattern may need to be re-examined. Some weak, semi-forbidden lines were detected recently which may be able to resolve some of the ambiguities. Studies of single elements have shown that depletion of carbon and oxgyen are much smaller than previously determined. The high resolution ultraviolet spectral scans of sigma Sco were originally made in 1973, but have only recently been analyzed. All these stars are bright and moderately reddened. All four stars will be analyzed in detail, but sigma Sco is the first one completed. The data has broad coverage of ions, making these stars excellent candidates for determination of accurate depletions. A profile-fitting analysis was used rather than curves-of-growth in order to determine separate abundances and depletions in components separated by several km/sec.

Allen, M. M.↗

Photometric properties of lunar terrains derived from Hapke's equation

Hapke's bidirectional reflectance equations provide the most rigorous available description of photometric behavior in terms of physically meaningful parameters. The primary objective of this study was to derive Hapke parameters for the lunar surface from both disk-integrated and disk-resolved photometric data. Hapke's equation was fit to the disk-integrated phase curves and disk-resolved data for dark, average, and bright terrain classes using an iterative, nonlinear least squares algorithm described by Helfenstein. Parameters were initially determined from the disk-integrated data, and the result was applied as a first guess to the iterative solution of parameters for individual terrain classes. Plots are presented of the disk-resolved data normalized to corresponding brightnesses predicted from the disk-integrated solution under the same illumination and viewing geometries. Systematic trends in disk-resolved parameters can be identified. Values for single scattering albedo (w) of the dominantly anorthositic average and bright terrains are significantly larger than the value for the basaltic dark terrains (mare). Values for surge brightness parameters, particle phase function, and average topographic slope angle are also discussed.

Helfenstein, Paul↗