Search NASA⌕ Search

SEARCH · Search NASA

Results for “Finite deformation theory”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Global properties of nuclei at finite-temperature within the covariant energy density functional theory

In stellar environments nuclei appear at finite temperatures, becoming extremely hot in core-collapse supernovae and neutron-star mergers. However, due to theoretical and computational complexity, most model calculations of nuclear properties are performed at zero temperature, while those existing at finite temperatures are limited only to selected regions of the nuclide chart. In this study we perform the global calculation of nuclear properties for even-even 8 ≤ Z ≤ 104 nuclei at temperatures in range 0 ≤ T ≤ 2 MeV. Calculations are based on the finite-temperature relativistic Hartree-Bogoliubov model supplemented by the Bonche-Levit-Vautherin vapor subtraction procedure. We find that near the neutron-drip line the continuum states have significant contribution already at moderate temperature T ≈ 1 MeV, thus emphasizing the necessity of the vapor subtraction procedure. Results include neutron emission lifetimes, quadrupole deformations, neutron-skin thickness, proton and neutron pairing gaps, entropy and excitation energy. Up to the temperature T ≈ 1 MeV, the nuclear landscape is influenced only moderately by the finite-temperature effects, mainly by reducing the pairing correlations. Here, as the temperature increases further, the effects on nuclear structures become pronounced, reducing both the deformations and the shell effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

GW190814: Impact of a 2.6 solar mass neutron star on the nucleonic equations of state

Is the secondary component of GW190814 the lightest black hole or the heaviest neutron star ever discovered in a double compact-object system? This is the central question animating this paper. Covariant density functional theory provides a unique framework to investigate both the properties of finite nuclei and neutron stars, while enforcing causality at all densities. By tuning existing energy density functionals we were able to: (i) account for a 2.6 M ⊙ neutron star, (ii) satisfy the original constraint on the tidal deformability of a 1.4 M ⊙ neutron star, and (iii) reproduce ground-state properties of finite nuclei. Yet, for the class of models explored in this work, we find that the stiffening of the equation of state required to support supermassive neutron stars is inconsistent with either constraints obtained from energetic heavy-ion collisions or from the low deformability of medium-mass stars. Furthermore, we speculate that the maximum neutron star mass can not be significantly higher than the existing observational limit and that the 2.6 M ⊙ compact object is likely to be the lightest black hole ever discovered.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

DiscoFluxM version 1.x

DiscoFluxM is software for modeling the deformation behavior of crystalline materials within the third-party finite-element based software, MOOSE. The software enables full-field crystal plasticity finite element simulations of deformation processes of metals and other crystalline materials. DiscoFluxM is an implementation of the DiscoFlux theoretical and algorithmic framework for coupling single-crystal plasticity with dislocation transport. The continuum mesoscale modeling framework is founded upon a nonlocal dislocation-density based crystal plasticity theory. The nonlocal theory couples continuum dislocation transport with an otherwise spatially-local single-crystal model employing nonlinear thermoelasticity and crystallographic plasticity. Dislocation transport is modeled by enforcing dislocation conservation at a slip-system level through the solution of advection-diffusion equations. The configuration of geometrically necessary dislocation density gives rise to a back-stress that inhibits or accentuates the flow of dislocations.

Luscher, Darby↗
DOE OSTI2024

Noniterative finite amplitude methods for giant resonances and the application to the neutron radiative capture cross sections

We calculate the electric dipole (E1) and the magnetic dipole (M1) giant resonances with noniterative finite amplitude methods and demonstrate how the fully microscopic density functional theory predicts the giant resonances without any phenomenological parameters. Then, we calculate neutron capture reactions based on the statistical Hauser-Feshbach theory with the result of E1 and M1 transitions and find that the capture cross sections for deformed nuclei are enhanced due to the contribution from the low energy M1 scissors mode.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Strain localization criteria for viscoplastic geomaterials

