SEARCH · Engineering Papers

Results for “deviated”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Constraints on Dynamical Dark Energy from Multiple Probes in the Full Dark Energy Survey

We present results on dark energy evolution, assuming a time-dependent equation of state $w(a)=w_0+w_a(1-a)$, from growth and geometric probes using the full six-year Dark Energy Survey dataset: type Ia supernovae, baryon acoustic oscillations, and weak gravitational lensing and galaxy clustering (3$\times$2pt). The combination yields $w_0=-0.84^{+0.10}_{-0.10}$ and $w_a=-0.44^{+0.60}_{-0.55}$, the tightest constraints ever obtained from a single survey, with $2.2σ$ deviation from a cosmological constant. Adding the DESI DR2 BAO data yields $w_0=-0.84^{+0.06}_{-0.07}$ and $w_a=-0.53^{+0.33}_{-0.28}$, representing the most stringent low-redshift-only test of dynamical dark energy to date, with a $2.3σ$ deviation. In this combination, adding 3$\times$2pt doubles the constraining power. Finally, when combined with primary CMB information, we obtain $w_0=-0.82^{+0.05}_{-0.05}$, $w_a=-0.63^{+0.21}_{-0.18}$, with a $3.0σ$ deviation. We find that including 3$\times$2pt in the previously studied SN + DESI BAO + CMB combination leaves the significance essentially unchanged ($3.2 σ$ to $3.0σ$) while improving the figure of merit by $\sim$10%. We systematically investigate the impact of leaving out each one of the probes and find that the significance of the deviation from a cosmological constant ranges from 2.3 to 3.2$σ$, with best-fit parameters consistently in the region $w_0 >-1$ and $w_a <0$. Excluding SN from the all data combination yields a $2.6σ$ departure from $Λ$CDM, providing a cross-check independent of supernova photometric calibration. These results support the weak preference for evolving dark energy reported by several recent cosmological analyses. By combining growth and geometric probes from a single survey, this work realizes the multi-probe dark energy program envisioned at the inception of DES.

Abbott, T. M.C. [Cerro-Tololo InterAmerican Obs.] ↗

DOE OSTI2026

Reliable statistics-based detection and investigation of anomalies in a SMART valve system

Reliable anomaly detection and diagnosis are critical for the safe operation of complex engineered systems. This study presents a unified framework that integrates statistical, model-based, and data-driven techniques for anomaly detection and investigation, demonstrated on SMART valve systems in hybrid energy applications. Four detection methods—mean deviation, seasonal extreme studentized deviate, ARIMA forecasting, and matrix profiling—were implemented and compared. Matrix profiling was particularly effective in revealing subtle deviations and hidden relationships among variables. Anomaly investigation was performed by analyzing variable-level and grouped signal profiles, with system topology incorporated to distinguish primary faults from propagated effects. Grouping signals by type enhanced interpretability, enabling accurate localization of anomalies across multi-dimensional datasets. Experimental results confirmed the framework's capability to consistently detect and isolate anomalies while providing actionable insights into system interdependencies. The proposed methodology offers a robust, interpretable, and scalable solution for condition monitoring, with potential applications in safety-critical domains such as nuclear energy, aerospace, and process industries.

ARIMA models↗

DOE OSTI2023

Characterization of Noise in CMOS Ring Oscillators at Cryogenic Temperatures

Allan deviation provides a means to characterize the time-dependence of noise in oscillators and potentially identify the source characteristics. Measurements on a 130 nm, 7-stage ring oscillator show that the Allan deviation declines from 300 K to 150 K as expected, but surprisingly increases from 150 K to 11 K. At low temperatures, the measured Allan deviation can be well fit using a few random telegraph noise (RTN) sources over the range of a few kilohertz to a few gigahertz. Further, the RTN characteristics evolve to reveal an enhanced role in low-frequency noise at lower temperatures.

47 OTHER INSTRUMENTATION↗

NASA NTRS1994

Sensor biases effect on the estimation algorithm for performance-seeking controllers

