Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Equations of states”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

Use of Generalized Fluid System Simulation Program (GFSSP) for Teaching and Performing Senior Design Projects at the Educational Institutions

This paper describes the experience of the authors in using the Generalized Fluid System Simulation Program (GFSSP) in teaching Design of Thermal Systems class at University of Alabama in Huntsville. GFSSP is a finite volume based thermo-fluid system network analysis code, developed at NASA/Marshall Space Flight Center, and is extensively used in NASA, Department of Defense, and aerospace industries for propulsion system design, analysis, and performance evaluation. The educational version of GFSSP is freely available to all US higher education institutions. The main purpose of the paper is to illustrate the utilization of this user-friendly code for the thermal systems design and fluid engineering courses and to encourage the instructors to utilize the code for the class assignments as well as senior design projects. The need for a generalized computer program for thermofluid analysis in a flow network has been felt for a long time in aerospace industries. Designers of thermofluid systems often need to know pressures, temperatures, flow rates, concentrations, and heat transfer rates at different parts of a flow circuit for steady state or transient conditions. Such applications occur in propulsion systems for tank pressurization, internal flow analysis of rocket engine turbopumps, chilldown of cryogenic tanks and transfer lines, and many other applications of gas-liquid systems involving fluid transients and conjugate heat and mass transfer. Computer resource requirements to perform time-dependent, three-dimensional Navier-Stokes computational fluid dynamic (CFD) analysis of such systems are prohibitive and therefore are not practical. Available commercial codes are generally suitable for steady state, single-phase incompressible flow. Because of the proprietary nature of such codes, it is not possible to extend their capability to satisfy the above-mentioned needs. Therefore, the Generalized Fluid System Simulation Program (GFSSP1) has been developed at NASA Marshall Space Flight Center (MSFC) as a general fluid flow system solver capable of handling phase changes, compressibility, mixture thermodynamics and transient operations. It also includes the capability to model external body forces such as gravity and centrifugal effects in a complex flow network. The objectives of GFSSP development are: a) to develop a robust and efficient numerical algorithm to solve a system of equations describing a flow network containing phase changes, mixing, and rotation; and b) to implement the algorithm in a structured, easy-to-use computer program. The analysis of thermofluid dynamics in a complex network requires resolution of the system into fluid nodes and branches, and solid nodes and conductors as shown in Figure 1. Figure 1 shows a schematic and GFSSP flow circuit of a counter-flow heat exchanger. Hot nitrogen gas is flowing through a pipe, colder nitrogen is flowing counter to the hot stream in the annulus pipe and heat transfer occurs through metal tubes. The problem considered is to calculate flowrates and temperature distributions in both streams. GFSSP has a unique data structure, as shown in Figure 2, that allows constructing all possible arrangements of a flow network with no limit on the number of elements. The elements of a flow network are boundary nodes where pressure and temperature are specified, internal nodes where pressure and temperature are calculated, and branches where flowrates are calculated. For conjugate heat transfer problems, there are three additional elements: solid node, ambient node, and conductor. The solid and fluid nodes are connected with solid-fluid conductors. GFSSP solves the conservation equations of mass and energy, and equation of state in internal nodes to calculate pressure, temperature and resident mass. The momentum conservation equation is solved in branches to calculate flowrate. It also solves for energy conservation equations to calculate temperatures of solid nodes. The equations are coupled and nonlinear; therefore, they are solved by an iterative numerical scheme. GFSSP employs a unique numerical scheme known as simultaneous adjustment with successive substitution (SASS), which is a combination of Newton-Raphson and successive substitution methods. The mass and momentum conservation equations and the equation of state are solved by the Newton-Raphson method while the conservation of energy and species are solved by the successive substitution method. GFSSP is linked with two thermodynamic property programs, GASP2 and WASP3 and GASPAK4, that provide thermodynamic and thermophysical properties of selected fluids. Both programs cover a range of pressure and temperature that allows fluid properties to be evaluated for liquid, liquid-vapor (saturation), and vapor region. GASP and WASP provide properties of 12 fluids. GASPAK includes a library of 36 fluids. GFSSP has three major parts. The first part is the graphical user interface (GUI), visual thermofluid analyzer of systems and components (VTASC). VTASC allows users to create a flow circuit by a 'point and click' paradigm. It creates the GFSSP input file after the completion of the model building process. GFSSP's GUI provides the users a platform to build and run their models. It also allows post-processing of results. The network flow circuit is first built using three basic elements: boundary node, internal node, and branch.

