Search NASA⌕ Search

SEARCH · Search NASA

Results for “CLOSURES”

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,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

From Weibel seeds to dynamo beyond pair-plasmas

Bridging the spatiotemporal scales of magnetic seed field generation and subsequent dynamo amplification in the weakly collisional intracluster medium presents an extreme numerical challenge. We perform collisionless turbulence simulations with initially unmagnetised electrons that capture both magnetic seed generation via the electron Weibel instability and the ensuing dynamo amplification. Going beyond existing pair-plasma studies, we use an ion-to-electron mass ratio of 100 for which we find electron and ion dynamics are sufficiently decoupled. These simulations are enabled by the 10-moment collisionless fluid solver of Gkeyll, which evolves the full pressure tensor for all species. The electron heat-flux closure regulates pressure isotropisation and effectively sets the magnetic Reynolds number. We investigate how the strength of the closure influences the transition between a regime reminiscent of previous kinetic pair-plasma simulations and a regime exhibiting dynamo behaviour qualitatively similar to magnetohydrodynamics.

astrophysical plasmas↗
DOE OSTI2022

The Role of Methylaryl Radicals in the Growth of Polycyclic Aromatic Hydrocarbons: The Formation of Five-Membered Rings

The regions of the C 13 H 11 potential energy surface (PES) related to the unimolecular isomerization and decomposition of the 1-methylbiphenylyl radical and accessed by the 1-/2-methylnaphthyl + C 2 H 2 reactions have been explored by ab initio G3(MP2,CC)//B3LYP/6-311G(d,p) calculations. The kinetics of these reactions relevant to the growth of polycyclic aromatic hydrocarbons (PAH) under high-temperature conditions in circumstellar envelopes and in combustion flames has been studied employing the RRKM-Master Equation approach. The unimolecular reaction of 1-methylbiphenylyl proceeding via a five-membered ring closure followed by H elimination is predicted to be very fast, on a submicrosecond scale above 1000 K and to result in the formation of an embedded five-membered ring in the 9H-fluorene product. The 1-/2-methylnaphthyl + C 2 H 2 reaction mechanism involves acetylene addition to the radical on the methylene group followed by a six- or five-membered ring closure and aromatization via an H atom loss. Despite of the complexity of the C 13 H 11 PES, these straightforward pathways are dominant in the high-temperature regime (above ~1000 K), with the prevailing products being phenalene, with a significant contribution of 1H-cyclopenta(a)naphthalene, for 1-methylnaphthyl + C 2 H 2 , and 1H-cyclopenta(b)naphthalene and 3H-cyclopenta(a)naphthalene, for 2-methylnaphthyl + C 2 H 2 . The methylnaphthyl reactions with acetylene represent a clean source of the three-ring PAHs, but they are relatively slow owing to the high entrance barriers of ~10 kcal/mol, with the rate constants of about an order of magnitude lower as compared to those for naphthyl + allene and σ-aryl + C 2 H 2 . Here, the 1-methylnaphthyl + C 2 H 2 and 2-methylnaphthyl + C 2 H 2 reactions represent prototypes for PAH growth by an extra six- and five-membered ring on a zigzag edge or a corner of PAH and the generated modified Arrhenius expressions are recommended for kinetic modeling of PAH expansion by the mechanism of acetylene addition to methylaryl radicals.

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

Thermodynamics–Structure–Dynamics Correlations and Nonuniversal Effects in the Elastically Collective Activated Hopping Theory of Glass-Forming Liquids

