Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “confinement”

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 1,765 records · Page 98

Tracer particle transport dynamics in the diffusive sandpile cellular automaton

The confinement properties of the diffusive running sandpile are characterized by tracking the motion of a population of marked grains of sand. It is found that, as the relative strength of the avalanching to the diffusive transport channel is varied, a point is reached at which the particle global confinement time and the probability density functions of the jump-sizes and waiting-times of the tracked grains experience a sudden change, thus revealing a dynamical transition, that is consistent with previous studies (Newman DE et al., Phys Rev Lett 2002;88(20):204304). Across this transition, the sandpile moves from a regime characterized by self-similarity and memory, where avalanches of all possible sizes dominate transport across the system, to another regime where transport is taken over by near system-size, quasi-periodic avalanches. Values for the fractional transport exponents that quantify effective transport across the sandpile prior to the transition are also obtained.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Pore accessibility and trapping of methane in Marcellus shale

Accessibility of pore space in unconventional reservoirs is an important factor influencing both methane storage capacity and the kinetics of methane desorption. The determination of open (accessible) versus closed (inaccessible) porosity is therefore vital for the prediction of gas production potential. This study investigates accessibility of pores to methane in overmature middle Devonian Marcellus shale samples (cut parallel and perpendicular to bedding) using small and ultra-small neutron scattering (SANS and USANS) with contrast matching (CM), supplemented by other complementary techniques, such as mercury injection capillary pressure (MICP) and low pressure gas (N2 and CO2) adsorption. Our results demonstrate that for the samples studied, only about 6 % of pores with diameter 25-500 nm are accessible to methane. The accessibility fraction for pores larger than 500 nm is 35%. For nanopores smaller than 25 nm, pore accessibility could not be quantitatively determined due to increased methane density and condensation effects in confinement. Our observations indicate that methane penetrates the accessible small mesopores and micropores down to at least 1 nm in diameter, and the density of confined deuterated methane (CD4) is 0.68 g/cm3 for pores of diameter 25 nm and gradually increases with the decreased pore size. Moreover, elevated gas pressure causes formation of additional high-density methane nano-clusters. These clusters have a form of slightly anisotropic polydisperse discs oriented along bedding plane, about 1-12 nm in diameter and with average thickness of 3.6 nm. Utilizing samples cut parallel and perpendicular to the bedding, this study also briefly addresses anisotropy of pores. Based on the iso-size intensity ratio R, our SANS and USANS results demonstrate anisotropy in the out-of-bedding direction and suggest that the degree of anisotropy depends on the pore size. Specifically, for pore diameters ~2.5 to 250 nm, the extent of anisotropy is smaller than for pores ~500 nm to 6 µm in diameter. Finally, comparison of pore size distribution results calculated from SANS/USANS to those obtained using MICP shows good agreement at low pressures, but large difference at pressures above 1000 bar. Here, this discrepancy requires further testing; it is possible that the high mercury pressure used in MICP alters the mesopore and micropore structure of shales.

03 NATURAL GAS↗

An Eulerian multimaterial framework for simulating high-explosive aquarium tests

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

42 ENGINEERING↗

Improved stellarator permanent magnet designs through combined discrete and continuous optimizations

A common optimization problem in the areas of magnetized plasmas and fusion energy is the design of magnets to produce a given three-dimensional magnetic field distribution to high precision. When designing arrays of permanent magnets for stellarator plasma confinement, such problems have tens of thousands of degrees of freedom whose solutions, for practical reasons, should be constrained to discrete spaces. We perform a direct comparison between two algorithms that have been developed previously for this purpose, and demonstrate that composite procedures that apply both algorithms in sequence can produce substantially improved results. One approach uses a continuous, quasi-Newton procedure to optimize the dipole moments of a set of magnets and then projects the solution onto a discrete space. The second uses an inherently discrete greedy optimization procedure that has been enhanced and generalized for this work. Further, the approaches are both applied to design arrays cubic rare-Earth permanent magnets to confine a quasi-axisymmetric plasma with a magnetic field on axis of 0.5 T. The first approach tends to find solutions with higher field accuracy, whereas the second can find solutions with substantially (up to 30%) fewer magnets. When the approaches are combined, they can obtain solutions with magnet quantities comparable to the second approach while matching the field accuracy of the first.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental flow-through a single fracture with monolayer proppant at reservoir conditions: A case study on Caney Shale, Southwest Oklahoma, USA

