Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Cases”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Near-wall model for compressible turbulent boundary layers based on an inverse velocity transformation

In this work, a near-wall model, which couples the inverse of a recently developed compressible velocity transformation (Griffin et al., Proc. Natl Acad. Sci., vol. 118, 2021, p. 34) and an algebraic temperature–velocity relation, is developed for high-speed turbulent boundary layers. As input, the model requires the mean flow state at one wall-normal height in the inner layer of the boundary layer and at the boundary-layer edge. As output, the model can predict mean temperature and velocity profiles across the entire inner layer, as well as the wall shear stress and heat flux. The model is tested in an a priori sense using a wide database of direct numerical simulation high-Mach-number turbulent channel flows, pipe flows and boundary layers (48 cases, with edge Mach numbers in the range 0.77–11, and semi-local friction Reynolds numbers in the range 170–5700). The present model is significantly more accurate than the classical ordinary differential equation (ODE) model for all cases tested. The model is deployed as a wall model for large-eddy simulations in channel flows with bulk Mach numbers in the range 0.7–4 and friction Reynolds numbers in the range 320–1800. When compared to the classical framework, in the a posteriori sense, the present method greatly improves the predicted heat flux, wall stress, and temperature and velocity profiles, especially in cases with strong heat transfer. In addition, the present model solves one ODE instead of two, and has a computational cost and implementation complexity similar to that of the commonly used ODE model.

42 ENGINEERING↗
DOE OSTI2021

Impact of particle loading and phase coupling on gas–solid flow dynamics: A case study of a two-phase, gas–solid flow in an annular pipe

The present study is devoted to a two-phase, gas–solid flow in an annular pipe (hollow cylinder) at an elevated pressure of 15 bars and moderate Reynolds number of circa 6000. The influence of the particle loading, the interaction between the phases, and turbulence dispersion on the flow dynamics is systematically studied by means of computational fluid dynamics simulations, employing the Ansys FLUENT commercial package. The cases with a particle volumetric fraction of 1.2% are referred to as “high particle loading,” and those with 0.13% are denoted as “low particle loading.” The following cases are investigated: (1) pure gas flow; (2) low particle loading two-phase flow with one-way coupling and with turbulence dispersion; (3) low particle loading two-phase flow with two-way coupling but without turbulence dispersion; (4) low particle loading two-phase flow with two-way coupling and with turbulence dispersion; (5) high particle loading two-phase flow with one-way coupling and with turbulence dispersion; (6) high particle loading two-phase flow with two-way coupling but without turbulence dispersion; and (7) high particle loading two-phase flow with two-way coupling and with turbulence dispersion. Additionally, the boundary layer is found to grow without fluctuations of the turbulent kinetic energy (TKE) for cases 1, 2, and 5. For case 4, the TKE fluctuations have been identified, although they appear to be less substantial than those in cases 6 and 7. The authors attribute the semi-chaotic nature of the TKE fluctuations to the particle loading and two-way coupling. In addition, the onset and development of the flow instability have been observed at a random axial distance in cases 4, 6, and 7. Such instability is also attributed to the two-way coupling with turbulence dispersion in the flow. It is concluded that the particle loading, one-way, or two-way coupling between the phases, and the turbulence dispersion models significantly influence the development of the flow dynamics with the same inlet and boundary conditions. Consequently, it is not a trivial question, which result a user should trust. The present computational results inspire to perform verification as well as experimental validation of the simulations, so the simulation results can subsequently be used with confidence for design analysis.

42 ENGINEERING↗
DOE OSTI2022

Design and execution of a verification, validation, and uncertainty quantification plan for a numerical model of left ventricular flow after LVAD implantation