Majumdar, A. K.↗
DOE OSTI2023

Shock-induced melting of [100] lithium fluoride: Sound speed and Hugoniot measurements to 230 GPa

Although [100] lithium flouride (LiF) is the most widely used optical window material in dynamic compression experiments, its high stress (>100 GPa) shock compression response, including melting, is not well understood. To address this need, we measured wave profiles in plate impact experiments to determine the Hugoniot states and longitudinal sound speeds in [100] LiF crystals shock compressed to 231 GPa. The measured peak states are fitted well by a linear shock velocity-particle velocity relation, providing an accurate determination of the LiF Hugoniot curve to significantly higher stresses than previous experiments. The longitudinal sound speeds show a near linear increase with density compression to 182 GPa. Between 182 GPa and 195 GPa, the sound speed and the longitudinal modulus decrease abruptly, due to shock-induced melting. The increasing sound speeds and moduli at higher stresses suggest that shock compressed LiF is fully liquid at 195 GPa and above, allowing determination of the Gruneisen parameter for liquid LiF. Here, the melt stress determined here differs from that predicted by current multiphase equations of state for LiF. Our results provide important insight into the high stress solid and liquid states of shock compressed LiF and point to the need for an improved multiphase equation of state at high pressures and high temperatures.

36 MATERIALS SCIENCE↗
NASA NTRS1991

Fluid dynamics of multiphase flow in volcanic environments on the Earth, Io, and other planets

Data was compiled on H2O, CO2, S, and SO2 and it was analyzed by various techniques to prepare phase diagrams relevant to conditions on the Earth, Io, Mars, and other planets. The use of Cray models for volcanic eruptions was examined. Attempts are being made to develop equations of state of the compounds listed above to be compatible with the Cray code so that they can be called as subroutines for the volcanic models. This is a major project, as the current code uses only the equation of state for a perfect gas, and requires that the atmosphere be the same composition as the volcanic gas. The analysis of the thermodynamic data is time consuming, and was complicated by changes from a VAX/VMS system to Macintosh's, Crays using CTSS, Crays using Unicos, and Sun workstations.

Kieffer, Susan Werner↗
DOE OSTI2021

Density-Functional Tight-Binding Parameters for Bulk Zirconium: A Case Study for Repulsive Potentials

Density-functional tight-binding (DFTB) parameters are presented for the simulation of the bulk phases of zirconium. Electronic parameters were obtained using a band structure fitting strategy, while two-center repulsive potentials were created by particle swarm optimization. As objective functions for the repulsive potential fitting, we employed the Birch–Murnaghan equations of state for hexagonal close-packed (HCP), body-centered cubic (BCC) and ω phases of Zr from density-functional theory (DFT). When fractional atomic coordinates are not allowed to change in the generation of the equation-of-state curves, long-range repulsive DFTB potentials are able to almost perfectly reproduce equilibrium structures, relative DFT energies of the bulk phases, and bulk moduli. However, the same potentials lead to artifacts in the DFTB potential energy surfaces when atom positions in the unit cell are allowed to fully relax during the change of unit cell parameters. Finally, conventional short-range repulsive DFTB potentials, while inferior in their ability to reproduce DFT bulk energetics, are able to correctly reproduce the qualitative shape of the DFT potential energy surfaces, including the location of global minima, and can therefore be considered more transferable.

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

Validation of Reactive Flow Capabilities in PAGOSA

