Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite density”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Effect of XC functionals and dispersion corrections on the DFT‐computed structural and vibrational properties of SrCl 2 –NaCl and ZrF 4 –LiF

Density functional theory (DFT) calculations were performed to examine the impact of exchange–correlation (XC) functionals and van der Waals corrections (specifically the D3 method) on the structural and vibrational properties of the SrCl 2 –NaCl and ZrF 4 –LiF salt systems. Multiple XC functionals, including the local density approximation (LDA), the generalized gradient approximation using the Perdew–Burke–Ernzerhof (PBE) model, and its modified form suitable for solids (PBEsol), the dispersion-corrected PBE-D3 and PBEsol-D3, were considered. Of these functionals, LDA was found to exhibit the highest degree of error, while PBEsol and PBE-D3 displayed the least error. Underestimated lattice parameters compared with experimental values were observed to result in higher force constants, leading to an overprediction of vibrational frequencies. Conversely, an overestimation of lattice parameters was associated with lower vibrational frequencies. The methodology presented in this study yielded results that are in good agreement with experiment, irrespective of the method (finite differences vs. density functional perturbation theory) employed for calculating infrared and Raman spectra. It was further demonstrated that for alkali halides with weak Raman scattering, utilizing a supercell constructed from primitive cells better predicts Raman features than does the use of conventional cells.

Raman↗
NASA NTRS1973

The effect of grid transparency and finite collector size on determining ion temperature and density by the retarding potential analyzer

