Search NASA⌕ Search

SEARCH · Search NASA

Results for “partially”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

High-pressure elastic properties of dolomite melt supporting carbonate-induced melting in deep upper mantle

Deeply subducted carbonates likely cause low-degree melting of the upper mantle and thus play an important role in the deep carbon cycle. However, direct seismic detection of carbonate-induced partial melts in the Earth’s interior is hindered by our poor knowledge on the elastic properties of carbonate melts. Here we report the first experimentally determined sound velocity and density data on dolomite melt up to 5.9 GPa and 2046 K by in-situ ultrasonic and sink-float techniques, respectively, as well as first-principles molecular dynamics simulations of dolomite melt up to 16 GPa and 3000 K. Using our new elasticity data, the calculated V P /V S ratio of the deep upper mantle (~180–330 km) with a small amount of carbonate-rich melt provides a natural explanation for the elevated V P /V S ratio of the upper mantle from global seismic observations, supporting the pervasive presence of a low-degree carbonate-rich partial melt (~0.05%) that is consistent with the volatile-induced or redox-regulated initial melting in the upper mantle as argued by petrologic studies. This carbonate-rich partial melt region implies a global average carbon (C) concentration of 80–140 ppm. by weight in the deep upper mantle source region, consistent with the mantle carbon content determined from geochemical studies.

58 GEOSCIENCES↗
DOE OSTI2026

The Physics Imposed on a Streaming Operator by Spherical Transport Problems

The streaming operator, which generates a displacement of a particle on a straight line at a constant speed in transport theory, is derived algebraically from a spherical coordinate formulation of Newton’s second law. This derivation leads to an operator that has more partial derivatives than a Cartesian coordinate formulation of the operator. The additional partial derivatives, which are with respect to the normalized velocity variables of a particle, take into account the intrinsic curvature of a ball. Moreover, these partial derivatives mitigate ray effects, which arise when a finite number of normalized velocities (also called directions or discrete ordinates) are used to simulate a continuous S 2 sphere of directions, by rotating the polar axis of the S 2 sphere into the radial direction of the coordinate system. As a consequence of this rotation, the number of actual discrete ordinates is greatly amplified to an enormous number of effective discrete ordinates by a multiplier that is equal to the number of patches that partitions a spherical surface. In addition to the derivation of the streaming operator, we provide in closed form a solution to the system of characteristic equations that is equivalent to the streaming operator. Furthermore, the solution to the system of characteristic equations enables the construction of an integral operator that is the inverse to the streaming operator. Examples in which ray effects are immensely mitigated by spherical coordinates are presented.

integral operator↗
DOE OSTI2023

Uniqueness of relaxation times determined by dielectric spectroscopy

Dielectric spectroscopy is extremely powerful to study molecular dynamics, because of the very broad frequency range. Often multiple processes superimpose resulting in spectra that expand over several orders of magnitude, with some of the contributions partially hidden. For illustration, we selected two examples, (i) normal mode of high molar mass polymers partially hidden by conductivity and polarization and (ii) contour length fluctuations partially hidden by reptation using the well-studied polyisoprene melts as example. The intuitive approach to describe experimental spectra and to extract relaxation times is the addition of two or more model functions. Here, we use the empirical Havriliak-Negami function to illustrate the ambiguity of the extracted relaxation time, despite an excellent agreement of the fit with experimental data. We show that there are an infinite number of solutions for which a perfect description of experimental data can be achieved. However, a simple mathematical relationship indicates uniqueness of the pairs of the relaxation strength and relaxation time. Sacrificing the absolute value of the relaxation time enables to find the temperature dependence of the parameters with a high accuracy. For the specific cases studied here, the time temperature superposition (TTS) is very useful to confirm the principle. However, the derivation is not based on a specific temperature dependence, hence, independent from the TTS. We compare new and traditional approaches and find the same trend for the temperature dependence. The important advantage of the new technology is the knowledge of the accuracy of the relaxation times. Relaxation times determined from data for which the peak is clearly visible are the same within the experimental accuracy for traditional and new technology. However, for data where a dominant process hides the peak, substantial deviations can be observed. Finally, we conclude that the new approach is particularly helpful for cases in which relaxation times need to be determined without having access to the associated peak position.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Aperiodically ordered nano-graphene on the quasicrystalline substrate