Several reactive flow modeling capabilities are studied and applied to physics problems of interest in the PAGOSA hydrocode with an emphasis on validation. The Davis equation of state and the Arrhenius shock temperature dependent WSD (AWSD) reactive flow model have been implemented by the author in PAGOSA over the last year. These high-explosive modeling options along with the previously existing HE-JWL equation of state and the SURF/SURFplus reactive flow models are discussed, exercised and compared. Validation of the models analyze numerical results against experimental data from cylinder tests, two-stage gas gun experiments and proton radiography (pRad) experiments. Applications of the modeling options are explored involving both HMX based PBX 9501 and TATB based PBX 9502. The effect of mesh size is considered for prediction of phenomena such as corner turning and dead zones, and the relative code performance of the modeling options is discussed.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
NASA NTRS1987

Outer magnetospheric fluctuations and pulsar timing noise

The Cheng, Ho, and Ruderman (1986) outer-magnetosphere gap model was used to investigate the stability of Crab-type outer magnetosphere gaps for pulsars having the parameter (Omega-square B) similar to that of the Crab pulsar. The Lamb, Pines, and Shaham (1978) fluctuating magnetosphere noise model was applied to the Crab pulsar to examine the type of the equation of state that best describes the structure of the neutron star. The noise model was also applied to other pulsars, and the theoretical results were compared with observational data. The results of the comparison are consistent with the stiff equation of state, as suggested by the vortex creep model of the neutron star interior. The timing-noise observations also contribute to the evidence for the existence of superfluid in the core of the neutron star.

Cheng, K. S.↗
DOE OSTI2023

Dark Energy Survey Year 3 results: Constraints on extensions to Λ CDM with weak lensing and galaxy clustering

We constrain six possible extensions to the Λ cold dark matter (CDM) model using measurements from the Dark Energy Survey’s first three years of observations, alone and in combination with external cosmological probes. The DES data are the two-point correlation functions of weak gravitational lensing, galaxy clustering, and their cross-correlation. We use simulated data vectors and blind analyses of real data to validate the robustness of our results to astrophysical and modeling systematic errors. In many cases, constraining power is limited by the absence of theoretical predictions beyond the linear regime that are reliable at our required precision. The Λ CDM extensions are dark energy with a time-dependent equation of state, nonzero spatial curvature, additional relativistic degrees of freedom, sterile neutrinos with eV-scale mass, modifications of gravitational physics, and a binned σ 8 ( z ) model which serves as a phenomenological probe of structure growth. For the time-varying dark energy equation of state evaluated at the pivot redshift we find ( w p , w a ) = ( - 0.9 9 - 0.17 + 0.28 , - 0.9 ± 1.2 ) at 68% confidence with z p = 0.24 from the DES measurements alone, and ( w p , w a ) = ( - 1.0 3 - 0.03 + 0.04 , - 0. 4 - 0.3 + 0.4 ) with z p = 0.21 for the combination of all data considered. Curvature constraints of Ω k = 0.0009 ± 0.0017 and effective relativistic species N eff = 3.1 0 - 0.16 + 0.15 are dominated by external data, though adding DES information to external low-redshift probes tightens the Ω k constraints that can be made without cosmic microwave background observables by 20%. For massive sterile neutrinos, DES combined with external data improves the upper bound on the mass m eff by a factor of 3 compared to previous analyses, giving 95% limits of ( Δ N eff , m eff ) ≤ ( 0.28 , 0.20 eV ) when using priors matching a comparable Planck analysis. For modified gravity, we constrain changes to the lensing and Poisson equations controlled by functions Σ ( k , z ) = Σ 0 Ω Λ ( z ) / Ω Λ , 0 and μ ( k , z ) = μ 0 Ω Λ ( z ) / Ω Λ , 0 , respectively, to Σ 0 = 0. 6 - 0.5 + 0.4 from DES alone and ( Σ 0 , μ 0 ) = ( 0.04 ± 0.05 , 0.0 8 - 0.19 + 0.21 ) for the combination of all data, both at 68% confidence. Overall, we find no significant evidence for physics beyond Λ CDM .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