Global energy systems are undergoing a crucial transition to clean energy sources, thereby moving away from fossil-fuel based energy with the ultimate aim of reducing carbon dioxide emissions. It is essential to continue to advance in exploring and developing shale gas resources as this will underpin the global supply of natural gas that is needed to support the energy transition. The challenge hereby is that the slow fluid flow of the shale matrix due to its very low permeability requires large-scale volume fracturing to produce a conductive fracture network that enhances the connectivity between the shale formation and the wellbore. This work examines the influence of a thin proppant layer on a single fracture using a flow-through experiment, and also explores whether proppant embedment in the Caney shale is affected by the rock mineralogy, surface roughness, fluids, confining stress, time, temperature and bedding. The experiment was conducted for a duration of nine days (216 h), experimental temperature was varied from ambient temperature to reservoir temperature of 125 °C (257 °F) and confining stress was varied from 367 psi (2.53MPa) to a maximum of 4011.82 psi (27.66 MPa). Here, we find that the conductivity of the fracture is primarily influenced by the layer of proppant used, surface roughness, mineralogy, fluids, temperature and closure stress.

Caney shale↗

Nonlinear control of the minimum safety factor in tokamaks by optimal allocation of spatially moving electron cyclotron current drive

The minimum value of the safety factor profile is related to the magnetohydrodynamic (MHD) stability of the plasma confined in a tokamak. Therefore, active control of the minimum safety factor may mitigate MHD instabilities that can degrade or even terminate plasma confinement. Typically, in most tokamak scenarios, the minimum safety factor evolves spatially with time, i.e., the location at which the safety factor achieves the minimum value changes with time. In addition to the inherent nonlinearities in the minimum safety factor evolution, its spatial variation makes the control design challenging. In particular, complexity in control design may arise from the need for time-dependent nonlinear models that account for spatial variation of the minimum safety factor. Furthermore, the minimum safety factor may drift to locations where the actuator authority is low. The problem of minimum safety factor control with target location tracking and moving electron cyclotron current drive (ECCD) is addressed in this work. A nonlinear time-dependent model that incorporates the spatial variation of the minimum safety factor is presented. A nonlinear controller based on optimal feedback linearization is developed to track a target minimum safety factor. The proposed controller treats the ECCD position as a controllable variable. In other words, the controller prescribes the ECCD position (in addition to the non-inductive powers) in real time based on an optimal criterion that is defined a priori. This work also presents the steps necessary to integrate the minimum safety factor controller with a total energy controller to achieve multiple control objectives simultaneously. In conclusion, the proposed integrated control algorithm is tested using nonlinear simulations in the Control Oriented Transport SIMulator (COTSIM) for a DIII-D tokamak scenario.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Sodium diffusion in heterogeneous porous media: Connecting laboratory experiments and simulations