Designing exotic structures in low dimensions is key in today's quest to tailor novel quantum states in materials with unique symmetries. Particularly intriguing materials in this regard are low dimensional aperiodic structures with non-conventional symmetries that are otherwise forbidden in translation symmetric crystals. In our work, we focus on the link between the structural and electronic properties of aperiodically ordered aromatic molecules on a quasicrystalline surface, which has largely been neglected so far. As an exemplary case, we investigate the self-assembly and the interfacial electronic properties of the nano-graphene-like molecule coronene on the bulk truncated icosahedral (i) Al–Pd–Mn quasicrystalline surface using multiple surface sensitive techniques. We find an aperiodically ordered coronene monolayer (ML) film on the i-Al–Pd–Mn surface that is characterized by the same local motifs of the P1 Penrose tiling model as the bare i-Al–Pd–Mn surface. The electronic valence band structure of the coronene/i-Al–Pd–Mn system is characterized by the pseudogap of thebare i-Al–Pd–Mn, which persists the adsorption of coronene confirming the quasiperiodic nature of the interface. In addition, we find a newly formed interface state of partial molecular character that suggests an at least partial chemical interaction between the molecule and the quasicrystalline surface. We propose that this partial chemical molecule–surface interaction is responsible for imprinting the quasicrystalline order of the surface onto the molecular film.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Comparison of divertor behavior and plasma confinement between argon and neon seeding in EAST

The exhaust of excessively high heat and particle fluxes on the divertor target is crucial for EAST long-pulse operation. In the recent EAST experiments, stable partial energy detachment around the upper outer strike point with H 98,y2 ~ 1 was achieved with either Ne or Ar seeding from the upper outer divetor target in the upper single null configuration with ITER-like tungsten divertor. With either Ar or Ne seeding, the electron temperature around the upper outer strike point (T et,UOSP ) was maintained at around 5 eV, the peak temperature of divertor target surface around the upper outer strike point (T div,UO ) decreased significantly, and material sputtering was well suppressed. Here, it was observed that there was less Ar seeding needed for partial energy detachment onset than Ne seeding, which shows that Ar is more efficient in the cooling of T et on the upper outer divertor than Ne. However, there was no detachment on the upper inner divertor with T et around strike point (T et,UISP ) remaining >10 eV with either Ar or Ne seeding from the upper outer divertor. Accompanied with the disappearance of double peak phenomenon of ion flux density on the upper inner divertor target (j s,UI ), the peak T div,UI around the strike point increased to around 300 °C. Although the heat flux on the upper inner divertor target (q t,UI ) is still in the acceptable level, either Ar or Ne seeding only from the upper outer divertor target is not enough to protect the upper inner divertor target from sputtering under current EAST conditions. On the other hand, Ar seeding always causes confinement degradation in the partial energy detachment state. It was observed that there is a slight confinement improvement (~10%) with Ne seeding, which may be due to density peaking, dilution effects and stabilization of the ion temperature gradient mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Sparse Cholesky factorization for solving nonlinear PDEs via Gaussian processes

