Search NASA⌕ Search

SEARCH · Search NASA

Results for “VELOCITY ERROR”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Dynamic population balance in molecular-level simulations of hypersonic flows

This report summarizes the work towards developing stochastic weighted particle methods (SWPM) for future application in hypersonic flows. Extensive changes to Sandia’s direct simulation Monte Carlo (DSMC) solver, SPARTA (Stochastic Particle Real Time Analyzer), were made to enable the necessary particle splitting and reduction capabilities for SWPM. The results from one-dimensional Couette and Fourier flows suggest that SWPM can reproduce the correct transport for a large range of Knudsen numbers with adequate accuracy. The associated velocity and temperature profiles are in good agreement with DSMC. An issue with particle placement during particle number reduction, is identified, to which, a simple but effective solution based on minimizing the center of mass error is proposed. High Mach wheel flows are simulated using the SWPM and DSMC methods. SWPM is capable of providing nearly an order of magnitude increase in efficiency over DSMC while retaining high accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A combined ensemble-volume average homogenization method for lattice structures with defects under dynamic and static loading

In the study of lattices structures, both experiments and numerical simulations are often conducted with small samples. Using combined ensemble and volume averaging, this work introduces a method to extract a macroscopic constitutive response of a lattice material from numerical simulations performed in periodic domains. The domain size needed to obtain statistically accurate results is investigated. Similar to molecular dynamics, the concept of the virial stress is introduced after homogenized equations are derived using the ensemble averaging method. Under static conditions, the virial stress is shown to agree with the volume averaged solid stress. Using the homogenization method, constitutive relations for this stress can be obtained from systems with uniform strains. Application of such obtained constitutive relations to more general cases results in an error proportional to the square of the ratio between the lattice length scale and the macroscopic length scale. Taking advantage of this property, numerical simulations are performed in systems with a uniform gradient of the average velocity. The volume average method is then used to accelerate convergence when studying lattices with defects. To avoid the artificial numerical time scale from the size of a representative volume element divided by the wave speed, a numerical scheme is developed to enforce a spatially uniform velocity gradient within the computational domain while allowing fluctuations of the velocity or displacement to develop naturally. To account for probability distribution of lattice defects, the stress is calculated as the ensemble-volume averaged value. For dynamic systems, energy dissipation properties are also studied.

36 MATERIALS SCIENCE↗
DOE OSTI2020

An LES-based airborne Doppler lidar simulator and its application to wind profiling in inhomogeneous flow conditions

Wind profiling by Doppler lidar is common practice and highly useful in a wide range of applications. Airborne Doppler lidar can provide additional insights relative to ground-based systems by allowing for spatially distributed and targeted measurements. Providing a link between theory and measurement, a first large eddy simulation (LES)-based airborne Doppler lidar simulator (ADLS) has been developed. Simulated measurements are conducted based on LES wind fields, considering the coordinate and geometric transformations applicable to real-world measurements. The ADLS provides added value as the input truth used to create the measurements is known exactly, which is nearly impossible in real-world situations. Thus, valuable insight can be gained into measurement system characteristics as well as retrieval strategies. As an example application, airborne Doppler lidar wind profiling is investigated using the ADLS. For commonly used airborne velocity azimuth display (AVAD) techniques, flow homogeneity is assumed throughout the retrieval volume, a condition which is violated in turbulent boundary layer flow. Assuming an ideal measurement system, the ADLS allows to isolate and evaluate the error in wind profiling which occurs due to the violation of the flow homogeneity assumption. Overall, the ADLS demonstrates that wind profiling is possible in turbulent wind field conditions with reasonable errors (root mean squared error of 0.36 m s –1 for wind speed when using a commonly used system setup and retrieval strategy for the conditions investigated). Nevertheless, flow inhomogeneity, e.g., due to boundary layer turbulence, can cause an important contribution to wind profiling error and is non-negligible. Results suggest that airborne Doppler lidar wind profiling at low wind speeds (<5 m s –1 ) can be biased, if conducted in regions of inhomogeneous flow conditions.

17 WIND ENERGY↗
DOE OSTI2020

Scaling Turbulent Combustion Fields in Explosions

