Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “structural collapse”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Wellbore Stability and Mud Loss Management in Geothermal Drilling: Optimizing Mud Weight to Mitigate Tensile Wellbore Fracturing at The Geysers, California

As part of a U.S. Department of Energy (DOE) Geothermal Technologies Office-funded initiative, Geysers Power Company, LLC, a subsidiary of Calpine Corporation, has been working to enhance drilling performance at the world’s largest geothermal field, The Geysers, in northern California. In a recent drilling operation of the GDC-36 well, excessive mud losses were encountered, initially addressed through repeated but largely ineffective cement plugging. Ultimately, the most effective strategy was to drill blind through the loss zones, made feasible by the high rate of penetration (ROP) achieved with PDC bits, allowing significant progress before the mud tanks were depleted and water-sensitive argillic formation layers could collapse. In response to these challenges, the project team explored alternative methods to minimize downtime and risks associated with cement plugging and continuous mud loss and to contemplate the driving mechanisms for the losses. Wellbore imaging using Formation MicroImager (FMI) and Ultrasonic Borehole Imager (UBI) tools revealed longitudinal tensile fractures, which were attributed to mud weights exceeding the minimum circumferential stress resulting from the native stress field and formation pressure. This study examines the mud losses encountered and leverages wellbore imaging data to understand the mechanisms behind mud induced tensile fracturing in specific rock facies. Understanding fracture behavior across different lithologies is crucial, as fractures within the reservoir can enhance steam migration throughout the system. The reservoir at The Geysers lies within the Mesozoic Franciscan Assemblage, a tectonic mélange formed by subduction. It consists of metamorphosed turbidite sandstone (greywacke) and mudstone (argillite), oceanic upper crust (including greenstone and chert), and serpentinized ultramafic rocks - each exhibiting distinct geomechanical fracturing properties. The structural fabric of the Franciscan Assemblage was shaped by low-angle Mesozoic thrust faulting and later overprinted by sub-vertical strike-slip structures related to the Pacific-North American plate boundary. A wellbore stability model was developed using core measurements and logs to simulate fracturing scenarios during drilling under varying stress conditions. These simulations guided the development of an optimized mud weight management strategy that should enable adaptive adjustments during drilling, reducing the likelihood of tensile fracturing and mud losses, ultimately improving operational efficiency.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Constraints on the normal branch of DGP gravity from SPT galaxy clusters with DES and HST weak-lensing mass calibration and from P l a n c k PR4 CMB anisotropies

We present constraints on the normal branch of the Dvali-Gabadadze-Porrati (nDGP) braneworld gravity model from the abundance of massive galaxy clusters. On scales below the nDGP crossover scale r c , the nDGP model features an effective gravitylike fifth force that alters the growth of structure, leading to an enhancement of the halo mass function (HMF) on cluster scales. The enhanced cluster abundance allows for constraints on the nDGP model using cluster samples. We employ the South Pole Telescope (SPT) cluster sample, selected through the thermal Sunyaev-Zel’dovich effect with the SPT and with mass calibration using weak-lensing data from the Dark Energy Survey (DES) and the Hubble Space Telescope (HST). The cluster sample contains 1,005 clusters with redshifts 0.25 < z < 1.78 , which are confirmed with the multicomponent matched filter algorithm using optical and near-infrared data. Weak-lensing data from DES and HST enable a robust mass measurement of the cluster sample. We use DES Year 3 data for 688 clusters with redshifts z < 0.95 , and HST data for 39 clusters with redshifts 0.6 < z < 1.7 . We account for the enhancement in the HMF through a semi-analytic correction factor to the standard cosmology HMF derived from the spherical collapse model in the nDGP model. We then further calibrate this model using N -body simulations. In addition, for the first time, we analyze the primary CMB temperature and polarization anisotropy measurements from Planck PR4 within the nDGP model. We obtain a competitive constraint from the joint analysis of the SPT cluster abundance with the Planck PR4 data, and report an upper bound of 1 / H 0 r c < 1.41 at 95% when assuming a cosmology with massive neutrinos.