In recent years, there has been widespread adoption of machine learning-based approaches to automate the solving of partial differential equations (PDEs). Among these approaches, Gaussian processes (GPs) and kernel methods have garnered considerable interest due to their flexibility, robust theoretical guarantees, and close ties to traditional methods. They can transform the solving of general nonlinear PDEs into solving quadratic optimization problems with nonlinear, PDE-induced constraints. However, the complexity bottleneck lies in computing with dense kernel matrices obtained from pointwise evaluations of the covariance kernel, and its partial derivatives, a result of the PDE constraint and for which fast algorithms are scarce. The primary goal of this paper is to provide a near-linear complexity algorithm for working with such kernel matrices. We present a sparse Cholesky factorization algorithm for these matrices based on the near-sparsity of the Cholesky factor under a novel ordering of pointwise and derivative measurements. The near-sparsity is rigorously justified by directly connecting the factor to GP regression and exponential decay of basis functions in numerical homogenization. We then employ the Vecchia approximation of GPs, which is optimal in the Kullback-Leibler divergence, to compute the approximate factor. This enables us to compute ϵ-approximate inverse Cholesky factors of the kernel matrices with complexity O(N log d (N/ϵ)) in space and O(N log 2d (N/ϵ)) in time. We integrate sparse Cholesky factorizations into optimization algorithms to obtain fast solvers of the nonlinear PDE. We numerically illustrate our algorithm’s near-linear space/time complexity for a broad class of nonlinear PDEs such as the nonlinear elliptic, Burgers, and Monge-Ampère equations. In summary, we provide a fast, scalable, and accurate method for solving general PDEs with GPs and kernel methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Computational properties of finite PNmatrices

Abstract Recent compositionality results in logic have highlighted the advantages of enlarging the traditional notion of logical matrix semantics, namely by incorporating non-determinism and partiality. Still, several important properties which are known to be computable for finite logical matrices have not been studied in the wider context of partial non-deterministic matrices (PNmatrices). In this paper, we study how incorporating non-determinism and/or partiality in logical matrices impacts on the computational properties of some natural problems regarding their induced logics and concretely their sets of theorems. We show that, while for some of these problems there is no relevant computational impact, there are problems whose computational complexity increases and still other problems that simply become undecidable. In particular, we show that the problem of checking whether the logics characterized by two finite PNmatrices have the same set of theorems is not decidable. This undecidability result explores the connection between PNmatrices and term-DAG-automata, where the universality problem is known to be undecidable. This link also motivates a final contribution, in the form of a pumping-like lemma, which can be used, in some cases, to show that a given logic cannot be characterized by a finite PNmatrix.

Filipe, Pedro↗
DOE OSTI2021

Singular value decomposition and similarity renormalization group evolution of nuclear interactions

One of the main challenges for ab initio nuclear many-body theory is the growth of computational and storage costs as calculations are extended to heavy, exotic, and structurally complex nuclei. Here, we investigate the factorization of nuclear interactions as a means to address this issue. We perform Singular Value Decompositions of nucleon-nucleon interactions in partial wave representation and study the dependence of the singular value spectrum on interaction characteristics like regularization scheme and resolution scales. We develop and implement the Similarity Renormalization Group (SRG) evolution of the factorized interaction, and demonstrate that this SVD-SRG approach accurately preserves two-nucleon observables. We find that low-resolution interactions allow the truncation of the SVD at low rank, and that a small number of relevant components is sufficient to capture the nuclear interaction and perform an accurate SRG evolution, while the Coulomb interaction requires special consideration. The rank is uniform across all partial waves, and almost independent of the basis choice in the tested cases. This suggests an interpretation of the relevant singular components as mere representations of a small set of abstract operators that can describe the interaction and its SRG flow. Following the traditional workflow for nuclear interactions, we discuss how the transformation between the center-of-mass and laboratory frames creates redundant copies of the partial wave components when implemented in matrix representation, and we discuss strategies for mitigation. Lastly, we test the low-rank approximation to the SRG-evolved interactions in many-body calculations using the In-Medium SRG. By including nuclear radii in our analysis, we verify that the implementation of the SRG using the singular vectors of the interaction does not spoil the evolution of other observables.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Mirror nucleon removal reactions in p-shell nuclei