Sodium has been suggested as a tracer for brine in reservoir formations where a significant amount of sodium ion is found mainly in the aqueous phase. Fortunately, the most abundant sodium isotope, 23 Na, is Nuclear Magnetic Resonance (NMR) active, making it possible to study the structure and dynamical properties of both bulk and pore confined fluid systems. We investigated the diffusion of key dissolved sodium species in bulk solution and porous media as a function of composition, ionic strength, and permeability of the carbonate samples using NMR relaxometry and pulsed-field gradient (PFG) NMR techniques. We use Indiana limestone as an example of natural porous media and water as a freely diffusible tracer and carrier. We demonstrate apparent diffusion measurement of sodium based on changes in spin–spin relaxation time (T 2 ) signal in pore confinement. The diffusion rate of sodium decreases with increasing counter-ion size. This effect is greater at higher ionic strengths and lower chemical potentials in porous media. The reactive transport code, CrunchFlow, was used to complement the NMR experiments to simulate diffusion behavior in porous media. Combining these two methods provides a powerful approach to estimating effective diffusion coefficients in heterogeneous matrices. The modeling considers the influence of physical properties (porosity and tortuosity) and chemical properties (geochemical composition and chemical gradients). The results presented in this work highlight the advantage of measuring apparent diffusivity using NMR T 2 relaxometry in conjunction with numerical simulation to derive effective diffusivity and the corresponding matrix properties (i.e., tortuosity) of the system.

58 GEOSCIENCES↗

Using high pressure solutions of polyfluoroacrylate and CO 2 to Seal cement cracks for improved wellbore integrity

Polyfluoroacrylate (PFA) is a hydrophobic and oleophobic polymer that is soluble in high pressure carbon dioxide (CO 2 ). In this study, the ability of PFA-CO 2 solutions to greatly reduce the apparent permeability of split or cracked Portland cement cylindrical samples is assessed. The apparent permeability values of confined samples were determined before and after treatment with PFA-CO 2 solutions. In four tests, PFA-CO 2 solutions were continuously displacing pure CO 2 from the cracked cement and the decrease in apparent permeability due to PFA adsorption and wettability alteration was monitored. The lowest apparent permeability cracked cement sample (81 nD) was completely sealed with a very small amount of solution. Samples with initial apparent permeabilities of 89 μD and 29.4 mD exhibited 92% and 99% reductions in permeability, respectively, before the experiments had to be stopped because of the excessively large increase in pressure drop. A 50% reduction in apparent permeability was observed with a 3.80 mD sample. Four other split cement samples (bound together with tape) with an initial apparent permeability in the 9.0–70 mD range were removed from the core holder and immersed in a PFA-CO 2 solution for 24 h to allow for PFA adsorption. Then the PFA-CO 2 solution was depressurized, allowing for the deposition of additional PFA from the solution within the crack as the pressure fell below the cloud point pressure of the PFA-CO 2 solution. These four samples were then confined again in a core holder and apparent permeability reductions of 29–93% were observed. Finally, results from these eight experiments indicates that the more substantial reductions in the nD – mD apparent permeability of the cracked cement correlated to lower initial crack permeability, higher PFA concentration, and slower injection rate of the PFA-CO 2 solution into the crack.

58 GEOSCIENCES↗

Experimental study of fracture slip due to stress perturbation in fractured geo-resources

Anthropogenic activities such as underground fluid injection/extraction, well stimulation, and dam impoundment will perturb the state-of-stress in the subsurface, which could activate critically stressed pre-existing fractures and faults. Fault and fracture activation can induce seismicity, change flow behavior, and damage asperities on the fracture surface to alter reservoir permeability. Here, in this study, triaxial shear tests were conducted on saw cut and rough fracture Barre granite specimens to better characterize how a decrease in normal stress and increase in shear stress will influence crystalline rock fracture hydromechanical properties. In a triaxial set up, confining pressure (CP) and differential stress (DS) were used to induce normal and shear stress on a 45° oriented fracture. Fracture slip and the corresponding stress relaxation was observed when the confining stress decreased. Fracture permeability increased when the fracture slipped. Moreover, new micro-cracks observed by x-ray images indicated the creation of new flow paths for fluid movement. The experimental outcomes of this study indicate that decrease in normal stress due to stress perturbations along the faults/fractures in the extensional stress regime can induce the slip along the fracture. This could have a significant influence on slope stability in quarries, stimulation and injection strategy in deep geo-resources, and potentially minimize the risk of induced seismicity. The observed stress dependency of hydraulic aperture in rough fracture during pre- and post-slip stages could provide key inputs for field-scale simulation of geo-resources.