Left ventricular assist devices (LVADs) are implantable pumps that act as a life support therapy for patients with severe heart failure. Despite improving the survival rate, LVAD therapy can carry major complications. Particularly, the flow distortion introduced by the LVAD in the left ventricle (LV) may induce thrombus formation. While previous works have used numerical models to study the impact of multiple variables in the intra-LV stagnation regions, a comprehensive validation analysis has never been executed. The main goal of this work is to present a model of the LV-LVAD system and to design and follow a verification, validation and uncertainty quantification (VVUQ) plan based on the ASME V&V40 and V&V20 standards to ensure credible predictions. The experiment used to validate the simulation is the SDSU cardiac simulator, a bench mock-up of the cardiovascular system that allows mimicking multiple operation conditions for the heart-LVAD system. The numerical model is based on Alya, the BSC’s in-house platform for numerical modelling. Alya solves the Navier-Stokes equation with an Arbitrary Lagrangian-Eulerian (ALE) formulation in a deformable ventricle and includes pressure-driven valves, a 0D Windkessel model for the arterial output and a LVAD boundary condition modeled through a dynamic pressure-flow performance curve. The designed VVUQ plan involves: (a) a risk analysis and the associated credibility goals; (b) a verification stage to ensure correctness in the numerical solution procedure; (c) a sensitivity analysis to quantify the impact of the inputs on the four quantities of interest (QoIs) (average aortic root flow $Q^{avg}_{Ao}$, maximum aortic root flow $Q^{max}_{Ao}$, average LVAD flow $Q^{avg}_{VAD}$, and maximum LVAD flow $Q^{max}_{VAD}$; (d) an uncertainty quantification using six validation experiments that include extreme operating conditions. Numerical code verification tests ensured correctness of the solution procedure and numerical calculation verification showed a grid convergence index (GCI)95% <3.3%. The total Sobol indices obtained during the sensitivity analysis demonstrated that the ejection fraction, the heart rate, and the pump performance curve coefficients are the most impactful inputs for the analysed QoIs. The Minkowski norm is used as validation metric for the uncertainty quantification. It shows that the midpoint cases have more accurate results when compared to the extreme cases. The total computational cost of the simulations was above 100 [core-years] executed in around three weeks time span in Marenostrum IV supercomputer. This work details a novel numerical model for the LV-LVAD system, that is supported by the design and execution of a VVUQ plan created following recognised international standards. We present a methodology demonstrating that stringent VVUQ according to ASME standards is feasible but computationally expensive.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Supo Thermal Model Development IV

This report describes the continuation of the Computational Fluid Dynamics (CFD) model of the Supo cooling system described in the reports, Supo Thermal Model Development, II, and III. In the previous reports, the solution natural convection of Supo was simulated using a 2D, axisymmetric model that altered the actual geometry. A 3D geometry of Supo will more accurately estimate the natural convection heat transfer coefficient (HTC) of the system since Supo contains several non-symmetric internal components that can only be modeled in 3D. This, in turn, will improve the natural convection HTC results for comparison purposes with past operational data and previous 2D CFD results. Unfortunately, the complexities of performing multiphase 3D CFD simulations has proven too large of a task for our current capabilities. Instead, this report will focus on the 3D model development for the 2.39 kW/L (31.3 kW) case and limitations and challenges associated with the 3D simulations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Divergence Reduction in Monte Carlo Neutron Transport with On-GPU Asynchronous Scheduling

While Monte Carlo Neutron Transport (MCNT) is near-embarrasingly parallel, the effectively unpredictable lifetime of neutrons can lead to divergence when MCNT is evaluated on GPUs. Divergence is the phenomenon of adjacent threads in a warp executing different control flow paths; on GPUS, it reduces performance because each work group may only execute one path at a time. The process of Thread Data Remapping (TDR) resolves these discrepancies by moving data across hardware such that data in the same warp will be processed through similar paths. A common issue among prior implementations of TDR is the synchronous nature of its remapping and processing cycles, which exhaustively sort data produced by prior processing passes and exhaustively evaluate the sorted data. In another work, we defined a method of remapping data through an asynchronous scheduler which allows for work to be stored in shared memory and deferred arbitrarily until that work is a viable option for low-divergence evaluation. This article surveys a wider set of cases, with the goal of characterizing performance trends across a more comprehensive set of parameters. These parameters include cross sections of scattering/capturing/fission, use of implicit capture, source neutron counts, simulation time spans, and tuned memory allocations. Across these cases, we have recorded minimum and average execution times, as well as a heuristically tuned near-optimal memory allocation size for both synchronous and asynchronous scheduling. Across the collected data, it is shown that the asynchronous method is faster and more memory efficient in the majority of cases, and that it requires less tuning to achieve competitive performance.

Computer Science↗
DOE OSTI2021

Tritium Diffusion in Fe-Al Aluminide Coating Bulk Phases

The theoretical investigations of interstitial tritium diffusion pathways in various Al-rich (>60% at. Al) iron aluminide coating phases found that in FeNiAl 5 tritium prefers to diffuse in between layers, while in Fe 4 Al 13 , tritium prefers to diffuse across layers. In the case of Fe 2 Al x phases, the simulations show that the variable Al content in the channels affects the diffusion behavior of tritium in the channels. For 100% of vacancy, tritium prefers to diffuse at the center of the channel, while for 0% Al vacancy, tritium prefers to diffuse along the channel’s wall. In the case where 50% Al vacancy are present in the channels, tritium can diffuse in a mixed way, with positions oscillating between the channel’s center and wall. However, the large energy barrier involved in this case suggests that tritium diffusion could be correlated with Al diffusion in the channels. While we found that at 600 K interstitial tritium diffuses faster in Fe 2 Al x phases (D Τ ≤10 -11 m 2 .s -1 ), then in Fe 4 Al 13 (D Τ ≈10 -12 m 2 .s -1 ), and finally in FeNiAl 5 (D Τ ≈10 -13 m 2 .s -1 ), we also find that interstitial tritium generally diffuses faster in Fe-Al coating phases than in γ-LiAlO 2 pellets (D Τ ≈10 -14 m 2 .s -1 ).

36 MATERIALS SCIENCE↗
DOE OSTI2024

Aerosol-induced closure of marine cloud cells: enhanced effects in the presence of precipitation

Abstract. The Weather Research Forecasting (WRF) version 4.3 model is configured within a Lagrangian framework to quantify the impact of aerosols on evolving cloud fields. Kilometer-scale simulations utilizing meteorological boundary conditions are based on 10 case study days offering diverse meteorology during the Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA). Measurements from aircraft, the ground-based Atmosphere Radiation Measurement (ARM) site at Graciosa Island in the Azores, and A-Train and geostationary satellites are utilized for validation, demonstrating good agreement with the WRF-simulated cloud and aerosol properties. Higher aerosol concentration leads to suppressed drizzle and increased cloud water content in all case study days. These changes lead to larger radiative cooling rates at cloud top, enhanced vertical velocity variance, and increased vertical and horizontal wind speed near the base of the lower-tropospheric inversion. As a result, marine cloud cell area expands, narrowing the gap between shallow clouds and increasing cloud optical thickness, liquid water content, and the top-of-atmosphere outgoing shortwave flux. While similar aerosol effects are observed in lightly to non-raining clouds, they tend to be smaller by comparison. These simulations show a relationship between cloud cell area expansion and the radiative adjustments caused by liquid water path and cloud fraction changes. The adjustments are positive and scale as 74 % and 51 %, respectively, relative to the Twomey effect. While higher-resolution large-eddy simulations may provide improved representation of cloud-top mixing processes, these results emphasize the importance of addressing mesoscale cloud-state transitions in the quantification of aerosol radiative forcing that cannot be attained from traditional climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Emulation of quantum correlations by classical dynamics in a spin-$\frac{1}{2}$ Heisenberg chain