An Eulerian multimaterial framework for simulating high-explosive aquarium tests

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

42 ENGINEERING↗
NASA NTRS1973

Rotating neutron star models and pulsars

Recently developed equations of state are used to construct fully relativistic models of rotating neutron stars. All of these stable neutron star models are bound. Consideration of binding energies suggests that very-low-mass neutron stars will not be formed in supernova explosions. Some criteria for equations of state are derived from a comparison with observations. Limits for the physical parameters of the Crab pulsar are obtained. In addition, we find that cosmic-ray protons cannot be produced by pulsars.

Boerner, G.↗
NASA NTRS1991

Efficient real gas upwinded Navier-Stokes computations of high speed flows

An efficient method to account for the chemically frozen thermodynamic and transport properties of air in three-dimensional Navier-Stokes calculations has been demonstrated. This approach uses an explicitly specified equation of state (EOS) so that the fluid pressure, temperature, and transport properties are directly related to the flow variables. The method is efficient since no subiterations are required to deduce the pressure and temperature from the flux variables and is modular by allowing different equations of state to be easily supplied to the code. The flexibility of the EOS approach is shown by its implementation into a high-order total variation diminishing upwinding scheme as well as a standard central-differencing scheme. The EOS approach is then demonstrated by computing the hypersonic flow through the corner region of two mutually perpendicular flat plates using both the upwind and central-differencing schemes.

Coirier, William J.↗
NASA NTRS1961

Experimental Determination of the Recovery Factor and Analytical Solution of the Conical Flow Field for a 20 deg Included Angle Cone at Mach Numbers of 4.6 and 6.0 and Stagnation Temperatures to 2600 degree R

The local recovery factor was determined experimentally along the surface of a thin-walled 20 deg included angle cone for Mach numbers near 6.0 at stagnation temperatures between 1200 deg R and 2600 deg R. In addition, a similar cone configuration was tested at Mach numbers near 4.5 at stagnation temperatures of approximately 612 deg R. The local Reynolds number based on flow properties at the edge of the boundary layer ranged between 0.1 x 10(exp 4) and 3.5 x 10(exp 4) for tests at temperatures above 1200 deg R and between 6 x 10(exp 4) and 25 x 10(exp 4) for tests at temperatures near 612 deg R. The results indicated, generally, that the recovery factor can be predicted satisfactorily using the square root of the Prandtl number. No conclusion could be made as to the necessity of evaluating the Prandtl number at a reference temperature given by an empirical equation, as opposed to evaluating the Prandtl number at the wall temperature or static temperature of the gas at the cone surface. For the tests at temperatures above 1200 deg R (indicated herein as the tests conducted in the slip-flow region), two definite trends in the recovery data were observed - one of increasing recovery factor with decreasing stagnation pressure, which was associated with slip-flow effects and one of decreasing recovery factor with increasing temperature. The true cause of the latter trend could not be ascertained, but it was shown that this trend was not appreciably altered by the sources of error of the magnitude considered herein. The real-gas equations of state were used to determine accurately the local stream properties at the outer edge of the boundary layer of the cone. Included in the report, therefore, is a general solution for the conical flow of a real gas using the Beattie-Bridgeman equation of state. The largest effect of temperature was seen to be in the terms which were dependent upon the internal energy of the gas. The pressure and hence the pressure drag terms were unaffected.

Pfyl, Frank A.↗
DOE OSTI2026

Long-Range Transverse-Momentum Correlations and Radial Flow in Pb-Pb Collisions at the LHC