This report presents a viscoplastic localization criterion to detect quasi-instantaneous (i.e., load-induced) and delayed (creep-induced) strain localization in rate-dependent solids. The study is based on the theory of controllability and a viscoplastic description of the mechanical response. Analytical precursors of unstable states are defined through systems of ordinary differential equations (OEDs). The use of the proposed criteria is illustrated at the material point level through a set of strain localization analyses simulating active strain localization of a porous rock. In addition, full-field finite element simulations of compression tests conducted under various pressures are reported to demonstrate the role of local unstable viscoplasticity in the spontaneous propagation of deformation bands under stationary boundary conditions. The study shows that the viscoplastic localization criterion maintains a negative sign as long as the behavior is unstable, that is, the rate of deformation is accelerating. The sign switch coincides with the transition to decelerating deformation. The analyses revealed that pulses of overstress always emerge in correspondence with the growth of unstable behavior, and the peak matches the transition to stable behavior. The local responses recovered from full-field analyses were consistent with those observed in analyses at material point level and the predictions of the presented theory.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Three-dimensional Skyrme Hartree-Fock-Bogoliubov solver in coordinate-space representation

The coordinate-space representation of the Hartree-Fock-Bogoliubov theory is the method of choice to study weakly bound nuclei whose properties are affected by the quasiparticle continuum space. To describe such systems, we developed a three-dimensional Skyrme-Hartree-Fock-Bogoliubov solver HFBFFT based on the existing, highly optimized and parallelized Skyrme-Hartree-Fock code Sky3D. The code does not impose any self-consistent spatial symmetries such as mirror inversions or parity. The underlying equations are solved in HFBFFT directly in the canonical basis using the fast Fourier transform. To remedy the problems with pairing collapse, we implemented the soft energy cutoff and pairing annealing. The convergence of HFB solutions was improved by a sub-iteration method. The Hermiticity violation of differential operators brought by Fourier-transform-based differentiation has also been solved. Furthermore, the accuracy and performance of HFBFFT were tested by benchmarking it against other HFB codes, both spherical and deformed, for a set of nuclei, both well-bound and weakly-bound.

3D coordinate-space representation↗
DOE OSTI2024

Coulomb Branch Amplitudes from a Deformed Amplituhedron Geometry

The amplituhedron provides, via geometric means, the all-loop integrand of scattering amplitudes in maximally supersymmetric Yang-Mills theory. Unfortunately, dimensional regularization, used conventionally for integration, breaks the beautiful geometric picture. This motivates us to propose a “deformed” amplituhedron. Focusing on the four-particle amplitude, we introduce two deformation parameters, which can be interpreted as particle masses. We provide evidence that the mass pattern corresponds to a specific choice of vacuum expectation values on the Coulomb branch. The deformed amplitude is infrared finite, making the answer well defined in four dimensions. Leveraging four-dimensional integration techniques based on differential equations, we compute the amplitude up to two loops. In the limit where the deformation parameters are taken to zero, we recover the known Bern-Dixon-Smirnov amplitude. In the limit where only one deformation parameter is taken to zero, we find a connection to the angle-dependent cusp anomalous dimension. Published by the American Physical Society 2024

Physics↗
DOE OSTI2020

Thermomechanical conversion in metals: dislocation plasticity model evaluation of the Taylor-Quinney coefficient

Using a partitioned-energy thermodynamic framework which assigns energy to that of atomic configurational stored energy of cold work and kinetic-vibrational, in this study we derive an important constraint on the Taylor-Quinney coefficient, which quantifies the fraction of plastic work that is converted into heat during plastic deformation. Associated with the two energy contributions are two separate temperatures – the ordinary temperature for the thermal energy and the effective temperature for the configurational energy. We show that the Taylor-Quinney coefficient is a function of the thermodynamically defined effective temperature that measures the atomic configurational disorder in the material. Finite-element analysis of recently published experiments on the aluminum alloy 6016-T4 [1], using the thermodynamic dislocation theory (TDT), shows good agreement between theory and experiment for both stress-strain behavior and temporal evolution of the temperature. The simulations include both conductive and convective thermal energy loss during the experiments, and significant thermal gradients exist within the simulation results. Computed values of the differential Taylor-Quinney coefficient are also presented and suggest a value which differs between materials and increases with increasing strain.