We considered the topic of explosions from spherical high-explosive (HE) charges. We studied how the turbulent combustion fields scale. On the basis of theories of dimensional analysis by Bridgman and similarity theories of Sedov and Barenblatt, we found that all fields scaled with the explosion length scale r0. This included the blast wave, the mean and root mean squared (RMS) profiles of thermodynamic variables, combustion variables, velocities, vorticity, and turbulent Reynolds stresses. This was a consequence of the formulation of the problem and our numerical method, which both satisfied the similarity conditions of Sedov. We performed numerical simulations of 1 g charges and 1 kg charges; the solutions were identical (within roundoff error) when plotted in scaled variables. We also explored scaling laws related to three-phase pyrotechnic explosions. We show that although the scaling formally broke down, the fireball still essentially scaled with the explosion length scale r0. However, the discrete Lagrange particles (DLP) (phase 2) and the heterogeneous continuum model (HCM) of the DLP wakes (phase 3) did not scale with r0, and mean and RMS profiles could differ by a factor of 10 in some regions. This was because the DLP particles and wakes introduced an additional scale that broke the similarity conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Online Learning of Effective Turbine Wind Speed in Wind Farms

To develop better wind farm controllers that can meet more complex objectives, methods of modeling the wind turbine wakes at low computational expense are needed. Gaussian process (GP) regression offers a computationally inexpensive framework for learning complex functions from noisy measurements with very few datapoints. In this work, an online learning approach is presented to learn the rotor-averaged wind velocity at downstream wind turbines with GPs, using the available datastream of wind field measurements and wind turbine control set-points. This framework can readily be integrated into model-based controls methods because the model a) is updated online at low computational expense, b) assumes a mathematically favorable Gaussian form, and c) explicitly quantifies the stochastic nature of the wake field so that the trade-off between exploration and exploitation, and the uncertainty in the prediction, can be utilized. We show that a GP-learned model can match true values with errors within 0.5% on average, with as few as 5 training data points.

Gaussian process↗
DOE OSTI2020

EXPRES. I. HD 3651 as an Ideal RV Benchmark

The next generation of exoplanet-hunting spectrographs should deliver up to an order of magnitude improvement in radial velocity (RV) precision over the standard 1 m s{sup −1} state-of-the-art spectrographs. This advance is critical for enabling the detection of Earth-mass planets around Sun-like stars. New calibration techniques such as laser frequency combs and stabilized etalons ensure that the instrumental stability is well characterized. However, additional sources of error include stellar noise, undetected short-period planets, and telluric contamination. To understand and ultimately mitigate error sources, the contributing terms in the error budget must be isolated to the greatest extent possible. Here, we introduce a new high-cadence RV program, the Extreme Precision Spectrograph (EXPRES) 100 Earths Survey, which aims to identify rocky planets around bright, nearby G and K dwarfs. We also present a benchmark case: the 62 day orbit of a Saturn-mass planet orbiting the chromospherically quiet star, HD 3651. The combination of high eccentricity (0.6) and a moderately long orbital period ensures significant dynamical clearing of any inner planets. Our Keplerian model for this planetary orbit has a residual rms of 58 cm s{sup −1} over a ∼6 month time baseline. By eliminating significant contributors to the RV error budget, HD 3651 serves as a standard for evaluating the long-term precision of extreme precision RV programs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

On the Use of Field RR Lyrae as Galactic Probes. V. Optical and Radial Velocity Curve Templates

We collected the largest spectroscopic catalog of RR Lyrae (RRLs) including ≈20,000 high-, medium-, and low-resolution spectra for ≈10,000 RRLs. We provide the analytical forms of radial velocity curve (RVC) templates. These were built using 36 RRLs (31 fundamental—split into three period bins—and five first-overtone pulsators) with well-sampled RVCs based on three groups of metallic lines (Fe, Mg, Na) and four Balmer lines (H$α$, H$β$, H$γ$, H$δ$). We tackled the long-standing problem of the reference epoch to anchor light-curve and RVC templates. For the V-band, we found that the residuals of the templates anchored to the phase of the mean magnitude along the rising branch are ~35% to ~45% smaller than those anchored to the phase of maximum light. For the RVC, we used two independent reference epochs for metallic and Balmer lines and we verified that the residuals of the RVC templates anchored to the phase of mean RV are from 30% (metallic lines) up to 45% (Balmer lines) smaller than those anchored to the phase of minimum RV. We validated our RVC templates by using both the single-point and the three phase point approaches. We found that barycentric velocities based on our RVC templates are two to three times more accurate than those available in the literature. We applied the current RVC templates to Balmer lines RVs of RRLs in the globular NGC 3201 collected with MUSE at VLT. We found the cluster barycentric RV of $V$$γ$ = 496.89 ± 8.37(error) ± 3.43 (standard deviation) km s -1 , which agrees well with literature estimates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Redshift-space streaming velocity effects on the Lyman- α forest baryon acoustic oscillation scale

