Search NASA⌕ Search

SEARCH · Search NASA

Results for “radius equation”

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 109 records · Page 6

Design optimization of gap distance for the capillary limitation of a heat pipe with annular-type wick structure

In this study, an experimental investigation was conducted on the rising height and contact angle of fluid in an annular wick-type heat pipe. The annular wick-type heat pipe was characterized by a small gap between the wick structure and tube wall, which compensated for the pressure drop along the porous media and created additional capillary force. To describe and model the advantage of this gap, the rising of a wetting liquid in the gap between a vertical solid plate and a mesh (with a small angle between them) was experimentally measured and analyzed. An additional experiment was performed to investigate the effect of curvature on the capillary rise using tubes and meshes of varying radii. Resultantly, we confirmed that the linear combination of the contact angles of the solid plate and mesh could be applied to calculate the rising height from the Laplace–Young equation. Furthermore, the effect of curvature on the rising height of the liquid was negligible. These results were extended to the investigation of finding the optimal gap distance for the annular wick-type heat pipe by referring to previous studies. We observed that a gap distance of 1.27 mm provided the largest permeability (K) over the effective pore radius (r eff ) value for a heat pipe with ethanol, which in turn resulted in the highest capillary limitation. For a sodium heat pipe, a gap distance of 0.84 mm resulted in the highest capillary limitation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Interpreting Mass and Radius Measurements of Neutron Stars with Dark Matter Halos

Abstract The high densities of neutron stars (NSs) could provide astrophysical locations for dark matter (DM) to accumulate. Depending on the DM model, these DM admixed NSs (DANSs) could have significantly different properties than pure baryonic NSs, accessible through X-ray observations of rotation-powered pulsars. We adopt the two-fluid formalism in general relativity to numerically simulate stable configurations of DANSs, assuming a fermionic equation of state (EOS) for the DM with repulsive self-interaction. The distribution of DM in the DANS as a halo affects the path of X-rays emitted from hot spots on the visible baryonic surface, causing notable changes in the pulse profile observed by telescopes such as NICER, compared to pure baryonic NSs. We explore how various DM models affect the DM mass distribution, leading to different types of dark halos. We quantify the deviation in observed X-ray flux from stars with each of these halos. We identify the pitfalls in interpreting mass and radius measurements of NSs inferred from electromagnetic radiation and constraining the baryonic matter EOS if these dark halos exist.

Shawqi, Shafayat (ORCID:0000000210956183)↗

Role of fast ions in spontaneous neoclassical tearing mode instabilities in NSTX

Spontaneous neoclassical tearing modes (spontaneous NTMs) have been observed in the National Spherical Torus Experiment (NSTX). Recently, a TRANSP based model for the modified Rutherford equation analysis was developed to accurately predict the effect of fast ions on the growth of magnetic islands. The goal of this paper is to utilize the TRANSP NTM model to understand the role of fast ions in NSTX spontaneous NTMs. Our analysis shows that when the passing fast ion driven kinetic neoclassical polarization current term is larger than 1% of the bootstrap current term in the modified Rutherford equation, fast ions play a decisive role in the early growth phase of spontaneous NTM. This result is applicable to any tokamak plasmas with a fast ion drift orbit width, or the fast ion poloidal Larmor radius, comparable to the magnetic island width.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Comparison of Eulerian Bin and Lagrangian Particle-Based Schemes in Simulations of Pi Chamber Dynamics and Microphysics