This Letter presents measurements of long-range transverse-momentum correlations using a new observable, 𝑣 0⁡ (𝑝 T ), serving as a probe of event-by-event radial-flow fluctuations, the underlying radial expansion, and the medium’s properties in heavy-ion collisions. Results are reported for inclusive charged particles, pions, kaons, and protons across various centrality intervals in Pb-Pb collisions at $\sqrt{𝑠_{\textrm{NN}}}$ = 5.02 TeV, recorded by the ALICE detector. A pseudorapidity-gap technique, similar to that used in anisotropic-flow studies, is employed to suppress short-range correlations. At low 𝑝 T , a characteristic mass ordering consistent with hydrodynamic collective flow is observed. At higher 𝑝 T (>3 GeV/𝑐), protons exhibit larger 𝑣 0 ⁡(𝑝 T ) than pions and kaons, in agreement with expectations from quark-recombination models. Comparisons to viscous hydrodynamic calculations with varying bulk viscosity and equation of state demonstrate the sensitivity of the 𝑣 0 ⁡(𝑝 T ) observable to these key medium properties. The findings establish 𝑣 0 ⁡(𝑝 T ) as a valuable addition to the set of observables used in Bayesian analyses for extracting the transport properties and constraining the equation of state of strongly interacting matter, while also helping to systematically explore its sensitivity and impact within such global studies.

Acharya, S. [INFN] (ORCID:0000000292135329)↗
DOE OSTI2020

Topaz, a Potential Volatile-Carrier in Cold Subduction Zone: Constraint from Synchrotron X-ray Diffraction and Raman Spectroscopy at High Temperature and High Pressure

The equation of state and stability of topaz at high-pressure/high-temperature conditions have been investigated by in situ synchrotron X-ray diffraction (XRD) and Raman spectroscopy in this study. No phase transition occurs on topaz over the experimental pressure–temperature ($P-T$) range. The pressure–volume data were fitted by the third-order Birch–Murnaghan equation of state (EoS) with the zero-pressure unit–cell volume V 0 = 343.86 (9) Å 3 , the zero-pressure bulk modulus K 0 = 172 (3) GPa, and its pressure derivative K’ 0 = 1.3 (4), while the obtained K 0 = 155 (2) GPa when fixed K’ 0 = 4. In the pressure range of 0–24.4 GPa, the vibration modes of in-plane bending OH-groups for topaz show non-linear changes with the increase in pressure, while the other vibration modes show linear changes. Moreover, the temperature–volume data were fitted by Fei’s thermal equation with the thermal expansion coefficient α 300 = 1.9 (1) × 10 –5 K –1 at 300 K. Finally, the $P-T$ stability of topaz was studied by a synchrotron-based single-crystal XRD at simultaneously high $P-T$ conditions up to ~10.9 GPa and 700 K, which shows that topaz may maintain a metastable state at depths above 370 km in the upper mantle along the coldest subducting slab geotherm. Thus, topaz may be a potential volatile-carrier in the cold subduction zone. It can carry hydrogen and fluorine elements into the deep upper mantle and further affect the geochemical behavior of the upper mantle.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Correlations of conserved quantities at finite baryon density

Correlations involving the seven conserved quantities, namely energy, baryon number, electric charge, strangeness, and the three components of momentum, give rise to correlations in heavy-ion collisions. Through the utilization of a simple one-dimensional hydrodynamic model, we calculate the evolution of the entire 7×7 matrix of correlations as a function of relative spatial rapidity. This comprehensive analysis accounts for finite baryon density, which results in off-diagonal correlations between the charge-related quantities and the energy-momentum quantities. These correlations in coordinate space are subsequently transformed into correlations in momentum space using statistical weighting. As a result, the entire matrix of correlations is revealed to be highly sensitive to the equation of state, viscosity, and diffusivity.

Equations of state of nuclear matter↗
DOE OSTI2020

Two-temperature warm dense hydrogen as a test of quantum protons driven by orbital-free density functional theory electronic forces