We simulate the dynamical spin structure factor (DSSF) 𝒮⁡(𝑞,𝜔) of the spin-1/2 Heisenberg antiferromagnetic chain using classical simulations. By employing Landau-Lifshitz Dynamics, we emulate quantum correlations through temperature-dependent corrections, including rescaling of magnetic dipoles and renormalization of exchange interactions. Here, our results closely match Quantum Monte-Carlo calculations for 𝑘 B⁢ 𝑇/𝐽≳1, extending the applicability of classical dynamics to the challenging case of gapless excitations. At higher temperatures, our simulations comply with general predictions for uncorrelated paramagnetic fluctuations in the infinite temperature limit. Entanglement witnesses derived from the quantum-equivalent DSSF act as sensitive diagnostics for the quantum-to-classical crossover. Their reliability stems from their dependence on spectral features alone, enabling classical dynamics to emulate quantum thresholds without genuine entanglement. This framework also reproduces transverse spin correlations in finite magnetic fields, in agreement with quantum simulations. Together, our results establish quantum-corrected classical dynamics as a scalable and predictive tool for interpreting scattering experiments and exploring quantum correlations in strongly correlated spin systems.

Inelastic neutron scattering↗
DOE OSTI2022

High gas throughput SOLPS-ITER simulations extending the ITER database to strong detachment

