Search NASA⌕ Search

SEARCH · Search NASA

Results for “explicit solution”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Fast explicit solutions for neutrino-electron scattering: Explicit asymptotic methods

Here, we present results of explicit asymptotic approximations applied to neutrino-electron scattering in a representative model of neutrino population evolution under conditions characteristic of core-collapse supernova explosions or binary neutron star mergers. It is shown that this approach provides stable solutions of these stiff systems of equations, with accuracy and time stepping comparable to that for standard implicit treatments such as backward Euler, fixed point iteration, and Anderson-accelerated fixed point iteration. Because each time step can be computed more rapidly with the explicit asymptotic approximation than with implicit methods, this suggests that algebraically stabilized explicit integration methods could be used to compute neutrino evolution coupled to hydrodynamics more efficiently in stellar explosions and mergers than the methods currently in use.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Explicit kink solutions in several one-parameter families of higher-order field theory models

Here, we present several one-parameter families of higher-order field theory models some of which admit explicit kink solutions with an exponential tail while others admit explicit kink solutions with a power-law tail. Various properties of these families of kink solutions are examined in detail. Further, by applying the recent Manton formalism, we provide estimates for the kink–kink (KK) and antikink–kink (AK-K) acceleration and hence the ratio of the corresponding AK-K and KK forces.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Osmotic Pressure Gradient Effects on Water Diffusion in Porous Rock: Can They Pervert Permeability Tests?

Abstract Generation of a large network of hydraulic cracks is of key importance not only for the success of fracking of shale but also for the recent scheme of sequestration of CO2 in deep formations of basalt and peridotite, which are mafic and ultramafic rocks that combine chemically with CO2. In numerical simulation of the creation of a fracture network in porous rock, an important goal is to enhance the rock permeability. The objective of this article is to calculate the effect of osmotic pressure gradients caused by gradients of concentration of the ions of Ca, Mg, Na, etc. on the effective permeability of the rock. The basic differential equations are formulated, and their explicit solutions for appropriate initial and boundary conditions are obtained under certain plausible simplifications. The main result is explicit approximate formulas for the critical time before which no water permeation through a test specimen can be observed. Depending on various parameters, this time can be unacceptably long, which is manifested as a zero water outflow. The solution may also explain the unreasonably small permeability values reported for some shales.

Mechanics↗
DOE OSTI2023

SCOMAP-XD : atomistic deuterium contrast matching for small-angle neutron scattering in biology

The contrast-variation method in small-angle neutron scattering (SANS) is a uniquely powerful technique for determining the structure of individual components in biomolecular systems containing regions of different neutron scattering length density ρ . By altering the ρ of the target solute and the solvent through judicious incorporation of deuterium, the scattering of desired solute features can be highlighted. Most contrast-variation methods focus on highlighting specific bulk solute elements, but not on how the scattering at specific scattering vectors q , which are associated with specific structural distances, changes with contrast. Indeed, many systems exhibit q -dependent contrast effects. Here, a method is presented for calculating both bulk contrast-match points and q -dependent contrast using 3D models with explicit solute and solvent atoms and SASSENA , an explicit-atom SANS calculator. The method calculates the bulk contrast-match points within 2.4% solvent D 2 O accuracy for test protein–nucleic acid and lipid nanodisc systems. The method incorporates a general model for the incorporation of deuterium at non-exchangeable sites that was derived by performing mass spectrometry on green fluorescent protein. The method also decomposes the scattering profile into its component parts and identifies structural features that change with contrast. The method is readily applicable to a variety of systems, will expand the understanding of q -dependent contrast matching and will aid in the optimization of next-generation neutron scattering experiments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Surface defects in holographic 5d SCFTs

We use holography to study codimension-2 surface defects in 5d SCFTs engineered by (p, q) 5-brane webs. The three-dimensional defects are realized by D3-branes ending on the brane web. We identify the holographic representation of the defects in Type IIB AdS 6 solutions as probe D3-branes, and study conformal and non-conformal defects which, respectively, preserve one half and one quarter of the supersymmetry. For a sample of 5d SCFTs, including the T N theories, we provide explicit solutions for conformal and non-conformal defects. For the conformal defects we obtain their contribution to the free energy on S 5 .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A Multicomponent Reactive Transport Model for Integrated Surface‐Subsurface Hydrology Problems