Nucleon removal reactions have been shown to be an effective tool for studying the single particle structure of nuclei. This work continues efforts to experimentally probe and benchmark the reaction and structure models used to calculate the removal reaction cross sections when using microscopic nuclear structure inputs. Here, three different single nucleon removal reactions were performed, from p -shell nuclei with masses A = 7, 9, and 10. The residual nuclei from the reactions were detected in coincidence with γ rays to determine partial cross sections to individual final states. The eikonal direct-reaction model is combined with overlap functions and residual nucleus densities from microscopic, variational Monte Carlo calculations to provide consistent nuclear structure input to the partial cross section calculations. Comparisons of measured and calculated cross sections, including for mirror reactions, are presented. The analysis of the partial cross sections leading to the ground states shows a similar behavior to the one observed from analyses of inclusive cross sections using shell model nuclear structure input: the theoretical description of the removal process is in better agreement with the data when removing weakly bound nucleons, than when removing well-bound ones. The two mirror reaction pairs presented here show consistent results between the respective members of the pairs. The results obtained for the population of the excited states, however, show a systematically different trend that appears connected to the structure part of the calculation. Additional cases are needed to better understand the respective roles of structure and dynamical effects in the deviations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Two-pion exchange as a leading-order contribution in chiral effective field theory

Pion exchange is the central ingredient to nucleon-nucleon interactions used in nuclear structure calculations, and one-pion exchange (OPE) enters at leading order in chiral effective field theory. In the 2S+1 L J = 1 S 0 partial wave, however, OPE and a contact term needed for proper renormalization fail to produce the qualitative, and quantitative, features of the scattering phase shifts. Cutoff variation also revealed a surprisingly low breakdown momentum Λ b ≈ 330 MeV in this partial wave. Here we show that potentials consisting of OPE, two-pion exchange (TPE), and a single contact address these problems and yield accurate and renormalization group (RG) invariant phase shifts in the 1 S 0 partial wave. Here, we demonstrate that a leading-order potential with TPE can be systematically improved by adding a contact quadratic in momenta. For momentum cutoffs Λ ≲ 500 MeV, the removal of relevant physics from TPE loops needs to be compensated by additional contacts to keep RG invariance. Inclusion of the Δ isobar degree of freedom in the potential does not change the strong contributions of TPE.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

First near-threshold measurements of the 13 C($α, n$ 1 ) 16 O reaction for low-background-environment characterization

The modeling of background sources in large volume detection systems requires accurate nuclear cross sections for a variety of reactions. Among the most important are (α,n) on light nuclei, where α particles, up to ≈9 MeV, are produced from the decay of actinides present in trace amounts in detection and structural material. In order to model the neutron energy spectra and production of other secondary particles, the partial cross sections are needed. Yet very little experimental data exists for these partial cross sections because past measurements, hampered by the experimental challenges of neutron detection, have focused mostly on total reaction cross section measurements. Here, in this paper, the partial cross section of the 13 C(α,n 1 ) 16 O reaction is reported for the first time. The measurements were made near the reaction threshold using a high energy resolution helium spectrometer. The measurements show a rapidly increasing cross section, which quickly becomes a significant fraction of the total reaction cross section. Measurements are compared with previous theory estimates and differences of more than an order of magnitude are found. A comparison is made with new total cross section measurements and the current level of consistency is demonstrated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Generalized susceptibilities and the properties of charm degrees of freedom across the QCD crossover temperature

We study the generalized charm susceptibilities in 2+1 flavor QCD on the lattice at several lattice spacings. We show that, below the chiral crossover, these susceptibilities are well described by the hadron resonance gas (HRG) model if charmed hadrons not listed in tables of the Particle Data Group are included. However, the HRG description abruptly breaks down just above the chiral crossover. To understand this, we use a model for the charm pressure in which it is expressed as the sum of partial pressures from charmed baryons, charmed mesons, and charm quarks. We present continuum estimates of these partial pressures and find that, while the partial pressures of charmed mesons and baryons drop below their respective HRG predictions, the charm quark pressure becomes nonzero above the chiral crossover.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Efficient end-to-end simulation of time-dependent coherent X-ray scattering experiments