Abstract SOLPS-ITER simulations performed for Q DT = 10, P SOL = 100 MW burning plasmas on ITER extend the existing database to high values of separatrix averaged neon impurity concentration (⟨ c Ne ⟩ ≈ 6%) and divertor neutral pressure (⟨ p div ⟩ > 25 Pa) in order to determine the heat flux mitigation capability of these scenarios and whether strongly detached states are accessible. In the existing database of ITER simulations, the level of detachment was limited to cases where the integral ion flux to the outer target was greater than 80% of the value at rollover, with the impurity radiation localized near the target. With the possibility of narrow heat flux channels and increased deposited power due to tile shaping, it is important to explore operation at a higher degree of detachment. Two series of simulations were explored to extend the database of SOLPS simulations. By increasing the deuterium and neon puff rates proportionally, the peak divertor energy flux ( q ⊥,max ) is decreased from 5 to 3 MW m −2 while ⟨ p div ⟩ increased from 11 to 27 Pa. By increasing only the neon puff, q ⊥, max can be reduced to <1MW m −2 while ⟨ p div ⟩ is maintained at ∼ 11 Pa. As the neon puff level is increased, the position of the impurity radiation peak is shifted towards the X-point. At the highest neon puff levels with steady-state solutions, the electron temperature is reduced below 1 eV across 50 cm of each divertor target. The new cases extend previously observed tight relationships in power and momentum loss factors to low electron temperature improving their utility for highly detached regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Experimental and computational investigations of ethane and ethylene kinetics with copper oxide particles for Chemical Looping Combustion

In this work, reaction pathways for the oxidation of methane, ethane, and ethylene with CuO was obtained by ReaxFF Molecular Dynamics (MD) simulations between temperatures of 1000 K and 2000 K. Experiments in a fixed-bed flow reactor were preformed with methane, ethane, and ethylene at temperatures ranging from 500 K to 1000 K with time-dependent species measurements from an Electron-Ionization Molecular Beam Mass Spectrometer (MBMS), and species validation with Gas Chromatography (GC) for detection of complete and intermediate combustion products. The MBMS and GC allow for the detection of oxygenated species and larger species produced from radical reformation. The simulation and experiment agree on the production of such species as CH 3 CHO, CH 2 O, CO, and H 2 O, which allow for the creation of simple C1 and C2 reaction pathways, which can be used in kinetic models of C2 species and larger fuels such as biofuels, which inherently depend on C1 and C2 kinetics and reaction pathway. The simulation and experiment disagree on the formation of C 2 H 2 , CH 3 OH, and CO 2 with large amounts of C 2 H 2 being measured in the ethylene oxidation simulations and CH 3 OH being formed in methane oxidation simulations, while neither species were experimentally found. In the case of CO 2 large amounts of CO 2 are rapidly produced in experiments with C2 fuels at 800 K, while little-to-no CO 2 was observed in simulations. This is believed to be resulting from the extremely short timescale of the simulations, preventing total oxidation of the fuel. Here, the differences in products produced between simulation and experiment allow for the potential to modify the ReaxFF potential functions to more accurately model the experimental products of Cu–H–O–C reaction kinetics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Impacts of Stochastic Mixing in Idealized Convection-Permitting Simulations of Squall Lines

This study investigates impacts of altering subgrid-scale mixing in “convection-permitting” kilometer-scale horizontal-grid-spacing (Δ h ) simulations by applying either constant or stochastic multiplicative factors to the horizontal mixing coefficients within the Weather Research and Forecasting Model. In quasi-idealized 1-km Δ h simulations of two observationally based squall-line cases, constant enhanced mixing produces larger updraft cores that are more dilute at upper levels, weakens the cold pool, rear-inflow jet, and front-to-rear flow of the squall line, and degrades the model’s effective resolution. Reducing mixing by a constant multiplicative factor has the opposite effect on all metrics. Completely turning off parameterized horizontal mixing produces bulk updraft statistics and squall-line mesoscale structure closest to an LES “benchmark” among all 1-km simulations, although the updraft cores are too undilute. The stochastic mixing scheme, which applies a multiplicative factor to the mixing coefficients that varies stochastically in time and space, is employed at 0.5-, 1-, and 2-km Δ h . It generally reduces midlevel vertical velocities and enhances upper-level vertical velocities compared to simulations using the standard mixing scheme, with more substantial impacts at 1- and 2-km Δ h compared to 0.5-km Δ h . Further, the stochastic scheme also increases updraft dilution to better agree with the LES for one case, but has less impact on the other case. Stochastic mixing acts to weaken the cold pool but without a significant impact on squall-line propagation. It also does not affect the model’s overall effective resolution unlike applying constant multiplicative factors to the mixing coefficients.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Effect of Ice Number Concentration on the Evolution of Boundary Layer Clouds During Arctic Marine Cold‐Air Outbreaks

