Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Compact model”

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 865 records · Page 48

DNS/LES Simulations of Separated Flows at High Reynolds Numbers

Direct numerical simulations (DNS) and large-eddy simulations (LES) simulations of flow through a periodic channel with a constriction are performed using the dynamic Smagorinsky model at two Reynolds numbers of 2800 and 10595. The LES equations are solved using higher order compact schemes. DNS are performed for the lower Reynolds number case using a fine grid and the data are used to validate the LES results obtained with a coarse and a medium size grid. LES simulations are also performed for the higher Reynolds number case using a coarse and a medium size grid. The results are compared with an existing reference data set. The DNS and LES results agreed well with the reference data. Reynolds stresses, sub-grid eddy viscosity, and the budgets for the turbulent kinetic energy are also presented. It is found that the turbulent fluctuations in the normal and spanwise directions have the same magnitude. The turbulent kinetic energy budget shows that the production peaks near the separation point region and the production to dissipation ratio is very high on the order of five in this region. It is also observed that the production is balanced by the advection, diffusion, and dissipation in the shear layer region. The dominant term is the turbulent diffusion that is about two times the molecular dissipation.

Balakumar, P.↗

Planet Hunters TESS IV: A Massive, Compact Hierarchical Triple Star System TIC 470710327

We report the discovery and analysis of a massive, compact, hierarchical triple system (TIC 470710327) initially identified by citizen scientists in data obtained by NASA’s Transiting Exoplanet Survey Satellite (TESS). Spectroscopic follow-up observations obtained with the HERMES spectrograph, combined with eclipse-timing variations (ETVs), confirm that the system is comprised of three OB stars, with a compact 1.10 d eclipsing binary and a non-eclipsing tertiary on a 52.04 d orbit. Dynamical modelling of the system (from radial velocity and ETVs) reveal a rare configuration wherein the tertiary star (O9.5-B0.5V; 14–17 M ꙩ ) is more massive than the combined mass of the inner binary (10.9–13.2 M ꙩ ). Given the high mass of the tertiary, we predict that this system will undergo multiple phases of mass transfer in the future, and likely end up as a double neutron star gravitational wave progenitor or an exotic Thorne–Żytkow object. Further observational characterization of this system promises constraints on both formation scenarios of massive stars as well as their exotic evolutionary end-products.

NASA's Transiting Exoplanet Survey Satellite (TESS↗

AGR-5/6/7 Post-Irradiation Examination Plan

This plan describes the PIE activities to be carried out on AGR-5/6/7 fuel and the non-fuel components of the AGR-5/6/7 irradiation test train. PIE of AGR-5/6/7 fuel and test train components will provide data on fuel performance over a range of irradiation conditions (e.g., burnup, neutron fluence, and temperature), test fuel under postulated accident conditions, and support development of fuel performance and fission product transport models. Examination of the non-fuel test train components will give information on TRISO fuel and compact performance and information on releases of radioactive fission products from the fuel compacts. PIE of the fuel compacts is used for many purposes including, but not limited to: updating the irradiation thermal calculations, measuring burnup, determining retention/release of fission products from the compacts and the TRISO particles, determining TRISO layer failure rates, determining fuel performance at elevated accident temperatures, determining fuel performance under oxidizing atmospheres, determining the distribution of fission products within the compact and fuel particles, and observing compact and particle morphology to explain mechanisms of TRISO layer failure and/or fission product release/retention.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Statistically-informed deep learning for gravitational wave parameter estimation

We introduce deep learning models to estimate the masses of the binary components of black hole mergers, $(m_1,m_2)$, and three astrophysical properties of the post-merger compact remnant, namely, the final spin, $a_\mathrm f$, and the frequency and damping time of the ringdown oscillations of the fundamental $\ell = m = 2$ bar mode, $(\omega_\mathrm R, \omega_\mathrm I)$. Our neural networks combine a modified WaveNet architecture with contrastive learning and normalizing flow. We validate these models against a Gaussian conjugate prior family whose posterior distribution is described by a closed analytical expression. Upon confirming that our models produce statistically consistent results, we used them to estimate the astrophysical parameters $(m_1,m_2, a_\mathrm f, \omega_\mathrm R, \omega_\mathrm I)$ of five binary black holes: GW150914, GW170104, GW170814, GW190521 and GW190630. We use PyCBC Inference to directly compare traditional Bayesian methodologies for parameter estimation with our deep learning based posterior distributions. Our results show that our neural network models predict posterior distributions that encode physical correlations, and that our data-driven median results and 90% confidence intervals are similar to those produced with gravitational wave Bayesian analyses. This methodology requires a single V100 NVIDIA GPU to produce median values and posterior distributions within two milliseconds for each event. Furthermore, this neural network, and a tutorial for its use, are available at the Data and Learning Hub for Science.