The performance-seeking-control algorithm (PSC) is designed to continuously optimize the performance of propulsion systems. The PSC uses a nominal model of the propulsion system and estimates, in flight, the engine deviation parameters (EDPs) characterizing engine deviations with respect to nominal conditions. In practice, the measurement biases (or model uncertainties) may prevent the estimated EDPs from reflecting the engine's actual off-nominal condition. This factor has a direct impact on the PSC scheme exacerbated by the open-loop character of the algorithm. An observability analysis shows that the biases cannot be estimated together with the EDPs. Moreover, biases and EDPs turn out to have equivalent effects on the measurements, leaving it undecided whether the estimated EDPs represent the actual engine deviation or whether they simply reflect the measurement biases. In this article, the effects produced by unknown measurement biases over the estimation algorithm are evaluated. This evaluation allows for identification of the most critical measurements for application of the PSC algorithm to an F100 engine.

F100 ENGINE↗

DOE OSTI2022

Side Chain Dynamics of Poly(norbornene)- g -Poly(propylene oxide) Bottlebrush Polymers

The segmental dynamics of the side chains of poly(norbornene)- g -poly(propylene oxide) (PNB- g -PPO) bottlebrush polymer in comparison to PPO is studied by quasi-elastic neutron scattering. Having experimental time and length scale information simultaneously allows to extract spatial information in addition to relaxation time. Tethering one end of the PPO side chain onto a stiff PNB backbone slows down the segmental relaxation, over the length and time scales investigated. The power law dependence of the relaxation time on the momentum transfer, Q , indicates a more heterogeneous relaxation pattern for the bottlebrush polymer, whereas the linear PPO has less deviations from a homogenous relaxation. Similar conclusions can be drawn from the time dependent mean square displacement, $\langle$ r 2 ( t )$\rangle$, and the non-Gaussian parameter, α 2 ( t ). Herein, the bottlebrush polymer shows a more restricted dynamics, whereas the linear PPO reaches $\langle$ r 2 ( t )$\rangle$∝ t 0.5 at the highest temperature. The deviations from Gaussian behavior are evident at the α 2 ( t ). Both samples show a decaying α 2 ( t ). Finally, the non-Gaussian parameter supports the results from the power law dependence of the relaxation times, with lower α 2 ( t ) values for PPO compared to those for PNB-g-PPO, pointing to less deviations from Gaussian behavior.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Measurement of the cross section for electroweak production of a Z boson, a photon and two jets in proton-proton collisions at $\sqrt{s} =$ 13 TeV and constraints on anomalous quartic couplings

A measurement is presented of the cross section for electroweak production of a Z boson and a photon in association with two jets (Zγjj) in proton-proton collisions. The Z boson candidates are selected through their decay into a pair of electrons or muons. The process of interest, electroweak Zγjj production, is isolated by selecting events with a large dijet mass and a large pseudorapidity gap between the two jets. The measurement is based on data collected at the CMS experiment at $ \sqrt{s} $ = 13 TeV, corresponding to an integrated luminosity of 35.9 fb$^{−1}$. The observed significance of the signal is 3.9 standard deviations, where a significance of 5.2 standard deviations is expected in the standard model. These results are combined with published results by CMS at $ \sqrt{s} $ = 8 TeV, which leads to observed and expected respective significances of 4.7 and 5.5 standard deviations. From the 13 TeV data, a value is obtained for the signal strength of electroweak Zγjj production and bounds are given on quartic vector boson interactions in the framework of dimension-eight effective field theory operators.[graphic not available: see fulltext]

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Measurement of the $t\bar{t} t\bar{t}$ production cross section in pp collisions at $ \sqrt{s} $ = 13 TeV with the ATLAS detector

A measurement of four-top-quark production using proton-proton collision data at a centre-of-mass energy of 13 TeV collected by the ATLAS detector at the Large Hadron Collider corresponding to an integrated luminosity of 139 fb -1 is presented. Events are selected if they contain a single lepton (electron or muon) or an opposite-sign lepton pair, in association with multiple jets. The events are categorised according to the number of jets and how likely these are to contain b-hadrons. A multivariate technique is then used to discriminate between signal and background events. The measured four-top-quark production cross section is found to be 26$^{+17}_{-15}$ fb, with a corresponding observed (expected) significance of 1.9 (1.0) standard deviations over the background-only hypothesis. The result is combined with the previous measurement performed by the ATLAS Collaboration in the multilepton final state. The combined four-top-quark production cross section is measured to be 24$^{+7}_{-6}$ fb, with a corresponding observed (expected) signal significance of 4.7 (2.6) standard deviations over the background-only predictions. It is consistent within 2.0 standard deviations with the Standard Model expectation of 12.0 ± 2.4 fb.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Investigating the Relationship Between Bolide Entry Angle and Apparent Direction of Infrasound Signal Arrivals