The baryon acoustic oscillation (BAO) scale acts as a standard ruler for measuring cosmological distances and has therefore emerged as a leading probe of cosmic expansion history. However, any physical effect that alters the length of the ruler can lead to a bias in our determination of distance and expansion rate. One of these physical effects is the streaming velocity, the relative velocity between baryons and dark matter in the early Universe, which couples to the BAO scale due to their common origin in acoustic waves at recombination. In this work, we investigate the impact of streaming velocity on the BAO feature of the Lyman- α forest auto-power spectrum, one of the main tracers being used by the recently commissioned Dark Energy Spectroscopic Instrument (DESI). To do this, we develop a perturbative model for Lyman- α flux fluctuations which is complete to second order for a certain set of fields, and applicable to any redshift-space tracer of structure since it is based only on symmetry considerations. We find that there are 8 biasing coefficients through second order. We find streaming velocity-induced shifts in the BAO scale of 0.081%–0.149% (transverse direction) and 0.053%–0.058% (radial direction), depending on the model for the biasing coefficients used. These are smaller than, but not negligible compared to, the DESI Lyman- α BAO error budget, which is 0.46% on the overall scale. Overall, the sensitivity of these results to our choice of bias parameters underscores the need for future work to measure the higher-order biasing coefficients from simulations, especially for future experiments beyond DESI.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Notes on Regression Analysis for Radar Parameter Estimation

A fundamental task of radar, beyond merely detecting a target, is to estimate some parameters associated with it. For example, this might include range, direction, velocity, etc. In any case, multiple measurements, often noisy, need to be processed to yield a ‘best estimate’ of the parameter. A common mathematical method for doing so is called “Regression” analysis. The goal is to minimize the expected squared error in the estimate. Even when alternate algorithms are considered, the least squared-error regression analysis is the benchmark against which alternatives are compared.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

An Eulerian multimaterial framework for simulating high-explosive aquarium tests

Aquarium tests of cylindrical high-explosive charges provide optical data of the detonation front velocity and shape, propagation of the shock wave in the surrounding water, and expansion rates of the detonation products behind the front. Data from aquarium experiments is often used for calibration of reactive burn models based on phenomenological equations of state (EOS) and reaction rate laws. This paper presents a multimaterial numerical modeling framework to solve the 2D axisymmetric reactive Euler equations for high-explosive aquarium tests, in particular for ammonium nitrate - fuel oil (ANFO) explosives. An extension of the Ghost Fluid Method (GFM) is used to handle the dynamic material interfaces for the ANFO explosion products, the charge-confining material (polymethyl methacrylate PMMA), and the surrounding water. This study analyzes the sensitivity of calculations (both computational efficiency and numerical accuracy) to different algorithms for the material interface models including the original GFM versus Riemann solver-based strategies. A novel method for defining the left and right states in the interfacial Riemann problem eliminates the need for sorting or nodal interpolation during the projection along the material interface. Numerical tests indicate that populating the interface node values using the Riemann solution mitigate the overheating error observed in steady-state calculations. Solution convergence and computational efficiency are explored as a function of the spatial and temporal order of the schemes. Results from the computational model with analytical equations of state and fitted reaction rate parameters show very good quantitative agreement with experimentally observed detonation front velocity, reaction products expansion, and shock wave propagation in the surrounding water for a cylindrical ANFO charge encased in PMMA. Finally, the proposed modeling framework, in conjunction with experimental tests, provides a reliable tool to assess equations of state and reaction rate expressions for reactive burn models of confined high explosives.

42 ENGINEERING↗
DOE OSTI2022

