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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

The shock physics of giant impacts: Key requirements for the equations of state

We discuss major challenges in modeling giant impacts between planetary bodies, focusing on the equations of state (EOS). During the giant impact stage of planet formation, rocky planets are melted and partially vaporized. However, most EOS models fail to reproduce experimental constraints on the thermodynamic properties of the major minerals over the required phase space. Here, we present an updated version of the widely-used ANEOS model that includes a user-defined heat capacity limit in the thermal free energy term. Our revised model for forsterite (Mg 2 SiO 4 ), a common proxy for the mantles of rocky planets, provides a better fit to material data over most of the phase space of giant impacts. We discuss the limitations of this model and the Tillotson equation of state, a commonly used alternative model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Determination of acoustic nonlinearity parameters using thermal modulation of ultrasonic waves

This study presents a test method and its theoretical framework to determine the acoustic nonlinearity parameters (α, β, δ) of material using thermal modulation of ultrasonic waves. Temperature change-induced thermal strain excites the nonlinear response of the material and modulates the ultrasonic wave propagating in it. Experimental results showed a strong correlation between the relative wave velocity change and the temperature change. With a quadratic polynomial model, the acoustic nonlinearity parameters were obtained from the polynomial coefficients by curve fitting the experimental data. Their effects on thermal-induced velocity change were discussed. The parameters α, β, and δ govern the hysteretic gap, average slope, and curvature of the correlation curve, respectively. The proposed theory was validated on aluminum, steel, intact and damaged concrete samples. The obtained nonlinear parameters show reasonable agreement with values reported in the literature. Compared to other nonlinear acoustic methods using vibration or acoustic excitation, the thermal modulation method generates more uniform, slow changing, and larger strain field in the test sample. Employing the thermal effect as the driving force for nonlinearity instead of an undesired influencing factor, this method can measure the absolute values of α, β, and δ with good accuracy using a simple ultrasonic test setup.

Physics↗
DOE OSTI2022

Transducerless time domain reflectance measurement of semiconductor thermal properties

The dominant techniques for laser-based thermal measurement, time-domain thermoreflectance (TDTR) and frequency-domain thermoreflectance, generally require that samples be coated with a metal transducer layer. The properties of the transducer and its interface to the sample are then important to the interpretation of results, and heat conduction in the transducer layer can limit measurements of ultra-thin samples or in-plane conductivity. Here, we demonstrate the feasibility of measuring semiconductor samples using standard TDTR instrumentation, without the use of a metal transducer. A mathematical model accounting for non-zero laser penetration depth as well as both thermal and charge carrier effects is presented, and the ability of the model to describe and fit to experimental data for silicon and germanium samples is demonstrated. Though at this stage measurement uncertainties are larger than for traditional TDTR, we anticipate that the use of this technique will expand the range of samples that can be measured with existing TDTR setups.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Optical readout of singlet fission biexcitons in a heteroacene with photoluminescence detected magnetic resonance

Molecular spin systems based on photoexcited triplet pairs formed via singlet fission (SF) are attractive as carriers of quantum information because of their potentially pure and controllable spin polarization, but developing systems that offer optical routes to readout as well as initialization is challenging. Herein, we characterize the electron spin magnetic resonance change in the photoluminescence intensity for a tailored organic molecular crystal while sweeping a microwave drive up to 10 GHz in a broadband loop structure. We observe resonant transitions for both triplet and quintet spin sublevel populations showing their optical sensitivity and revealing the zero-field parameters for each. We map the evolution of these spectra in both microwave frequency and magnetic field, producing a pattern of optically detected magnetic resonance (ODMR) peaks. Fits to these data using a suitable model suggest significant spin polarization in this system with orientation selectivity. Unusual excitation intensity dependence is also observed, which inverts the sign of the ODMR signal for the triplet features, but not for the quintet. These observations demonstrate optical detection of the spin sublevel population dictated by SF and intermolecular geometry, and highlight anisotropic and multi-scale dynamics of triplet pairs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Acoustic properties, elasticity, and equation of state of glycerol under pressure

Herein we employed high-pressure Brillouin scattering to study the pressure dependencies of acoustic modes of glycerol up to 14 GPa at 300 K. We observed longitudinal acoustic velocities and transverse acoustic velocities for the first time from 5 to 14 GPa. The results allow the determination of a complete set of elastic properties and an accurate determination of the pressure–volume (P–V) equation of state (EOS). EOS parameters, K 0 = 14.9 ± 1.8 GPa and K' 0 = 5.6 ± 0.5, were determined from fits to the data from ambient pressure to 14 GPa. Direct volume measurements of the P–V EOS are consistent with those determined by Brillouin scattering. A deviation from a Cauchy-like relationship for elastic properties was observed, and the pressure dependencies of the photoelastic constants and relaxation times were documented from 5 to 14 GPa. These results have broad implications for glass-forming liquids, viscoelastic theory, and mode coupling theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Calibration and validation of the foundation for a multiphase strength model for tin