We employ the microscopic Elastically Collective Nonlinear Langevin Equation (ECNLE) theory of activated dynamics in combination with crystal-avoiding simulations to study four inter-related questions for metastable monodisperse hard sphere fluids. The first is how significantly improved integral equation theory structural input (Modified-Verlet (MV) closure) changes the dynamical predictions of ECNLE theory. The main consequence is a modest enhancement of the importance of the collective elastic barrier relative to its local cage contribution, which increases the alpha relaxation time and fragility relative to prior results based on the Percus–Yevick closure. Second, ECNLE-MV theory predictions for the alpha time and self-diffusion constant in the metastable regime are quantitatively compared to our new simulations. The small adjustment of a numerical prefactor that enters the collective elastic barrier leads to quantitative agreement over three decades. Third, using the more accurate MV structural input, ECNLE theory is shown to predict thermodynamics–structure–dynamics “correlations” based on various long and short wavelength scalar properties all related to static two-point collective density fluctuations. The logarithm of the alpha relaxation time scales as a power law with these scalar metrics with an exponent that is significantly lower in the less dense noncooperative activated regime compared to the very dense highly cooperative regime. However, the discovered correlation of activated relaxation with a thermodynamic property (dimensionless compressibility) is not causal in ECNLE theory, but rather reflects a strong connection between the local structural quantities that quantify kinetic constraints in the theory with the amplitude of long wavelength density fluctuations. Fourth, the consequences of chemically specific nonuniversalities associated with the onset condition and relative importance of collective elasticity are studied. In conclusion, the predicted thermodynamics–structure–dynamics correlations are found to be robust, albeit with nontrivial shifts of the onset condition.

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

Geometrical Frustration and Defect Formation in Growth of Colloidal Nanoparticle Crystals on a Cylinder: Implications for Assembly of Chiral Nanomaterials

Using a combination of experiment and simulation, we study how two-dimensional (2D) crystals of colloidal nanoparticles grow on cylindrical substrates. The cylindrical geometry allows us to examine growth in the absence of Gaussian curvature but in the presence of a closure constraint - the requirement that a crystal loops back onto itself. In some cases, this constraint results in structures that have been observed previously in theory and nonequilibrium packing experiments: chiral crystals and crystals with linear defects known as "line slips". More generally, though, the structures we see differ from those that have been observed: the line slips are kinked and contain partial vacancies. We show that these structures arise because the cylinder changes how the crystal grows. After a crystal wraps around the cylinder and touches itself, it must grow preferentially along the cylinder axis. As a result, crystals with a chiral line slip tend to trap partial vacancies. Indeed, we find that line slips that are less aligned with the cylinder axis incorporate more partial vacancies on average than the ones that are more aligned. Here, these results show that crystal growth on a cylinder is frustrated by the closure requirement, a finding that may shed some light on the assembly of biological nanosystems such as tobacco mosaic virus and might inform ways to fabricate chiral optical nanomaterials.

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

Drift Instabilities in Thin Current Sheets Using a Two–Fluid Model With Pressure Tensor Effects

The integration of kinetic effects in fluid models is important for global simulations of the Earth's magnetosphere. We use a two–fluid 10–moment model, which includes the pressure tensor and has been used to study reconnection, to study the drift kink and lower hybrid drift instabilities. Using a nonlocal linear eigenmode analysis, we find that for the kink mode, the 10–moment model shows good agreement with kinetic calculations with the same closure model used in reconnection simulations, while the electromagnetic and electrostatic lower hybrid instabilities require modeling the effects of the ion resonance using a Landau fluid closure. As a result, comparisons with kinetic simulations and the implications of the results for global magnetospheric simulations are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

The Impact of Pore-Scale Flow Regimes on Upscaling of Immiscible Two-Phase Flow in Porous Media

Empirical or theoretical extensions of Darcy’s law for immiscible two-phase flow have shown significant limitations in properly modeling the flow at the continuum scale. We tackle this problem by proposing a set of upscaled equations based on pore-scale flow regimes, that is, the topology of flowing phases. The incompressible Navier-Stokes equation is upscaled by means of multiple-scale expansions and its closures derived from the mechanical energy balance for different flow regimes at the pore scale. We also derive the applicability conditions of the upscaled equations based on the order of magnitude of relevant dimensionless numbers, that is, Eotvos, Reynolds, capillary, Froude numbers, and the viscosity and density ratio of the system, as well as a set of closures valid for the basic flow regimes of low Eotvos number systems, that is, core-annular and plug and drop traffic flows. We provide analytical expressions for the relative permeability of the wetting and nonwetting phases in different flow regimes and demonstrate that the effect of the flowing-phases topology on the relative permeabilities is significant. Finally, we show that the classical two-phase Darcy law is recovered for a limited range of operative conditions, while specific terms accounting for interfacial and wall interactions should be incorporated to accurately model ganglia or drop traffic flow.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