Estimation of the Surface Fluxes for Heat and Momentum in Unstable Conditions with Machine Learning and Similarity Approaches for the LAFE Data Set

Abstract Measurements of three flux towers operated during the land atmosphere feedback experiment (LAFE) are used to investigate relationships between surface fluxes and variables of the land–atmosphere system. We study these relations by means of two machine learning (ML) techniques: multilayer perceptrons (MLP) and extreme gradient boosting (XGB). We compare their flux derivation performance with Monin–Obukhov similarity theory (MOST) and a similarity relationship using the bulk Richardson number (BRN). The ML approaches outperform MOST and BRN. Best agreement with the observations is achieved for the friction velocity. For the sensible heat flux and even more so for the latent heat flux, MOST and BRN deviate from the observations while MLP and XGB yield more accurate predictions. Using MOST and BRN for latent heat flux, the root mean square errors (RMSE) are 107 Wm $$^{-2}$$ - 2 and 121 Wm $$^{-2}$$ - 2 , respectively, as well as the intercepts of the regression lines are $$\approx 110$$ ≈ 110 Wm $$^{-2}$$ - 2 . For the ML methods, the RMSEs reduce to 31 Wm $$^{-2}$$ - 2 for MLP and 33 Wm $$^{-2}$$ - 2 for XGB as well as the intercepts to just 4 Wm $$^{-2}$$ - 2 for MLP and $$-1$$ - 1 Wm $$^{-2}$$ - 2 for XGB with slopes of the regression lines close to 1, respectively. These results indicate significant deficiencies of MOST and BRN, particularly for the derivation of the latent heat flux. In fact, in contrast to the established theories, feature importance weighting demonstrates that the ML methods base their improved derivations on net radiation, the incoming and outgoing shortwave radiations, the air temperature gradient, and the available water contents, but not on the water vapor gradient. The results imply that further studies of surface fluxes and other turbulent variables with ML techniques provide great promise for deriving advanced flux parameterizations and their implementation in land–atmosphere system models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Interpreting Sunyaev–Zel’dovich observations with MillenniumTNG: mass and environment scaling relations

ABSTRACT Sunyaev–Zel’dovich (SZ) measurements can dramatically improve our understanding of the intergalactic medium and the role of feedback processes in galaxy formation, allowing us to calibrate important astrophysical systematics in cosmological constraints from weak lensing galaxy clustering surveys. However, the signal is only measured in a two-dimensional projection, and its correct interpretation relies on understanding the connection between observable quantities and the underlying intrinsic properties of the gas, in addition to the relation between the gas and the underlying matter distribution. One way to address these challenges is through the use of hydrodynamical simulations such as the high-resolution, large-volume MillenniumTNG suite. We find that measurements of the optical depth, τ, and the Compton-y parameter, Y, receive large line-of-sight contributions that can be removed effectively by applying a compensated aperture photometry filter. In contrast with other τ probes (e.g. X-rays and fast radio bursts), the kinematic SZ-inferred τ receives most of its signal from a confined cylindrical region around the halo due to the velocity decorrelation along the line of sight. Additionally, we perform fits to the Y–M and τ–M scaling relations and report best-fitting parameters adopting the smoothly broken power law formalism. We note that subgrid physics modelling can broaden the error bar on these by 30 per cent for intermediate-mass haloes (${\sim }10^{13} \, {\rm M}_{\odot }$). The scatter of the scaling relations can be captured by an intrinsic dependence on concentration and an extrinsic dependence on tidal shear. Finally, we comment on the effect of using galaxies rather than haloes in observations, which can bias the inferred profiles by ∼20 per cent for L* galaxies.

Hadzhiyska, Boryana (ORCID:0000000223123121)↗
DOE OSTI2022

Eddy current flow meter model validation with a moving solid rod*

We report that the eddy current flow meter (ECFM) has been used to measure velocities and temperatures of conductive flows, such as liquid metal flows in a nuclear fission reactor. The goal of this paper is to develop a finite element electromagnetic model that can characterize the ECFM sensor performance and validate this finite element model with detailed velocity measurements of a controlled, well-characterized moving conductive solid rod. Both measurements and modeling were performed for various parameters that are important for ECFM performance such as rod velocity, rod material, ECFM sensor coil length, number of sensor coils, applied alternating current (AC) current amplitude, and applied ECFM AC frequency. For all parametric scans, the measurement and modeling agree well in both magnitude and trend. The normalized root-mean-square error between measurement and modeling is less than 10% for all cases. These results suggest that electromagnetic modeling could eventually be used to cost-effectively design future ECFM sensors in arbitrary geometry for more challenging applications such as liquid metal nuclear fission reactors.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

