Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “quantum probability”

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

Angle-dependent strong-field ionization and fragmentation of carbon dioxide measured using rotational wave packets

In this work, we experimentally study the angle-dependent single ionization of carbon dioxide (CO 2 ) by linearly and circularly polarized pulses. The angle dependence of the ionization probability by linearly polarized pulses extracted from time-domain measurements on an impulsively excited rotational wave packet is compared with data obtained from a direct angle-scan measurement. The results from the measurement with linear and circular polarization are consistent with the adiabatic ionization approximation. Further, we extend the time-domain method to extract the dependence of the asymptotic momentum distribution of fragment ions on the orientation of the molecular axis, and apply it to investigate dissociative double ionization of CO 2 . We show that such measurements can directly test the validity of the axial recoil approximation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2013

Galactic Cosmic Ray Event-Based Risk Model (GERM) Code

This software describes the transport and energy deposition of the passage of galactic cosmic rays in astronaut tissues during space travel, or heavy ion beams in patients in cancer therapy. Space radiation risk is a probability distribution, and time-dependent biological events must be accounted for physical description of space radiation transport in tissues and cells. A stochastic model can calculate the probability density directly without unverified assumptions about shape of probability density function. The prior art of transport codes calculates the average flux and dose of particles behind spacecraft and tissue shielding. Because of the signaling times for activation and relaxation in the cell and tissue, transport code must describe temporal and microspatial density of functions to correlate DNA and oxidative damage with non-targeted effects of signals, bystander, etc. These are absolutely ignored or impossible in the prior art. The GERM code provides scientists data interpretation of experiments; modeling of beam line, shielding of target samples, and sample holders; and estimation of basic physical and biological outputs of their experiments. For mono-energetic ion beams, basic physical and biological properties are calculated for a selected ion type, such as kinetic energy, mass, charge number, absorbed dose, or fluence. Evaluated quantities are linear energy transfer (LET), range (R), absorption and fragmentation cross-sections, and the probability of nuclear interactions after 1 or 5 cm of water equivalent material. In addition, a set of biophysical properties is evaluated, such as the Poisson distribution for a specified cellular area, cell survival curves, and DNA damage yields per cell. Also, the GERM code calculates the radiation transport of the beam line for either a fixed number of user-specified depths or at multiple positions along the Bragg curve of the particle in a selected material. The GERM code makes the numerical estimates of basic physical and biophysical quantities of high-energy protons and heavy ions that have been studied at the NASA Space Radiation Laboratory (NSRL) for the purpose of simulating space radiation biological effects. In the first option, properties of monoenergetic beams are treated. In the second option, the transport of beams in different materials is treated. Similar biophysical properties as in the first option are evaluated for the primary ion and its secondary particles. Additional properties related to the nuclear fragmentation of the beam are evaluated. The GERM code is a computationally efficient Monte-Carlo heavy-ion-beam model. It includes accurate models of LET, range, residual energy, and straggling, and the quantum multiple scattering fragmentation (QMSGRG) nuclear database.

Cucinotta, Francis A.↗
DOE OSTI2024

Oak Ridge National Laboratory Neutron Spin Echo Beamline on NB-2