Vogt, S. M.L. [Munich U. Observ.; LMU Munich (main↗
DOE OSTI2022

Concurrent Pressure-Induced Spin-State Transitions and Jahn–Teller Distortions in MnTe

Cooperative Jahn-Teller (JT) distrotion accompanied by electromagnetic-lattice coupling plays a crucial role in numerous important phenomena in materials science. The JT distortion occurring in high-spin 3d 4 and 3d 9 configurations in octahedral complexes is common, but similarly large distortion due to lifting the degeneracy of t 2g orbits has so far seldomly observed. Here we report the discovery of pressure-induced cooperative JT effect coupled with a large volume-collapse of ~12.8% and a magnetic moment collapse from S = 5/2 to S = 1/2 at ~20.9 GPa in MnTe. The first-principles calculation indicates that weakened p-d hybridization and enhanced d states localization result in lifting t 2g orbital degeneracy in Mn 2+ (d 5 ) which induced a cooperative JT effect-triggered electromagnetic-lattice coupling. Furthermore, these findings suggest a new mechanism for the cooperative JT distortion via pressure induced ligand distortion in antiferromagnetic materials.

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

Clumpy structures within the turbulent primordial cloud

ABSTRACT The primordial clouds in the mini-haloes hatch the first generation stars of the Universe, which play a crucial role in cosmic evolution. In this paper, we investigate how turbulence impacts the structure of primordial star-forming clouds. Previous cosmological simulations of the first star formation predicted a typical mass of around $\mathrm{ 100 \, M_\odot }$. This conflicts with recent observations of extremely metal-poor stars, suggesting a lower mass scale of about $\mathrm{25 \, M_\odot }$. The discrepancy may arise from unresolved turbulence in the star-forming cloud, driven by primordial gas accretion during mini-halo formation in the previous simulations. To quantitatively examine the turbulence effect on the primordial cloud formation, we employ the adaptive mesh refinement code Enzo to model the gas cloud with primordial composition, including artificially driven turbulence on the cloud scale and relevant gas physics. This artificially driven turbulence utilizes a stochastic forcing model to mimic the unresolved turbulence inside mini-haloes. Our results show that the turbulence with high Mach number and compressional mode effectively fragments the cloud into several clumps, each with dense cores of $\mathrm{22.7 - 174.9 \, M_\odot }$ that undergo Jeans instability to form stars. Fragmentation caused by intense and compressive turbulence prevents a runaway collapse of the cloud. The self-bound clumps with smaller masses in the turbulent primordial clouds suggest a possible pathway to decrease the theoretical mass scale of the first stars, further reconciling the mass discrepancy between simulations and observations.

Astronomy & Astrophysics↗
DOE OSTI2024

Global properties of nuclei at finite-temperature within the covariant energy density functional theory

In stellar environments nuclei appear at finite temperatures, becoming extremely hot in core-collapse supernovae and neutron-star mergers. However, due to theoretical and computational complexity, most model calculations of nuclear properties are performed at zero temperature, while those existing at finite temperatures are limited only to selected regions of the nuclide chart. In this study we perform the global calculation of nuclear properties for even-even 8 ≤ Z ≤ 104 nuclei at temperatures in range 0 ≤ T ≤ 2 MeV. Calculations are based on the finite-temperature relativistic Hartree-Bogoliubov model supplemented by the Bonche-Levit-Vautherin vapor subtraction procedure. We find that near the neutron-drip line the continuum states have significant contribution already at moderate temperature T ≈ 1 MeV, thus emphasizing the necessity of the vapor subtraction procedure. Results include neutron emission lifetimes, quadrupole deformations, neutron-skin thickness, proton and neutron pairing gaps, entropy and excitation energy. Up to the temperature T ≈ 1 MeV, the nuclear landscape is influenced only moderately by the finite-temperature effects, mainly by reducing the pairing correlations. Here, as the temperature increases further, the effects on nuclear structures become pronounced, reducing both the deformations and the shell effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Novel challenges in tracking self-interacting dark matter subhalos

Cosmological N-body simulations are among the primary tools for studying structure formation in the Universe. Analyses of these simulations critically depend on accurately identifying and tracking dark matter subhalos over time. In recent years, several new algorithms have been developed to improve the accuracy and consistency of subhalo tracking in cold dark matter (CDM) simulations. These algorithms should be revisited in the context of new physics beyond gravity, which can modify the evolution and final properties of subhalo populations. In this work, we apply the particle-tracking-based subhalo finder Symfind to velocity-dependent self-interacting dark matter (SIDM) simulations with large cross section amplitudes to assess the performance of particle-tracking methods beyond the CDM paradigm. We find that the core-particle-tracking technique, which is key to the success of these algorithms in CDM, does not always yield accurate results in SIDM. In particular, the interplay between dark matter self-interactions and tidal stripping can cause the diffusion of core particles to larger radii, leading particle-tracking-based algorithms to prematurely lose track of SIDM subhalos. For massive core-expansion subhalos and core-collapse subhalos that experience close or repeated pericentric passages, a significant fraction of core particles can be lost, and particle-tracking-based finders such as Symfind offer no clear advantage over traditional methods that rely on identifying phase-space overdensities. On the other hand, for subhalos with large pericentric distances or fewer, more distant passages, Symfind tends to outperform. These differences depend sensitively on the cross section amplitude and turnover velocity of the underlying SIDM model. We therefore recommend a hybrid approach that leverages the strengths of both techniques to produce complete and robust catalogs of core-expansion and core-collapse SIDM subhalos.

cosmological simulations↗
DOE OSTI2026

Pre-Equilibrium De-Excitations for Neutrino-Nucleus Interactions

The Deep Underground Neutrino Experiment (DUNE) is sensitive to MeV-scale energy depositions from low-energy astrophysical neutrinos, including those from core-collapse supernovae. Interpreting these detector signals requires accurate modeling of the nuclear de-excitation that follows from the neutrino-nucleus interaction. The MARLEY (Model of Argon Reaction Low-Energy Yields) event generator specializes in the low-energy regime. MARLEY currently assumes the residual nucleus equilibrates immediately after the primary interaction. This omits the intermediate pre-equilibrium stage in which energy redistributes among nucleons until statistical equilibrium is reached. While pre-equilibrium effects are well established for nucleon-induced reactions, they have not previously been studied for neutrino-nucleus interactions. This work addresses that gap by implementing a two-component exciton model, which is the first dedicated treatment of pre-equilibrium de-excitation for neutrino-nucleus interactions, restructured around MARLEY's existing class hierarchy to prepare for direct integration, including particle-hole state densities, internal transition rates, and pre-equilibrium particle emission. The calculations show encouraging agreement with the TALYS-2.2 nuclear reaction code for neutron-nucleus interactions. We further propose a concrete integration path into the full MARLEY event generator, including derived class structure and an extended event record for pre-equilibrium vertices in support of future reweighting. Remaining work focuses on refining the emission width calculation, adding $\gamma$-ray emission, and completing this integration to quantify the impact of pre-equilibrium effects on the expected low-energy neutrino signals in DUNE and similar experiments.

Visser, Erin [Michigan State U., East Lansing (mai↗
DOE OSTI2026

Nutrient limitation shapes functional traits of mycorrhizal fungi and phosphorus-cycling bacteria across an elevation gradient

In nutrient-limited high-elevation ecosystems, plants rely on arbuscular mycorrhizal (AM) fungi to provide mineral phosphorus (P) in the form of phosphate (PO43-). AM fungi gather these nutrients from phosphorus-cycling bacteria (PCBs) that can mineralize PO43- from organic matter and solubilize mineral-bound P. How climate, soil factors, and nutrient limitation influence AM fungi and PCB assembly remains unclear. We collected soil from montane meadows across a 1,000-m elevation gradient on three replicate mountainsides and analyzed AM fungal marker genes, P-cycling genes from shotgun metagenomes, and edaphic measurements. High-elevation soils had nearly 50-fold less soil PO₄³⁻ and 60% more AM fungal hyphae than low-elevation soils. AM fungal turnover was linked to changes in pH, organic carbon, and PO₄³-. The composition of 198 P-cycling genes was influenced by the AM fungal community structure. Drivers of individual PCB functional genes, including pH and organic carbon, varied with gene phylogeny. We found a trade-off in P-cycling strategies across elevation: P-rich, low-elevation soils supported root-colonizing AM fungi and organic P-mineralizing bacteria. P-poor, high-elevation soils were dominated by stress-tolerant AM fungi and mineral P-solubilizing bacteria. Our results suggest that AM fungi and PCB community turnover across elevation are both shaped by pH, organic carbon, and P availability. With continued climate warming, the structure and function of mountaintop ecosystems might shift to resemble lower elevations, disrupting long-established and specialized microbial assemblages, with consequences for P-cycling dynamics and the total P available to plant communities.IMPORTANCEPhosphorus (P) limits plant productivity in high-elevation ecosystems, yet the microbial networks that mobilize P, including arbuscular mycorrhizal (AM) fungi and phosphorus-cycling bacteria (PCBs), remain under-characterized in these nutrient-poor soils. We show that across a 10,00-m elevation gradient, AM fungi and P-cycling gene assemblages shift predictably with pH, organic carbon, and phosphate availability. Higher elevations, with less available P, select for stress-tolerant AM fungal taxa and PCB strategies geared toward mineral solubilization, while low-elevation sites favor root colonization by AM fungi and organic P mineralization. These results suggest that nutrient limitation can constrain microbial community assembly in consistent ways across landscapes. High mountain soils are low in P and rely on a network of underground AM fungi and PCB to deliver nutrients to plants. This study shows how those underground relationships reorganize with elevation and how climate change could collapse long-standing microbial strategies by pushing high-elevation ecosystems toward lowland conditions. As soils warm and dry, the microbial scaffolding that supports alpine plant life may become increasingly unstable.

arbuscular mycorrhizal fungi↗
DOE OSTI2023

Neutralizing topological obstructions to bubbles of nothing

Theories with compact extra dimensions can exhibit a vacuum instability known as a bubble of nothing. These decay modes can be obstructed if the internal manifold is stabilized by fluxes, or if it carries Wilson lines for background gauge fields, or if the instanton is incompatible with the spin structure. In each of these cases the decay can proceed by adding dynamical charged membranes or gauge fields. We give a general, bottom-up procedure for constructing approximate bubble of nothing solutions in models with internal spheres stabilized by flux and study the influence of the brane tension on the tunneling exponent, finding two branches of solutions that merge at a minimal superextremal value of the tension. In the case of Wilson operators and incompatible fermions, the relevant bubble is shown to be the Euclidean Reissner-Nordstrom black hole, and the ordinary decay exponent is modified by 1/g 2 effects. We examine the Dirac operator on this background and comment on the relevance for models of supergravity with gauged R-symmetry.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Evidence for X-Ray Emission in Excess to the Jet-afterglow Decay 3.5 yr after the Binary Neutron Star Merger GW 170817: A New Emission Component

For the first ~3 yrs after the binary neutron star merger event GW 170817, the radio and X-ray radiation has been dominated by emission from a structured relativistic off-axis jet propagating into a low-density medium with n < 0.01 cm -3 . We report on observational evidence for an excess of X-ray emission at δt > 900 days after the merger. With L x ≈ 5 × 10 38 erg s -1 at 1234 days, the recently detected X-ray emission represents a ≥ 3.2σ (Gaussian equivalent) deviation from the universal post-jet-break model that best fits the multiwavelength afterglow at earlier times. In the context of JetFit afterglow models, current data represent a departure with statistical significance ≥ 3.1σ, depending on the fireball collimation, with the most realistic models showing excesses at the level of ≥ 3.7σ. A lack of detectable 3 GHz radio emission suggests a harder broadband spectrum than the jet afterglow. These properties are consistent with the emergence of a new emission component such as synchrotron radiation from a mildly relativistic shock generated by the expanding merger ejecta, i.e., a kilonova afterglow. In this context, we present a set of ab initio numerical relativity binary neutron star (BNS) merger simulations that show that an X-ray excess supports the presence of a high-velocity tail in the merger ejecta, and argues against the prompt collapse of the merger remnant into a black hole. Radiation from accretion processes on the compact-object remnant represents a viable alternative. Neither a kilonova afterglow nor accretion-powered emission have been observed before, as detections of BNS mergers at this phase of evolution are unprecedented

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Insights into the origin of halo mass profiles from machine learning

ABSTRACT The mass distribution of dark matter haloes is the result of the hierarchical growth of initial density perturbations through mass accretion and mergers. We use an interpretable machine-learning framework to provide physical insights into the origin of the spherically-averaged mass profile of dark matter haloes. We train a gradient-boosted-trees algorithm to predict the final mass profiles of cluster-sized haloes, and measure the importance of the different inputs provided to the algorithm. We find two primary scales in the initial conditions (ICs) that impact the final mass profile: the density at approximately the scale of the haloes’ Lagrangian patch RL ($R\sim 0.7\, R_L$) and that in the large-scale environment (R ∼ 1.7 RL). The model also identifies three primary time-scales in the halo assembly history that affect the final profile: (i) the formation time of the virialized, collapsed material inside the halo, (ii) the dynamical time, which captures the dynamically unrelaxed, infalling component of the halo over its first orbit, (iii) a third, most recent time-scale, which captures the impact on the outer profile of recent massive merger events. While the inner profile retains memory of the ICs, this information alone is insufficient to yield accurate predictions for the outer profile. As we add information about the haloes’ mass accretion history, we find a significant improvement in the predicted profiles at all radii. Our machine-learning framework provides novel insights into the role of the ICs and the mass assembly history in determining the final mass profile of cluster-sized haloes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Temperature controlled transformations of giant unilamellar vesicles of amphiphilic triblock copolymers synthesized via microfluidic mixing

We report on a simple approach for synthesis of temperature-responsive giant unilamellar vesicles (GUVs) from poly(N-vinylcaprolactam) 15 -block-poly(dimethylsiloxane) 65 -block-poly(N-vinylcaprolactam) 15 (PVCL 15 -PDMS65-PVCL 15 ) triblock copolymer and non-temperature responsive small and giant vesicles from novel poly(N-vinylpyrrolidone)-block-poly(dimethylsiloxane)-block-poly(N-vinylpyrrolidone) (PVPON 15 -PDMS 65 -PVPON 15 and PVPON 6- PDMS 30- PVPON 6 ) triblock copolymers using microfluidic mixing at 25 °C. We show that temperature-responsive PVCL 15 -PDMS 65- PVCL 15 GUVs with the average diameter of 1.4 ± 0.2 µm while being stable at room temperature for at least 14 days, transformed irreversibly into small vesicles of 168 ± 40 nm after incubation of their aqueous solution at 42 °C for 24 h. We hypothesized that this transformation is induced by local compressive stresses of the vesicle membrane due to the collapse of PVCL blocks above the copolymer lower critical solution temperature (LCST) leading to the decrease of the vesicle membrane thickness. Consequently, we found that the temperature-induced size transformation of the PVCL-based GUVs at 42°C can be suppressed by substituting PVCL with its hydrophilic homologue PVPON, or by suppressing the PVCL's LCST behavior through hydrogen-bonding with tannic acid molecules. In the former case, novel PVPON n -PDMS m -PVPON n triblock copolymers (n = 15, m = 65 and n = 6, m = 30) assemble into vesicles stable from 25 °C to 55 °C as confirmed by optical and electron microscopy, dynamic light scattering (DLS) and small-angle neutron scattering (SANS). In the latter, hydrogen bonding interactions of PVCL with the polyphenol tannic acid (TA) at room temperature resulted in stable PVCL-based GUVs at 42°C as confirmed by optical, electron, and atomic force microscopies. We also found that physical crosslinking of the PVCL corona through hydrogen bonding with TA in PVCL 15- PDMS 65- PVCL 15 GUVs will delay their low pH-induced degradation at 37°C by 48 h compared to non-modified GUVs. Our findings open opportunities for the development of temperature-regulated stable micro-vehicles that would change their structural characteristics in the physiologically relevant temperature range from 25 to 42°C and can be utilized for cell mimicking studies. The developed GUVs also have potential in theranostic drug delivery as substitutes for polymer microcapsules and lipid microbubbles as well as for stimuli-triggered sensing, protection, and rapid response in an aqueous environment.

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

Quantum Theory of the Classical: Einselection, Envariance, Quantum Darwinism and Extantons

Core quantum postulates including the superposition principle and the unitarity of evolutions are natural and strikingly simple. I show that—when supplemented with a limited version of predictability (captured in the textbook accounts by the repeatability postulate)—these core postulates can account for all the symptoms of classicality. In particular, both objective classical reality and elusive information about reality arise, via quantum Darwinism, from the quantum substrate. This approach shares with the Relative State Interpretation of Everett the view that collapse of the wavepacket reflects perception of the state of the rest of the Universe relative to the state of observer’s records. However, our “let quantum be quantum” approach poses questions absent in Bohr’s Copenhagen Interpretation that relied on the preexisting classical domain. Thus, one is now forced to seek preferred, predictable, hence effectively classical but ultimately quantum states that allow observers keep reliable records. Without such (i) preferred basis relative states are simply “too relative”, and the ensuing basis ambiguity makes it difficult to identify events (e.g., measurement outcomes). Moreover, universal validity of quantum theory raises the issue of (ii) the origin of Born’s rule, p k =|ψ k | 2 , relating probabilities and amplitudes (that is simply postulated in textbooks). Last not least, even preferred pointer states (defined by einselection—environment—induced superselection)—are still quantum. Therefore, unlike classical states that exist objectively, quantum states of an individual system cannot be found out by an initially ignorant observer through direct measurement without being disrupted. So, to complete the ‘quantum theory of the classical’ one must identify (iii) quantum origin of objective existence and explain how the information about objectively existing states can appear to be essentially inconsequential for them (as it does for states in Newtonian physics) and yet matter in other settings (e.g., thermodynamics). I show how the mathematical structure of quantum theory supplemented by the only uncontroversial measurement postulate (that demands immediate repeatability—hence, predictability) leads to preferred states. These (i) pointer states correspond to measurement outcomes. Their stability is a prerequisite for objective existence of effectively classical states and for events such as quantum jumps. Events at hand, one can now enquire about their probability—the probability of a pointer state (or of a measurement record). I show that the symmetry of entangled states—(ii) entanglement—assisted invariance or envariance—implies Born’s rule. Envariance also accounts for the loss of phase coherence between pointer states. Thus, decoherence can be traced to symmetries of entanglement and understood without its usual tool—reduced density matrices. A simple and manifestly noncircular derivation of p k =|ψ k | 2 follows. Monitoring of the system by its environment in course of decoherence typically leaves behind multiple copies of its pointer states in the environment. Only pointer states can survive decoherence and can spawn such plentiful information-theoretic progeny. This (iii) quantum Darwinism allows observers to use environment as a witness—to find out pointer states indirectly, leaving systems of interest untouched. Quantum Darwinism shows how epistemic and ontic (coexisting in epiontic quantum state) separate into robust objective existence of pointer states and detached information about them, giving rise to extantons—composite objects with system of interest in the core and multiple records of its pointer states in the halo comprising of environment subsystems (e.g., photons) which disseminates that information throughout the Universe.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Quantification of Mono- and Multivalent Counterion-Mediated Bridging in Polyelectrolyte Brushes

Multivalent counterion-induced bridging interactions have been identified as the key mechanism of drastic collapse of the height of polyelectrolyte (PE) brushes. In this article, we employ all-atom molecular dynamics (MD) simulations to quantify the bridging interactions in PE brushes for counterions of different sizes and valences. We identify that unlike the current notion, bridging interactions are not the sole function of the counterion valence. Rather the bridging interactions depend on the fraction of counterions (of a given type) that get physically condensed on the PE backbone as well as the size of the counterion solvation shell. These mechanisms ensure that certain monovalent counterions demonstrate much stronger bridging interactions than that witnessed for certain divalent and trivalent counterions, while certain counterions of identical valences show drastically different bridging. Here, we argue that these counterion-specific bridging interactions eventually enable not only the significant reduction of the PE brush height in presence of certain multivalent screening counterions, but may also give rise to scenarios where the brush height reduction for certain monovalent counterions is larger than certain divalent and trivalent counterions. This latter observation contradicts the experimental findings where the multivalent counterions invariably led to a larger decrease in the height of the PE brushes: we argue that this discrepancy stems from the fact that in our simulations we only consider densely grafted and short (and hence less flexible) PE brushes that hinder the formation of different laterally inhomogeneous structures (like pinned micelles and cylindrical bundles) that would have led to a larger brush height reduction (in experiments, which invariably consider longer and less densely grafted brushes, the formation of such inhomogeneous structures are primarily responsible for larger brush height reduction in presence of multivalent counterions). Finally, we also probe the dynamic properties of the counterions (i.e., their time-dependent displacements) and their bridging interactions (i.e., lifetime of bridging interactions).

42 ENGINEERING↗
DOE OSTI2021

Numerical Investigation of Enhanced Dehumidification Processes By Using Dielectrophoresis Principles in Moist Airflows

Dispersed particle-laden flows are encountered in many building and industrial applications, such as flow in a fluidized bed, hydrocarbon transportation in pipelines, and the fouling of air-cooled heat exchangers (Kuruneru et al., 2016; Ray et al., 2019; Wang et al., 2019). Computational fluid dynamic (CFD) models have been developed in recent years to depict particle-fluid and particle-particle interactions in laminar or turbulent flows with increasing accuracy and stability. One particular particle-laden system of interest for moisture control is electrically-enhanced condensation in air and water droplet flows. Electrically-enhanced condensation consists of the use of highly charged water droplets injected in the moist air. The droplets become electric seeds that attract polar water vapor molecules to their surfaces and promote condensation. The nucleation and growth of the charged droplets deplete the vapor phase near a droplet, which is compensated for by the dielectrophoresis flow and diffusion. Dielectrophoresis flow involves surrounding vapor at a distance of about 10 to 100 nm for droplets charged by an electrospray compared to ~2 nm for a single electron charge in a droplet. As the vapor molecules collapse on the surface of the droplets, their initial electrical charge decreases with time due to the neutralization of the ions. While the physics of this phenomena is well known, engineering models for predicting the condensation rates are not available. This work computationally investigates dehumidification of moist airflow in a converging rectangular duct. The objective is to develop an engineering model that predicts water vapor condensation by employing dielectrophoresis principles. We construct a Computational Fluid Dynamics (CFD) model of the duct with electrically-enhanced condensation. The model is implemented in the open-source software OpenFOAM. We utilize the Multi-Phase Particle-In-Cell (MP-PIC) method coupled with a Population Balance Equation (PBE) approach to simulate the particle-laden system. This methodology is an Eulerian-Lagrangian approach used to simulate the droplets' behavior in the humid air. The MP-PIC approach (Andrews and O'Rourke, 1996) mitigates the computational cost by parceling several fundamental particles with similar properties (such as types, sizes, and temperature) into one computational particle. Thus, the billions of particles can be substituted by millions of computational particles without significant loss of information. The PBE was considered with the Lagrangian frame to combine the particle distribution function used in MP-PIC (Kim et al., 2020). This approach preserves mass and energy conservation between the phases in the Eulerian and Lagrangian structures. The PBE in this procedure was directly linked to the discrete parcels, making the simulation of the particle distribution computationally efficient and robust. The MP-PIC-PBE approach used in the present work was applied to the dehumidification of air. Water droplets were injected in the air stream and forced to grow according to experimentally derived correlation. The experiments were conducted on a converging duct with the same geometry and boundary conditions used to build the CFD model. This approach enabled us to approximate the effect of dielectrophoresis phenomena on the droplet and air interface. This presentation will discuss the details of the new CFD model built for the duct, the implementation of the model in OpenFOAM CFD programming language, and the experimental validation of the newly developed model. The results revealed a moderate yet measurable increase in droplet diameter due to water vapor condensation at the vapor-liquid interface of the electrically charged droplets' surface. The seed water droplet particles grew in size by capturing the water vapor in the surrounding air. The OpenFOAM model predicted reductions of humidity in the air from 5 to 10 percent.

Yel Mahi, Maliha↗
DOE OSTI2024

Classical field approximation of ultralight dark matter: Quantum break times, corrections, and decoherence

The classical field approximation is widely used to better understand the predictions of ultralight dark matter. Here, in this work, we use the truncated Wigner approximation method to test the classical field approximation of ultralight dark matter. This method approximates a quantum state as an ensemble of independently evolving realizations drawn from its Wigner function. The method is highly parallelizable and allows the direct simulation of quantum corrections and decoherence times in systems many times larger than have been previously studied in reference to ultralight dark matter. Our study involves simulation of systems in 1, 2, and 3 spatial dimensions. We simulate three systems, the condensation of a Gaussian random field in three spatial dimensions, a stable collapsed object in three spatial dimensions, and the merging of two stable objects in two spatial dimensions. We study the quantum corrections to the classical field theory in each case. We find that quantum corrections grow exponentially during nonlinear growth with the timescale being approximately equal to the system dynamical time. In stable systems the corrections grow quadratically. We also find that the primary effect of quantum corrections is to reduce the amplitude of fluctuations on the de Broglie scale in the spatial density. Finally, we find that the timescale associated with decoherence due to gravitational coupling to baryonic matter is at least as fast as the quantum corrections due to gravitational interactions. These results are consistent with the predictions of the classical field theory being accurate.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Many-Body Collective Neutrino Oscillations: Recent Developments

Neutrino flavor transformations in core-collapse supernovae and binary neutron star mergers represent a complex and unsolved problem that is integral to our understanding of the dynamics and nucleosynthesis in these environments. The high number densities of neutrinos present in these environments can engender various collective effects in neutrino flavor transformations, driven either by neutrino-neutrino coherent scattering, or in some cases, through collisional (incoherent) interactions. An ensemble of neutrinos undergoing coherent scattering among themselves is an interacting quantum many-body system -- as such, there is a tantalising prospect of quantum entanglement developing between the neutrinos, which can leave imprints on their flavor evolution histories. Here, in this study, we seek to summarize recent progress that has been made towards understanding this phenomenon.

Hypernuclei↗
DOE OSTI2023

Intergranular Hotspots: A Molecular Dynamics Study on the Influence of Compressive and Shear Work

Numerous crystal- and microstructural-level mechanisms are at play in the formation of hotspots, which are known to govern high explosives initiation behavior. Most of these mechanisms, including pore collapse, interfacial friction, and shear banding, involve both compressive and shear work done within the material and have thus far remained difficult to separate. Here, we assess hotspots formed at shocked crystal–crystal interfaces using quasi-1D molecular dynamics simulations that isolate effects due to compression and shear. Two high explosive materials are considered (TATB and PETN) that exhibit distinctly different levels of molecular conformational flexibility and crystal packing anisotropy. Temperature and intramolecular strain energy localization in the hotspot are assessed through parametric variation of the crystal orientation and two velocity components that respectively modulate compression and shear work. The resulting hotspots are found to be highly localized to a region within 5–20 nm of the crystal–crystal interface. Compressive work plays a considerably larger role in localizing temperature and intramolecular strain energy for both materials and all crystal orientations considered. Shear induces a moderate increase in energy localization relative to unsheared cases only for relatively weak compressive shock pressures of approximately 10 GPa. These results help isolate and rank the relative importance of hotspot generation mechanisms and are anticipated to guide the treatment of crystal–crystal interfaces in coarse-grained models of polycrystalline high-explosive materials.

36 MATERIALS SCIENCE↗