High-fidelity retrieval from instantaneous line-of-sight returns of nacelle-mounted lidar including supervised machine learning

Abstract. Wind turbine applications that leverage nacelle-mounted Doppler lidar are hampered by several sources of uncertainty in the lidar measurement, affecting both bias and random errors. Two problems encountered especially for nacelle-mounted lidar are solid interference due to intersection of the line of sight with solid objects behind, within, or in front of the measurement volume and spectral noise due primarily to limited photon capture. These two uncertainties, especially that due to solid interference, can be reduced with high-fidelity retrieval techniques (i.e., including both quality assurance/quality control and subsequent parameter estimation). Our work compares three such techniques, including conventional thresholding, advanced filtering, and a novel application of supervised machine learning with ensemble neural networks, based on their ability to reduce uncertainty introduced by the two observed nonideal spectral features while keeping data availability high. The approach leverages data from a field experiment involving a continuous-wave (CW) SpinnerLidar from the Technical University of Denmark (DTU) that provided scans of a wide range of flows both unwaked and waked by a field turbine. Independent measurements from an adjacent meteorological tower within the sampling volume permit experimental validation of the instantaneous velocity uncertainty remaining after retrieval that stems from solid interference and strong spectral noise, which is a validation that has not been performed previously. All three methods perform similarly for non-interfered returns, but the advanced filtering and machine learning techniques perform better when solid interference is present, which allows them to produce overall standard deviations of error between 0.2 and 0.3 m s−1, or a 1 %–22 % improvement versus the conventional thresholding technique, over the rotor height for the unwaked cases. Between the two improved techniques, the advanced filtering produces 3.5 % higher overall data availability, while the machine learning offers a faster runtime (i.e., ∼ 1 s to evaluate) that is therefore more commensurate with the requirements of real-time turbine control. The retrieval techniques are described in terms of application to CW lidar, though they are also relevant to pulsed lidar. Previous work by the authors (Brown and Herges, 2020) explored a novel attempt to quantify uncertainty in the output of a high-fidelity lidar retrieval technique using simulated lidar returns; this article provides true uncertainty quantification versus independent measurement and does so for three techniques rather than one.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Blood flow imaging by optimal matching of computational fluid dynamics to 4D-flow data

Three-dimensional, time-resolved blood flow measurement (4D-flow) is a powerful research and clinical tool, but improved resolution and scan times are needed. Therefore, this study aims to (1) present a postprocessing framework for optimization-driven simulation-based flow imaging, called 4D-flow High-resolution Imaging with a priori Knowledge Incorporating the Navier-Stokes equations and the discontinuous Galerkin method (4D-flow HIKING), (2) investigate the framework in synthetic tests, (3) perform phantom validation using laser particle imaging velocimetry, and (4) demonstrate the use of the framework in vivo. An optimizing computational fluid dynamics solver including adjoint-based optimization was developed to fit computational fluid dynamics solutions to 4D-flow data. Synthetic tests were performed in 2D, and phantom validation was performed with pulsatile flow. Reference velocity data were acquired using particle imaging velocimetry, and 4D-flow data were acquired at 1.5 T. In vivo testing was performed on intracranial arteries in a healthy volunteer at 7 T, with 2D flow as the reference. Results Synthetic tests showed low error (0.4%-0.7%). Phantom validation showed improved agreement with laser particle imaging velocimetry compared with input 4D-flow in the horizontal (mean -0.05 vs -1.11 cm/s, P < .001; SD 1.86 vs 4.26 cm/s, P < .001) and vertical directions (mean 0.05 vs -0.04 cm/s, P = .29; SD 1.36 vs 3.95 cm/s, P < .001). In vivo data show a reduction in flow rate error from 14% to 3.5%. Phantom and in vivo results from 4D-flow HIKING show promise for future applications with higher resolution, shorter scan times, and accurate quantification of physiological parameters.