Neutron Spin Echo (NSE) spectroscopy uniquely measures the Q-dependence of slow relaxation dynamics, and having such a machine at HFIR will advance polarized neutron spectroscopy and promote the study of biophysical and energy materials using neutrons. The 2018 Instrument Advisory Board (IAB) advises Neutron Sciences Directorate, ORNL to upgrade the cold neutron delivery system at HFIR, addressing geometrical challenges in neutron transfer from the bright cold source. This is being optimized now: the entrance to the guide system is moving closer to the source, and more guides are being added. However, there are still significant challenges to making the cold flux at the sample world-class. With the proposed guide system, NB-2 can deliver more than 10 6 polarized neutrons per cm 2 /sec below 8Å, and probably 10Å with further optimization. With no potential for running user experiments with neutrons longer than 15Å, the utility of such a machine for biological studies is limited. However, there is a large demand for this class of spectrometer in materials science, chemical engineering, and nanotechnology. Workshops held over the past decade and the three-source vision have highlighted the community's need for a low-angle, high-energy resolution spectrometer at HFIR. A detailed design study of this spectrometer is needed, focusing on optimizing the neutron delivery system for SANS-type studies. This neutron spectrometer would complement studies already performed on both SNSNSE and BASIS by overlapping and extend the dynamic range accessible, while acknowledging the proposed machine, EXPANSE, at the Second Target Station. Two technologies currently exist for such a spectrometer: traditional DC solenoids like those on BL-15 at SNS and IN15 at the ILL, and RF-flipper-based machines similar to RESEDA at FRMII. Either instrument would satisfy the scientific justification. However, there are strong business and scientific cases to utilize the thermal flux, the resonant expertise developed at ORNL, and the potential to utilize an entangled beam of neutrons to probe quantum matter by adopting the neutron resonant spin-echo configuration.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Flow-driven spectral chaos (FSC) method for simulating long-time dynamics of arbitrary-order non-linear stochastic dynamical systems

Uncertainty quantification techniques such as the time-dependent generalized polynomial chaos (TD-gPC) use an adaptive orthogonal basis to better represent the stochastic part of the solution space (aka random function space) in time. However, because the random function space is constructed using tensor products, TD-gPC-based methods are known to suffer from the curse of dimensionality. Here, we introduce a new numerical method called the flow-driven spectral chaos (FSC) which overcomes this curse of dimensionality at the random-function-space level. The proposed method is not only computationally more efficient than existing TD-gPC-based methods but is also far more accurate. The FSC method uses the concept of enriched stochastic flow maps to track the evolution of a finite-dimensional random function space efficiently in time. To transfer the probability information from one random function space to another, two approaches are developed and studied herein. In the first approach, the probability information is transferred in the mean-square sense, whereas in the second approach the transfer is done exactly using a new theorem that was developed for this purpose. The FSC method can quantify uncertainties with high fidelity, especially for the long-time response of stochastic dynamical systems governed by ODEs of arbitrary order. Six representative numerical examples, including a nonlinear problem (the Van-der-Pol oscillator), are presented to demonstrate the performance of the FSC method and corroborate the claims of its superior numerical properties. Finally, a parametric, high-dimensional stochastic problem is used to demonstrate that when the FSC method is used in conjunction with Monte Carlo integration, the curse of dimensionality can be overcome altogether.

(nonlinear) stochastic dynamical systems↗
DOE OSTI2025

Langevin Dynamics modeling of gas-phase ion-ion recombination (Final Technical Report)