In this work, the Common Model of Multi-phase Strength and Equation of State (CMMP) model was applied to tin. Specifically, calibrations of the strength-specific elements of the CMMP foundation were developed with a combination of experiments and theory, and then the model was validated experimentally. The first element of the foundation is a multi-phase analytic treatment of the melt temperature and the shear modulus for the solid phases. These models were parameterized for each phase based on ab initio calculations using the software VASP (Vienna Ab initio Simulations Package) based on density functional theory. The shear modulus model for the ambient phase was validated with ultrasonic sound speed measurements as a function of pressure and temperature. The second element of the foundation is a viscoplastic strength model for the phase, upon which strength for inaccessible higher-pressure phases can be scaled as necessary. The stress–strain response of tin was measured at strain rates of 10 -3 to 3 x 10 3 s -1 and temperatures ranging from 87 to 373 K. The Preston–Tonks–Wallace (PTW) strength model was fit to that data using Bayesian model calibration. For validation, six forward and two reverse Taylor impact experiments were performed at different velocities to measure large plastic deformation of tin at strain rates up to ⁠10 5 s -1 . The PTW model accurately predicted the deformed shapes of the cylinders, with modest discrepancies attributed to the inability of PTW to capture the effects of twinning and dynamic recrystallization. Some material in the simulations of higher velocity Taylor cylinders reached the melting temperature, thus testing the multiphase model because of the presence of a second phase, the liquid. In simulations using a traditional modeling approach, the abrupt reduction of strength upon melt resulted in poor predictions of the deformed shape and non-physical temperatures. With CMMP, the most deformed material points evolved gradually to a mixed solid–liquid but never a fully liquid state, never fully lost strength, stayed at the melt temperature as the latent heat of fusion was absorbed, and predicted the deformed shape well.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Validation of predictive performance models for supersonic gas-jet nozzles at the Laboratory for Laser Energetics

We present results characterizing the neutral-density distributions produced by the supersonic nozzles used in experiments on the OMEGA-60 and OMEGA-EP laser systems at the University of Rochester’s Laboratory for Laser Energetics (LLE). Axisymmetric Fluent® simulations using LLE nozzle specifications capture the viscous effects, gas expansion, and shock waves that complicate flow predictions for offsets above the nozzle exit. Here, these simulations show good agreement with neutral-density measurements obtained using a four-wave shearing interferometer. An analytical form is given for the plateau length. Fits to simulation data for boundary layer thickness, mean plateau density, and density ramps are given as functions of nozzle offset and nozzle backing pressure for a number of nozzles and gases.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Predicting the viscoplastic response of a crystallizing fluoropolymer using transient network theory

We employ a molecular theory of dynamic polymer networks to describe the viscoplastic response of rubbery FK-800, a thermoplastic copolymer of chlorotrifluoroethylene and vinylidene fluoride, over a broad range of thermal histories. The kinetics of crystallization at different annealing temperatures was modeled using a modified Avrami equation, whose parameters were found to evolve through simple relationships over the full temperature range of the rubbery state. By fitting experimental compression data, we discovered predictable trends for the physical parameters in our mechanical model over its full range of crystallinities (up to ≈20%) and provided insights based on molecular-level physics to justify them. Using this, an end-to-end model was developed to predict the yielding and post-yield behavior of rubbery FK-800 for arbitrary thermal histories. The model successfully predicted the highly nonlinear evolution of characteristic mechanical signatures (stiffness, yield point, post-yield drop) throughout the crystallization process. A statistical analysis of variance test was employed to determine that the measured variations in the mechanical behavior of rubbery FK-800 are primarily dictated by its fractional crystallinity, regardless of its exact thermal history.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Deriving reliable nucleation rates from metadynamics simulations: Application to Yukawa fluids

In order to solidify the usefulness of metadynamics in studying nucleation of crystals from supercooled liquids, this work provides a specific procedure to calculate nucleation free energy barriers. After a pedagogical review of the important elements of classical nucleation theory and how metadynamics is used to find nucleation free energy barriers, we explain the benefits of local collective variables over more common global collective variables. We show how a metadynamics free energy barrier must be carefully post-processed so that classical nucleation theory can be applied to calculate nucleation rates. We apply our procedure to a Yukawa plasma and show that a particular physically motivated fit to metadynamics data reproduces low-temperature reference data, justifying the usefulness of metadynamics to predict nucleation rates.