58 GEOSCIENCES↗

The role of annealing and grain boundary controls on the mechanical properties of limestones and marbles

Chemical and mechanical processes are coupled in many geological and geochemical environments. Reactive processes anneal defects and restructure grain boundaries, modifying their elastic properties, levels of internal friction, wave propagation rates and fracture behaviors. The nature of these changes is, however, contingent on the initial state of the rock. Here, in this study, impulse excitation was used to measure changes in mechanical properties as a function of dry and steam heating time at 300 °C for three carbonate rocks: Carrara marble, Carthage marble (Burlington Limestone), and Texas Cream limestone (Edwards limestone), with initial porosities of about 1 %, 3 %, and 27 %, respectively. Frequency-dependent phenomena along with mineral recrystallization were observed. This was coupled with small-angle X-ray and neutron scattering analysis to determine the relationship between changes in bulk mechanical and microstructural properties. An observed decrease in both the Young's and shear moduli in the experimental limestones as a function of heating time reflects and quantifies a reduction in the stiffness of the rock due to annealing. The internal friction of the samples first increases then decreases with reaction time, reflecting defect annealing, but this was balanced against grain boundary dissolution apparently driven by condensation of steam in the confined grain-boundary environment. This suggests an increase in the boiling point under confinement leading to dissolution, increased porosity and widening of the grain boundaries. In addition, comparison of small-angle X-ray and neutron scattering results suggests that it is inappropriate to assume that pores are empty for quantitative analysis of typical small-angle scattering samples.

58 GEOSCIENCES↗

A non-isothermal breakage-damage model for plastic-bonded granular materials incorporating temperature, pressure, and rate dependencies

Plastic-bonded granular materials (PBM) are widely used in industrial sectors, including building construction, abrasive applications, and defense applications such as plastic-bonded explosives. The mechanical behavior of PBM is highly nonlinear, irreversible, rate dependent, and temperature sensitive governed by various micromechanical attributions such as grain crushing and binder damage. This paper presents a thermodynamically consistent, microstructure-informed constitutive model to capture these characteristic behaviors of PBM. Key features of the model include a breakage internal variable to upscale the grain-scale information to the continuum level and to predict grain size evolution under mechanical loading. In addition, a damage internal state variable is introduced to account for the damage, deterioration, and debonding of the binder matrix upon loading. Temperature is taken as a fundamental external state variable to handle non-isothermal loading paths. The proposed model is able to capture with good accuracy several important aspects of the mechanical properties of PBM, such as pressure-dependent elasticity, pressure-dependent yield strength, brittle-to-ductile transition, temperature dependency, and rate dependency in the post-yielding regime. Furthermore, the model is validated against multiple published datasets obtained from confined and unconfined compression tests, covering various PBM compositions, confining pressures, temperatures, and strain rates.

Breakage↗

Pentaerythritol structure-directing agents bias aluminum siting in MFI zeolites to increase para -xylene selectivity during toluene methylation

MFI zeolites contain ten-membered ring channels (∼0.55 nm diameter) that intersect to form larger voids (∼0.70 nm diameter). The distribution of framework Al sites dictates the void environments that confine active H + sites. MFI synthesized with pentaerythritol (PET) and Na + as structure-directing agents are interrogated using kinetically controlled toluene methylation rates and xylene isomer selectivities. High para-xylene selectivity and low toluene methylation rates indicate that MFI-PET samples contain H + sites predominantly in smaller channels. As-synthesized samples contain more framework Al than Na + , indicating that some charge compensation is provided by protonated PET-solvent complexes that bias Al siting towards these smaller channels.

36 MATERIALS SCIENCE↗

Hybrid bilayered vanadium oxide electrodes with large and tunable interlayer distances in lithium-ion batteries