A self-consistent trajectory simulation approach to model MN reactions (Fig. 1) which incorporates the probability of electron transfer as a Monte Carlo operator (Fig. 2) was developed and published as Liu et al. J. Chem. Phys. 159, 114111 (2023). The electron transfer probability p ET estimated using the two-state Landau-Zener (LZ) theory was incorporated into classical trajectory simulations to elicit predictions of MN reaction cross-section σ (vacuum) or rate constant β (finite pressure). Electronic structure calculations with multireference configuration interaction (MRCI) and large correlation consistent basis sets were used to derive inputs to the LZ theory. The key advance of our trajectory simulation approach is the incorporation of electron transfer probability and the inclusion of the effect of ion-neutral interactions on MN using a Langevin representation of the effect of neutral gas on ions. For H + – H - and Li + – H(D) - pairs, our approach quantitatively agrees with measured speed-dependent cross-sections for up to ~10 5 m/s. For the ion pair Ne + – Cl - , our predictions of the MN rate constant at ~1 torr are a factor of ~2 – 3 higher than the experimentally measured value. Similarly, for Xe + – F - in the pressure range of ~20000 – 80000 Pa, our predictions of the MN rate constant are ~20% lower but are in excellent qualitative agreement with experimental data. The paradigm of using trajectory simulations to self-consistently model MN reactions is the basis for inclusion of additional non-classical, and static magnetic and electric field effects. Subsequent work, published as Roy et al. focused on modeling recombination rate constant for three ion pairs (rare gas Ar + cation and halide anions): Ar + – Cl - , Ar + – Br - , Ar + – I - , 2) considering spin-orbit couplings in the electronic structure calculations to obtain high-fidelity estimates of the electron transfer probability and incorporated within the classical trajectory simulations to elicit predictions. In addition to calculations of ion-ion recombination rate constants, a classical trajectory simulation technique (published as Roy et al. J. Chem. Phys. 162(9), 094104 (2023)) that uses quaternions to represent orientation of non-spherical particles (ions or aerosol particles) was developed to simulate the recombination of diatomic or more generally, polyatomic molecules. Finally, several other ion pairs such as Ne + – Cl - , Kr + – Cl - , were explored using the developed semi-classical trajectory simulations to understand various challenges in tackling electronic structure calculations. Using empirical approaches to parameterize the electron transfer radius, trajectory simulations were also used to probe the effect of ion number density on MN rate constant.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Quantum Choreography of the Nucleus: Rotations, Vibrations, and Emergent Structure

Abstract Nuclei are complex many-body quantum systems where interactions of the neutrons and protons via the strong, the weak, and the electromagnetic forces lead to the emergence of simple patterns of energy states that have been described by various theoretical approaches. One of the goals of all the theoretical models is the development of a universal theory that can be applied across the entire chart of nuclides. Significant progress has been made by experiments as well as the increasing sophistication of models, but a universal theory has yet to be established. A recent review of nuclei in the Z = 50–82 region of the chart of nuclides has analyzed all the available compiled data from several decades of studies towards a clarification of the low-lying structure of nuclei. Other reviews have reported and explained the emergence of multiple different shapes in nuclei at somewhat higher excitation energies than the ground state. Some have challenged the interpretation of the first excited K π = 0 + band as a vibration of ground state. This work attempts to provide a guide to determining the nature of the first excited K π = 0 + band in nuclei by the combined use of nuclear lifetimes, energy level evolutions, dynamic moments of inertia, and intrinsic quadrupole moments extracted from transition probabilities. The result is that for a subset of the nuclei in this region, the K π = 0 + band is consistent with the traditional 𝛽-vibration description of an oscillation built on the ground state.

CHFB+5DCH↗
DOE OSTI2025

Nonperturbative and perturbative dynamics of a light QCD axion: Dark matter and the strong 𝐶⁢𝑃 problem

Considerable theoretical efforts have gone into expanding the reach of the quantum chromodynamics (QCD) axion beyond its canonical mass–decay-constant relation. The 𝑍 𝒩 QCD axion model reduces the QCD axion mass naturally, by invoking a discrete 𝑍 𝒩 symmetry through which the axion field is coupled to 𝒩 copies of the Standard Model. Before the QCD phase transition at temperature 𝑇 QCD , the 𝑍 𝒩 potential has a minimum at misalignment angle 𝜃 = 𝜋. At 𝑇 QCD , 𝜃 = 𝜋 becomes a maximum; the axion potential becomes exponentially suppressed and develops 𝒩 minima—only one of which actually solves the strong 𝐶⁢𝑃 problem. Before 𝑇 QCD , 𝜃 relaxes toward 𝜋. After 𝑇 QCD , the axion field starts from around the hilltop and may have sufficient kinetic energy to overcome the newly suppressed potential barriers. Such a field evolution leads to nonperturbative effects via the self-interactions near the hilltop, which can cause the exponential growth of fluctuations and backreaction on the coherent motion. This behavior can influence the relic density of the field and the minimum in which it settles. We conduct the first lattice simulations of the 𝑍 𝒩 QCD axion using 𝒞osmoℒattice to accurately calculate dark matter abundances and find nonperturbative dynamics reduce the abundance by up to a factor of two. We furthermore find that the probability of solving the strong 𝐶⁢𝑃 problem tends to diverge considerably from the naïve expectation of 1/𝒩.