Abstract Despite the widespread use of integrated hydrology models in a variety of applications, consideration of multicomponent reactive transport is still not common. The implementation of these processes requires coupling transport at the surface‐subsurface interface and efficient solution of the non‐linear geochemical model that is consistent with the integrated hydrology solution. The Advanced Terrestrial Simulator provides a flexible multiphysics framework that facilitated this process. In this work, the integrated reactive transport process kernel (PK) was weakly coupled to the integrated hydrology PK. In turn, integrated transport and reactions were coupled using an operator splitting approach. This splitting enabled an explicit solution of the integrated transport problem, including a novel algorithm to calculate exchange fluxes across the surface‐subsurface interface and a point‐by‐point solution of the geochemical problem. Geochemical capabilities were added using well‐established external codes, but rather than using a custom interface to each, a generic interface was used that clearly specifies the variables and operations used by the chemistry PK. The implementation is demonstrated with two example simulations: transport of a tracer in a soil column as it saturates over time and water ponds on the surface and reactive transport in a hillslope driven by successive wet‐dry cycles that result in infiltration, runoff and exfiltration processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Analytic bounds on late-time axion-scalar cosmologies

The cosmological dynamics of multiple scalar/pseudoscalar fields are difficult to solve, especially when the field-space metric is curved. This presents a challenge in determining whether a given model can support cosmic acceleration, without solving for the on-shell solution. In this work, we present bounds on late-time FLRW-cosmologies in classes of theories that involve arbitrary numbers of scalar and pseudoscalar fields coupled both kinetically (leading to a curved field space metric) and through scalar potentials. Such bounds are proven analytically, independently of initial conditions, with no approximation in the field equations and without referring to explicit solutions. Besides their broad applications to cosmological model building, our bounds can be applied to studying asymptotic cosmologies of certain classes of string compactifications.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Negative complexity of formation: the compact dimensions strike back

We show that the vacuum-subtracted maximal volume, the proposed holographic dual to complexity of formation, can be negative when contributions from compact directions are included. We construct explicit solutions with arbitrarily negative complexity of formation in asymptotically AdS 4 × S 7 SUGRA. These examples rely critically on the compact directions, specifically the fact that the full eleven-dimensional spacetime is not asymptotically AdS 11 . While there is some ambiguity in the extension of the holographic complexity proposal to the compact directions, we show that the two natural candidates can both have arbitrarily negative complexity of formation in SUGRA solutions. We further find examples in which complexity can even decrease at late times, including cases of both single-sided geometries and two-sided wormholes. In particular, we construct a cosmological wormhole with simultaneously negative and decreasing complexity of formation (as computed by volume) at late times. We find a distinguished role for relevant primaries in these constructions and comment on possible interpretations.

2D Gravity↗
DOE OSTI2021

Current-conserving relativistic linear response for collisional plasmas

Here we investigate the response of a relativistic plasma to electromagnetic fields in the framework of the Boltzmann equation incorporating a collision term in the relaxation rate approximation selected in a form assuring current conservation. We obtain an explicit solution for the linearized perturbation of the Fermi–Dirac equilibrium distribution in terms of the average relaxation rate k. We study the resulting covariant, gauge invariant, and current conserving form of the polarization tensor in the ultrarelativistic and non-relativistic limits. We evaluate the susceptibility in the ultrarelativistic limit and explore their dependence on k. Finally, we study the dispersion relations for the longitudinal and transverse poles of the propagator. We show that for k >2ω p , where ω p is the plasma frequency, the plasma wave modes are overdamped. In the opposite case, k$\ll$ ω p , the propagating plasma modes are weakly damped.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

On Quantum Rainbows: Density Operator in the Frequency-Bin Representation for Entangled Twin-Photons Generated With Sub-Threshold Microcombs

Kerr optical frequency combs are generated by pumping a high-Q integrated microresonator with a resonant laser. Below threshold, the pump laser field mediates the phenomenon of spontaneous four-wave mixing, where two pump photons are symmetrically up- and down-converted as twin photons that can be entangled across up to tens of eigenmodes in the spectral domain. While these room-temperature integrated photonic circuits are expected to play a central role in quantum technology, their high dimensionality and dissipative nature are a challenge for their theoretical description, therefore hindering the understanding of their properties and potential of performance. In this article, we develop a framework that permits to obtain an explicit solution for the density operator of quantum microcombs below threshold. Furthermore, this self-consistent theoretical description allows for their complete characterization, as well as for the analytical determination of various performance metrics such as fidelity, purity, and entropy.

42 ENGINEERING↗
DOE OSTI2023

Explicit solvent machine-learned coarse-grained model of sodium polystyrene sulfonate to capture polymer structure and dynamics