36 MATERIALS SCIENCE↗
DOE OSTI2026

ACT University Collaboration Proposal: Role of manufacturing defects on material failure under dynamic loading for developing enhanced failure models and theories (Final Report)

Research on the tensile behavior of additively manufactured 316L stainless steel coupon specimens at increasing strain rates was conducted over the past year at Penn State. Dynamic loading rates into the 1000 1/s loading rates were performed on specimens with nearly 1.32” gage lengths. Stress-strain plots show ductile and plastic behavior well beyond 12% strain with necked specimens, having smaller effective gage lengths showing up to 56% ultimate strain. Additional tests performed on compact tension specimens helped with simulation work to understand the deformation behavior. Using finite elements, it was possible to determine the feasibility of a comprehensive experimental validation study towards an improved model for AM material failure – specifically, the Bai-Wierzbicki approach which accounts for lode angle and triaxiality. A non-significant number of tests are projected for eight specific failure nodes, with additional replicates to provide strain-rate capabilities to the existing formulation. One such correction factor is explored for quasi-static, notched specimens. Finally, stress intensity factor was explored using a set of compact-tension specimens which would also provide useful validation data for any simulation work.

36 MATERIALS SCIENCE↗
DOE OSTI2022

(U) Marker Lagrangian Simulations using Geometry from CAD Assembly Files [Slides]

The Eulerian code PAGOSA now contains an option for the forward time integration of any user selected material in the Lagrangian frame of reference. This option is accomplished using a Multifield version of Brackbill’s FLIP plus Sulsky’s MPM. FLIP+MPM is a Marker Lagrangian scheme that is not to be confused with well known Grid Lagrangian schemes. This Marker Lagrangian method allows materials to undergo completely arbitrary deformation, up to and including fracture and separation. Using this option makes PAGOSA into a mixed frame simulation tool such that fluids can be integrated in the Eulerian frame and solids integrated in the Lagrangian frame. Boundaries between the fields are said to be ’immersed’ in the computational domain, and their nature is formulated using Multifield Theory. Each Marker in a field of FLIP+MPM material tracks the thermodynamic state of a small (but finite) piece of mass. Initiation of the Markers using the geometry and material information from CAD Assembly files is the subject of the talk. The ease and speed with which the dynamics of the assembly can be studied is demonstrated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Critical Scaling of Solid Fragmentation at Quasistatic and Finite Strain Rates

Here, using two-dimensional simulations of sheared, brittle solids, we characterize the resulting fragmentation and explore its underlying critical nature. Under quasistatic loading, a power-law distribution of fragment masses emerges after fracture which grows with increasing strain. With increasing strain rate, the maximum size of a grain decreases and a shallower distribution is produced. We propose a scaling theory for distributions based on a fractal scaling of the largest mass with system size in the quasistatic limit or with a correlation length that diverges as a power of rate in the finite-rate limit. Critical exponents are measured using finite-size scaling techniques.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Quasiparticle random-phase approximation calculations for M 1 transitions with the noniterative finite-amplitude method and application to neutron radiative capture cross sections

Here, we derive the equations of quasiparticle random-phase approximation (QRPA) based on the finite amplitude method with the Hartree-Fock+Bardeen–Cooper–Schrieffer (HF+BCS) single-particle states, and calculate the magnetic dipole (M⁢1) transition for deformed gadolinium isotopes. Our QRPA calculation shows both large spin-flip transitions in the 5 to 10 MeV excitation energy and the low energy orbital transition that would correspond to the M⁢1 scissors mode observed in nuclear experiments. Then, we calculate neutron capture reactions based on the statistical Hauser-Feshbach theory with the photoabsorption cross sections of even-even nuclei given by QRPA. We find that the capture cross section is enhanced due to the contribution from the low energy M⁢1 transition although the calculated capture cross section still underestimates the experimental data. This issue in the calculated capture cross section could be improved by uncertainties of the low energy E⁢1 transition neglected in our QRPA calculation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Effect of static local distortions vs. dynamic motions on the stability and band gaps of cubic oxide and halide perovskites