Axions↗
DOE OSTI2020

Weak scaling of the contact distance between two fluctuating interfaces with system size

A pair of flat parallel surfaces, each freely diffusing along the direction of their separation, will eventually come into contact. If the shapes of these surfaces also fluctuate, then contact will occur when their centers-of-mass remain separated by a nonzero distance ℓ. An example of such a situation is the motion of interfaces between two phases at conditions of thermodynamic coexistence, and in particular the annihilation of domain wall pairs under periodic boundary conditions. Here we present a general approach to calculate the probability distribution of the contact distance ℓ and determine how its most likely value ℓ* depends on the surfaces' lateral size L. Using the Edward-Wilkinson equation as a model for interfaces, we demonstrate that ℓ* scales weakly with system size, i.e., the dependence of ℓ* on L for both (1+1)- and (2+1)-dimensional interfaces is such that lim L→∞ (ℓ*/L) = 0. In particular, for (2+1)-dimensional interfaces ℓ* is an algebraic function of logL, a result that is confirmed by computer simulations of slab-shaped domains formed under periodic boundary conditions. Overall, this weak scaling implies that such domains remain topologically intact until ℓ becomes very small compared to the lateral size of the interface, contradicting expectations from equilibrium thermodynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2011

Quantum-Classical Hybrid for Information Processing

Based upon quantum-inspired entanglement in quantum-classical hybrids, a simple algorithm for instantaneous transmissions of non-intentional messages (chosen at random) to remote distances is proposed. The idea is to implement instantaneous transmission of conditional information on remote distances via a quantum-classical hybrid that preserves superposition of random solutions, while allowing one to measure its state variables using classical methods. Such a hybrid system reinforces the advantages, and minimizes the limitations, of both quantum and classical characteristics. Consider n observers, and assume that each of them gets a copy of the system and runs it separately. Although they run identical systems, the outcomes of even synchronized runs may be different because the solutions of these systems are random. However, the global constrain must be satisfied. Therefore, if the observer #1 (the sender) made a measurement of the acceleration v(sub 1) at t =T, then the receiver, by measuring the corresponding acceleration v(sub 1) at t =T, may get a wrong value because the accelerations are random, and only their ratios are deterministic. Obviously, the transmission of this knowledge is instantaneous as soon as the measurements have been performed. In addition to that, the distance between the observers is irrelevant because the x-coordinate does not enter the governing equations. However, the Shannon information transmitted is zero. None of the senders can control the outcomes of their measurements because they are random. The senders cannot transmit intentional messages. Nevertheless, based on the transmitted knowledge, they can coordinate their actions based on conditional information. If the observer #1 knows his own measurements, the measurements of the others can be fully determined. It is important to emphasize that the origin of entanglement of all the observers is the joint probability density that couples their actions. There is no centralized source, or a sender of the signal, because each receiver can become a sender as well. An observer receives a signal by performing certain measurements synchronized with the measurements of the others. This means that the signal is uniformly and simultaneously distributed over the observers in a decentralized way. The signals transmit no intentional information that would favor one agent over another. All the sequence of signals received by different observers are not only statistically equivalent, but are also point-by-point identical. It is important to assume that each agent knows that the other agent simultaneously receives the identical signals. The sequences of the signals are true random, so that no agent could predict the next step with the probability different from those described by the density. Under these quite general assumptions, the entangled observers-agents can perform non-trivial tasks that include transmission of conditional information from one agent to another, simple paradigm of cooperation, etc. The problem of behavior of intelligent agents correlated by identical random messages in a decentralized way has its own significance: it simulates evolutionary behavior of biological and social systems correlated only via simultaneous sensoring sequences of unexpected events.

Zak, Michail↗
DOE OSTI2022