Infrasound sensing offers critical capabilities for detecting and geolocating bolide events globally. However, the observed back azimuths, directions from which infrasound signals arrive at stations, often differ from the theoretical expectations based on the bolide’s peak brightness location. For objects with shallow entry angles, which traverse longer atmospheric paths, acoustic energy may be emitted from multiple points along the trajectory, leading to substantial variability in back azimuth residuals. This study investigates how the entry angle of energetic bolides affects the back azimuth deviations, independent of extrinsic factors such as atmospheric propagation, station noise, and signal processing methodologies. A theoretical framework, the Bolide Infrasound Back-Azimuth EXplorer Model (BIBEX-M), was developed to compute predicted back azimuths solely from geometric considerations. The model quantifies how these residuals vary as a function of source-to-receiver distance, revealing that bolides entering at shallow angles, e.g., 10°, can produce average residuals of 20°, with deviations reaching up to 46° at distances below 1000 km, and remaining significant even at 5000 km (up to 8°). In contrast, bolides with steeper entry angles, e.g., > 60°, show smaller deviations, typically under 5° at 1000 km and diminishing to less than ~1° beyond 5000 km. These findings attest to the need for careful interpretation when evaluating signal detections and estimating bolide locations. This work is not only pertinent to bolides but also to other high-energy, extended-duration atmospheric phenomena such as space debris and reentry events, where similar geometric considerations can influence infrasound arrival directions.

Acoustics↗

DOE OSTI2025

Part distortion monitoring in additive manufacturing using machining

In additive manufacturing, accumulation of residual stresses can result in severe part distortion from the desired preform shape. Current methods for in-situ part distortion monitoring in additive manufacturing typically require expensive sensors, or capital equipment, and require time-consuming post-processing to understand the shape deviation. This paper presents an in-situ method, in the context of hybrid manufacturing, for part distortion detection using machining of additively manufactured parts. As a surrogate, three test artifacts were used to represent different distorted geometries. The tool axis positions from the machine tool controller and the cutting power were monitored during a facing operation. Cutting power data was used to detect the tool entry and exit in the workpiece using a novel approach with power standard deviation metric. The workpiece geometry and distorted configuration was subsequently predicted for positional and rotational deviations to within 2 mm accuracy using synchronized tool position data with cutting power. The proposed method can be used in a hybrid (additive and subtractive) machine tool to periodically check part distortion in the additive build. The method is applicable for any additive process and is low-cost and computationally inexpensive.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Techno-economic analysis of advanced small modular nuclear reactors

Here, small modular nuclear reactors (SMRs) represent a robust opportunity to develop low-carbon and reliable power with the potential to meet cost parity with conventional power systems. This study presents a detailed, bottom-up economic evaluation of a 12 × 77 MW e (924 MW e total) light-water SMR (LW-SMR) plant, a 4 × 262 MW e (1,048 MW e ) gas-cooled SMR (GC-SMR) plant, and a 5 × 200 MW e (1,000 MW e total) molten salt SMR (MS-SMR) plant. Cost estimates are derived from equipment costs, labor hours, material inputs, and process-engineering models. The advanced SMRs are compared to natural gas combined cycle plants with and without post-combustion carbon capture and a conventional large nuclear reactor. Overnight capital cost (OCC) and levelized cost of energy (LCOE) estimates are developed. The OCC of the LW-SMR, GC-SMR, and MS-SMR are found to be $\$4,844$/kW, $\$4,355$/kW, and $\$3,985$/kW respectively. The LCOE of the LW-SMR, GC-SMR, and MS-SMR are found to be $\$89.6$/MWh, $\$81.5$/MWh, and $\$80.6$/MWh respectively. A Monte Carlo analysis is performed, for which the OCC and construction time of the LW-SMR is found to have a lower mean and standard deviation than a conventional large reactor. The LW-SMR OCC is found to have a mean of $\$5,233$/kW with a standard deviation of $\$658$/kW and a 90 % probability of remaining between $\$4,254$/kW and $\$6,399$/kW, while the construction duration is found to have a mean of 4.5 years with a standard deviation of 0.8 years and a 90 % probability of remaining between 3.4 and 6.0 years. The economic impact of economies of scale, simplification, modularization, and construction time for SMRs are discussed. Additionally, policy implications for direct SMR capital subsidies and the impact of a carbon tax on natural gas emissions are explored.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

Condensation heat transfer and pressure drop of low-global warming potential refrigerants in smooth aluminum tubes