42 ENGINEERING↗
DOE OSTI2021

Pseudorapidity distributions of charged particles in pp($\bar p$), p(d)A and AA collisions using Tsallis thermodynamics

The pseudorapidity distributions of charged particles measured in p + p($\overline{\text{p}}$) collisions for energies ranging from $\sqrt{{s}_{\mathrm{N}\mathrm{N}}}=23.6$ GeV to 13 TeV and A + A collisions at RHIC and LHC are investigated in the fireball model with Tsallis thermodynamics. We assume that the rapidity axis is populated with fireballs following q-Gaussian distribution and the charged particles follow the Tsallis distribution in the fireball. We also extend the fireball model to asymmetric collision systems, i.e. d + Au collisions at $\sqrt{{s}_{\mathrm{N}\mathrm{N}}}=200$ GeV and p + Pb collisions at $\sqrt{{s}_{\mathrm{N}\mathrm{N}}}=5.02$ TeV, by taking into account the asymmetric geometry configuration. The model can fit the experimental data well for all the collision systems and centralities investigated. The collision energy and centrality dependence of the model parameters for the symmetric (asymmetric) collision system, i.e. the central position y 0 (y 0a , y 0A ) and its width σ (σ a , σ A ) of the fireball distribution, are also investigated and discussed. Furthermore, our results suggest that the fireball model with Tsallis thermodynamics can be used as a universal framework for the pseudorapidity distributions of charged particles in high energy collisions at RHIC and LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Uniqueness of relaxation times determined by dielectric spectroscopy

Dielectric spectroscopy is extremely powerful to study molecular dynamics, because of the very broad frequency range. Often multiple processes superimpose resulting in spectra that expand over several orders of magnitude, with some of the contributions partially hidden. For illustration, we selected two examples, (i) normal mode of high molar mass polymers partially hidden by conductivity and polarization and (ii) contour length fluctuations partially hidden by reptation using the well-studied polyisoprene melts as example. The intuitive approach to describe experimental spectra and to extract relaxation times is the addition of two or more model functions. Here, we use the empirical Havriliak-Negami function to illustrate the ambiguity of the extracted relaxation time, despite an excellent agreement of the fit with experimental data. We show that there are an infinite number of solutions for which a perfect description of experimental data can be achieved. However, a simple mathematical relationship indicates uniqueness of the pairs of the relaxation strength and relaxation time. Sacrificing the absolute value of the relaxation time enables to find the temperature dependence of the parameters with a high accuracy. For the specific cases studied here, the time temperature superposition (TTS) is very useful to confirm the principle. However, the derivation is not based on a specific temperature dependence, hence, independent from the TTS. We compare new and traditional approaches and find the same trend for the temperature dependence. The important advantage of the new technology is the knowledge of the accuracy of the relaxation times. Relaxation times determined from data for which the peak is clearly visible are the same within the experimental accuracy for traditional and new technology. However, for data where a dominant process hides the peak, substantial deviations can be observed. Finally, we conclude that the new approach is particularly helpful for cases in which relaxation times need to be determined without having access to the associated peak position.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Signatures of superconducting triplet pairing in Ni–Ga-bilayer junctions

Ni–Ga bilayers are a versatile platform for exploring the competition between strongly antagonistic ferromagnetic and superconducting phases. We characterize the impact of this competition on the transport properties of highly-ballistic Al/Al 2 O 3 (/EuS)/Ni–Ga tunnel junctions from both experimental and theoretical points of view. While the conductance spectra of junctions comprising Ni (3 nm)–Ga (60 nm) bilayers can be well understood within the framework of earlier results, which associate the emerging main conductance maxima with the junction films’ superconducting gaps, thinner Ni (1.6 nm)–Ga (30 nm) bilayers entail completely different physics, and give rise to novel large-bias (when compared to the superconducting gap of the thin Al film as a reference) conductance-peak subseries that we term conductance shoulders. These conductance shoulders might attract considerable attention also in similar magnetic superconducting bilayer junctions, as we predict them to offer an experimentally well-accessible transport signature of superconducting triplet pairings that are induced around the interface of the Ni–Ga bilayer. We further substantiate this claim performing complementary polarized neutron reflectometry measurements on the bilayers, from which we deduce (1) a nonuniform magnetization structure in Ga in a several nanometer-thick area around the Ni–Ga boundary and can simultaneously (2) satisfactorily fit the obtained data only considering the paramagnetic Meissner response scenario. While the latter provides independent experimental evidence of induced triplet superconductivity inside the Ni–Ga bilayer, the former might serve as the first experimental hint of its potential microscopic physical origin. Finally, we introduce a simple phenomenological toy model to confirm also from the theoretical standpoint that superconducting triplet pairings around the Ni–Ga interface can indeed lead to the experimentally observed conductance shoulders, which convinces that our claims are robust and physically justified. Arranging our work in a broader context, we expect that Ni–Ga-bilayer junctions could have a strong potential for future superconducting-spintronics applications whenever an efficient engineering of triplet-pairing superconductivity is required.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Dark matter cores and cusps in spiral galaxies and their explanations