Improved Heralded Single-Photon Source with a Photon-Number-Resolving Superconducting Nanowire Detector

Deterministic generation of single photons is essential for many quantum information technologies. A bulk optical nonlinearity emitting a photon pair, where the measurement of one of the photons heralds the presence of the other, is commonly used with the caveat that the single-photon emission rate is constrained due to a trade-off between multiphoton events and pair emission rate. Using an efficient and low noise photon-number-resolving superconducting nanowire detector we herald, in real time, a single photon at telecommunication wavelength. We perform a second-order photon correlation g 2 ( 0 ) measurement of the signal mode conditioned on the measured photon number of the idler mode for various pump powers and demonstrate an improvement of a heralded single-photon source. We develop an analytical model using a phase-space formalism that encompasses all multiphoton effects and relevant imperfections, such as loss and multiple Schmidt modes. We perform a maximum-likelihood fit to test the agreement of the model to the data and extract the best-fit mean photon number μ of the pair source for each pump power. A maximum reduction of 0.118 ± 0.012 in the photon g 2 ( 0 ) correlation function at μ = 0.327 ± 0.007 is obtained, indicating a strong suppression of multiphoton emissions. For a fixed g 2 ( 0 ) = 7 × 10 − 3 , we increase the single pair generation probability by 25%. Our experiment, built using fiber-coupled and off-the-shelf components, delineates a path to engineering ideal sources of single photons.

Davis, Samantha I.↗
DOE OSTI2026

Improved energies and local energies with weighted variational Monte Carlo

Neural network parametrizations have increasingly been used to represent the ground and excited states in variational Monte Carlo (VMC) with promising results. However, traditional VMC methods only optimize the wave function in regions of peak probability. The wave function is uncontrolled in the tails of the probability distribution, which can limit the accuracy of the trained wave function. To improve the approximation accuracy in the probability tails, this paper interprets VMC as a gradient flow in the space of wave functions, followed by a projection step. From this perspective, arbitrary probability distributions can be used in the projection step, allowing the user to prioritize accuracy in different regions of state space. Motivated by this theoretical perspective, the paper tests a weighted VMC method on the antiferromagnetic Heisenberg model for a periodic spin chain. Compared to traditional VMC, weighted VMC reduces the error in the ground state energy by a factor of 2, and it reduces the errors in the local energies away from the mode by large factors of 10 2 –10 4 .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Silica nanocluster binding rate coefficients from molecular dynamics trajectory calculations

Oxide nanoparticle growth from vapor phase precursors occurs in high temperature aerosol reactors first via the formation of nanoclusters, which are nanometer-scale condensed-phase species composed of 10 1 -10 2 atoms. The binding rate for nanoclusters, defined as the rate at which two nanoclusters collide with and stick to one another, is hence of interest in predicting nanoparticle size distribution functions in an aerosol. We have utilized molecular dynamics (MD) trajectory calculations to determine the homogeneous (equal-sized) and heterogeneous (disparate-sized) binding rate coefficients of SiO2 (silica) nanoclusters composed of 18, 144, and 333 atoms. MD trajectory calculations incorporated all-atom models of nanoclusters using a combined Born–Huggins–Mayer-Lennard-Jones potential model, which accounts for both short range interactions and electrostatic interactions. MD calculations were utilized to determine the binding probability as a function of both initial relative velocity and impact parameter; integration of the binding probability across impact parameter and relative velocity space yields the binding rate coefficient. MD trajectory calculations reveal that the most common type of encounters between nanoclusters leading to binding are grazing collisions, i.e. instances where collision would not occur without induced dipole potential influences. The resulting binding rate coefficients are found to be extremely weakly dependent on system temperature, which is in contrast to the use of rate coefficient models which are the product of a hard-sphere collision rate coefficient and a constant enhancement factor (leading to a rate coefficient proportional to the square root of temperature). Enhancement factors defined from MD trajectory calculations fall in the range of 3–9 as system temperature decreases from 1500 K to 300 K. Furthermore, such large values suggest that potential interactions need to be considered when modeling oxide nanocluster growth in the gas phase.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Mixed quantum-classical methods for polaron spectral functions