Physical optics simulations for beamlines and experiments allow users to test experiment feasibility and optimize beamline settings ahead of beam time in order to optimize valuable beam time at synchrotron light sources like NSLS-II. Further, such simulations also help to develop and test experimental data processing methods and software in advance. The Synchrotron Radiation Workshop ( SRW ) software package supports such complex simulations. We demonstrate how recent developments in SRW significantly improve the efficiency of physical optics simulations, such as end-to-end simulations of time-dependent X-ray photon correlation spectroscopy experiments with partially coherent undulator radiation (UR). The molecular dynamics simulation code LAMMPS was chosen to model the sample: a solution of silica nanoparticles in water at room temperature. Real-space distributions of nanoparticles produced by LAMMPS were imported into SRW and used to simulate scattering patterns of partially coherent hard X-ray UR from such a sample at the detector. The partially coherent UR illuminating the sample can be represented by a set of orthogonal coherent modes obtained by simulation of emission and propagation of this radiation through the coherent hard X-ray (CHX) scattering beamline followed by a coherent-mode decomposition. GPU acceleration is added for several key functions of SRW used in propagation from sample to detector, further improving the speed of the calculations. The accuracy of this simulation is benchmarked by comparison with experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The three cellulose synthase isoforms for secondary cell wall make specific contributions to microfibril synthesis

Cellulose is synthesized at the plasma membrane by the cellulose synthase complex, a structure that contains three distinct isoforms of the catalytic subunit, cellulose synthase A (CESA). The division into three subunits appears early in land plant evolution and is highly conserved, particularly for the secondary cell wall. However, what if any unique roles each isoform plays in the complex remain unclear. Here, we assessed the contributions of specific isoforms to microfibril synthesis. First, we expressed CESA isoforms of the primary cell wall or the moss Physcomitrium patens in Arabidopsis thaliana backgrounds missing a secondary cell wall CESA. While the primary cell wall isoforms rescued the cesa knockout phenotype with partial isoform specificity, those from the moss rescued with fewer restrictions. Then, we recreated various CESA missense mutations in all three of the secondary cell wall isoforms; while results are consistent with isoform specificity, they are difficult to interpret further without molecular structures. Finally, we show that catalytically inactive CESA isoforms restore growth and cellulose content in the corresponding knockout in an isoform-specific manner; along with partial rescue of the growth and cellulose content of the inflorescence stem, the replacement lines have fiber cells with partially disorganized microfibrils and secondary cell wall cellulose with narrow crystal width. Generally, effects were more pronounced in lines where CESA8 was inactivated compared with inactivating CESA4 or 7, which tended to have similar phenotypes to each other. Here, we account for these results with a model for cellulose synthase structure with the isoforms assigned specific localization within the cellulose synthase complex.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Probing Thermal Transport in Fluidized Bed Using Modulated Photothermal Radiometry

Abstract In concentrated solar power (CSP) applications, fluidized bed is a promising approach for high heat transfer coefficient (HTC) solar receivers and heat exchangers. However, the complexity of multiphase mixing has made it difficult to characterize and analyze the heat transfer mechanism. This paper presents an experimental study on simultaneously characterizing heat transfer in both the near-wall and the bulk regions of a fluidized bed using modulated photothermal radiometry (MPR). The MPR is a non-contact frequency-domain technique using an intensity-modulated laser as the heat source and surface infrared emission as thermometry. The thermal penetration depth of the laser heating is varied by controlling its modulation frequency, and thus the measurement can resolve the near-wall and the bulk thermal resistances. With the MPR technique, we measured fluidized silica sands with a mean size of 164 μm in a vertical channel of 6 mm depth. Our results show that the near-wall thermal resistance is substantially increased with increasing gas velocity, which partially offsets the benefit of higher HTC brought by stronger particle mixing during the fluidization. We also used the MPR to quantify the improvement in particle-wall heat transfer in an inclined channel. We found that an 8° inclination towards the heat exchanging side led to a lower near-wall thermal resistance and a higher HTC at high gas velocities. This work demonstrates that the MPR technique is a useful tool to quantify the important near-wall thermal resistance from a bulk particle bed, which not only advances our understanding of heat transfer in fluidized beds, but may also contribute to the design of fluidized bed heat exchangers with higher HTC.