We compare proposed solutions to the core vs cusp issue of spiral galaxies, which has also been framed as a diversity problem, and demonstrate that the cuspiness of dark matter halos is correlated with the stellar surface brightness. We compare the rotation curve fits to the SPARC sample from a self-interacting dark matter (SIDM) model, which self-consistently includes the impact of baryons on the halo profile, and hydrodynamical N-body simulations with cold dark matter (CDM). The SIDM model predicts a strong correlation between the core size and the stellar surface density, and it provides the best global fit to the data. The CDM simulations without strong baryonic feedback effects fail to explain the large dark matter cores seen in low surface brightness galaxies. On the other hand, with strong feedback, CDM simulations do not produce galaxy analogs with high stellar and dark matter densities, and therefore they have trouble in explaining the rotation curves of high surface brightness galaxies. This implies that current feedback implementations need to be modified. We also explicitly show how the concentration-mass and stellar-to-halo mass relations together lead to a radial acceleration relation (RAR) in an averaged sense, and reiterate the point that the RAR does not capture the diversity of galaxy rotation curves in the inner regions. These results make a strong case for SIDM as the explanation for the cores and cusps of field galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Cosmological constraints from unWISE and Planck CMB lensing tomography

A number of recent, low-redshift, lensing measurements hint at a universe in which the amplitude of lensing is lower than that predicted from the ΛCDM model fit to the data of the Planck CMB mission. Here we use the auto- and cross-correlation signal of unWISE galaxies and Planck CMB lensing maps to infer cosmological parameters at low redshift. In particular, we consider three unWISE samples (denoted as "blue", "green" and "red") at median redshifts z ~ 0.6, 1.1 and 1.5, which fully cover the Dark Energy dominated era. Our cross-correlation measurements, with combined significance S/N ~ 80, are used to infer the amplitude of low-redshift fluctuations, σ 8 ; the fraction of matter in the Universe, Ω m ; and the combination S 8 ≡ σ 8 (Ω m /0.3) 0.5 to which these low-redshift lensing measurements are most sensitive. The combination of blue, green and red samples gives a value S m = 0.784 ± 0.015, that is fully consistent with other low-redshift lensing measurements and in 2.4σ tension with the CMB predictions from Planck. Finally, this is noteworthy, because CMB lensing probes the same physics as previous galaxy lensing measurements, but with very different systematics, thus providing an excellent complement to previous measurements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Towards fully predictive gyrokinetic full- f simulations: validation and triangularity studies in TCV

Designing economical magnetic confinement fusion power plants motivates computational tools that can estimate plasma behavior from engineering parameters without direct reliance on experimental measurement of the plasma profiles. In this work, we present full-f global long-wavelength gyrokinetic simulations of edge and scrape-off layer turbulence in tokamaks that use only magnetic geometry, heating power, and particle inventory as inputs. Unlike many modeling approaches that employ free parameters fitted to experimental data, raising uncertainties when extrapolating to reactor scales. This approach directly simulates turbulence and resulting profiles through gyrokinetics without such empirical adjustments. This is achieved via an adaptive sourcing algorithm in Gkeyll that strictly controls energy injection and emulates particle sourcing due to neutral recycling. We show that the simulated kinetic profiles compare reasonably well with Thomson scattering and Langmuir probe data for Tokamak á Configuration Variable (TCV) discharge #65125, and that the simulations reproduce characteristic features such as blob transport and self-organized electric fields. Applying the same framework to study triangularity effects suggests mechanisms contributing to the improved confinement reported for negative triangularity (NT). Simulations of TCV discharges #65125 and #65130 indicate that NT increases the E x B flow shear (by about 20% in these cases), which correlates with reduced turbulent losses and a modest change in the distribution of power exhaust to the vessel wall. While the physical models contain approximations that can be refined in future work, the predictive capability demonstrated here, evolving multiple profile relaxation times with kinetic electron and ion models in hundreds of GPU hours, indicates the feasibility of using Gkeyll to support design studies of fusion devices.