An Extended MHD Study of the 16 October 2015 MMS Diffusion Region Crossing

The Magnetospheric Multiscale (MMS) mission has given us unprecedented access to high cadence particle and field data of magnetic reconnection at Earth's magnetopause. MMS first passed very near an X-line on 16 October 2015, the Burch event, and has since observed multiple X-line crossings. Subsequent 3-D particle-in-cell (PIC) modeling efforts of and comparison with the Burch event have revealed a host of novel physical insights concerning magnetic reconnection, turbulence-induced particle mixing, and secondary instabilities. In this study, we employ the Gkeyll simulation framework to study the Burch event with different classes of extended, multifluid magnetohydrodynamics (MHD), including models that incorporate important kinetic effects, such as the electron pressure tensor, with physics-based closure relations designed to capture linear Landau damping. Such fluid modeling approaches are able to capture different levels of kinetic physics in global simulations and are generally less costly than fully kinetic PIC. Herein we focus on the additional physics one can capture with increasing levels of fluid closure refinement via comparison with MMS data and existing PIC simulations. In particular, we find that the ten-moment model well captures the agyrotropic structure of the pressure tensor in the vicinity of the X-line and the magnitude of anisotropic electron heating observed in MMS and PIC simulations. However, the ten-moment model is found to have difficulty resolving the lower hybrid drift instability, which plays a fundamental role in heating and mixing electrons in the current layer.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Improving the Physical Basis for Updraft Dynamics in Deep Convection Parameterizations

This article presents a new deep convective parameterization that determines cloud characteristics based on a specified cloud size distribution. The vertical profiles of cloud properties are determined by analytic equations, which formulate entrainment with an inverse relationship to cloud width. In line with recent studies of large eddy simulations (LES), cloud widths are assumed to be constant with height and vertical mass flux (M) characteristics of the clouds are therefore regulated by the vertical velocity profile. The parameterization is configured to work with existing cloud base M closure formulations, with the closure predicting the total cloud area rather than the cloud base M directly. Analytic formula are also used to connect the vertical wind shear magnitude to the cloud size distribution, wherein larger shear magnitudes result in more numerous large updrafts than weaker shear magnitudes, which is in line with recent research results. The parameterization is compared against 10 deep convective LES with varying thermodynamic and vertical wind shear profiles. Results show dramatic improvements in the prediction of normalized M, detrainment, and the properties of detrained air over the existing Zhang and McFarlane (1995) scheme. In particular, the new model is able to correctly portray the transition from a bottom-heavy M profile in weakly sheared environments, to a top-heavy M profile in strongly sheared environments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Linking Deep and Shallow Convective Mass Fluxes via an Assumed Entrainment Distribution in CAM5‐CLUBB: Parameterization and Simulated Precipitation Variability

Abstract We modify the Zhang‐McFarlane (ZM) deep convection scheme in the Community Atmosphere Model version 5 to couple it with a unified parameterization for boundary‐layer turbulence and shallow convection, that is, Cloud Layers Unified by Binormals (CLUBB). By assuming a lognormal distribution of entrainment rate across the entire moist convective regimes, we link mass fluxes between shallow and deep convection, which are partitioned by the entrainment rate of the shallowest deep convective plume. Hence, a new deep convective closure is established which is coupled to the sub‐grid vertical motion variability in CLUBB. The convection feedback (or memory) effects are also considered to decrease the entrainment spectrum width and enhance the vertical velocity variability that further affect deep convection. Results show that the revised scheme improves the precipitation simulations in terms of the mean state and variability at various timescales, such as the alleviated double‐intertropical convergence zone and more realistic simulations of the seasonal variation of monsoon precipitation over East Asia, Madden‐Julian Oscillation, and precipitation diurnal phase propagations downstream of large terrains. The improvements are still seen in many aspects such as the mean‐state precipitation when turning off the convection feedback impacts in the revised scheme, emphasizing the benefits of using the modified mass‐flux closure. However, the convection feedbacks have considerable effects on the precipitation diurnal cycle simulations over regions with late‐afternoon precipitation peaks. Overall, the revised scheme provides a unified treatment for sub‐grid vertical motions across regimes of boundary‐layer turbulence, shallow convection, and deep convection, leading to better‐simulated precipitation at various timescales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Automated Symbolic Upscaling: 1. Model Generation for Extended Applicability Regimes