In this work we report condensation heat transfer coefficient and pressure drop data for R-454B, R-32, and R-454C. These refrigerants are low Global Warming Potential (GWP) replacements for R-410A. We obtained these data in a smooth 7.2 mm outer diameter (5.94 mm inner diameter) aluminum tube. In these experiments, the key variables were condensation temperature (35 °C ≤ T cond ≤ 50 °C) and refrigerant mass flux (151.1 ≤ G ≤ 347.3 kg m –2 s –1 ). In general, for similar nominal operating conditions, we measured the highest heat transfer coefficients for R-32 and the lowest for R-454C. Among the correlations we used to predict the heat transfer data for all three refrigerants, the one developed by Thome et al. (2003) predicted the data with the greatest accuracy. Deviations between experimental data and predictions of this correlation were 21.7% for R-454B (with SBG correction), 22.0% for R-32, and 18.8% for R-454C (with SBG correction). Experimental pressure gradient data were predicted by the Friedel correlation Friedel (1979) with deviations of 13.8% for R-454B, 9.8% for R-32, and 10.2% for R-454C. In this work, we also compare the predictions of the Cavallini et al. correlation Cavallini et al. (2006) for R-454B and R-32 in smooth copper versus smooth aluminum tubes. To improve this correlation’s predictive capability for aluminum tubes, we have proposed some modifications. After incorporating these modifications, the deviations between the correlation’s predictions and our data were 10% for R-454B (with SBG correction), 8.1% for R-32, and 5.2% for R-454C (with SBG correction). Furthermore, these data and findings will assist the HVAC and refrigeration industry to implement low-GWP refrigerants in condensers using smooth aluminum tubes.

Aluminum tube↗

DOE OSTI2023

A general predictive methodology for fuel-mixture properties up to supercritical conditions

A predictive thermodynamic model is utilized for the calculation of fuel properties of oxymethylene dimethyl ethers (OME 3–4 ), surrogates for gasoline, diesel and aviation fuel, as well as alcohol blends with gasoline and diesel. The alcohols used for these blends are methanol, ethanol, propanol, butanol and pentanol; their mixing ratio ranges from 10 to 50% by volume. The model is based on the Perturbed-Chain Statistical Association Fluid Theory (PC-SAFT) equation of state (EoS) and Vapor Liquid Equilibrium (VLE) calculations at constant temperature, density and composition. The model includes the association term, with the assumption of two association sites (2B scheme), to enable the modeling of alcohols. The pure-component parameters are estimated based on the Group Contribution (GC) method of various sources, as well as a parametrization model specifically designed for the case of OME 3–4 . The results of the computational model for the density, vapor pressure and distillation curves at various conditions, including high-pressure, high-temperature (HPHT), are compared to experimental and computational data available in the literature. In the cases where no measurements are available for the surrogates, experimental data for the corresponding target fuel are used, taking into consideration the inherent deviation in properties between real and surrogate fuel. Overall, the results are in good agreement with the data from the literature, with the average deviation not exceeding 12% for temperature (Kelvin) on the distillation curves, 10% for density and 46% for vapor pressure and the general trend being captured successfully. The use of different pure component parameter estimation techniques can further improve the prediction quality in the cases of OME3–4 and the aviation fuel surrogate, especially for the vapor pressure, leading to an average deviation lower than 18%. These results demonstrate the predictive capabilities of the model, which extend to a wide range of fuel types and pressure/temperature conditions. Through this investigation, the present work aims to establish the limits of applicability of this thermodynamic property prediction methodology.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2020

A deep insight on the coal ash-to-slag transformation behavior during the entrained flow gasification process