Hoffmann, Antoine Cyril David [Princeton Plasma Ph↗
DOE OSTI2020

The mystery of photometric twins DES17X1boj and DES16E2bjy

In this paper, we present an analysis of DES17X1boj and DES16E2bjy, two peculiar transients discovered by the Dark Energy Survey (DES). They exhibit nearly identical double-peaked light curves that reach very different maximum luminosities ($M_r$ = -15.4 and -17.9, respectively). The light-curve evolution of these events is highly atypical and has not been reported before. The transients are found in different host environments: DES17X1boj was found near the nucleus of a spiral galaxy, while DES16E2bjy is located in the outskirts of a passive red galaxy. Early photometric data are well fitted with a blackbody and the resulting moderate photospheric expansion velocities (1800 km s –1 for DES17X1boj and 4800 km s –1 for DES16E2bjy) suggest an explosive or eruptive origin. Additionally, a feature identified as high-velocity Ca II absorption ($v ≈$ 9400 km s –1 ) in the near-peak spectrum of DES17X1boj may imply that it is a supernova. While similar light-curve evolution suggests a similar physical origin for these two transients, we are not able to identify or characterize the progenitors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Signatures of Light Massive Relics on non-linear structure formation

ABSTRACT Cosmologies with Light Massive Relics (LiMRs) as a subdominant component of the dark sector are well-motivated from a particle physics perspective, and can also have implications for the σ8 tension between early and late time probes of clustering. The effects of LiMRs on the cosmic microwave background (CMB) and structure formation on large (linear) scales have been investigated extensively. In this paper, we initiate a systematic study of the effects of LiMRs on smaller, non-linear scales using cosmological N-body simulations; focusing on quantities relevant for photometric galaxy surveys. For most of our study, we use a particular model of non-thermal LiMRs but the methods developed generalizing to a large class of LiMR models – we explicitly demonstrate this by considering the Dodelson–Widrow velocity distribution. We find that, in general, the effects of LiMR on small scales are distinct from those of a ΛCDM universe, even when the value of σ8 is matched between the models. We show that weak lensing measurements around massive clusters, between ∼0.1 h−1Mpc and ∼10 h−1Mpc, should have sufficient signal-to-noise in future surveys to distinguish between ΛCDM and LiMR models that are tuned to fit both CMB data and linear scale clustering data at late times. Furthermore, we find that different LiMR cosmologies indistinguishable by conventional linear probes can be distinguished by non-linear probes if their velocity distributions are sufficiently different. LiMR models can, therefore, be best tested by jointly analyzing the CMB and late-time structure formation on both large and small scales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The metabolic origins of non-photorespiratory CO 2 release during photosynthesis: a metabolic flux analysis

Respiration in the light (R L ) releases CO 2 in photosynthesizing leaves and is a phenomenon that occurs independently from photorespiration. Since R L lowers net carbon fixation, understanding R L could help improve plant carbon-use efficiency and models of crop photosynthesis. Although R L was identified more than 75 years ago, its biochemical mechanisms remain unclear. To identify reactions contributing to R L , we mapped metabolic fluxes in photosynthesizing source leaves of the oilseed crop and model plant camelina (Camelina sativa). We performed a flux analysis using isotopic labeling patterns of central metabolites during 13 CO 2 labeling time course, gas exchange, and carbohydrate production rate experiments. To quantify the contributions of multiple potential CO 2 sources with statistical and biological confidence, we increased the number of metabolites measured and reduced biological and technical heterogeneity by using single mature source leaves and quickly quenching metabolism by directly injecting liquid N 2 ; we then compared the goodness-of-fit between these data and data from models with alternative metabolic network structures and constraints. Our analysis predicted that R L releases 5.2 μmol CO2 g –1 FW h –1 of CO 2 , which is relatively consistent with a value of 9.3 μmol CO 2 g –1 FW h –1 measured by CO 2 gas exchange. The results indicated that ≤10% of R L results from TCA cycle reactions, which are widely considered to dominate R L . Further analysis of the results indicated that oxidation of glucose-6-phosphate to pentose phosphate via 6-phosphogluconate (the G6P/OPP shunt) can account for >% of CO 2 released by R L .

59 BASIC BIOLOGICAL SCIENCES↗