Abstract In porous media theory, upscaling techniques are fundamental to deriving rigorous Darcy‐scale models for flow and reactive transport in subsurface systems. Due to limitations in classical techniques, a number of ad hoc approaches have been proposed to address physical regimes in which reactive time scales are similar to, or faster than, diffusive time scales. In Part 1 of this two part series, we present a strategy for expanding the applicability of classical homogenization theory by generalizing the assumed closure form. We detail the implementation of this strategy on two reactive mass transport problems with moderately reactive physics. The strategy produces nontrivial homogenized models with emergent terms and effective parameters that couple reactive, diffusive, and advective transport. The differences in equation forms between the macroscopic and pore‐scale descriptions advise caution to further studies where the forms of macroscopic equations are assumed, as opposed to rigorously derived. Numerical validation is provided for each problem to show that the error estimates of homogenization theory are satisfied, and to justify the implemented strategy. In Part 2, the presented strategy is automated using symbolic computing to expedite its implementation.

Pietrzyk, Kyle↗
DOE OSTI2018

Evidence for prevalent Z = 6 magic number in neutron-rich carbon isotopes

Here, the nuclear shell structure, which originates in the nearly independent motion of nucleons in an average potential, provides an important guide for our understanding of nuclear structure and the underlying nuclear forces. Its most remarkable fingerprint is the existence of the so-called magic numbers of protons and neutrons associated with extra stability. Although the introduction of a phenomenological spin–orbit (SO) coupling force in 1949 helped in explaining the magic numbers, its origins are still open questions. Here, we present experimental evidence for the smallest SO-originated magic number (subshell closure) at the proton number six in C obtained from systematic analysis of point-proton distribution radii, electromagnetic transition rates and atomic masses of light nuclei. Performing ab initio calculations on C, we show that the observed proton distribution radii and subshell closure can be explained by the state-of-the-art nuclear theory with chiral nucleon–nucleon and three-nucleon forces, which are rooted in the quantum chromodynamics.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2013

Resistive magnetohydrodynamic simulations of helicity-injected startup plasmas in National Spherical Torus eXperiment