Abstract Marine cold‐air outbreaks (MCAOs) are crucial for Arctic Ocean heat loss, featuring convective cloud rolls that transition into convection cells downstream. Understanding factors controlling this transformation is the key for improving MCAO cloud representation in climate models. This study employs large‐eddy simulations to investigate how cloud ice number concentrations () affect cloud evolution using a case from the Cold‐Air Outbreaks in the Marine Boundary Layer Experiment (COMBLE) campaign. The simulations, performed in a Lagrangian framework following an air mass trajectory, are driven by ERA5 reanalysis data. Initially, all simulations produce similar cloud patterns, but higher leads to earlier breakup of cloud rolls. Between 4 and 10 hr, surface precipitation rates are similar across simulations, but precipitation initiates earlier, and the cloud‐base precipitation rates are higher when is higher. The stronger precipitation evaporation leads to increased stability of the boundary layer and reduced intensity of vertical mixing between the surface and cloud layer. An increased sink of cloud layer moisture via precipitation and decreased source through diminished vertical transport result in earlier cloud breakup in higher conditions. Simulations with different sea surface temperatures (SST) indicate that this cloud breakup mechanism remains valid for MCAOs of different strengths, although the cloud organization is more sensitive to SST changes in low environments. This work highlights the importance of accurate representations of ice processes in simulating MCAO clouds and suggests the need for observational constraints of ice nucleating particles and over the mixed‐phase cloud regimes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Fully compressible simulations of waves and core convection in main-sequence stars

Context: Recent, nonlinear simulations of wave generation and propagation in full-star models have been carried out in the anelastic approximation using spectral methods. Although it makes long time steps possible, this approach excludes the physics of sound waves completely and requires rather high artificial viscosity and thermal diffusivity for numerical stability. A direct comparison with observations is thus limited. Aims: In this study, we explore the capabilities of our compressible multidimensional Seven-League Hydro (SLH) code to simulate stellar oscillations. Methods: We compare some fundamental properties of internal gravity and pressure waves in 2D SLH simulations to linear wave theory using two test cases: (1) an interval gravity wave packet in the Boussinesq limit and (2) a realistic 3 M ⊙ stellar model with a convective core and a radiative envelope. Oscillation properties of the stellar model are also discussed in the context of observations. Results: Our tests show that specialized low-Mach techniques are necessary when simulating oscillations in stellar interiors. Basic properties of internal gravity and pressure waves in our simulations are in good agreement with linear wave theory. As compared to anelastic simulations of the same stellar model, we can follow internal gravity waves of much lower frequencies. The temporal frequency spectra of velocity and temperature are flat and compatible with the observed spectra of massive stars. Conclusion: The low-Mach compressible approach to hydrodynamical simulations of stellar oscillations is promising. Our simulations are less dissipative and require less luminosity boosting than comparable spectral simulations. The fully-compressible approach allows for the coupling of gravity and pressure waves in the outer convective envelopes of evolved stars to be studied in the future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Analysis of Screening Current Effects in a Hybrid Nb$_{3}$Sn/REBCO Superconducting Accelerator Magnet Using a T - A Formulation

To explore the feasibility of using high-temperature superconducting (HTS) REBCO coated conductors in future accelerator magnets, two REBCO flat racetrack coils were fabricated using 4-mm wide EuBCO tapes at the High Energy Accelerator Research Organization (KEK). These coils were tested as an insert inside a Nb$_{3}$Sn common-coil dipole magnet, which provides a background field of up to $\sim$ 9.5T, at the Brookhaven National Laboratory (BNL). REBCO tapes offer exceptionally high critical current density under strong magnetic fields; however, they also exhibit significant magnetization due to screening currents, leading to magnetic field errors. Here, this study presents a 2D finite element model of screening current-induced fields (SCIF) in REBCO coils using the T-A formulation, along with the results obtained. Simulations were then performed for two KEK test cases: one where the REBCO conductors were oriented with the HTS tapes parallel to the background field, and another where the tapes were perpendicular to it. Since screening currents also influence the stress distribution and increase the peak stress in the coils, the mechanical effects of these currents were analyzed. The implications of these simulation and test results for the design of Nb$_{3}$Sn/REBCO superconducting accelerator magnets are discussed.