Recent research provided deep insight on the coal ash-to-slag transformation characterization during the entrained flow gasification process, with experimentation on a 40 kg/hr (dry basis) coal-fed opposed multi-burner (OMB) entrained flow gasifier and simulation via FactSage™ software. A general mechanism is presented to relate the gasifier design temperature, ash fluid temperature, and operating temperature with the degree of the slag polymerization. The change of the high temperature zone, the corresponding particle residence time in the high temperature zone, and syngas composition have obvious effects on the slag mineral transformation behavior. Mineral types formed on the wall of the gasifier chamber were mainly anorthite (CaAl 2 Si 2 O 8 ), aluminum oxide (Al 2 O 3 ), and calcium sulfide (CaS). These minerals transformed to anorthite and diopside (CaMgSi 2 O 6 ) at the slag hole zone, while the minerals at the lock hopper were anorthite, orthoclase (KAlSi 3 O 8 ), quartz (SiO 2 ), gypsum (CaSO 4 ), calcite (CaCO 3 ), and halite (NaCl). FactSage™ predicted minerals as anorthite, diopside, orthoclase, and albite (NaAlSi 3 O 8 ), etc., where the slag temperature was below the ash fluid temperature and when the ratios of CO/CO 2 and (CO + H 2 )/CO 2 were lower than 1.0 and 2.0, respectively. By simulation, residual carbon was found to be the dominant factor over syngas composition to cause mineral transformation, and this was verified experimentally. The Ca-based crystals, typically anorthite, was shifted to diopside, near the slag hole zone, and a linear relationship was found between the content ratios of diopside/(anorthite + diopside), CaO/SiO 2 , and (CaO + MgO)/SiO 2 . A dimensionless number, θ, was defined to characterize the changing chemical composition and the degree of slag polymerization, with temperature deviation from the design condition. Three zones of θ were identified and related to the deviation between the actual gasification condition from the design condition. Finally, a low slag polymerization degree corresponded with a higher temperature deviation between the actual condition and design condition, and this proved that increased residual carbon content and changing iron valence state increased the mineral types when the slag temperature was below the ash fluid temperature.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2025

Modeling distributed energy resource aggregations in security constrained unit commitment and economic dispatch

The Federal Energy Regulatory Commission (FERC) recently issued Order 2222, which requires all wholesale electricity markets in the US to allow distributed energy resources (DERs) to participate in the market as aggregated resources. These DER aggregations may be composed of many individual resources that are offered and dispatched by the market as a single entity. We present here a model of a distributed energy resource aggregator (DERA) that is scheduled by a market operator’s security constrained unit commitment (SCUC) and security constrained economic dispatch (SCED). The DERA model includes constraints for battery energy storage systems (BESSs), demand response resources (DRRs), and a simple distributed energy resource (DER). This paper describes a model for each resource type and presents two methods for the DERA to generate market offer curves: a profit-maximizing optimization to compute cost curves and a direct cost algorithm to determine dispatch costs for each resource and combine into cost curves. Once all participating DERAs are scheduled in SCUC/SCED, the model is then modified to dispatch individual DERs to maximize profit or minimize schedule deviation of the DERAs. A simulation of a representative day illustrates the DERA offers, the scheduled generation, and the DERA dispatch. Findings show the potential for unavoidable schedule deviations due to internal DER constraints and due to economic incentives to deviate from the SCUC/SCED schedules. This highlights the importance of DERA offer construction on market efficiency and system reliability. Novel aspects of our approach include: (1) We consider the asymmetry of price incentives impacting DERAs from the wholesale market compared to those impacting consumers from the retail market, as imposed by current regulations and laws. (2) We model aggregate consumer response through statistically parameterizable utility functions rather than a potentially impractical approach of modeling each individual consumer. (3) We show how to use the DERA operational dispatch model to create offers into the wholesale electricity market. (4) We show how DERAs may fail to meet their scheduled dispatch because the market offer format may not permit them to fully express their operational features such as intertemporal costs and constraints to the market.

aggregations↗

DOE OSTI2024

Quantifying replication through repeated analysis of UVM-A, a liquid reference material for cosmogenic 10 Be and 26 Al studies

In this study, in-situ produced cosmogenic nuclide sample preparation and analysis is a complicated, multi-step process with numerous possible sources of error, many of which have not been robustly quantified. Here, we use a liquid reference material (UVM-A) to test whether the analytic precision associated with individual isotopic ratios (n = 96) measured at two Accelerator Mass Spectrometry facilities (Lawrence Livermore National Laboratory and Purdue Rare Isotope Measurement Laboratory) accurately represents the range of values determined for aliquots prepared in a single laboratory over a period of five years. We find that 10 Be/ 9 Be ratios measured at the different accelerator facilities (n = 26 and n = 70) have statistically indistinguishable central tendencies, matching to within 0.1 %. Based on this large dataset, we suggest preliminary consensus values for UVM-A of 1.45 ± 0.06 x 10 -13 for 10 Be/ 9 Be (mean, one standard deviation, n = 96) and 4.47 ± 0.22 x 10 -13 for 26 Al/ 27 Al (mean, one standard deviation, n = 27). For both 10 Be/ 9 Be and 26 Al/ 27 Al, the relative standard deviations of the measured ratios are similar to the reported analytic uncertainties of sample measurement, suggesting that sample preparation introduced little if any additional scatter beyond the uncertainty of isotopic analysis. This dataset demonstrates that robust sample preparation and measurement can generate isotope ratio data reproducible at the level of counting statistics.