4D-flow MRI↗
DOE OSTI2022

Mind-the-gap - Part III: Doppler velocity measurements from space

Convective motions and hydrometeors microphysical properties are highly sought-after parameters for evaluating atmospheric numerical models. With most of the Earth's surface covered by water, space-borne Doppler radars are ideal for acquiring such measurements at a global scale. While these systems have proven to be useful tools for retrieving cloud microphysical and dynamical properties from the ground, their adequacy, and specific requirements for spaceborne operation still need to be evaluated. Comprehensive forward simulations enable us to assess the advantages and drawbacks of six different Doppler radar architectures currently planned or under consideration by space agencies for the study of cloud dynamics. Radar performance is examined against state-of-the-art numerical model simulations of well-characterized shallow and deep, continental, and oceanic convective cases. Mean Doppler velocity (MDV) measurements collected at multiple frequencies (13, 35 and 94 GHz) provide complementary information in deep convective cloud systems. The high penetration capability of the 13-GHz radar enables to obtain a complete, albeit horizontally under-sampled, view of deep convective storms. The smaller instantaneous field of view (IFOV) of the 35-GHz radar captures more precise information about the location and size of convective updrafts above 5-8 km height of most systems which was determined is the portion of storms where the mass flux peak is typically located. Finally, the lower mean Doppler velocity uncertainty of displaced phase center antenna (DPCA) radars makes them an ideal system for studying microphysics in shallow convection and frontal systems, as well as ice and mixed-phase clouds. It is demonstrated that a 94-GHz DCPA system can achieve retrieval errors as low as 0.05-0.15 mm for raindrop volume-weighted mean diameter and 25% for rime fraction (for a -10 dBZ echo).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Heavy-Duty Nonroad Material Handler Electrification Part 1: Real-World Drive Cycle Development

Knowing a detailed operating cycle is critical for developing and testing equipment. Operating cycles can be separated by two clear distinctions: (1) regulatory or non-regulatory and (2) application at the engine-only or full machine level. The Environmental Protection Agency’s (EPA) Nonroad Transient Cycle (NRTC) may be a good representation of engine use in many types of equipment, but there is a gap in standardized and validated drive cycles specifically for nonroad material handlers. Lacking a standardized drive cycle makes it difficult to accurately benchmark machine performance and validate new powertrain technologies. The objective of this investigation is to illustrate the development of a custom drive cycle augmented with real-world customer use data that serves multiple purposes: (1) understand the range of operation and utilization that formulated inputs for electrified architecture analysis and (2) develop a repetitive and consistent maneuver to establish baseline energy consumption enabling equivalent comparison to future electrified prototype builds. This article presents a solution specifically for a 23-ton nonroad material handler in which material handling, machine transport, and extended idle were homologated to form representative short cycles defined by machine velocity and hydraulic cylinder position. The most intensive material handling short cycles had a load factor of 40% and an average fuel rate of 16 L/h. Combined with a visual aid, the short cycles exhibited low variability, having less than 5% root mean square (RMS) error in lift and reach position with respect to the average. The machine’s performance on these short cycles at the Advanced Power Systems Research Center (APSRC) was compared to results from two real-world customer locations operating the instrumented test machine in a cyclical manner, and for similar ground conditions were found to be comparable in fuel consumption.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Modeling athermal phonons in novel materials using the G4CMP simulation toolkit

Understanding phonon and charge propagation in superconducting devices plays an important role in both performing low-threshold dark matter searches and limiting correlated errors in superconducting qubits. The Geant4 Condensed Matter Physics (G4CMP) package, originally developed for the Cryogenic Dark Matter Search (CDMS) experiment, models charge and phonon transport within silicon and germanium detectors and has been validated by experimental measurements of phonon caustics, mean charge-carrier drift velocities, and heat pulse propagation times. Here, in this work, we present a concise framework for expanding the capabilities for phonon transport to a number of other novel substrate materials of interest to the dark matter and quantum computing communities, including sapphire (Al 2 O 3 ), gallium arsenide (GaAs), lithium fluoride (LiF), calcium tungstate (CaWO 4 ), and calcium fluoride (CaF 2 ). We demonstrate the use of this framework in generating phonon transport properties of these materials and compare these properties with experimentally-determined values where available.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