14 SOLAR ENERGY↗
DOE OSTI2024

Extracting and Interpreting Electrochemical Impedance Spectra (EIS) from Physics-Based Models of Lithium-Ion Batteries

This paper implements a highly efficient algorithm to extract electrochemical impedance spectra (EIS) from physics-based battery models (e.g., a P2D model). The mathematical approach is different from how EIS is practiced experimentally. Experimentally, the voltage (current) is harmonically perturbed over a wide range of frequencies and the amplitude and phase shift of the corresponding current (voltage) is measured. The experimental approach can be implemented in simulation software, but is computationally expensive. The approach here is to determine locally linear state-space models from the full physical model. The four Jacobian matrices that are the basis of the state-space models can be derived by numerical differentiation of the physical model. The EIS is then extracted from the state-space model using computationally efficient matrix-manipulation techniques. The algorithm can evaluate the full EIS at an instant in time during a transient, independent of whether the battery is in a stationary state. The approach is also able to separate the full-cell impedance to evaluate partial EIS, such as for a battery anode alone. Although such partial EIS is difficult to measure experimentally, the partial EIS provides valuable insights in interpreting the full-cell EIS.

25 ENERGY STORAGE↗
DOE OSTI2021

HyPerPy (Hydrogen Extraction and Parabolic Trough Plant Performance Models using Python) [SWR-21-53]

The HyPerPy package consists of three scripts: 1) Model for Hydrogen Tracking in Parabolic Trough Power Plants, 2) Model for Receiver Performance, and 3) Model for Hydrogen Extraction Process. The power plant model (1) tracks hydrogen generation and transport within the circulating heat transfer fluid (HTF) of the power plant. This script is a transient, initial value simulation, in which the hydrogen concentration in the circulating HTF is 0 moles per cubic meter everywhere at time 0 seconds. Hydrogen concentration is calculated at discrete locations within the circulating HTF with 4-second time resolution. During each time step, the change in hydrogen concentration due to hydrogen generation and permeation is calculated for each location according to the local HTF temperature, vessel or piping properties, hydrogen concentration and partial pressure. This model predicts hourly hydrogen concentrations for typical operating days in the spring, summer, fall, and winter seasons. Data is used to create hourly mappings of hydrogen concentrations for a typical operating year. The receiver model (2) uses a 1-hour time step to estimate getter loading and annulus pressure. For each time step, the model uses HTF and ambient temperature data to estimate absorber tube, bellows, and getter temperatures for the time step. In addition, the model uses hourly HTF hydrogen concentrations that are generated by the power plant model (1), and annulus hydrogen pressure from the previous time step. With these data, the model calculates the moles of hydrogen permeating across the absorber tube and bellows during the time step. The net change in moles hydrogen is added or subtracted from the getter loading for the previous time step, and the hydrogen pressure is re-calculated based on getter loading and temperature. The model uses this algorithm to simulate hydrogen permeation and loading 24 hours per day, 365 days per year using seasonal temperature data. The model repeats these calculations for 25 years to create a mapping of receiver getter loading and annulus hydrogen pressure for four seasons of each year. The model for hydrogen extraction (3) estimates hydrogen extraction rates for a specific separation module configuration. The rate depends primarily on membrane area, vacuum pump performance, headspace gas flowrate to the membrane, and headspace gas hydrogen partial pressure. This model has two versions. The steady-state version predicts hydrogen extraction rates when the module is operating in separation mode. The dynamic version predicts hydrogen transfer through the membrane when the module is operating in sensor mode. The steady-state version is used with the plant model (1) to predict hydrogen partial pressures in the power plant when the extraction process is operating.

Glatzmaier, Gregory↗