The interlayer distances in layered electrode materials, influenced by the chemical composition of the confined interlayer regions, have a significant impact on their electrochemical performance. Chemical preintercalation of inorganic metal ions affects the interlayer spacing, yet expansion is limited by the hydrated ion radii. Herein, we demonstrate that using varying concentrations of decyltrimethylammonium (DTA + ) and cetyltrimethylammonium (CTA + ) cations in chemical preintercalation synthesis followed by hydrothermal treatment, the interlayer distance of hybrid bilayered vanadium oxides (BVOs) can be tuned between 11.1 Å and 35.6 Å. Our analyses reveal that these variations in interlayer spacing are due to different amounts of structural water and alkylammonium cations confined within the interlayer regions. Increased concentrations of alkylammonium cations not only expand the interlayer spacing but also induce local bending and disordering of the V-O bilayers. Electrochemical cycling of hybrid BVO electrodes in non-aqueous lithium-ion cells show that specific capacities decrease as interlayer regions expand, suggesting that the densely packed alkylammonium cations obstruct intercalation sites and hinder Li + ion transport. Furthermore, we found that greater layer separation facilitates the dissolution of active material into the electrolyte, resulting in rapid capacity decay during extended cycling. In conclusion, this study emphasizes that layered electrode materials require both spacious interlayer regions as well as high structural and chemical stabilities, providing guidelines for structural engineering of organic–inorganic hybrids.

25 ENERGY STORAGE↗

Nanoconfinement-induced shift in photooxidative degradation pathway of polystyrene

Polymer nanocomposites with high concentrations of nanoparticles (NPs) possess exceptional mechanical, transport, and thermal properties. To enable their widespread use in structural applications and functional coatings, it is crucial to understand how nanoconfinement and the polymer-NP interface influence polymer degradation under various environmental conditions, including prolonged UV exposure. In this study, we investigate the photooxidative degradation of polystyrene (PS)-confined in the interstices of SiO 2 NP films. These nanocomposite films are prepared by the capillary rise infiltration (CaRI) of PS into interstices of SiO 2 NP packings, and subsequently subjected to UV irradiation. Our investigation reveals that PS degradation progresses uniformly across the thickness, with degradation initiating from the center of the NP interstitial pores and extending towards the NP surface. We rationalize this degradation mechanism based on the disparity in the surface energies of PS and the NP surface, as well as slow oxygen diffusion through confined PS. Here, we also demonstrate that packing smaller NPs at a given thickness or thicker packing of a specific NP size facilitates photooxidative degradation, highlighting the critical role of the number of interstitial pores in influencing degradation processes.

36 MATERIALS SCIENCE↗

On numerical errors to the fields surrounding a relativistically moving particle in PIC codes

The particle-in-cell (PIC) method is widely used to model the self-consistent interaction between discrete particles and electromagnetic fields. It has been successfully applied to problems across plasma physics including plasma based acceleration, inertial confinement fusion, magnetically confined fusion, space physics, astrophysics, high energy density plasmas. In many cases the physics involves how relativistic particles (those with high relativistic γ factors) are generated and interact with plasmas. However, when relativistic particles stream across the grid, both in “vacuum” and in plasma, many numerical issues may arise which can lead to unphysical results. We present a detailed analysis of how discretized Maxwell solvers used in PIC codes can lead to numerical errors to the fields that surround particles that move at relativistic speeds across the grid. Expressions for the axial electric field as integrals in k space are presented that reveal two types of errors. The first arises from errors to the numerator of the integrand and leads to unphysical fields that are antisymmetric about the particle. Furthermore, the second arises from errors to the denominator of the integrand and lead to Cherenkov like radiation in “vacuum”. These fields are not anti-symmetric, extend behind the particle, and cause the particle to accelerate or decelerate depending on the solver and parameters. The unphysical fields are studied in detail for two representative solvers - the Yee solver and the FFT based solver. Although the Cherenkov fields are absent, the space charge fields are still present in the fundamental Brillouin zone for the FFT based solvers. In addition, the Cherenkov fields are present in higher order zones for the FFT based solvers. Comparison between the analytical solutions and PIC simulation results are presented. A solution for eliminating these unphysical fields by modifying the k operator in the axial direction is also presented. Using a customized finite difference solver, this solution was successfully implemented into OSIRIS [1]. Results from the customized solver are also presented. Additionally, this solution will be useful for a beam of particles that all move in one direction with a small angular divergence.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Surrogate approximation of the Grad–Shafranov free boundary problem via stochastic collocation on sparse grids

