Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical Relativity”

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,693 records · Page 94

Warming and snow loss increase reliance on old groundwater in a Colorado River headwater

Atmospheric warming is reducing snowpack, with uncertain effects on mountainous streamflow, a crucial water resource. Despite limited historical observations of groundwater–streamflow interactions above 2,500 m, new measurements in the Upper Colorado River headwaters indicate declining groundwater storage that is dated decades to millennia old. Here we use integrated hydrologic modelling spanning water years 2015–2021 to determine whether the loss of old-age groundwater buffers streamflow during low-snow years and whether that loss is exacerbated with warming. Results show that old-groundwater contributions to streams remain relatively steady through time, unlike the more variable contributions from young groundwater. Numerical experiments of increased surface air temperatures (+2.5 °C and +4 °C) increase rain–snow fractions and evapotranspiration and decrease runoff ratio by 2–3% per degree Celsius increase. As streamflow declines with warming, the age of groundwater supporting it gets older, in part owing to intermediate-aged (1–3 year) groundwater declining twice as fast. Simulations show that water table depths at higher elevations (>3,700 m) decline disproportionately and fail to recover even during wet years. These findings suggest altered groundwater–streamflow interactions with warming and snow loss, with implications for water resources.

Siirila-Woodburn, Erica R. [Lawrence Berkeley Nati↗

Analytical nuclear gradients for the range-separated many-body dispersion model of noncovalent interactions

An accurate treatment of the long-range electron correlation energy, including van der Waals (vdW) or dispersion interactions, is essential for describing the structure, dynamics, and function of a wide variety of systems. Among the most accurate models for including dispersion into density functional theory (DFT) is the range-separated many-body dispersion (MBD) method [A. Ambrosetti et al., J. Chem. Phys., 2014, 140, 18A508], in which the correlation energy is modeled at short-range by a semi-local density functional and at long-range by a model system of coupled quantum harmonic oscillators. In this work, we develop analytical gradients of the MBD energy with respect to nuclear coordinates, including all implicit coordinate dependencies arising from the partitioning of the charge density into Hirshfeld effective volumes. To demonstrate the efficiency and accuracy of these MBD gradients for geometry optimizations of systems with intermolecular and intramolecular interactions, we optimized conformers of the benzene dimer and isolated small peptides with aromatic side-chains. We find excellent agreement with the wavefunction theory reference geometries of these systems (at a fraction of the computational cost) and find that MBD consistently outperforms the popular TS and D3(BJ) dispersion corrections. To demonstrate the performance of the MBD model on a larger system with supramolecular interactions, we optimized the C 60 @C 60 H 28 buckyball catcher host–guest complex. In our analysis, we also find that neglecting the implicit nuclear coordinate dependence arising from the charge density partitioning, as has been done in prior numerical treatments, leads to an unacceptable error in the MBD forces, with relative errors of ~20% (on average) that can extend well beyond 100%.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structural, electronic, vibrational, and dielectric properties of LaBGeO{sub 5} from first principles

Structural, electronic, vibrational, and dielectric properties of LaBGeO{sub 5} with the stillwellite structure are determined based on ab initio density functional theory. The theoretically relaxed structure is found to agree well with the existing experimental data with a deviation of less than 0.2%. Both the density of states and the electronic band structure are calculated, showing five distinct groups of valence bands. Furthermore, the Born effective charge, the dielectric permittivity tensors, and the vibrational frequencies at the center of the Brillouin zone are all obtained. Compared to existing model calculations, the vibrational frequencies are found in much better agreement with the published experimental infrared and Raman data, with absolute and relative rms values of 6.04 cm{sup −1}, and 1.81%, respectively. Consequently, numerical values for both the parallel and perpendicular components of the permittivity tensor are established as 3.55 and 3.71 (10.34 and 12.28), respectively, for the high-(low-)frequency limit.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Note: Precise radial distribution of charged particles in a magnetic guiding field

Current high precision beta decay experiments of polarized neutrons, employing magnetic guiding fields in combination with position sensitive and energy dispersive detectors, resulted in a detailed study of the mono-energetic point spread function (PSF) for a homogeneous magnetic field. A PSF describes the radial probability distribution of mono-energetic electrons at the detector plane emitted from a point-like source. With regard to accuracy considerations, unwanted singularities occur as a function of the radial detector coordinate which have recently been investigated by subdividing the radial coordinate into small bins or employing analytical approximations. In this note, a series expansion of the PSF is presented which can numerically be evaluated with arbitrary precision.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Degenerate variational integrators for magnetic field line flow and guiding center trajectories

Symplectic integrators offer many benefits for numerically approximating solutions to Hamiltonian differential equations, including bounded energy error and the preservation of invariant sets. Two important Hamiltonian systems encountered in plasma physics—the flow of magnetic field lines and the guiding center motion of magnetized charged particles—resist symplectic integration by conventional means because the dynamics are most naturally formulated in non-canonical coordinates. New algorithms were recently developed using the variational integration formalism; however, those integrators were found to admit parasitic mode instabilities due to their multistep character. This work eliminates the multistep character, and therefore the parasitic mode instabilities via an adaptation of the variational integration formalism that we deem “degenerate variational integration.” Both the magnetic field line and guiding center Lagrangians are degenerate in the sense that the resultant Euler-Lagrange equations are systems of first-order ordinary differential equations. We show that retaining the same degree of degeneracy when constructing discrete Lagrangians yields one-step variational integrators preserving a non-canonical symplectic structure. Lastly, numerical examples demonstrate the benefits of the new algorithms, including superior stability relative to the existing variational integrators for these systems and superior qualitative behavior relative to non-conservative algorithms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Probing the nature of dark matter with accreted globular cluster streams

ABSTRACT The steepness of the central density profiles of dark matter (DM) in low-mass galaxy haloes (e.g. dwarf galaxies) is a powerful probe of the nature of DM. We propose a novel scheme to probe the inner profiles of galaxy subhaloes using stellar streams. We show that the present-day morphological and dynamical properties of accreted globular cluster (GC) streams – those produced from tidal stripping of GCs that initially evolved within satellite galaxies and later merged with the Milky Way (MW) – are sensitive to the central DM density profile and mass of their parent satellites. GCs that accrete within cuspy cold dark matter (CDM) subhaloes produce streams that are physically wider and dynamically hotter than streams that accrete inside cored subhaloes. A first comparison of MW streams ‘GD-1’ and ‘Jhelum’ (likely of accreted GC origin) with our simulations indicates a preference for cored subhaloes. If these results hold up in future data, the implication is that either the DM cusps were erased by baryonic feedback, or their subhaloes naturally possessed cored density profiles implying particle physics models beyond CDM. Moreover, accreted GC streams are highly structured and exhibit complex morphological features (e.g. parallel structures and ‘spurs’). This implies that the accretion scenario can naturally explain the recently observed peculiarities in some of the MW streams. We also propose a novel mechanism for forming ‘gaps’ in stellar streams when the remnant of the parent subhalo (which hosted the GC) later passes through the GC stream. This encounter can last a longer time (and have more of an impact) than the random encounters with DM subhaloes previously considered, because the GC stream and its parent subhalo are on similar orbits with small relative velocities. Current and future surveys of the MW halo will uncover numerous faint stellar streams and provide the data needed to substantiate our preliminary tests with this new probe of DM.

(Galaxy:) globular clusters: individual↗

Novel Probes Project: Tests of gravity on astrophysical scales

The Novel Probes Project, an initiative to advance the field of astrophysical tests of the dark sector by creating a forum that connects observers and theorists, is introduced. This review focuses on tests of gravity and is intended to be of use primarily to observers, as well as theorists with an interest in the development of experimental tests. It is twinned with a separate upcoming review on dark matter self-interactions. The review focuses on astrophysical tests of gravity in the weak-field regime, ranging from stars to quasilinear cosmological scales. This regime is complementary to both strong-field tests of gravity and background and linear probes in cosmology. In particular, the nonlinear screening mechanisms that are an integral part of viable modified-gravity models lead to characteristic signatures, specifically on astrophysical scales. The potential of these probes is not limited by cosmic variance but comes with the challenge of building robust theoretical models of the nonlinear dynamics of stars, galaxies, and large-scale structure. The groundwork is laid for a thorough exploration of the weak-field, nonlinear regime, with an eye to using the current and next generation of observations for tests of gravity. The scene is set by showing how gravitational theories beyond general relativity are expected to behave, focusing primarily on screening mechanisms. Analytic and numerical techniques for exploring the relevant astrophysical regime are described, as are the pertinent observational signals. With these in hand a range of astrophysical tests of gravity are presented, and prospects for future measurements and theoretical developments are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

Assessing hybrid vigour using the thermal sensitivity of physiological trade-offs in tiger salamanders

Hybridization between species affects biodiversity and population sustainability in numerous ways, many of which depend on the fitness of the hybrid relative to the parental species. Hybrids can exhibit fitter phenotypes compared to the parental lineages, and this ‘hybrid vigour’ can then lead to the extinction of one or both parental lines. In this study, we analysed the relationship between water loss and gas exchange to compare physiological performance among three tiger salamander genotypes—the native California tiger salamander (CTS), the invasive barred tiger salamanders (BTS) and CTS × BTS hybrids across multiple temperatures (13.5°C, 20.5°C and 23.5°C). We developed a new index of performance, the water-gas exchange ratio (WGER), which we define as the ratio of gas exchange to evaporative water loss (μL V O 2 /μL H 2 O). The ratio describes the ability of an organism to support energetically costly activities with high levels of gas exchange while simultaneously limiting water loss to lower desiccation risk. We used flow through respirometry to measure the thermal sensitivity of metabolic rate and resistance to water loss of each salamander genotype to compare indices of physiological performance. We found that temperature had a significant effect on metabolic rate and resistance to water loss, with both traits increasing as temperatures warmed. Across genotypes, we found that hybrids have a higher WGER than the native CTS, owing to a higher metabolic rate despite having a lower resistance to water loss. These results provide a greater insight into the physiological mechanisms driving hybrid vigour and offer a potential explanation for the rapid spread of salamander hybrids. More broadly, our introduction of the WGER may allow for species- or lineage-wide comparisons of physiological performance across changing environmental conditions, highlighting the insight that can be gleaned from multitrait analysis of organism performance.

59 BASIC BIOLOGICAL SCIENCES↗

Verification and Validation of the FLAG Hydrocode for Impact Cratering Simulations

Verification and validation (V&V) are necessary processes to ensure accuracy of the computational methods used to solve problems key to vast numbers of applications and industries. Simulations are essential for addressing impact cratering problems, because these problems often exceed experimental capabilities. We show that the free Lagrange (FLAG) hydrocode, developed at Los Alamos National Laboratory (Los Alamos, NM), can be used for impact cratering simulations by verifying FLAG against two analytical models of aluminum-on-aluminum impacts at different impact velocities and validating FLAG against a glass-into-water laboratory impact experiment. Our verification results show good agreement with the theoretical maximum pressures, with relative errors as low in magnitude as 1.00%. Our validation results demonstrate FLAG's ability to model various stages of impact cratering, with crater radius relative errors as low as 3.48% and crater depth relative errors as low as 0.79%. Our mesh resolution study shows that FLAG converges at resolutions low enough to reduce the required computation time from about 28 h to about 25 min. We anticipate that FLAG can be used to model larger impact cratering problems with increased accuracy and decreased computational cost on current systems relative to other hydrocodes tested by Pierazzo et al. (2008, “Validation of Numerical Codes for Impact and Explosion Cratering: Impacts on Strengthless and Metal Targets,” MAPS, 43(12), pp. 1917–1938).

97 MATHEMATICS AND COMPUTING↗

Acarbose impairs gut Bacteroides growth by targeting intracellular glucosidases

ABSTRACT Acarbose is a type 2 diabetes medicine that prevents dietary starch breakdown into glucose by inhibiting host amylase and glucosidase enzymes. Numerous gut species in theBacteroidesgenus enzymatically break down starch and change in relative abundance within the gut microbiome in acarbose-treated individuals. To mechanistically explain this observation, we used two model starch-degradingBacteroides,Bacteroides ovatus(Bo), andBacteroides thetaiotaomicron(Bt). Bt growth on starch polysaccharides is severely impaired by acarbose, whereas Bo growth is much less affected by the drug. TheBacteroidesuse a starch utilization system (Sus) to grow on starch. We hypothesized that Bo and Bt Sus enzymes are differentially inhibited by acarbose. Instead, we discovered that although acarbose primarily targets the Sus periplasmic GH97 enzymes in both organisms, the drug affects starch processing at multiple other points. Acarbose competes for transport through the TonB-dependent SusC proteins and binds to the Sus transcriptional regulators. Furthermore, Bo expresses a non-Sus GH97 (BoGH97D) when grown in starch with acarbose. The Bt homolog, BtGH97H, is not expressed in the same conditions, nor can overexpression of BoGH97D complement the Bt growth inhibition in the presence of acarbose. This work informs us about unexpected complexities of Sus function and regulation inBacteroides, including variation between related species. Furthermore, this indicates that the gut microbiome may be a source of variable response to acarbose treatment for diabetes. IMPORTANCE Acarbose is a type 2 diabetes medication that works primarily by stopping starch breakdown into glucose in the small intestine. This is accomplished by the inhibition of host enzymes, leading to better blood sugar control via reduced ability to derive glucose from dietary starches. The drug and undigested starch travel to the large intestine where acarbose interferes with the ability of some bacteria to grow on starch. However, little is known about how gut bacteria interact with acarbose, including microbes that can use starch as a carbon source. Here, we show that two gut species,Bacteroides ovatus(Bo) andBacteroides thetaiotaomicron(Bt), respond differently to acarbose: Bt growth is inhibited by acarbose, while Bo growth is less affected. We reveal a complex set of mechanisms involving differences in starch import and sensing behind the different Bo and Bt responses. This indicates the gut microbiome may be a source of variable response to acarbose treatment for diabetes via complex mechanisms in common gut microbes.

Microbiology↗

Nonadiabatic effects in (1~2){sup 1}Π–X{sup 1}Σ{sup +} rovibronic transitions of KRb molecules

Numerical simulation of the effect of intramolecular electrostatic interactions on redistribution of relative intensities in the vibrational structure of (1 ~ 2){sup 1}Π–X{sup 1}Σ{sup +} rotationally resolved transitions of the KRb molecule is performed within the precision nonadiabatic model of coupled vibrational channels. It is established that mutual perturbation of electronically excited states modifies in a nontrivial way a nodal structure of nonadiabatic wavefunction of the (1 ~ 2){sup 1}Π complex, which is possible to use for rising efficiency of twostep laser synthesis and stabilization of ultracold ensembles of KRb molecules in the ground electronic state.

74 ATOMIC AND MOLECULAR PHYSICS↗

Conductance distribution near the Anderson transition

Using a modification of the Shapiro approach, we introduce the two-parameter family of conductance distributions W(g), defined by simple differential equations, which are in the one-to-one correspondence with conductance distributions for quasi-one-dimensional systems of size L{sup d–1} × L{sub z}, characterizing by parameters L/ξ and L{sub z}/L (ξ is the correlation length, d is the dimension of space). This family contains the Gaussian and log-normal distributions, typical for the metallic and localized phases. For a certain choice of parameters, we reproduce the results for the cumulants of conductance in the space dimension d = 2 + ϵ obtained in the framework of the σ-model approach. The universal property of distributions is existence of two asymptotic regimes, log-normal for small g and exponential for large g. In the metallic phase they refer to remote tails, in the critical region they determine practically all distribution, in the localized phase the former asymptotics forces out the latter. A singularity at g = 1, discovered in numerical experiments, is admissible in the framework of their calculational scheme, but related with a deficient definition of conductance. Apart of this singularity, the critical distribution for d = 3 is well described by the present theory. One-parameter scaling for the whole distribution takes place under condition, that two independent parameters characterizing this distribution are functions of the ratio L/ξ.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

project_0vbb v1.0-beta

This software is designed for the sole purpose of analyzing numerical lattice QCD results computed for one of the developers' research projects related to neutrinoless double beta decay, and will be distributed with their research publication. This software is provided to maximize the dissemination of the research results, to enable other scientists and the general public to easily check the results of their analysis, and also to be in compliance with the public availability of "digital data" as required by the DOE.

Walker-Loud, Andree↗

Atmospheric Flow Development and Associated Changes in Turbulent Sensible Heat Flux over a Patchy Mountain Snow Cover

In this study, the small-scale boundary layer dynamics and the energy balance over a fractional snow cover are numerically investigated. The atmospheric boundary layer flows over a patchy snow cover were calculated with an atmospheric model (Advanced Regional Prediction System) on a very high spatial resolution of 5 m. The numerical results revealed that the development of local flow patterns and the relative importance of boundary layer processes depend on the snow patch size distribution and the synoptic wind forcing. Energy balance calculations for quiescent wind situations demonstrated that well-developed katabatic winds exerted a major control on the energy balance over the patchy snow cover, leading to a maximum in the mean downward sensible heat flux over snow for high snow-cover fractions. This implies that if katabatic winds develop, total melt of snow patches may decrease for low snow-cover fractions despite an increasing ambient air temperature, which would not be predicted by most hydrological models. In contrast, stronger synoptic winds increased the effect of heat advection on the catchment’s melt behavior by enhancing the mean sensible heat flux over snow for lower snow-cover fractions. A sensitivity analysis to grid resolution suggested that the grid size is a critical factor for modeling the energy balance of a patchy snow cover. The comparison of simulation results from coarse (50 m) and fine (5 m) horizontal resolutions revealed a difference in the spatially averaged turbulent heat flux over snow of 40%–70% for synoptic cases and 95% for quiescent cases.

17 WIND ENERGY↗

Deformed Materials: Towards a Theory of Materials Morphology Dynamics

This grant supported work on the response of crystals to external stress. Our primary work described how disordered structural materials break in two (statistical models of fracture in disordered materials), studied models of deformation bursts (avalanches) that mediate deformation on the microscale, and developed continuum dislocation dynamics models for plastic deformation (as when scooping ice cream bends a spoon, Fig. 9). Glass is brittle -- it breaks with almost atomically smooth fracture surfaces. Many metals are ductile -- when they break, the fracture surface is locally sheared and stretched, and it is this damage that makes them hard to break. Bone and seashells are made of brittle material, but they are strong because they are disordered -- lots of little cracks form as they are sheared and near the fracture surface, diluting the external force. We have studied materials like bone and seashells using simulations, mathematical tools, and statistical mechanics models from physics. In particular, we studied the extreme values of fracture strengths (how likely will a beam in a bridge break far below its design strength), and found that the traditional engineering tools could be improved greatly. We also studied fascinating crackling-noise precursors -- systems which formed microcracks of a broad range of sizes before they broke. Ductile metals under stress undergo irreversible plastic deformation -- the planes of atoms must slide across one another (through the motion of dislocations) to change the overall shape in response to the external force. Microscopically, the dislocations in crystals move in bursts of a broad range of sizes (termed 'avalanches' in the statistical mechanics community, whose motion is deemed 'crackling noise'). In this grant period, we resolved a longstanding mystery about the average shape of avalanches of fixed duration (using tools related to an emergent scale invariance), we developed the fundamental theory describing the shapes of avalanches and how they are affected by the edges of the microscope viewing window, we found that slow creep of dislocations can trigger an oscillating response explaining recent experiments, we explained avalanches under external voltage, and we have studied how avalanches in experiments on the microscale relate to deformation of large samples. Inside the crystals forming the metal, the dislocations arrange into mysterious cellular structures, usually ignored in theories of plasticity. Writing a natural continuum theory for dislocation dynamics, we found that it spontaneously formed walls -- much like models of traffic jams and sonic booms. These walls formed rather realistic cellular structures, which we examined in great detail -- our walls formed fractal structures with fascinating scaling properties, related to those found in turbulent fluids. We found, however, that the numerical and mathematical tools available to solve our equations were not flexible enough to incorporate materials-specific information, and our models did not show the dislocation avalanches seen experimentally. In the last year of this grant, we wrote an invited review article, explaining how plastic flow in metals shares features with other stressed materials, and how tools of statistical physics used in these other systems might be crucial for understanding plasticity.

36 MATERIALS SCIENCE↗

Pronghorn Theory Manual

This manual is intended to provide a thorough description of the governing equations, the relations used for material, heat transfer, and fluid flow parameters, and the numerical method utilized in Pronghorn.

42 ENGINEERING↗

Collaborative Research: The Evolution of Magnetic Skeletons During 3D Reconnection. Final report

We investigated how the plasmoid instability affects the evolution of topological features of a magnetic field and the process of magnetic reconnection with three sub-projects. (i) We extensively studied the criterion for plasmoid-mediated current sheet disruption and the onset of fast reconnection. Scaling relations with respect to plasma parameters have been established, using both direct numerical simulations and analytic theories, for the disruption time, the current sheet width, the instability linear growth rate, and the dominant mode wavenumber. (ii) We applied the theories of plasmoid-mediated current sheet disruption to magnetohydrodynamic turbulence and predicted a steepening of the turbulence energy spectrum. This steepening has been confirmed by direct numerical simulations. (iii) We performed numerical simulations of ultraviolet burst events in the solar atmosphere with a spine-dome magnetic topology. The plasmoid instability leads to unsteady reconnection, which is consistent with solar observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