The generation of helicity-injected startup plasmas in National Spherical Torus eXperiment (NSTX), including flux surface closure, is studied using resistive-magnetohydrodynamic simulations with plasma flows, currents, ohmic heating and anisotropic thermal conduction. An injection-voltage pulse shape is used that separates the injection and closure phases allowing elucidation of the physics. The formation of an X-point near the helicity-injection gap is triggered as the injector voltage drops to zero. Near the forming X-point, magnetic pressure due to toroidal field entrained in the E × B plasma flow from the helicity-injection gap drops, allowing resistive magnetic reconnection even though the total injected current is almost constant. Where appropriate, the simulations are compared with Transient Coaxial Helicity Injection experiments in the NSTX spherical tokamak, which have demonstrated the formation of a promising candidate for non-inductive startup plasmas [Raman et al., Phys. Rev. Lett. 90, 075005 (2003)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Shell model nuclear matrix elements for competing mechanisms contributing to double beta decay

Recent progress in the shell model approach to the nuclear matrix elements for the double beta decay process are presented. This includes nuclear matrix elements for competing mechanisms to neutrionless double beta decay, a comparison between closure and non-closure approximation for {sup 48}Ca, and an updated shell model analysis of nuclear matrix elements for the double beta decay of {sup 136}Xe.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Accurate small and wide angle x-ray scattering profiles from atomic models of proteins and nucleic acids

A new method is introduced to compute X-ray solution scattering profiles from atomic models of macromolecules. The three-dimensional version of the Reference Interaction Site Model (RISM) from liquid-state statistical mechanics is employed to compute the solvent distribution around the solute, including both water and ions. X-ray scattering profiles are computed from this distribution together with the solute geometry. We describe an efficient procedure for performing this calculation employing a Lebedev grid for the angular averaging. The intensity profiles (which involve no adjustable parameters) match experiment and molecular dynamics simulations up to wide angle for two proteins (lysozyme and myoglobin) in water, as well as the small-angle profiles for a dozen biomolecules taken from the BioIsis.net database. The RISM model is especially well-suited for studies of nucleic acids in salt solution. Use of fiber-diffraction models for the structure of duplex DNA in solution yields close agreement with the observed scattering profiles in both the small and wide angle scattering (SAXS and WAXS) regimes. In addition, computed profiles of anomalous SAXS signals (for Rb{sup +} and Sr{sup 2+}) emphasize the ionic contribution to scattering and are in reasonable agreement with experiment. In cases where an absolute calibration of the experimental data at q = 0 is available, one can extract a count of the excess number of waters and ions; computed values depend on the closure that is assumed in the solution of the Ornstein–Zernike equations, with results from the Kovalenko–Hirata closure being closest to experiment for the cases studied here.

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

Systematic and simulation-free coarse graining of homopolymer melts: A structure-based study

We propose a systematic and simulation-free strategy for coarse graining of homopolymer melts, where each chain of N{sub m} monomers is uniformly divided into N segments, with the spatial position of each segment corresponding to the center-of-mass of its monomers. We use integral-equation theories suitable for the study of equilibrium properties of polymers, instead of many-chain molecular simulations, to obtain the structural and thermodynamic properties of both original and coarse-grained (CG) systems, and quantitatively examine how the effective pair potentials between CG segments and the thermodynamic properties of CG systems vary with N. Our systematic and simulation-free strategy is much faster than those using many-chain simulations, thus effectively solving the transferability problem in coarse graining, and provides the quantitative basis for choosing the appropriate N-values. It also avoids the problems caused by finite-size effects and statistical uncertainties in many-chain simulations. Taking the simple hard-core Gaussian thread model [K. S. Schweizer and J. G. Curro, Chem. Phys. 149, 105 (1990)] as the original system, we demonstrate our strategy applied to structure-based coarse graining, which is quite general and versatile, and compare in detail the various integral-equation theories and closures for coarse graining. Our numerical results show that the effective CG potentials for various N and closures can be collapsed approximately onto the same curve, and that structure-based coarse graining cannot give thermodynamic consistency between original and CG systems at any N < N{sub m}.

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

Vulnerability of Karangkates dams area by means of zero crossing analysis of data magnetic

Study with entitled Vulnerability Karangkates Dam Area By Means of Zero Crossing Analysis of Data Magnetic has been done. The study was aimed to obtain information on the vulnerability of two parts area of Karangkates dams, i.e. Lahor dam which was inaugurated in 1977 and Sutami dam inaugurated in 1981. Three important things reasons for this study are: 1). The dam age was 36 years old for Lahor dam and 32 years old for Sutami dam, 2). Geologically, the location of the dams are closed together to the Pohgajih local shear fault, Selorejo local fault, and Selorejo limestone-andesite rocks contact plane, and 3). Karangkates dams is one of the important Hydro Power Plant PLTA with the generating power of about 400 million KWH per year from a total of about 29.373MW installed in Indonesia. Geographically, the magnetic data acquisition was conducted at coordinates (112.4149oE;-8.2028oS) to (112.4839oE;-8.0989oS) by using Proton Precession Magnetometer G-856. Magnetic Data acquisition was conducted in the radial direction from the dams with diameter of about 10 km and the distance between the measurements about 500m. The magnetic data acquisition obtained the distribution of total magnetic field value in the range of 45800 nT to 44450 nT. Residual anomalies obtained by doing some corrections, including diurnal correction, International Geomagnetic Reference Field (IGRF) correction, and reductions so carried out the distribution of the total magnetic field value in the range of -650 nT to 700 nT. Based on the residual anomalies, indicate the presence of 2 zones of closed closures dipole pairs at located in the west of the Sutami dam and the northwest of the Lahor dam from 5 total zones. Overlapping on the local geological map indicated the lineament of zero crossing patterns in the contour of residual anomaly contour with the Pohgajih shear fault where located at about 4 km to the west of the Sutami dam approximately and andesite-limestone rocks contact where located at about 6 km to the west of the Lahor dam approximately. These shown a possible of vulnerability on geohazards at the west zone of the Karangkates (Lahor-Sutami) dams area if there are triggers by the vibration (earthquake) on the Pohgajih shear fault, andesite-limestone contact plane, and instability rocks on two zones of closed closure dipole pairs area. Reality, on the location of the study shown some local landslide at the several locations and the main road that need considering for disaster mitigation.