79 ASTRONOMY AND ASTROPHYSICS↗

Graph neural networks predict energetic and mechanical properties for models of solid solution metal alloy phases

Here, we developed a PyTorch-based architecture called HydraGNN that implements graph convolutional neural networks (GCNNs) to predict the formation energy and the bulk modulus for models of solid solution alloys for various atomic crystal structures and relaxed volumes. We trained the GCNN surrogate model on a dataset for nickel–niobium (NiNb) generated by the embedded atom model (EAM) empirical interatomic potential for demonstration purposes. The dataset was generated by calculating the formation energy and the bulk modulus as a prototypical elastic property for optimized geometries starting from initial body-centered cubic (BCC), face-centered cubic (FCC), and hexagonal compact packed (HCP) crystal structures, with configurations spanning the possible compositional range for each of the three types of initial crystal structures. Numerical results show that the GCNN model effectively predicts both the formation energy and the bulk modulus as function of the optimized crystal structure, relaxed volume, and configurational entropy of the model structures for solid solution alloys.

36 MATERIALS SCIENCE↗

Non-nebular Origin of Dark Mantles Around Chondrules and Inclusions in CM Chondrites

Our examination of nine CM chondrites that span the aqueous alteration sequence leads us to conclude that compact dark fine mantles surrounding chondrules and inclusions in CM chondrites are not discrete fine-grained rims acquired in the solar nebula as modeled by Metzler et al. [Accretionary dust mantles in CM chondrites: evidence for solar nebula processes. Geochim. Cosmochim. Acta 56, 1992, 2873-28971. Nebular processes that lead to agglomeration produce materials with porosities far higher than those in the dark mantles. We infer that the mantles were produced from porous nebular materials on the CM parent asteroid by impact-compaction (a process that produces the lowest porosity adjacent to chondrules and inclusions). Compaction was followed by aqueous alteration that formed tochilinite, serpentine, Ni-bearing sulfide, and other secondary products in voids in the interchondrule regions. Metzler et al. reported a correlation between mantle thickness and the radius of the enclosed object. In Yamato 791 198 we find no correlation when all sizes of central objects and dark lumps are included but a significant correlation (r(sup 2) = 0.44) if we limit consideration to central objects with radii >35 microns; a moderate correlation is also found in QUE 97990. We suggest that impact-induced shear of a plum-pudding-like precursor produced the observed "mantles"; these were shielded from comminution during impact events by the adjacent stronger chondrules and inclusions. Some mantles in CM chondrites with low degrees of alteration show distinct layers that may largely reflect differences in porosity. Typically, a gray, uniform inner layer is surrounded by an outer layer consisting of darker silicates with BSE-bright speckles. The CM-chondrite objects characterized as "primary accretionary rocks" by Metzler et al. did not form in the nebula, but rather on the parent body. The absence of solar-flare particle tracks and solar-wind-implanted rare gases in these clasts reflect their lithified nature and low surface/volume ratios during the period when they resided in the regolith and were subject to irradiation by solar particles. The clasts are analogous to the light-colored metamorphosed clasts in ordinary-chondrite regolith breccias (which also lack solar-flare particle tracks and solar-wind gas).

Trigo-Rodriquez, Josep M.↗

A Coupled Aeroelastic Model for Launch Vehicle Stability Analysis

A technique for incorporating distributed aerodynamic normal forces and aeroelastic coupling effects into a stability analysis model of a launch vehicle is presented. The formulation augments the linear state-space launch vehicle plant dynamics that are compactly derived as a system of coupled linear differential equations representing small angular and translational perturbations of the rigid body, nozzle, and sloshing propellant coupled with normal vibration of a set of orthogonal modes. The interaction of generalized forces due to aeroelastic coupling and thrust can be expressed as a set of augmenting non-diagonal stiffness and damping matrices in modal coordinates with no penalty on system order. While the eigenvalues of the structural response in the presence of thrust and aeroelastic forcing can be predicted at a given flight condition independent of the remaining degrees of freedom, the coupled model provides confidence in closed-loop stability in the presence of rigid-body, slosh, and actuator dynamics. Simulation results are presented that characterize the coupled dynamic response of the Ares I launch vehicle and the impact of aeroelasticity on control system stability margins.

Orr, Jeb S.↗

The Muon Collider

Muons offer a unique opportunity to build a compact high-energy electroweak collider at the 10 TeV scale. A Muon Collider enables direct access to the underlying simplicity of the Standard Model and unparalleled reach beyond it. It will be a paradigm-shifting tool for particle physics representing the first collider to combine the high-energy reach of a proton collider and the high precision of an electron-positron collider, yielding a physics potential significantly greater than the sum of its individual parts. A high-energy muon collider is the natural next step in the exploration of fundamental physics after the HL-LHC and a natural complement to a future low-energy Higgs factory. Such a facility would significantly broaden the scope of particle colliders, engaging the many frontiers of the high energy community. The last European Strategy for Particle Physics Update and later the Particle Physics Project Prioritisation Panel in the US requested a study of the muon collider, which is being carried on by the International Muon Collider Collaboration. In this comprehensive document we present the physics case, the state of the work on accelerator design and technology, and propose an R\&D project that can make the muon collider a reality.

Accelerator Physics (physics.acc-ph)↗

Prediction of strain, inter-layer interaction and critical current in CORC ® wires under axial strain by T-A modeling

Superconducting conductors on round core (CORC ® ) cables and wires can meet the needs of large high-field magnets, such as particle accelerators and compact nuclear fusion machines, due to their simple cabling process, high current-carrying capacity and reliable operation under high mechanical stresses. Many high-field magnets require CORC ® cables to carry a current of thousands of amperes in a background magnetic field exceeding 20 T. As a result, the large electromagnetic forces will deform the cable in the axial direction due to hoop stress and in the transverse direction by compressive stress. Therefore, it is essential to determine the irreversible deformation limit of the CORC ® cable under axial tensile load and optimize the cabling parameters to potentially extend this limit. Analytical and numerical methods are developed to assess the performance degradation of CORC ® wires under axial tensile load. The strain level, interlayer contact pressure and friction and their impact on the critical current are calculated by combining the mechanical response and the T-A method. Analyzing the results shows that the winding angle of the tape and the Poisson’s ratio of the inner core are key factors affecting the irreversible tensile strain limit of CORC ® wires. The smaller the winding angle and the higher the Poisson’s ratio of the inner core, the higher the irreversible tensile strain limit. For multi-layer CORC ® wires, the initial contact pressure caused by the cabling process must also be considered. The inter-layer interaction is coupled with the tape strain of each layer. The results of this research can serve as a basis for optimizing and designing CORC ® wires with extended irreversible strain limits.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Computations of Complex Three-Dimensional Turbulent Free Jets

Three-dimensional, incompressible turbulent jets with rectangular and elliptical cross-sections are simulated with a finite-difference numerical method. The full Navier- Stokes equations are solved at low Reynolds numbers, whereas at high Reynolds numbers filtered forms of the equations are solved along with a sub-grid scale model to approximate the effects of the unresolved scales. A 2-N storage, third-order Runge-Kutta scheme is used for temporary discretization and a fourth-order compact scheme is used for spatial discretization. Although such methods are widely used in the simulation of compressible flows, the lack of an evolution equation for pressure or density presents particular difficulty in incompressible flows. The pressure-velocity coupling must be established indirectly. It is achieved, in this study, through a Poisson equation which is solved by a compact scheme of the same order of accuracy. The numerical formulation is validated and the dispersion and dissipation errors are documented by the solution of a wide range of benchmark problems. Three-dimensional computations are performed for different inlet conditions which model the naturally developing and forced jets. The experimentally observed phenomenon of axis-switching is captured in the numerical simulation, and it is confirmed through flow visualization that this is based on self-induction of the vorticity field. Statistical quantities such as mean velocity, mean pressure, two-point velocity spatial correlations and Reynolds stresses are presented. Detailed budgets of the mean momentum and Reynolds stresses are presented. Detailed budgets of the mean momentum and Reynolds stress equations are presented to aid in the turbulence modeling of complex jets. Simulations of circular jets are used to quantify the effect of the non-uniform curvature of the non-circular jets.

Wilson, Robert V.↗

Relativistic fireballs and their impact on external matter - Models for cosmological gamma-ray bursts

We discuss the production of cosmological gamma-ray bursts intense enough to be detected at cosmological distances. Events such as the coalescence of compact binaries can create sufficient energy on time scales much less than 1 s. A short 'primary' burst is expected when the resultant fireball becomes optically thin, but this may be weak because the bulk of the radiative energy has been converted into kinetic energy while still trapped within the fireball. But when this expanding material impacts on an external medium, its bulk kinetic energy can be rerandomized. The requirements on the composition of the fireball itself are less stringent than for other interpretations of cosmological gamma-ray bursts. Moreover, our model suggests that the spectra and time structure of the bursts may depend in interesting ways on the environment in which the energy-generating event occurs.

Meszaros, P.↗

An Isolated Mass-gap Black Hole or Neutron Star Detected with Astrometric Microlensing

We present the analysis of five black hole candidates identified from gravitational microlensing surveys. Hubble Space Telescope astrometric data and densely sampled light curves from ground-based microlensing surveys are fit with a single-source, single-lens microlensing model in order to measure the mass and luminosity of each lens and determine if it is a black hole. One of the five targets (OGLE-2011-BLG-0462/MOA-2011-BLG-191 or OB110462 for short) shows a significant >1 mas coherent astrometric shift, little to no lens flux, and has an inferred lens mass of 1.6–4.4 M ⨀ . This makes OB110462 the first definitive discovery of a compact object through astrometric microlensing and it is most likely either a neutron star or a low-mass black hole. This compact-object lens is relatively nearby (0.70–1.92 kpc) and has a slow transverse motion of <30 kms -1 . OB110462 shows significant tension between models well fit to photometry versus astrometry, making it currently difficult to distinguish between a neutron star and a black hole. Additional observations and modeling with more complex system geometries, such as binary sources, are needed to resolve the puzzling nature of this object. For the remaining four candidates, the lens masses are <2M ⨀ , and they are unlikely to be black holes; two of the four are likely white dwarfs or neutron stars. We compare the full sample of five candidates to theoretical expectations on the number of black holes in the Milky Way (∼10 8 ) and find reasonable agreement given the small sample size.

Astrophysical black holes↗

Testing and Characterization to Develop a Mechanistic Explanation for Unsaturated Drift of Fiber Optic Sensors during High-Dose Irradiation

The primary limitation for any optical fiber-based sensor for nuclear reactor applications is radiation-induced attenuation (RIA) of the transmitted and/or reflected signals. Based on several recent studies, RIA is tolerable for some fused silica optical fibers with the proper choice of sensing wavelength and fiber dopants. For extreme temperature applications (> 1000°C), sapphire optical fibers have been proposed; however, recent optical transmission measurements performed on bulk sapphire samples showed prohibitively large RIA. For some sensors, radiation-induced dimensional changes in the fiber materials can also cause significant drift. Moreover, the drift that was observed in numerous experiments performed in the High Flux Isotope Reactor (HFIR), the Advanced Test Reactor, the Massachusetts Institute of Technology Reactor, and other international facilities far exceeded what would be expected based on compaction of fused silica glass. Clearly, additional work is needed to better understand the origins of both RIA and radiation-induced drift in both silica and sapphire optical fiber-based sensors before these sensors can be reliably deployed for nuclear applications. This work evaluated the underlying mechanisms that may be responsible for RIA and drift in silica and sapphire materials. First, detailed characterization was performed on bulk fused silica glass samples that were previously irradiated to different neutron fluences at different temperatures to better understand the structural changes that drive radiation-induced drift in the absence of coating effects that are discussed later. Results show that the non-monotonic compaction that occurs with increasing neutron fluence continues up to fast neutron fluences approaching 10 22 n/cm 2 , which has important implications for physics-based models that may be used to compensate for the sensor drift. Initial Raman spectroscopy and synchrotron x-ray diffraction provide insights into the nature of the structural changes. Next, detailed characterizations were performed on silica fibers with various coatings that were subjected to several different thermal treatments. The hypothesis is that the coatings convert to carbon-rich materials that compact under irradiation, putting a large compressive strain on the fiber. Out-of-pile testing confirms that both polyimide and acrylate fiber coatings convert to glassy carbon (GC) materials when heated under inert conditions, and the degree of order (i.e., graphitization) increases with increasing temperature. The results provide increasingly strong evidence that the combination of polymeric coatings and inert (or vacuum) conditions render fiber optic sensors susceptible to significant radiation-induced drift that would not otherwise exist in uncoated fibers. Finally, transmission electron microscopy was performed on bulk sapphire samples that were irradiated to two neutron fluences at different temperatures to gain insights into the potential mechanisms driving the prohibitive RIA at higher neutron fluences and temperatures. Contrary to previous hypotheses, results show that scattering from radiation-induced voids cannot explain the observed RIA. Similarly, models for scattering losses from dislocation loops also do not agree with the experimental results. Instead, fitting to the experimental data shows that increased absorption from aluminum vacancy centers is the most likely explanation for the the prohibitively large RIA that was observed at high irradiation temperature and dose. In addition, the voids that formed in these single-crystal samples were found to align along the basal plane (a-axis) as opposed to that seen in previous observations of c-axis alignment in polycrystalline samples, which could have important implications for anisotropic swelling and other phenomena that could affect sensor performance at high neutron fluence.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Testing and Characterization to Develop a Mechanistic Explanation for Unsaturated Drift of Fiber Optic Sensors during High-Dose Irradiation

The primary limitation for any optical fiber-based sensor for nuclear reactor applications is radiation-induced attenuation (RIA) of the transmitted and/or reflected signals. Based on several recent studies, RIA is tolerable for some fused silica optical fibers with the proper choice of sensing wavelength and fiber dopants. For extreme temperature applications (> 1000°C), sapphire optical fibers have been proposed; however, recent optical transmission measurements performed on bulk sapphire samples showed prohibitively large RIA. For some sensors, radiation-induced dimensional changes in the fiber materials can also cause significant drift. Moreover, the drift that was observed in numerous experiments performed in the High Flux Isotope Reactor (HFIR), the Advanced Test Reactor, the Massachusetts Institute of Technology Reactor, and other international facilities far exceeded what would be expected based on compaction of fused silica glass. Clearly, additional work is needed to better understand the origins of both RIA and radiation-induced drift in both silica and sapphire optical fiber-based sensors before these sensors can be reliably deployed for nuclear applications. This work evaluated the underlying mechanisms that may be responsible for RIA and drift in silica and sapphire materials. First, detailed characterization was performed on bulk fused silica glass samples that were previously irradiated to different neutron fluences at different temperatures to better understand the structural changes that drive radiation-induced drift in the absence of coating effects that are discussed later. Results show that the non-monotonic compaction that occurs with increasing neutron fluence continues up to fast neutron fluences approaching 10 22 n/cm 2 , which has important implications for physics-based models that may be used to compensate for the sensor drift. Initial Raman spectroscopy and synchrotron x-ray diffraction provide insights into the nature of the structural changes. Next, detailed characterizations were performed on silica fibers with various coatings that were subjected to several different thermal treatments. The hypothesis is that the coatings convert to carbon-rich materials that compact under irradiation, putting a large compressive strain on the fiber. Out-of-pile testing confirms that both polyimide and acrylate fiber coatings convert to glassy carbon (GC) materials when heated under inert conditions, and the degree of order (i.e., graphitization) increases with increasing temperature. The results provide increasingly strong evidence that the combination of polymeric coatings and inert (or vacuum) conditions render fiber optic sensors susceptible to significant radiation-induced drift that would not otherwise exist in uncoated fibers. Finally, transmission electron microscopy was performed on bulk sapphire samples that were irradiated to two neutron fluences at different temperatures to gain insights into the potential mechanisms driving the prohibitive RIA at higher neutron fluences and temperatures. Contrary to previous hypotheses, results show that scattering from radiation-induced voids cannot explain the observed RIA. Similarly, models for scattering losses from dislocation loops also do not agree with the experimental results. Instead, fitting to the experimental data shows that increased absorption from aluminum vacancy centers is the most likely explanation for the the prohibitively large RIA that was observed at high irradiation temperature and dose. In addition, the voids that formed in these single-crystal samples were found to align along the basal plane (a-axis) as opposed to that seen in previous observations of c-axis alignment in polycrystalline samples, which could have important implications for anisotropic swelling and other phenomena that could affect sensor performance at high neutron fluence.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Plasma environment effects on K lines of astrophysical interest IV. IPs, K thresholds, radiative rates, and Auger widths in Fe ii–Fe viii

Aims. Within the framework of compact-object accretion disks, we calculate plasma environment e ects on the atomic structure and decay parameters used in the modeling of K lines in lowly charged iron ions, namely Fe ii–Fe viii. Methods. For this study, we used the fully relativistic multiconfiguration Dirac–Fock method approximating the plasma electron–nucleus and electron-electron screenings with a time-averaged Debye-Hückel potential. Results. We report modified ionization potentials, K-threshold energies, wavelengths, radiative emission rates, and Auger widths for plasmas characterized by electron temperatures and densities in the ranges 10(exp 5)-10(exp 7) K and 10(exp 18)-10(exp 22) per cu.cm. In addition, we propose two universal fitting formulae to predict the IP and K-threshold lowerings in any elemental ion. Conclusions. We conclude that the high-resolution X-ray spectrometers onboard the future XRISM and ATHENA space missions will be able to detect the lowering of the K edges of these Fe ions due to the extreme plasma conditions occurring in the accretion disks around compact objects.

J. Deprince↗

Influence of Shear Strength Assumptions on BISON Debonding Simulations

Accurately predicting the thermomechanical response of buffer–IPyC debonding in TRISO fuel particles requires reliable mechanical property inputs for each coating layer, particularly the normal and shear strengths that influence interlayer delamination and stress concentrations. Micro tensile testing of AGR-2 fuel particles provided experimentally measured normal strengths for the buffer, IPyC, and buffer–IPyC interface; however, shear strength was not measured. As a result, BISON simulations of interface debonding must rely on assumed shear strength values, typically estimated as 20–40% of the measured ultimate tensile strength. This study evaluates how these assumed shear strength values influence cohesive zone model (CZM) predictions of buffer–IPyC separation in AGR 2 TRISO particles. Using micro tensile data from three AGR 2 compacts (2 1 3, 5 1 3, and 6 3 3), BISON simulations were performed with multiple shear strength assumptions to quantify their effect on radial and tangential stress evolution, debonding, and gap propagation. The results show that shear strength is a high sensitivity parameter: increasing the assumed shear strength significantly alters the stress distribution at the buffer–IPyC junction, shifts the predicted debonding location, and changes the extent of partial gap formation. While normal strength controls the initiation of interface separation, shear strength strongly influences the mode mixity of the failure process and the resulting stress concentrations transmitted to the IPyC and SiC layers. These findings highlight a critical gap in current TRISO mechanical characterization. Without experimentally measured shear strength, BISON simulations must rely on approximations that introduce uncertainty into predictions of coating layer integrity and fission product barrier performance. Future fuel qualification campaigns should therefore consider measurement of shear strength at the interlayer interfaces to reduce model uncertainty and improve the fidelity of TRISO fuel performance simulations.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Quantum simulation of bosons with the contracted quantum eigensolver

Abstract Quantum computers are promising tools for simulating many-body quantum systems due to their potential scaling advantage over classical computers. While significant effort has been expended on many-fermion systems, here we simulate a model entangled many-boson system with the contracted quantum eigensolver (CQE). We generalize the CQE to many-boson systems by encoding the bosonic wavefunction on qubits. The CQE provides a compact ansatz for the bosonic wave function whose gradient is proportional to the residual of a contracted Schrödinger equation. We apply the CQE to a bosonic system, where N quantum harmonic oscillators are coupled through a pairwise quadratic repulsion. The model is relevant to the study of coupled vibrations in molecular systems on quantum devices. Results demonstrate the potential efficiency of the CQE in simulating bosonic processes such as molecular vibrations with good accuracy and convergence even in the presence of noise.

Physics↗

Pore‐Scale Modeling of Reactive Transport with Coupled Mineral Dissolution and Precipitation

Abstract We present a new pore‐scale model for multicomponent advective‐diffusive transport with coupled mineral dissolution and precipitation. Both dissolution and precipitation are captured simultaneously by introducing a phase transformation vector field representing the direction and magnitude of the overall phase change. An effective viscosity model is adopted in simulating fluid flow during mineral dissolution‐precipitation that can accurately capture the velocity field without introducing any empirical parameters. The proposed approach is validated against analytical solutions and interface tracking simulations in simplified structures. After validation, the proposed approach is employed in modeling realistic rocks where mineral dissolution and precipitation are dominant at different locations. We have identified three regimes for mineral dissolution‐precipitation coupling: (a) compact dissolution‐precipitation where dissolution is dominant near the inlet and precipitation is dominant near the outlet, (b) wormhole dissolution with clustered precipitation where dissolution generates wormholes in the main flow paths and precipitation clogs the secondary flow paths, and (c) dissolution dominant where all solid grains are gradually dissolved. In the three regimes, the proposed approach provides reliable porosity‐permeability relationships that cannot be described well by traditional macroscale models. We find that the permeability can increase while the overall porosity decreases when the main flow paths are expanded by dissolution and adjacent pore spaces are clogged by precipitation.

58 GEOSCIENCES↗