58 GEOSCIENCES↗

DOE OSTI2024

2-D soil zymography: Accounting for the spatial variation of pH

Soil zymography is commonly used to quantify spatial distribution of hydrolytic enzyme activities on soil and plant root surfaces. It is recommended to adjust pH in zymography substrates and calibration solutions with respect to soil/root pH. However, pH values may vary greatly within a few mm of plant rhizosphere, potentially altering the distribution of pH in zymography membranes. Despite the fact that the effect of pH on the calibration of zymography membranes is generally known, its potential impact on zymography results is unaccounted for in processing zymography images and calculations of enzyme activity. In this study we assessed the effect of pH variations on the persistency of the methylumbelliferone (MUF) calibration. The studied pH values ranged from 4.5 to 7.5. The MUF calibration curves greatly deviated from that at a reference pH of 6.5, with a marked nonlinear increase of deviation with greater membrane brightness. We suggest that the problem can be partially alleviated by reducing the membrane incubation time. However, such deviations suggest the need for a more comprehensive resolution via mapping pH and using pH-specific calibrations to process zymography images. Here, we developed a MATLAB code to implement a pixel-based correction of enzyme activity for pH in processing time-lapse zymography images.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Deriving Material Properties from Feedback Error Signals in Scanning Tunneling Microscopy

Imperfections in measurements, e.g., deviations and broadening, are not devoid of information; rather, they can reveal valuable physical properties and processes. Scanning tunneling microscopy leverages a negative feedback loop to regulate the tunneling current. Practically, current fluctuates around its set point, and such deviations are considered insignificant and ignored. Here, in this study, we investigate the information embedded in these deviations. In the constant-current mode with an active feedback loop, we observe an unexpected persistent DC current offset from its set point when the tunneling junction is periodically perturbed. We demonstrate both experimentally and theoretically that such error signals encode local tunneling barrier heights and the square of local differential conductance as a consequence of the interplay between rectification and active feedback compensation. We provide evidence on the generalizability of this phenomenology to other negative feedback systems. This new approach has the potential to broadly impact physical sciences by allowing rapid measurements without lock-in amplifiers.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Interface and Thermophysical Properties of R 32 Refrigerant

Driven by the urgent demand for efficient cooling in microelectronics and advanced thermal management systems, difluoromethane (R32/CH 2 F 2 ) has emerged as a promising candidate owing to its favorable thermophysical properties, including high heat transfer efficiency and low viscosity. While bulk properties such as density, viscosity, and thermal conductivity have been widely studied, interfacial properties, including surface tension and interfacial thickness, remain comparatively underexplored, despite their importance in phase-transition dynamics. Here, we perform molecular dynamics (MD) simulations from 180 to 300 K using an optimized transferable force field for fluoropropenes with enhanced electrostatics to assess both bulk and interfacial behavior of R32. Simulations reproduced density within ±2.1%, viscosity within 3.05%, and thermal conductivity within 7.41% of NIST reference data. Heat capacities (C p and C v ) were predicted within 5%. For interfacial properties, surface tension trends were reproduced within 13.58% deviation, and the vapor–liquid coexistence curve closely matched reference data, yielding a critical temperature of 345.7 K (1.6% deviation) and a critical density of 0.397 g/cm 3 (6.4% deviation). Importantly, the vapor–liquid interface exhibited pronounced temperature-dependent broadening across the 180–290 K range. This behavior correlates with increasing molecular kinetic energy, reduction in intermolecular cohesive interactions, and a progressive loss of preferential dipole alignment, which collectively enhance thermal fluctuation amplitudes at elevated temperatures. These validated results provide predictive molecular-level insights, particularly for interfacial properties that remain less characterized. By reducing property prediction errors in key parameters such as critical temperature, this work provides reliable inputs for heat-exchanger and system models. Such correlations can support optimized component sizing, improved performance, and reduced refrigerant charge. Beyond R32, the methodology offers a transferable framework for blended and next-generation low-GWP refrigerants, contributing to sustainable thermal management aligned with the 2027 EU F-Gas regulation and 2030 Kigali Amendment.

Fluids↗