Ternary ABX 3 perovskites made of corner-sharing BX 6 octahedra have long featured prominently in solid-state chemistry and condensed matter physics. Still, the joint understanding of their two main subgroups—halides and oxides—has not been fully developed. Indeed, unlike the case in simpler compounds having a single, robust repeated motif (“monomorphous”), certain cubic perovskites can manifest a non-thermal (= intrinsic) distribution of local motifs (“polymorphous networks”). Such static deformations can include positional degrees of freedom (e.g., atomic displacements and octahedral tilting) or magnetic moment degrees of freedom in paramagnets. Unlike thermal motion, such static distortions do not time-average to zero, being an expression of the intrinsic symmetry breaking preference of the chemical bonding. Here, the present study compares electronic structure features of oxide and halide perovskites starting from the static polymorphous distribution of motifs described by Density Functional Theory (DFT) minimization of the internal energy, continuing to finite temperature thermal disorder modeled via finite temperature DFT molecular dynamics. We find that (i) different oxide vs. halide ABX 3 compounds adopt different energy-lowering symmetry-breaking modes. The calculated pair distribution function (PDF) of SrTiO 3 from the first-principles agrees with recently measured PDF. (ii) In both oxides and halides, such static distortions lead to band gap blueshifts with respect to undistorted cubic Pm-3m structure. (iii) For oxide perovskites, high-temperature molecular dynamics simulations initiated from the statically distorted polymorphous structures reveal that the thermally-induced distortions can lead to a band gap redshift. (iv) In contrast, for cubic halide perovskite CsPbI 3 , both the intrinsic distortions and the thermal distortions contribute in tandem to band gap blueshift, the former, intrinsic effect being dominant. (v) In the oxide SrTiO 3 and CaTiO 3 (but not in halide) perovskites, octahedral tilting leads to the emergence of a distinct Γ–Γ direct band gap component as a secondary valley minimum to the well-known indirect R–Γ gap. Understanding such intrinsic vs. thermal effects on oxide vs. halide perovskites holds the potential for designing target electronic properties.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Chiral limit of 2d QCD revisited with lightcone conformal truncation

We study the chiral limit of 2d QCD with a single quark flavor at finite Nc using LCT. By modifying the LCT basis according to the quark mass in a manner motivated by ’t Hooft’s analysis, we are able to restore convergence for quark masses much smaller than the QCD strong coupling scale. For such small quark masses, the IR of the theory is expected to be well described by the Sine-Gordon model. We verify that LCT numerics are able to capture in detail the spectrum and correlation functions of the Sine-Gordon model. This opens up the possibility for studying deformations of various integrable CFTs using LCT by considering the chiral limit of QCD like theories.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Magnetic dipole γ-ray strength functions in the crossover from spherical to deformed neodymium isotopes

We calculate the magnetic dipole $\gamma$-ray strength functions in a chain of even-mass neodymium isotopes $^{144-152}$Nd in the framework of the configuration-interaction (CI) shell model. We infer the strength function by applying the maximum entropy method (MEM) to the exact imaginary-time response function calculated with the shell-model Monte Carlo (SMMC) method. The success of the MEM depends on the choice of a good strength function as a prior distribution. We investigate two choices for the prior strength function: the static path approximation (SPA) and the quasiparticle random-phase approximation (QRPA). We find that the QRPA is a better approximation at low temperatures (i.e., near the ground state), while the SPA is a better choice at finite temperatures. We identify a low-energy enhancement (LEE) in the MEM deexcitation $M1$ strength functions of the even-mass neodymium isotopes and compare with recent experimental results for the total deexcitation $\gamma$-ray strength functions. The LEE is already seen in the SPA strength function but not in the QRPA strength function, indicating the importance of large-amplitude static fluctuations around the mean field in reproducing the LEE. Our method is currently the only one which can reproduce LEE in heavy open-shell nuclei where conventional CI shell model calculations are prohibited. With the onset of deformation as number of neutrons increases along the chain of neodymium isotopes, we observe that some of the LEE strength transfers to a low-energy excitation, which we interpret as a finite-temperature ``scissors'' mode. Here, we also observe a finite-temperature spin-flip mode.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Comments on the quantum field theory of the Coulomb gas formalism