58 GEOSCIENCES↗
DOE OSTI2015

A validation of the fibre orientation and fibre length attrition prediction for long fibre-reinforced thermoplastics

To improve the mechanical performance of polymeric parts, fibre reinforcement has established in industrial applications during the last decades. Next to the widely used Short Fibre-reinforced Thermoplastics (SFT) the use of Long Fibre-reinforced Thermoplastics (LFT) is increasingly growing. Especially for non-polar polymeric matrices like polypropylene (PP), longer fibres can significantly improve the mechanical performance. As with every kind of discontinuous fibre reinforcement the fibre orientations (FO) show a high impact on the mechanical properties. On the contrary to SFT where the local fibre length distribution (FLD) can be often neglected, for LFT the FLD show a high impact on the material’s properties and has to be taken into account in equal measure to the FOD. Recently numerical models are available in commercial filling simulation software and allow predicting both the local FOD and FLD in LFT parts. The aim of this paper is to compare i.) the FOD results and ii) the FLD results from available orientation- and fibre length attrition-models to those obtained from experimental data. The investigations are conducted by the use of different injection moulded specimens made from long glass fibre reinforced PP. In order to determine the FOD, selected part sections are examined by means of Computed Tomographic (CT) analyses. The fully three dimensional measurement of the FOD is then performed by digital image processing using grey scale correlation. The FLD results are also obtained by using digital image processing after a thermal pyrolytic separation of the polymeric matrix from the fibres. Further the FOD and the FLD are predicted by using a reduced strain closure (RSC) as well as an anisotropic rotary diffusion - reduced strain closure model (ARD-RSC) and Phelps-Tucker fibre length attrition model implemented in the commercial filling software Moldflow, Autodesk Inc., San Rafael, CA, USA.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Ion thermal and dispersion effects in Farley-Buneman instabilities

Farley-Buneman modes are an example of the collisional instability, which is thought to be the dominant mechanism for the irregularities in low ionosphere region. Despite high collisionality due to electron-neutral and ion-neutral collisions, the kinetic effects associated with finite temperature are important for determination of the mode frequencies and growth rate. This is especially important for ion component that is largely unmagnetized due to low ion cyclotron frequency. The ion thermal effects are strongly pronounced for shorter wavelengths and are crucial for the growth rate cut-off at high wavenumbers. We develop an extended fluid model for ion dynamics to incorporate the effects of ion thermal motion. The model is based on the extended MHD model that includes the evolution equations for higher order moments such as ion viscosity and ion heat flux. We also develop the generalized Chapman-Enskog closure model that provides exact linear closures based on the linearized kinetic equation. The results of these models are compared and tested against the linear kinetic model. The dispersion of Farley-Buneman modes and growth rate behavior are investigated in the short wavelength region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