This paper presents a comparison of simulations applying either a traditional Eulerian bin microphysics or a novel particle-based Lagrangian approach to represent CCN activation and cloud droplet growth. The Eulerian microphysics solve the evolution equation for the spectral density function, whereas the Lagrangian approach follows computational particles referred to as superdroplets. Each superdroplet represents a multiplicity of natural droplets that makes the Lagrangian approach computationally feasible. The two schemes apply identical representation of CCN activation and use the same droplet growth equation; these make direct comparison between the two schemes practical. The comparison, the first of its kind, applies an idealized simulation setup motivated by laboratory experiments with the Pi Chamber and previous model simulations of the Pi Chamber dynamics and microphysics. The Pi Chamber laboratory apparatus considers interactions between turbulence, CCN activation, and cloud droplet growth in moist Rayleigh–Bénard convection. Simulated steady-state droplet spectra averaged over the entire chamber are similar, with the mean droplet concentration, mean radius, and spectral width close in Eulerian and Lagrangian simulations. Small differences that do exist are explained by the inherent differences between the two schemes and their numerical implementation. The local droplet spectra differ substantially, again in agreement with the inherent limitations of the theoretical foundation behind each approach. There is a general agreement between simulations and Pi Chamber observations, with simplifications of the CCN activation and droplet growth equation used in the simulations likely explaining specific differences.

54 ENVIRONMENTAL SCIENCES↗

Leptodermic corrections to the TOV equations and nuclear astrophysics within the effective surface approximation

The macroscopic model for a neutron star (NS) as a liquid drop at the equilibrium is used to extend the Tolman-Oppenheimer-Volkoff (TOV) equations taking into account the gradient terms responsible for the system surface. The parameters of the Schwarzschild metric in the spherical case are found with these surface corrections to the known leading (zero) order of the leptodermic approximation a/R << 1, where a is the NS effective-surface (ES) thickness, and R is the effective NS radius. The energy density $\mathscr{E}$ is considered in a general form including the functions of the particle number density and of its gradient terms. The macroscopic gravitational component $Φ$(ρ) of the energy density is taken into account in the simplest form as expansion in powers of $ρ$ – $\overline{ρ}$, where $\overline{ρ}$ is the saturation density, up to second order, in terms of its contributions to the separation particle energy and incompressibility. Density distributions ρ across the NS ES in the normal direction to the ES, which are derived in the simple analytical form at the same leading approximation, was used for the derivation of the modified TOV (MTOV) equations by accounting for their NS surface corrections. As a result, the MTOV equations are analytically solved at first order and the results are compared with the standard TOV approach of the zero order.