In this article, strongly charged polyelectrolytes (PEs) demonstrate complex solution behavior as a function of chain length, concentrations, and ionic strength. The viscosity behavior is important to understand and is a core quantity for many applications, but aspects remain a challenge. Molecular dynamics simulations using implicit solvent coarse-grained (CG) models successfully reproduce structure, but are often inappropriate for calculating viscosities. To address the need for CG models which reproduce viscoelastic properties of one of the most studied PEs, sodium polystyrene sulfonate (NaPSS), we report our recent efforts in using Bayesian optimization to develop CG models of NaPSS which capture both polymer structure and dynamics in aqueous solutions with explicit solvent. We demonstrate that our explicit solvent CG NaPSS model with the ML-BOP water model [Chan et al. Nat Commun 10, 379 (2019)] quantitatively reproduces NaPSS chain statistics and solution structure. The new explicit solvent CG model is benchmarked against diffusivities from atomistic simulations and experimental specific viscosities for short chains. We also show that our Bayesian-optimized CG model is transferable to larger chain lengths across a range of concentrations. Overall, this work provides a machine-learned model to probe the structural, dynamic, and rheological properties of polyelectrolytes such as NaPSS and aids in the design of novel, strongly charged polymers with tunable structural and viscoelastic properties

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

First saturation correction in high energy proton-nucleus collisions. Part I. Time evolution of classical Yang-Mills fields beyond leading order

In high energy proton-nucleus collisions, the single- and double-inclusive soft gluon productions at the leading order have been calculated and phenomenologically studied in various approaches for many years. These studies do not take into account the saturation and multiple rescatterings in the field of the proton. The first saturation correction to these leading order results (the terms that are enhanced by the combination \( {\alpha}_s^2{\mu}^2 \) , where μ 2 is the proton’s color charge squared per unit transverse area) has not been completely derived despite recent attempts using a diagrammatic approach. This paper is the first in a series of papers towards analytically completing the first saturation correction to physical observables in high energy proton-nucleus collisions. Our approach is to analytically solve the classical Yang-Mills equations in the dilute-dense regime using the Color Glass Condensate effective theory and compute physical observables constructed from classical gluon fields. In the current paper, the Yang-Mills equations are solved perturbatively in the field of the dilute object (the proton). Next-to-leading order and next-to-next-to-leading order analytic solutions are explicitly constructed. A systematic way to obtain all higher order analytic solutions is outlined.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Kink solutions with power law tails

We present a comprehensive review about the various facets of kink solutions with a power law tail, which have received considerable attention during the last few years. This area of research is in its early stages; although several aspects have become clear by now, there are a number of issues which have only been partially understood or not understood at all. We first discuss the aspects which are reasonably well known and then address in some detail the issues which are only partially or not understood at all. We present a wide class of higher (than sixth) order field theory models admitting implicit kink as well as mirror kink solutions where the two tails facing each other have a power law or a power-tower type fall off, whereas the other two ends not facing each other could have either an exponential or a power law tail. The models admitting implicit kink solutions where the two ends facing each other have an exponential tail while the other two ends have a power law tail are also discussed. Moreover, we present several field theory models which admit explicit kink solutions with a power law fall off; we note that in all these polynomial models while the potential V ( ϕ ) is continuous, its derivative is discontinuous. We also discuss one of the most important and only partially understood issues of the kink–kink and the kink–antikink forces in case the tails facing each other have a power law fall off. Finally, we briefly discuss the kink–antikink collisions at finite velocity and present some open questions.

power-tower tail↗
DOE OSTI2024

Insights into the Assembly and Conformation of Nanoparticle Organic Hybrid Materials (NOHMs) in Solution with Varying Grafting Type

Recently, nanoparticle organic hybrid materials (NOHMs) have been considered promising electrolytes for energy storage applications due to their unique combination of properties, which include high thermal stability, negligible vapor pressure, and easy synthesis. However, the structural features of NOHMs in electrochemical systems are not well understood. Hence, we investigate the impact of grafting type (ionic vs covalent) on the structure and performance of NOHMs that consist of a silica core and grafted HPE polymer (Jeffamine M2070) in aqueous solution with and without the presence of the supporting electrolyte (0.1 M KHCO 3 ) using small-angle neutron scattering. Careful analyses of the results indicate that even at low grafting density (0.8 chains nm –2 ) and concentration (1–3 wt %), ionic NOHMs solutions contain a free polymer that may interact with the grafted polymer, altering the assembly of NOHMs in solution. Our investigation also elucidates that the covalent NOHM solutions consist of both aggregated and dispersed NOHMs in the solution. Moreover, the addition of the supporting electrolyte affects the assembly and structure of the polymer in the ionic NOHMs solution significantly, whereas only a slight change in the conformation of the grafted polymer is observed for covalent NOHMs. These conformational changes alter the performance of the NOHMs solution, resulting in a dramatic change in viscosity. Careful analysis shows that this significant alteration of solution viscosity in the ionic NOHMs can be traced back to the presence of the free polymer that interacts with the grafted layer and the attenuation of this interaction between the free and grafted polymers with the addition of salt. Furthermore, these studies elucidate the impact of grafting type on the overall structure, assembly, and performance of functionalized NOHMs in aqueous solution and explicitly correlate the structure and assembly of NOHMs to their transport properties and must be considered when designing functional systems that incorporate NOHMs.