Accelerator magnets↗
DOE OSTI2024

Advanced Sensor Deployment for Distribution System State Estimation and Fault Identification

Distribution systems are currently facing steep operational challenges as a result of the rapidly increasing integration of renewables and other distributed energy resources (DERs) at both the primary and secondary circuit levels. Distribution utilities and system operators have traditionally had some visibility of their primary circuits using low-frequency supervisory control and data acquisition systems, and they have had very poor if not zero visibility of the secondary circuits where the presence of DERs is constantly increasing. Therefore, this paper presents simulation studies to demonstrate the benefits of an advanced, high-fidelity sensor technology, called as the Meta-Alert System (MAS), developed by Electrical Grid Monitoring, Ltd. (EGM), on the distribution grid. First, a reliable model of the EGM sensors is developed, and then two use cases, distribution system state estimation (DSSE) and fault identification are simulated to evaluate the performance of the MAS technology. Simulation results on the Electric Power Research Institute J1 feeder demonstrate that the MAS can effectively participate in system-level DSSE programs and can detect and locate faults faster than traditional distribution protection schemes.

distribution system↗
DOE OSTI2022

Orographic enhancement of rainfall over the Congo Basin

Abstract The Congo rainforest located in central equatorial Africa is an important, yet understudied part of the globe surrounded by complex orographic features. A primitive understanding of precipitation processes such as mesoscale convective dynamics magnifies uncertainties in the future climate projections of the hydrological cycle over the Congo. Furthermore, the effects of orography, which is an important forcing for convection and precipitation, are poorly resolved by climate models, and ill‐conceptualized over the Congo. To address this knowledge gap, perturbed orographic forcing experiments are conducted using the Weather Research and Forecasting (WRF) mesoscale numerical model in a high‐resolution convection‐permitting model setup. The model simulated selected dates in November 2014. The thunderstorms and rainfall simulated in the control run for the case study analyzed in this article compared reasonably well to satellite‐derived brightness temperature and rainfall data. The results from this case study show that the dynamical impact of increasing the height of the East African Highlands is the blocking of the lower‐tropospheric tropical easterlies. This weakening of the lower‐tropospheric zonal winds increases the windshear over the Congo Basin resulting in slower propagating, more intense mesoscale convective systems with enhanced rainfall.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Measurement of temperature of a dusty plasma from its configuration

A new concept called “configurational temperature” is introduced in the context of dusty plasma, where the temperature of the dust particles submerged in the plasma can be measured directly from the positional information of the individual dust particles and the interaction potential between the dust grains. This method does not require the velocity information of individual particles, which is a key parameter to measure the dust temperature in the conventional method. The technique is initially tested using two-dimensional (2D) OpenMP parallel molecular dynamics and Monte Carlo simulation and then compared with the temperature evaluated from experimental data. The experiment have been carried out in the Dusty Plasma Experimental (DPEx) device, where a 2D stationary plasma crystal of melamine formaldehyde particles is formed in the cathode sheath of a DC glow discharge argon plasma. The kinetic temperature of the dust is calculated using the standard particle image velocimetry technique at different pressures. An extended simulation result for the three-dimensional case is also presented, which can be employed for the temperature measurement of a three-dimensional dust crystal in laboratory devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Testing the parametric model for self-interacting dark matter using matched halos in cosmological simulations

Here, we systemically evaluate the performance of the self-interacting dark matter (SIDM) halo model proposed in Ref.Yang et al. (2023) with matched halos from high-resolution cosmological CDM and SIDM simulations. The model incorporates SIDM effects along mass evolution histories of CDM halos and it is applicable to both isolated halos and subhalos. We focus on the accuracy of the model in predicting halo density profiles at z = 0 and the evolution of maximum circular velocity. We find the model predictions agree with the simulations within 10%–50% for most of the simulated (sub)halos, 50%–100% for extreme cases. This indicates that the model effectively captures the gravothermal evolution of the halos with very strong, velocity-dependent self-interactions. For an example application, we apply the model to study the impact of various SIDM scenarios on strong lensing perturber systems, demonstrating its utility in predicting SIDM effects for small-scale structure analyses. Our findings confirm that the model is an effective tool for mapping CDM halos into their SIDM counterparts.

Yang, Daneng [Chinese Academy of Sciences, Nanjing↗