Magner, A. G. [Institute for Nuclear Research, Kyi↗

Determination of the characteristic magnetic pre-sheath length at divertor surfaces using micro-engineered targets on DiMES at DIII-D

The magnetic pre-sheath (MPS) width, L MPS , is a critical parameter to define the sheath potential, which controls the ion trajectory of low-Z species (D, T, He, and C), as well as the prompt re-deposition of high-Z species. In this work, to determine L MPS , we fabricated micro-trenches (30×30×4 µm) via focused ion beam (FIB) milling on a silicon surface and exposed them to L-mode deuterium plasmas in DIII-D via the Divertor Material Evaluation System (DiMES) removable sample exposure probe. The areal distribution of impurity depositions, mainly consisting of carbon, was measured by energy-dispersive X-ray spectroscopy (EDS) to reveal the deuterium ion shadowing effect on the trench floors. The carbon deposition profiles showed that the erosion was maximized for the azimuthal direction of φ = -40° (referenced to the toroidal magnetic field direction) as well as the polar angle of θ = 80°. A Monte Carlo equation-of-motion model, based on a collisionless MPS, was used to calculate the azimuthal and polar deuterium ion angle distributions (IADs) for a range of L MPS = k × ρ i , where ρ i is the ion gyro radius and k = 0.5-4. Then, gross erosion profiles were calculated by a Monte Carlo micro-patterning and roughness (MPR) code for ion sputtering using as input the calculated azimuthal and polar IADs for each value of k . Good agreement with the experimental C deposition profiles was obtained for the case k = 2.5-3.5. This result is consistent with a previous kinetic modeling prediction of k ~ 3, as well as previous analytical investigations that predicted the L MPS to be several ion gyro radii. A validation of theoretical sheath models supports its applicability to ITER and pilot plant divertors to successfully predict plasma-materials interactions.

Erosion↗

Testing $f(R)$ gravity with scale dependent cosmic void velocity profiles

Here, we study the impact of cosmological scale modifications to general relativity on the dynamics of halos within voids by comparing N-body simulations incorporating Hu-Sawicki f(R) gravity, with |f R0 | = 10 -6 and 10 -5 , to those of Λ CDM. By examining the radial velocity statistics within voids classified based on their size and density-profile, as “rising” (R-type) or “shell” (S-type), we find that halo motions in small R-type voids, with effective radius <15 Mpc/h, reveal distinctive differences between f(R) and Λ CDM cosmologies. To understand this observed effect, we study the linear and nonlinear fifth forces, and develop an iterative algorithm to accurately solve the nonlinear fifth force equation. We use this to characterize the Chameleon screening mechanism in voids and contrast the behavior with that observed in gravitationally collapsed objects. The force analysis underscores how smaller R-type voids exhibit the highest ratios of fifth force to Newtonian force, which source distinguishable differences in the velocity profiles and thereby provide rich environments in which to constrain gravity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modelling of nanometer scale dust grains in tokamak

Dust poses a serious threat to tokamak operation and safety. It is important to study the behaviour of dust grains under tokamak's discharge conditions, which depends heavily on their size and charge. Existing simulations mainly address issues on dust grains with radii larger than 1 μm, in which case, the drift effect due to electromagnetic fields can be safely ignored. For nanometer scale dust grains, however, the drift effect becomes significant and a new model based on guiding-centre system needs to be established. In this work, the NDS has been done under BOUT++ framework. The simulation contains two parts. Part one, NDS evaluates the charging and ablation processes of the dust grains. In the second part, the guiding-centre orbits of dust particles are tracked in tokamak plasmas, whose parameters are obtained from BOUT++, a highly desirable C++ code package for performing parallel plasma fluid simulations with an arbitrary number of equations in 3D curvilinear coordinates. The orbit of nanodust dynamics is described by guiding centre equations for simplicity, and these equations are numerically solved by conventional fourth-order Runge Kutta method. Simulations provide results such as trajectories and evolutions of dust particles with different sizes and velocities for different tokamak geometries. Results show tungsten dust grains with a radius of a few nanometers launched from outer midplane will oscillate before totally ablated in C-Mod. The oscillation in this case is driven by the ion drag force. Larger Nanodust with a radius of 100 nm, on the contrary, cannot be completely constrained by the electromagnetic field. In conclusion, the high plasma temperature and density in the seperatrix region causes severe dust ablation, resulting in total ablation within several ms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Eigenmode analysis of the sheared-flow Z-pinch

Experiments have demonstrated that a Z-pinch can persist for thousands of times longer than the growth time of global magnetohydrodynamic (MHD) instabilities such as the m=0 sausage and m=1 kink modes. These modes have growth times on the order of ta=a/vi, where vi is the ion thermal speed and a is the pinch radius. Axial flows with duz/dr ≲ vi/a have been measured during the stable period, and the commonly accepted theory is that this amount of shear is sufficient to stabilize these modes as predicted by numerical studies using the ideal MHD equations. However, these studies only consider specific equilibrium profiles that typically have a modest magnitude for the logarithmic pressure gradient, qP≡d ln P/d ln r, and may not represent experimental conditions. Linear stability of the sheared-flow Z-pinch is studied here via a direct eigen-decomposition of the matrix operator obtained from the linear ideal MHD equations. Several equilibrium profiles with a large variation of qP are examined. Considering a practical range of k, 1/3 ≲ ka ≲ 10, it is shown that the shear required to stabilize m=0 modes can be expressed as duz/dr≥Cγ0/(ka)α. Here, γ0=γ0(ka) is the profile-specific growth rate in the absence of shear, which scales approximately with |qP|. Both C and α are profile-specific constants, but C is order unity and α≈1. It is further demonstrated that even a large value of shear, duz/dr=3vi/a, is not sufficient to provide linear stabilization of the m=1 kink mode for all profiles considered. This result is in contrast to the currently accepted theory predicting stabilization at much lower shear, duz/dr=0.1vi/a, and suggests that the experimentally observed stability cannot be explained within the linear ideal-MHD model.

Angus, J. R. (ORCID:0000000314740002)↗

Analysis of MHD stability and active mode control on KSTAR for high confinement, disruption-free plasma

Long-pulse plasma operation at high normalized beta, $β_N$, above the $\textit{n}$= 1 ideal MHD no-wall stability limit in KSTAR is presently limited by tearing instabilities rather than resistive wall modes. H-mode plasma operation during the recent KSTAR device campaign produced discharges having strong $\textit{m/n}$= 2/1 tearing instabilities at $β_N$ lower than the ideal MHD no-wall beta limit. The unstable tearing mode consequently reduced plasma confinement and toroidal plasma rotation significantly. We report the experiment confirmed that an extended duration of electron cyclotron heating (ECH) at the initial phase of the discharge plays a critical role in mode destabilization. To study destabilizing mechanisms that affect the mode growth, the stability of the observed tearing modes from plasmas with significantly different $β_N$ is computed by using the resistive DCON code and the M3D-C 1 code. Equilibrium reconstructions that include constraints from internal profile diagnostics, and computed fast particle pressure are used as input for reliable computation of stability. The classical tearing stability index, Δ', from resistive DCON is computed to be unstable when the island is fully saturated with large amplitude, while the unstable mode is not computed by M3D-C 1 . The modified Rutherford equation (MRE) describing the evolution of the island width has been constructed for KSTAR plasmas by using plasma parameters computed by the TRANSP code. The MRE model estimates a saturated island width corresponding to ~10% of the plasma minor radius for equilibrium at high $β_N$ 3 having a stable Δ' from the resistive DCON. In preparation for long-pulse plasma operation at higher beta utilizing increased plasma heating power, a resistive wall mode (RWM) active feedback control algorithm has been completed and enabled on KSTAR. To accurately determine the $\textit{n}$= 1component produced by RWMs, an algorithm has been developed that includes magnetic sensor compensation of the prompt applied field and the field from the induced current on the passive conductors. Use of multiple toroidal sensor arrays is enabled by modifying the sensor toroidal angles assumed in mode decomposition to include the effect of varied mode helicities in the outboard region where the mode measurement is made. This analysis on stability, transport, and control provides the required foundation for disruption prediction and avoidance research on KSTAR.

resistive wall mode (RWM)↗

Three-dimensional imaging of proton in basis light-front quantization

Here, we employ an effective Hamiltonian that includes the transverse and longitudinal confinement and the one-gluon exchange interaction with fixed coupling constant. By solving the eigenvalue equation in basis light-front quantization (BLFQ), we generate the light-front wavefunctions (LFWFs) for the nucleon in the valence quark Fock space. Fitting the model parameters, we obtain high quality descriptions of electromagnetic form factors and radius for proton while the results deviate somewhat from experimental data for neutron.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Cation Disorder of ${\text{Mg}}_{\mathbf{2}}{\text{SiO}}_{\mathbf{4}}$ in Super‐Earth Mantles

Understanding the mineralogy of exoplanets is essential for unraveling their interior structures, dynamics, and evolution. For large super-Earths, the post-post spinel ${\text{Mg}}_{\mathbf{2}}{\text{SiO}}_{\mathbf{4}}$, one of the major mantle phases, may undergo the order-disorder transition (ODT) at high temperatures. However, the ODT phase boundary of ${\text{Mg}}_{\mathbf{2}}{\text{SiO}}_{\mathbf{4}}$ has not been rigorously constrained. Additionally, fundamental thermodynamic properties of the disordered ${\text{Mg}}_{\mathbf{2}}{\text{SiO}}_{\mathbf{4}}$ remain poorly investigated. Here, we develop a unified machine learning potential (MLP) for ${\text{Mg}}_{\mathbf{2}}{\text{SiO}}_{\mathbf{4}}$ of ab initio accuracy under super-Earth mantle conditions. With the efficient MLP, we extensively calculate the free energy of post-post spinel ${\text{Mg}}_{\mathbf{2}}{\text{SiO}}_{\mathbf{4}}$ via the thermodynamic integration method. The results are used to constrain the ODT phase boundary. Furthermore, we report the P-V-T equation of state and Grüneisen parameters for post-post spinel ${\text{Mg}}_{\mathbf{2}}{\text{SiO}}_{\mathbf{4}}$ across various degrees of disorder. These thermodynamic properties are further applied to update the adiabatic thermal profiles and the mass-radius relation of super-Earths.

36 MATERIALS SCIENCE↗

What Controls the Entrainment Rate of Dry Buoyant Thermals with Varying Initial Aspect Ratio?

Here, this study uses theory and numerical simulations to analyze the nondimensional spreading rate a (change in radius with height) of buoyant thermals as they rise and entrain surrounding environmental fluid. A focus is on how a varies with initial thermal aspect ratio A r , defined as height divided by width of the initial buoyancy perturbation. An analytic equation for thermal ascent rate w t that depends on a is derived from the thermal-volume-averaged momentum bud get equation. The thermal top height when w t is maximum, defining a critical height z c , is inversely proportional to a. The height z c also corresponds to the thermal top height when buoyant fluid along the thermal’s vertical axis is fully replaced by entrained nonbuoyant environmental fluid rising from below the thermal. The time scale for this process is controlled by the vertical velocity of parcels rising upward through the thermal’s core. This parcel vertical velocity is approximated from Hill’s analytic spherical vortex, yielding an analytic inverse relation between a and A r . Physically, this a–A r relation is connected to changes in circulation as A r is modified. Numerical simulations of thermals with A r varied from 0.5 to 2 give a values close to the analytic theoretical relation, with a factor of ~3 decrease in a as A r is increased from 0.5 to 2. The theory also explains why a of initially spherical thermals from past laboratory and modeling studies is about 0.15. Overall, this study provides a theoretical underpinning for understanding the entrainment behavior of thermals, relevant to buoyantly driven atmospheric flows.

54 ENVIRONMENTAL SCIENCES↗

High-Order Multipole and Binary Love Number Universal Relations

Using a data set of approximately 2 million phenomenological equations of state consistent with observational constraints, we construct new equation-of-state-insensitive universal relations that exist between the multipolar tidal deformability parameters of neutron stars, Λ l , for several high-order multipoles (l = 5,6,7,8), and we consider finite-size effects of these high-order multipoles in waveform modeling. We also confirm the existence of a universal relation between the radius of the 1.4M ⊙ NS, R 1.4 and the reduced tidal parameter of the binary, Λ˜, and the chirp mass. We extend this relation to a large number of chirp masses and to the radii of isolated NSs of different mass M, R M . We find that there is an optimal value of M for every M such that the uncertainty in the estimate of R M is minimized when using the relation. We discuss the utility and implications of these relations for the upcoming LIGO O4 run and third-generation detectors.

97 MATHEMATICS AND COMPUTING↗

Solution of the linear wave-particle kinetic equation for global modes of arbitrary frequency in a tokamak

The linear response of a plasma to perturbations of arbitrary frequency and wavelength is derived for any axisymmetric magnetized toroidal plasma. An explicit transformation to action-angle coordinates is achieved using orthogonal magnetic coordinates and the Littlejohn Lagrangian, establishing the validity of this result to arbitrary order in normalized Larmor radius. The global resonance condition for compressional modes is clarified in more detail than in previous works, confirming that the poloidal orbit-average of the cyclotron frequency gives the desired result at lowest order in Larmor radius. The global plasma response to the perturbation at each resonance is captured by a poloidal and gyroaverage of the perturbing potential. A “global gyroaveraging” of the potential is a natural by-product of this analysis which takes into account the changing of the magnetic field over an orbit. The resonance condition depends on two arbitrary integers which completely separately capture the effects poloidal non-uniformity and finite Larmor radius in generating sidebands. We learn that poloidal sidebands generated for compressional modes are dominated by the change in gyrofrequency over the orbit, which is very different to shear modes where the gyrofrequency only contributes via a finite Larmor radius effect. This increases the number of bounce harmonics required to compute the linear drive, giving a more complicated resonance map. An example calculation is given comparing resonance of shear and compressional modes in a published DIII-D case.

Compressional↗

Properties of the neutron star crust: Quantifying and correlating uncertainties with improved nuclear physics

A compressible liquid-drop model (CLDM) is used to correlate uncertainties associated with the properties of the neutron star (NS) crust with theoretical estimates of the uncertainties associated with the equation of state (EOS) of homogeneous neutron and nuclear matter. For the latter, we employ recent calculations based on Hamiltonians constructed using chiral effective-field theory (χEFT). Fits to experimental nuclear masses are employed to constrain the CLDM further, and we find that they disfavor some of the χEFT Hamiltonians. The CLDM allows us to study the complex interplay between bulk, surface, curvature, and Coulomb contributions, and their impact on the NS crust. It also reveals how the curvature energy alters the correlation between the surface energy and the bulk symmetry energy. Here our analysis quantifies how the uncertainties associated with the EOS of homogeneous matter implies significant uncertainties for the composition of the crust, its proton fraction, and the volume fraction occupied by nuclei. We find that the finite-size effects impact the crust composition but have a negligible effect on the net isospin asymmetry of matter. The isospin asymmetry is largely determined by the bulk properties and the isospin dependence of the surface energy. The most significant uncertainties associated with matter properties in the densest regions of the crust, the precise location of the crust-core transition, are found to be strongly correlated with uncertainties associated with the Hamiltonians. By adopting a unified model to describe the crust and the core of NSs, we tighten the correlation between their global properties such as their mass-radius relationship, moment of inertia, crust thickness, and tidal deformability with uncertainties associated with the nuclear Hamiltonians.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Perpendicular magnetofluid theory for magnetically confined plasmas in the collisionless limit

Note that the magnetically confined plasmas in thermonuclear fusion turn to be collisionless. Some particles are trapped by the inhomogeneous equilibrium field or by the waves, while some others keep circulating. One cannot expect the charged particles to be collectively localized in space by collisions along the magnetic field lines. Consequently, the fluid parameters relating to the parallel motion, such as parallel fluid velocity, become inadequate. In this paper, the perpendic- ular magnetofluid theory for magnetically confined plasmas in the collisionless limit is developed. It revises Braginskii’s transport theory by releasing the collisional dominance assumption. In the perpendicular direction, the particle spatial localization is solo resulted from the strong magnetic field, while the finite Larmor radius (FLR) effects are taken into account in the next order. In the parallel direction, the particle mobility feature is fully retained by constructing only the per- pendicular moments. Therefore, besides the FLR effects till order O($k$ $^{2}_{⊥}$ $ρ$ $^{2}_{i}$ ), the Landau damping, the trapped particle effects, and other parallel kinetic effects are all kept in the current formalism. Here, k ⊥ is the perpendicular wave number and ρ i is the ion Larmor radius. Here, the theory can be used both for the analytical and numerical studies of magnetically confined plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Higher Ion Selectivity with Lower Energy Usage Promoted by Electro-osmotic Flow in the Transport through Conical Nanopores

The tradeoff between selectivity and throughput presents fundamental challenges to improve desalination and charge storage, salinity gradient-based energy harvesting, memory device/circuit development, and so forth. The well-known ion-current rectification and the recently resolved time-dependent transport hysteresis in conical nanopores or asymmetric nanointerfaces offer new opportunities for the selective transport of matter. This report shows that electro-osmotic flow (EOF) is an overlooked factor that increases ion selectivity while maintaining enhanced transport throughput in rectified nanoscale electrokinetic transport. The increased selectivity originates primarily from the suppression of anion flux by the fluid flow in the opposite direction under the applied electrical field. By solving the Poisson and Nernst–Planck (PNP) equations without and with coupled Navier–Stokes (PNP–NS) equations, the EOF effects on cation and anion transport are unequivocally revealed in asymmetric nanopipettes. The flux of cations and anions as well as the transference number and flow velocity are elucidated using the models and boundary conditions validated by previous experiments. A dimensionless parameter, radius over the Debye length, reveals optimal ion selectivity and energy cost at intermediate ion concentrations and nanopore sizes, up to hundreds of millimolars and tens of nanometers under time-dependent potential stimulus. Furthermore, the fundamental insights into EOF at nanointerfaces suggest new routes/strategies for better separation, analysis, and energy applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