The holomorphic Coulomb gas formalism, as developed by Feigin-Fuchs, Dotsenko-Fateev and Felder, is a set of rules for computing minimal model observables using free field techniques. We attempt to derive and clarify these rules using standard techniques of quantum field theory. We begin with a careful examination of the timelike linear dilaton. Although the background charge of the model breaks the scalar field’s continuous shift symmetry, the exponential of the action remains invariant under a discrete shift because the background charge is imaginary. Gauging this symmetry makes the dilaton compact and introduces winding modes into the spectrum. One of these winding operators corresponds to the anti-holomorphic completion of the BRST current first introduced by Felder, and the full left/right cohomology of this BRST charge isolates the irreducible representations of the Virasoro algebra within the degenerate Fock space of the linear dilaton. The “supertrace” in the BRST complex reproduces the minimal model partition function and exhibits delicate cancellations between states with both momentum and winding. The model at the radius $R=\sqrt{pp^{\prime }}$ has two marginal operators corresponding to the Dotsenko-Fateev “screening charges”. Deforming by them, we obtain a model that might be called a “BRST quotiented compact timelike Liouville theory”. The Hamiltonian of the zero-mode quantum mechanics of this model is not Hermitian, but it is PT -symmetric and exactly solvable. Its eigenfunctions have support on an infinite number of plane waves, suggesting an infinite reduction in the number of independent states in the full quantum field theory. Applying conformal perturbation theory to the exponential interactions reproduces the Coulomb gas calculations of minimal model correlation functions. In contrast to spacelike Liouville, these “resonance correlators” are finite because the zero mode is compact. We comment on subtleties regarding the reflection operator identification, as well as naive violations of truncation in correlators with multiple reflection operators inserted. This work is part of an attempt to understand the relationship between the JT model of two dimen- sional gravity and the worldsheet description of the (2 , p ) minimal string as suggested by Seiberg and Stanford.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

$T\overline{T}$-deformed free energy of the Airy model

Sharpening the correspondence of Jackiw-Teitelboim (JT) gravity and its dual matrix model description at a finite radial cutoff λ through the $T\overline{T}$ deformation is of interest. To proceed, we simplify the problem by considering the Airy model and deform Airy correlators in the same way as in $T\overline{T}$ -deformed JT gravity. We use those correlators to compute the annealed and quenched free energies for both λ > 0 and λ < 0 from an integral representation of the replica trick. At the leading order in λ and low temperatures, we confirm that the genus-zero quenched free energy monotonically decreases as a function of temperature when perturbation theory is valid. We then study the all-genus quenched free energy at low temperatures, where we discover and discuss subtleties due to non-perturbative effects in the Airy model, as well as the contributions from the non-perturbative branch under the $T\overline{T}$ deformation.

2D gravity↗
DOE OSTI2024

Approach and rotation of reconnecting topological defect lines in liquid crystal

Topological defects are a universal concept across many disciplines, such as crystallography, liquid-crystalline physics, low-temperature physics, cosmology, and even biology. In nematic liquid crystals, topological defects called disclinations have been widely studied. For their three-dimensional (3D) dynamics, however, only recently have theoretical approaches dealing with fully 3D configurations been reported. Further, recent experiments have observed 3D disclination line reconnections, a phenomenon characteristic of defect line dynamics, but detailed discussions were limited to the case of approximately parallel defects. In this paper, we focus on the case of two disclination lines that approach at finite angles and lie in separate planes, a more fundamentally 3D reconnection configuration. Observing and analyzing such reconnection events, we find the square-root law of the distance between disclinations and the decrease of the interdisclination angle over time. We compare the experimental results with theory and find qualitative agreement on the scaling of distance and angle with time, but quantitative disagreement on distance and angle relative mobilities. To probe this disagreement, we derive the equations of motion for systems with reduced twist constant and also carry out simulations for this case. These, together with the experimental results, suggest that deformations of disclinations may be responsible for the disagreement. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