The analysis of ion data from retarding potential analyzers (RPA's) is generally done under the planar approximation, which assumes that the grid transparency is constant with angle of incidence and that all ions reaching the plane of the collectors are collected. These approximations are not valid for situations in which the ion thermal velocity is comparable to the vehicle velocity, causing ions to enter the RPA with high average transverse velocity. To investigate these effects, the current-voltage curves for H+ at 4000 K were calculated, taking into account the finite collector size and the variation of grid transparency with angle. These curves are then analyzed under the planar approximation. The results show that only small errors in temperature and density are introduced for an RPA with typical dimensions; and that even when the density error is substantial for non-typical dimensions, the temperature error remains minimal.

Troy, B. E., Jr.↗
DOE OSTI2020

Cosmological decoherence from thermal gravitons

A bstract We study the effects of gravitationally-driven decoherence on tunneling processes associated with false vacuum decays, such as the Coleman-De Luccia instanton. We compute the thermal graviton-induced decoherence rate for a wave function describing a perfect fluid of nonzero energy density in a finite region. When the effective cosmological constant is positive, the thermal graviton background sourced by a de Sitter horizon provides an unavoidable decoherence effect, which may have important consequences for tunneling processes in cosmological history. We discuss generalizations and consequences of this effect and comment on its observability and applications to black hole physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Ensemble variational Fokker-Planck methods for data assimilation

Particle flow filters solve Bayesian inference problems by smoothly transforming a set of particles into samples from the posterior distribution. Particles move in state space under the flow of an McKean-Vlasov-Itˆo process. This work introduces the Variational Fokker-Planck (VFP) framework for data assimilation, a general approach that includes previously known particle flow filters as special cases. The McKean-Vlasov-Itˆo process that transforms particles is defined via an optimal drift that depends on the selected diffusion term. It is established that the underlying probability density - sampled by the ensemble of particles - converges to the Bayesian posterior probability density. For a finite number of particles the optimal drift contains a regularization term that nudges particles toward becoming independent random variables. Based on this analysis, we derive computationally-feasible approximate regularization approaches that penalize the mutual information between pairs of particles, and avoid particle collapse. Moreover, the diffusion plays a role akin to a particle rejuvenation approach that aims to alleviate particle collapse. The VFP framework is very flexible. Different assumptions on prior and intermediate probability distributions can be used to implement the optimal drift, and localization and covariance shrinkage can be applied to alleviate the curse of dimensionality. A robust implicit-explicit method is discussed for the efficient integration of stiff McKean- Vlasov-Itˆo processes. Here, the effectiveness of the VFP framework is demonstrated on three progressively more challenging test problems, namely the Lorenz ’63, Lorenz ’96 and the quasi-geostrophic equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Hysteresis between gas breakdown and plasma discharge

In direct-current (DC) discharge, it is well known that hysteresis is observed between the Townsend (gas breakdown) and glow regimes. Forward and backward voltage sweep is performed using a one-dimensional particle-in-cell Monte Carlo collision (PIC-MCC) model considering a ballast resistor. When increasing the applied voltage after reaching the breakdown voltage (V b ), transition from Townsend to glow discharges is observed. When decreasing the applied voltage from the glow regime, the discharge voltage (V d ) between the anode–cathode gap can be smaller than the breakdown voltage, resulting in a hysteresis, which is consistent with experimental observations. Next, the PIC-MCC model is used to investigate the self-sustaining voltage (V s ) in the presence of finite initial plasma densities between the anode and cathode gap. It is observed that the self-sustaining voltage coincides with the discharge voltage obtained from the backward voltage sweep. In addition, the self-sustaining voltage decreases with increased initial plasma density and saturates above a certain initial plasma density, which indicates a change in plasma resistivity. The decrease in self-sustaining voltage is associated with the electron heat loss at the anode for the low pd (rarefied) regime. In the high pd (collisional) regime, the ion energy loss toward the cathode due to the cathode fall and the inelastic collision loss of electrons in the bulk discharge balance out. Finally, it is demonstrated that the self-sustaining voltage collapses to a singular value, despite the presence of a initial plasma, for microgaps when field emission is dominant, which is also consistent with experimental observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Schwinger-Keldysh effective action for hydrodynamics with approximate symmetries

We study hydrodynamic theories with approximate symmetries in the recently developed effective action approach on the Schwinger-Keldysh contour. We employ the method of spurious symmetry transformation for small explicit symmetry-breaking parameters to systematically constrain symmetry-breaking effects in the nonequilibrium effective action for hydrodynamics. We apply our method to the hydrodynamic theory of chiral symmetry in quantum chromodynamics at finite temperature and density and its explicit breaking by quark masses. We show that the spurious symmetry and the Kubo-Martin-Schwinger relation dictate that the Ward-Takahashi identity for the axial symmetry, i.e., the partial conservation of axial vector current (PCAC) relation, contains a relaxational term proportional to the axial chemical potential, whose kinetic coefficient is at least of the second order in the quark mass. In the phase where the chiral symmetry is spontaneously broken, and the pseudo-Nambu-Goldstone pions appear as hydrodynamic variables, this relaxation effect is subleading compared to the conventional pion mass term in the PCAC relation, which is of the first order in the quark mass. On the other hand, in the chiral symmetry restored phase, we show that our relaxation term, which is of the second order in the quark mass, becomes the leading contribution to the axial charge relaxation. Therefore, the leading axial charge relaxation mechanism is parametrically different in the quark mass across a chiral phase transition.

Goldstone bosons↗
DOE OSTI2022

Possible quadrupole-order-driven commensurate-incommensurate phase transition in B20 CoGe

For this study, the B20-type cobalt germanide CoGe was investigated by measuring the specific heat, resistivity, and 59 Co nuclear magnetic resonance (NMR). We observed a phase transition at T Q = 13.7 K, evidenced by a very narrow peak of the specific heat and sharp changes of the nuclear spin-spin (T 2 –1 ) and spin-lattice ( T 1 – 1 ) relaxation rates. The fact that the entropy release is extremely small and the Knight shift is almost independent of temperature down to low temperatures as anticipated in a paramagnetic metal indicates that the T Q transition is of nonmagnetic origin. In addition, we detected a crossover scale T 0 ~ 30 K below which the resistivity and the NMR linewidth increase, and T 1 –1 is progressively distributed in space, that is, a static and dynamical spatial inhomogeneity develops. While the order parameter for the T Q transition remains an open question, a group-theoretical analysis suggests that the finite electric quadrupole density arising from the low local site symmetry at cobalt sites could drive the crystal symmetry lowering from the P2 1 3 symmetry that is commensurate to the R3 symmetry with an incommensurate wave vector, which fairly well accounts for the T Q transition. The quadrupole-order-driven commensurate-incommensurate phase transition may be another remarkable phenomenon arising from the structural chirality inherent in the noncentrosymmetric B20 family.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Extending a scaling equation of state to QCD

Whether quantum chromodynamics (QCD) exhibits a phase transition at finite temperature and density is an open question. It is important for hydrodynamic modeling of heavy ion collisions and neutron-star mergers. Lattice QCD simulations have definitively shown that the transition from hadrons to quarks and gluons is a crossover when the baryon chemical potential is zero or small. Here, we combine the parametric scaling equation of state, usually associated with the three-dimensional Ising model, with a background equation of state based on a smooth crossover from hadrons to quarks and gluons. Comparison to experimental data from the Beam Energy Scan II at the Relativistic Heavy Ion Collider or in heavy ion experiments at other accelerators may allow the critical exponents and amplitudes in the scaling equation of state to be determined for QCD if a critical point exists.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Fluctuations and correlations of baryonic chiral partners

The exploration of critical phenomena in phase transitions of strongly interacting matter governed by quantum chromodynamics (QCD) is one of the goals of present ultrarelativistic heavy-ion collision experiments at BNL and CERN. The key research direction is to locate the putative critical point on the phase diagram of QCD linked to the chiral symmetry restoration at finite temperature and/or density. One of the main theoretical tools used for this purpose is the fluctuations of conserved charges, such as the net-baryon number. However, due to experimental limitations, analyses of heavy-ion collision data suffer from a very doubtful basing of the net-proton number being a proxy for the total net-baryon number fluctuations. In this work, we use the parity doublet model to investigate the fluctuations of the net-baryon number density in hot and dense hadronic matter. The model accounts for chiral criticality within the mean-field approximation. We focus on the qualitative properties and systematics of the first- and second-order susceptibility of the net-baryon number density, and their ratios for nucleons of positive and negative parity, as well as their correlator. We show that the fluctuations of the positive-parity nucleon do not necessarily reflect the fluctuations of the total net-baryon number density at the phase boundary of the chiral phase transition. We also investigate the nontrivial structure of the correlator. Furthermore, we discuss and quantify the differences between the fluctuations of the net-baryon number density in the vicinity of the chiral and liquid-gas phase transition in nuclear matter. We indicate a possible relevance of our results with the interpretation of the experimental data on net-proton number fluctuations in heavy-ion collisions. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

QCD moat regime and its real-time properties

Dense quantum chromodynamics (QCD) matter may exhibit crystalline phases. Their existence is reflected in a moat regime, where mesonic correlations feature spatial modulations. We study the real-time properties of pions at finite temperature and density in QCD in order to elucidate the nature of this regime. We show that the moat regime arises from particle-hole-like fluctuations near the Fermi surface. This gives rise to a characteristic peak in the spectral function of the pion at nonzero spacelike momentum. This peak can be interpreted as a new quasi particle, the moaton. In addition, our framework also allows us to directly test the stability of the homogeneous chiral phase against the formation of an inhomogeneous condensate in QCD. We find that an inhomogeneous instability is highly unlikely for baryon chemical potentials μ B ≤ 630 MeV . Published by the American Physical Society 2025

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Ionization state and dielectric constant in cold rarefied hydrocarbon plasmas of inertial confinement fusion

A combined approach to study cold rarefied matter is introduced that includes a semianalytical method based on the free-energy minimization and ab initio calculations based on the finite-temperature density-functional theory. The approach is used to calculate the ionization state of hydrocarbon (CH) under the shock-release conditions in inertial confinement fusion. Furthermore, the dielectric constant of CH is calculated using the Kubo-Greenwood formulation and contribution from atomic polarizabilities is found to be as important as the free-electron contribution. Using the ionization state and dielectric constant, the electron density profile in the rarefaction wave of the shock-release plasma is obtained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1977

Calculation of sheath and wake structure about a pillbox-shaped spacecraft in a flowing plasma

A computer program was used for studies of the disturbed zones around bodies in flowing plasmas, particularly spacecraft and their associated sheaths and wakes. The program solved a coupled Poisson-Vlasov system of nonlinear partial differential integral equations to obtain distributions of electric potential and ion and electron density about a finite length cylinder in a plasma flow at arbitrary ion Mach numbers. The approach was applicable to a larger range of parameters than other available approaches. In sample calculations, bodies up to 100 Debye lengths in radius were treated, that is, larger than any previously treated realistically. Applications were made to in-situ satellite experiments.

Parker, L. W.↗
NASA NTRS1981

Plasma lasers (a strong source of coherent radiation in astrophysics)

The generation of electromagnetic radiation from the free energy available in electron streams is discussed. The fundamental principles involved in a particular class of coherent plasma radiation sources, i.e., plasma lasers, are reviewed, focusing on three wave coupling, nonlinear parametric instabilities, and negative energy waves. The simplest case of plasma lasers, that of an unmagnetized plasma containing a finite level of density fluctuations and electrons streaming with respect to the ions, is dealt with. A much more complicated application of plasma lasers to the case of auroral kilometric radiation is then examined. The concept of free electron lasers, including the role of relativistic scattering, is elucidated. Important problems involving the escape of the excited radiation from its generation region, effects due to plasma shielding and nonlinear limits, are brought out.

Papadopoulos, K.↗
NASA NTRS1983

Numerical analysis of plasma flows in an MPC duct with diverter

The two-dimensional MHD-flow of the ideal plasma in a channel of the magnetoplasma compressor (MPC) with an axial hole in the central electrode (divertor) is numerically simulated. The steady-state regime of the flow is obtained. The influence of finite and periodical density perturbation at the entry of the channel on the compressional flow properties is investigated. It is shown that the flow is stable under such perturbations.

Badin, V. V.↗
NASA NTRS1990

Self-similar perturbations of a Friedmann universe

The present analysis of spherically symmetric self-similar solutions to the Einstein equations gives attention to those solutions that are asymptotically k = 0 Friedmann at large z values, and possess finite but perturbed density at the origin. Such solutions represent nonlinear density fluctuations which grow at the same rate as the universe's particle horizon. The overdense solutions span only a narrow range of parameters, and resemble static isothermal gas spheres just within the sonic point; the underdense solutions may have arbitrarily low density at the origin while exhibiting a unique relationship between amplitude and scale. Their relevance to large-scale void formation is considered.

Carr, Bernard J.↗
NASA NTRS1996

Studies of Interstellar Pickup Ions in the Solar Wind

The work under this grant involves studies of the interaction of interstellar pickup ions with the solar wind, with the goal of a comprehensive model of the particle distributions and wave intensities to be expected throughout the heliosphere, as well as the interactions of those distributions with the solar wind termination shock. In the past year, we have completed a number of projects, including observations and modeling of the effects of a large scattering mean free path on the pickup He(+) seen at AMPTE, an analytical model of anisotropic pickup tons in a steady radial magnetic field, and a derivation of a reduced solar wind Mach number due to increased estimates on the inflowing hydrogen density allowing for a weak termination shock. In the next year, we plan to investigate in more detail the correspondence between our models of anisotropic pickup ions and the data on spectra, variations, and proton-He(+) correlation provided by AMPTE, Ulysses, and our instrument on SOHO. We will model the time-dependent pickup ion density resulting from finite periods of radial magnetic field. We will also incorporate the effects of a large mean free path into our analysis of the He(+) focusing cone, leading to more accurate parameter values for the interstellar helium gas. This progress report also includes a discussion of our Space Physics Educational Outreach activities in the past year and plans for the next year.

Isenberg, Philip A.↗
NASA NTRS2013

The Effect of Scale Dependent Discretization on the Progressive Failure of Composite Materials Using Multiscale Analyses

A multiscale modeling methodology, which incorporates a statistical distribution of fiber strengths into coupled micromechanics/ finite element analyses, is applied to unidirectional polymer matrix composites (PMCs) to analyze the effect of mesh discretization both at the micro- and macroscales on the predicted ultimate tensile (UTS) strength and failure behavior. The NASA code FEAMAC and the ABAQUS finite element solver were used to analyze the progressive failure of a PMC tensile specimen that initiates at the repeating unit cell (RUC) level. Three different finite element mesh densities were employed and each coupled with an appropriate RUC. Multiple simulations were performed in order to assess the effect of a statistical distribution of fiber strengths on the bulk composite failure and predicted strength. The coupled effects of both the micro- and macroscale discretizations were found to have a noticeable effect on the predicted UTS and computational efficiency of the simulations.

Ricks, Trenton M.↗