In this work, using two distinct semiclassical approaches—namely, the mean-field Ehrenfest method and the mapping approach to surface hopping—we investigate the spectral function of a single charge interacting with phonons on a lattice. This quantity is relevant for the description of angle-resolved photoemission experiments. Focusing on the one-dimensional Holstein model, we compare the performance of these approaches across a range of coupling strengths and lattice sizes, exposing the relative strengths and weaknesses of each. We demonstrate that these approaches can be efficiently applied with reasonable accuracy to ab initio polaron models. Furthermore, our work provides a route to the calculation of spectral properties in realistic electron–phonon-coupled systems in a computationally inexpensive manner with encouraging accuracy.

Atomic and molecular spectra↗
DOE OSTI2021

A quantitative model for a nanoscale switch accurately predicts thermal actuation behavior

Manipulation of temperature can be used to actuate DNA origami nano-hinges containing gold nanoparticles. We develop a physical model of this system that uses partition function analysis of the interaction between the nano-hinge and nanoparticle to predict the probability that the nano-hinge is open at a given temperature. The model agrees well with experimental data and predicts experimental conditions that allow the actuation temperature of the nano-hinge to be tuned over a range of temperatures from 30 °C to 45 °C. Additionally, the model identifies microscopic interactions that are important to the macroscopic behavior of the system, revealing surprising features of the system. This combination of physical insight and predictive potential is likely to inform future designs that integrate nanoparticles into dynamic DNA origami structures or use strand binding interactions to control dynamic DNA origami behavior. Furthermore, our modeling approach could be expanded to consider the incorporation, stability, and actuation of other types of functional elements or actuation mechanisms integrated into nucleic acid devices.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Forecasting cosmological constraints from the weak lensing magnification of type Ia supernovae measured by the Nancy Grace Roman Space Telescope

We report that the weak lensing magnification of type Ia supernovae (SNe Ia) is sensitive to the clustering of matter and provides an independent cosmological probe complementary to SN Ia distance measurements. The Nancy Grace Roman Space Telescope is uniquely sensitive to this measurement as it can discover high redshift SNe Ia and measure them with high precision. We present a methodology for reconstructing the probability distribution of the weak lensing magnification μ of SNe Ia, p(μ), from observational data, and using it to constrain cosmological parameters. We find that the reconstructed p(μ) can be fitted accurately by a stretched Gaussian distribution and used to measure the variance of μ, ξ μ , which can be compared to theoretical predictions in a likelihood analysis. Applying our methodology to a set of realistically simulated SNe Ia expected from the Roman Space Telescope, we find that using the weak lensing magnification of the SNe Ia constrains a combination of matter density Ω m and matter clustering amplitude σ 8 . SN Ia distances alone lead to a better than 1% measurement of Ωm. The combination of SN Ia weak lensing magnification and distance measurements result in a ~10 % measurement on σ 8 . The SNe Ia from Roman will be powerful in constraining the cosmological model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Short-period Variability in the Globular Cluster NGC 419 and the SMC Field