We consider a steady-state (but transient) situation in which a warm dense aggregate is a two-temperature system with equilibrium electrons at temperature Te, ions at Ti, and Te ≠ Ti. Such states are achievable by pump–probe experiments. For warm dense hydrogen in such a two-temperature situation, we investigate nuclear quantum effects (NQEs) on structure and thermodynamic properties, thereby delineating the limitations of ordinary ab initio molecular dynamics. We use path integral molecular dynamics (PIMD) simulations driven by orbital-free density functional theory (OFDFT) calculations with state-of-the-art noninteracting free-energy and exchange-correlation functionals for the explicit temperature dependence. We calibrate the OFDFT calculations against conventional (explicit orbitals) Kohn–Sham DFT. We find that when the ratio of the ionic thermal de Broglie wavelength to the mean interionic distance is larger than about 0.30, the ionic radial distribution function is meaningfully affected by the inclusion of NQEs. Moreover, NQEs induce a substantial increase in both the ionic and electronic pressures. This confirms the importance of NQEs for highly accurate equation-of-state data on highly driven hydrogen. For Te > 20 kK, increasing Te in the warm dense hydrogen has slight effects on the ionic radial distribution function and equation of state in the range of densities considered. In addition, we confirm that compared with thermostatted ring-polymer molecular dynamics, the primitive PIMD algorithm overestimates electronic pressures, a consequence of the overly localized ionic description from the primitive scheme.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1981

Theory of dense hydrogen

Equations of state for molecular and metallic hydrogen are calculated in an accurate and comparable manner by using the density-functional method applied to static lattices. By isolating the electronic energies in this way, attention can be drawn to the relative importance of the protonic degrees of freedom. From the equation-of-state results it can be concluded that a remnant molecular pairing is preferred in a band-overlap metallic state, complete dissociation occurring only at very high densities. An accurate determination of the pressure needed to achieve complete dissociation is shown to require a self-consistent treatment of electron and proton degrees of freedom.

Chakravarty, S.↗
DOE OSTI2026

Systematics of the chemical freeze-out line in the high baryon density regime explored at SIS100

The systematic uncertainties of chemical freeze-out fits at SIS100 energies (Au+Au reactions at $\sqrt{s_{NN}}$ = 3 – 5 GeV) are studied using UrQMD simulations. Although hadron production in UrQMD does not occur on a sharp chemical freeze-out hyper-surface, the extracted fit quality is shown to be very good. The extracted chemical parameters depend on the selected hadron species as well as the underlying equation of state (EoS) of the matter. Including light nuclei and anti-protons in the fit increases the expected freeze-out temperature, while a stiffer EoS increases the obtained chemical potential. Similarly, the baryon densities extracted by the thermal fits depend on the choice of hadrons as well as the underlying equation of state. These results are important for the upcoming CBM@FAIR physics program and highlight that a degree of caution is advised when one relates the chemical freeze-out curve to features on the QCD phase diagram like the critical endpoint or a possible phase transition.

Hofmann, Emma Lilith [GSI Helmholtzzentrum für Sch↗
DOE OSTI2020

Implementing Chiral Three-Body Forces in Terms of Medium-Dependent Two-Body Forces

Three-nucleon (3N) forces are an indispensable ingredient for accurate few-body and many- body nuclear structure and reaction theory calculations. While the direct implementation of chiral 3N forces can be technically very challenging, a simpler approach is given by employing instead a medium-dependent NN interaction V med that reflects the physics of three-body forces at the two-body normal-ordered approximation. We review the derivation and construction of Vmed from the chiral 3N interaction at next-to-next-to-leading order (N2LO), consisting of a long-range 2π-exchange term, a mid-range 1π-exchange component and a short-range contact-term. Several applications of V med to the equation of state of cold nuclear and neutron matter, the nucleon single-particle potential in nuclear matter, and the nuclear quasiparticle interaction are discussed. We also explore differences in using local vs. nonlocal regulating functions on 3N forces and make direct comparisons to exact results at low order in perturbation theory expansions for the equation of state and single-particle potential. We end with a discussion and numerical calculation of the in-medium NN potential V med from the next-to-next-to-next-to-leading order (N3LO) chiral 3N force, which consists of a series of long-range and short-range terms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