In magnetic confinement fusion devices, the equilibrium configuration of a plasma is determined by the balance between the hydrostatic pressure in the fluid and the magnetic forces generated by an array of external coils and the plasma itself. The location of the plasma is not known a priori and must be obtained as the solution to a free boundary problem. The partial differential equation that determines the behavior of the combined magnetic field depends on a set of physical parameters (location of the coils, intensity of the electric currents going through them, magnetic permeability, etc.) that are subject to uncertainty and variability. The confinement region is in turn a function of these stochastic parameters as well. In this work, we consider variations on the current intensities running through the external coils as the dominant source of uncertainty. This leads to a parameter space of dimension equal to the number of coils in the reactor. With the aid of a surrogate function built on a sparse grid in parameter space, a Monte Carlo strategy is used to explore the effect that stochasticity in the parameters has on important features of the plasma boundary such as the location of the x-point, the strike points, and shaping attributes such as triangularity and elongation. Furthermore, the use of the surrogate function reduces the time required for the Monte Carlo simulations by factors that range between 7 and over 30.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Identifying recharge processes into a vast “fossil” aquifer based on dynamic groundwater 81 Kr age evolution

Water in deep aquifers in arid regions is often considered to be "fossil" when modern recharge rates are negligible relative to the reservoir capacity. Over the past five decades, the Nubian Sandstone Aquifer (NSA) in the arid region of the Sinai Peninsula (Egypt) and the Negev Desert (Israel) has been considered to contain fossil water based on 14 C dating, which revealed 14 C ages of about 30 kyr over most of the aquifer. However, this relatively homogeneous age distribution contradicts the expected increase in groundwater age in the direction of decreasing piezometric head along the flow trajectories. Here, dating results with the longer-lived 81 Kr radioisotope (t 1/2 = 229 ± 11 kyr) are presented, highlighting a wide age range of 40 kyr to 630 kyr in the confined sections of the aquifer, all with very low 14 C activity (<1 pmC). Elevated 81Kr and 14 C activities were only observed within or close to the system's recharge areas. These findings support a new perception of groundwater replenishment during different epochs from the early mid- Pleistocene to the Holocene. By tracking the downstream age evolution, rejuvenation was identified in places where the confinement had been breached. At other locations, the existence of an older groundwater body contributing to the aquifer was detected by means of strongly depleted 81 Kr activity. High spatial heterogeneity in groundwater ages close to the discharge zone of the system is attributed to pronounced age stratification with depth. Calculated ages in the more isolated sections of the system were used to assess regional flow velocity, hydraulic conductivity, and their agreement with present recharge rates. We conclude that regional aquifers in which low 14 C activity suggests homogeneous ancient groundwater age should be reevaluated with 81Kr. With an effective age range beyond one million years, this may enable the restoration of recharge history and provide crucial information for the management of groundwater resources.

54 ENVIRONMENTAL SCIENCES↗

Quarkyonic or baryquark matter

Here, it has been proposed that at high densities nuclear matter will consist of a Fermi sea of quarks surrounded by a small shell of confined baryon at the large momenta, so called Quarkyonic matter. In this contribution we will discuss an alternative configuration, dubbed Baryquark matter, which in a sense is a complement of Quarkyonic matter. Baryquark matter consists of a Fermi sea of confined baryons surrounded by a shell of deconfined quarks. Following Koch and Vovchenko (2023) we will show that for certain (simplified) implementations Baryquark matter is energetically favored over Quarkyonic matter. We will then briefly discuss how the inclusion of the quark structure of nucleons will lead to a configuration which resembles the picture of Quarkyonic matter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