25 ENERGY STORAGE↗
DOE OSTI2024

Zero modes of massive fermions delocalize from axion strings

Massless chiral excitations can arise from the interactions between a fermion and an axion string, propagating along the string and allowing it to superconduct. The properties of these excitations, or zero modes, dictate how the string interacts with light and can thus have important phenomenological consequences. In this paper, we add a nowhere-vanishing Dirac mass for the fermion in the usual model of axion electrodynamics. We find that the zero modes exhibit an interesting phase structure in which they delocalize from the string’s core as the mass increases, up until a critical value past which they disappear. We study this structure from an analytic perspective, with explicit numerical solutions, and via anomaly inflow arguments. Finally, we derive the two-dimensional effective theory of the zero mode and its interactions with the four-dimensional gauge field and show how this effective theory breaks down as the zero modes delocalize.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

A fast reduced model for a shell-and-tube based latent heat thermal energy storage heat exchanger and its application for cost optimal design by nonlinear programming

Numerical simulation of latent heat thermal energy storage (LHTES) systems plays a fundamental role in studying the physical process and guiding the engineering design. Discretization of the PDEs describing the nonlinear solidification/melting process of phase change materials (PCMs) leads to a large-scale complex dynamics system, where the system behavior depends on a set of parameters. In a design setting, repeated model evaluations are required over the set of parameters results in significant computational burden. In this paper, an explicit analytic solution was built for the propagation of the solidification front in a cylindrical coordinate. The analytic solution approach is further employed to develop a low computational reduced model (RM) as a module for a shell-and-tube based LHTES heat exchanger. The levelized Cost of Energy (LCOE) is used as a design metric and the RM model is used to apply system-level constraints in the nonlinear programming formulation that facilitates efficient global optimal design of the PCM properties, flow conditions and tube geometries. The use of LCOE as the design metric prevents over design of the heat transfer rate and also establishes a fair ground for evaluation of different thermal storage technologies and their integrated applications with other systems. Optimal results showed that a higher effectiveness results in a higher LCOE; the velocity of the HTF and the length of the channel are highly correlated with each other; both larger PCM latent energy and conductivity result in lower LCOE.

25 ENERGY STORAGE↗
DOE OSTI2026

Mathematical modelling of the concave front in the adjacent high explosive detonation problem

This study presents an analysis of the transition-zone in adjacent high explosive (HE) detonation problems which uses a $D, 𝜅, \dot{D}$ relationship, where $D$ is the detonation front-normal velocity, 𝜅 is the detonation front curvature and $\dot{D}$ is the time derivative of detonation front-normal velocity. Our approach extends the traditional $(D, 𝜅)$ model to accurately predict the behaviour of both diverging and converging detonation shock fronts. Our findings affirm that a hyperbolic type of front evolution equation, enhanced with wave acceleration, provides a robust framework for modelling complex shock front dynamics in HE materials. This approach not only captures the natural effects of straightness and boundary slope jumps in the transition-zone but also bridges the gap between mathematical predictions and experimental observations, offering insights into the behaviour of both diverging and converging detonation propagations in a homogeneous HE.

acceleration↗
DOE OSTI2024

Simulating hindered grain boundary diffusion using the smoothed boundary method

Abstract Grain boundaries can greatly affect the transport properties of polycrystalline materials, particularly when the grain size approaches the nanoscale. While grain boundaries often enhance diffusion by providing a fast pathway for chemical transport, some material systems, such as those of solid oxide fuel cells and battery cathode particles, exhibit the opposite behavior, where grain boundaries act to hinder diffusion. To facilitate the study of systems with hindered grain boundary diffusion, we propose a model that utilizes the smoothed boundary method to simulate the dynamic concentration evolution in polycrystalline systems. The model employs domain parameters with diffuse interfaces to describe the grains, thereby enabling solutions with explicit consideration of their complex geometries. The intrinsic error arising from the diffuse interface approach employed in our proposed model is explored by comparing the results against a sharp interface model for a variety of parameter sets. Finally, two case studies are considered to demonstrate potential applications of the model. First, a nanocrystalline yttria-stabilized zirconia solid oxide fuel cell system is investigated, and the effective diffusivities are extracted from the simulation results and are compared to the values obtained through mean-field approximations. Second, the concentration evolution during lithiation of a polycrystalline battery cathode particle is simulated to demonstrate the method’s capability.

Materials Science↗