Delta Scuti (δ Sct) stars have been extensively studied in our Galaxy but far less in extragalactic systems. Here we study the population of δ Sct variables in NGC 419, an intermediate-age globular cluster of the Small Magellanic Cloud (SMC), using g, r, and i Gemini-S/Gemini Multi-Object Spectrograph (GMOS) time-series observations. Our goal is to study the role of such variables in the cluster extended main-sequence turnoff (MSTO). We report the discovery of 54 δ Sct stars and three eclipsing binaries in the NGC 419 field. We find only a handful of the δ Sct stars at the MSTO of NGC 419, while the majority are fainter, indicating that the cluster is younger (≲1.2 Gyr) than previously thought. Considering their radial distribution, we identify only six δ Sct stars as probable members of NGC 419, while the 48 remaining are likely δ Sct stars of the SMC field. Cluster δ Sct stars appear close to the red edge of the MSTO, supporting the idea that the extended MSTO has its origin in an age spread. The 48 field δ Sct stars represent the largest detection of δ Sct stars made in the SMC. The period distribution of these newly detected δ Sct stars (0.04 day ≲ P ≲ 0.15 day) is similar to that detected in other systems. The amplitude distribution (0.05 mag ≲ Δr ≲ 0.60 mag) is likely biased because of the lack of low-amplitude stars. We finally use the δ Sct stars to calculate distances using different period–luminosity relations. The average distance moduli obtained are 18.76 ± 0.14 mag for NGC 419 and 18.86 ± 0.11 mag for the SMC field, which agree with previous measurements.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Identifying and tracking bubbles and drops in simulations: A toolbox for obtaining sizes, lineages, and breakup and coalescence statistics

Knowledge of bubble and drop size distributions in two-phase flows is important for characterizing a wide range of phenomena, including combustor ignition, sonar communication, and cloud formation. The physical mechanisms driving the background flow also drive the time evolution of these distributions. Accurate and robust identification and tracking algorithms for the dispersed phase are necessary to reliably measure this evolution and thereby quantify the underlying mechanisms in interface-resolving flow simulations. The identification of individual bubbles and drops traditionally relies on an algorithm used to identify connected regions. This traditional algorithm can be sensitive to the presence of spurious structures. A cost-effective refinement is proposed to maximize volume accuracy while minimizing the identification of spurious bubbles and drops. An accurate identification scheme is crucial for distinguishing bubble and drop pairs with large size ratios. The identified bubbles and drops need to be tracked in time to obtain breakup and coalescence statistics that characterize the evolution of the size distribution, including breakup and coalescence frequencies, and the probability distributions of parent and child bubble and drop sizes. An algorithm based on mass conservation is proposed to construct bubble and drop lineages using simulation snapshots that are not necessarily from consecutive time steps. These lineages are then used to detect breakup and coalescence events, and obtain the desired statistics. Accurate identification of large-size-ratio bubble and drop pairs enables accurate detection of breakup and coalescence events over a large size range. Accurate detection of successive breakup and coalescence events requires that the snapshot interval be an order of magnitude smaller than the characteristic breakup and coalescence times to capture these successive events while minimizing the identification of repeated confounding events. Together, these algorithms serve as a toolbox for detailed analysis of two-phase simulations, and enable insights into the mechanisms behind bubble and drop formation and evolution in flows of practical importance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A Deep Study of an Intermediate-age Open Cluster SAI 35 (Juchert 20) Using Ground-based Imaging and Gaia EDR3 Astrometry

We present a CCD UBVI photometric study of poorly studied intermediate-age open cluster SAI 35 (Juchert 20) for the first time. To accomplish this study, we also used LAMOST DR5, Two Micron All Sky Survey, and Gaia EDR3 databases. We identified 214 most probable cluster members with membership probability higher than 50%. The mean proper motion of the cluster is found as μ{sub α}cos δ=1.10±0.01 and μ {sub δ} = −1.66 ± 0.01 mas yr{sup −1}. We find the normal interstellar extinction law using the various two-color diagrams. The age, distance, reddening, and radial velocity of the cluster are estimated to be 360 ± 40 Myr, 2.9 ± 0.15 kpc, 0.72 ± 0.05 mag, and −91.62 ± 6.39 km s{sup −1}, respectively. The overall mass function slope for main-sequence stars is found to be 1.49 ± 0.16 within the mass range 1.1–3.1 M {sub ⊙}, which is in agreement with Salpeter’s value within uncertainty. The present study demonstrates that SAI 35 is dynamically relaxed. Galactic orbital parameters are determined using Galactic potential models. We found that this object follows a circular path around the Galactic center.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