Lattice Spring Model & Simulation Based Parameter Inference
Presentation of mathematical content of SDRD project related to subsidence and shock physics
97 MATHEMATICS AND COMPUTING↗
SEARCH · Search NASA
Results for “Lattice Spring Model”
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.
A high-resolution pseudo-polygon discrete element model for regional sea ice
Here, this work presents a pseudo-polygon discrete element model for high-resolution sea ice simulations. A scale-invariant bonded particle contact model is proposed to model joints between sea ice floes based on the smeared fracture model and a lattice spring beam model, and the Mohr–Coulomb failure criterion is implemented to represent the shearing failure mechanism of sea ice packings under complex loadings. All mechanical parameters of the bond model can be directly determined from laboratory tests. Validations of the proposed model are made by investigations of mechanical response and failure criteria of field sea ice sheets. Compared with the field observations of sea ice from satellite radar and in situ stress sensors, the proposed model is capable of reproducing the typical constitutive behavior and the Coulomb friction envelope of field sea ice. Finally, the proposed discrete element sea ice model is used to study the effect of loading rates on mechanical behavior including failure strength of regional sea ice.
54 ENVIRONMENTAL SCIENCES↗
Strongly nonlinear wave propagation in elasto-plastic metamaterials: Low-order dynamic modeling
Nonlinear elastic metamaterials are known to support a variety of dynamic phenomena that enhance our capacity to manipulate elastic waves. Since these properties stem from complex, subwavelength geometry, full-scale dynamic simulations are often prohibitively expensive at scales of interest. Prior studies have therefore utilized low-order effective medium models, such as discrete mass-spring lattices, to capture essential properties in the long-wavelength limit. While models of this type have been successfully implemented for a wide variety of nonlinear elastic systems, they have predominantly considered dynamics depending only on the instantaneous kinematics of the lattice, neglecting history-dependent effects, such as wear and plasticity. Here, to address this limitation, the present study develops a lattice-based modeling framework for nonlinear elastic metamaterials undergoing plastic deformation. Due to the history- and rate-dependent nature of plasticity, the framework generally yields a system of differential-algebraic equations whose computational cost is significantly greater than an elastic system of comparable size. We demonstrate the method using several models inspired by classical lattice dynamics and continuum plasticity theory and explore means to obtain empirical plasticity models for general geometries, thereby gaining insight into the influence of microstructural plasticity on effective material performance, which can be used to improve the design of nonlinear mechanical metamaterials.
Dynamic simulation↗
Field-driven cluster formation in two-dimensional colloidal binary mixtures
Here, we study size- and charge-asymmetric oppositely charged colloids driven by an external electric field. The large particles are connected by harmonic springs, forming a hexagonal-lattice network, while the small particles are free of bonds and exhibit fluidlike motion. We show that this model exhibits a cluster formation pattern when the external driving force exceeds a critical value. The clustering is accompanied with stable wave packets in vibrational motions of the large particles.
42 ENGINEERING↗
Data-driven reduced-order models for port-Hamiltonian systems with operator inference
Hamiltonian operator inference has been developed in Sharma et al. (2022) to learn structure-preserving reduced-order models (ROMs) for Hamiltonian systems. The method constructs a low-dimensional model using only data and knowledge of the functional form of the Hamiltonian. The resulting ROMs preserve the intrinsic structure of the system, ensuring that the mechanical and physical properties of the system are maintained. In this work, we extend this approach to port-Hamiltonian systems, which generalize Hamiltonian systems by including energy dissipation, external input, and output. Based on snapshots of the system’s state and output, together with the information about the functional form of the Hamiltonian, reduced operators are inferred through optimization and are then used to construct data-driven ROMs. To further alleviate the complexity of evaluating nonlinear terms in the ROMs, a hyper-reduction method via discrete empirical interpolation is applied. Accordingly, we derive error estimates for the ROM approximations of the state and output. Lastly, we demonstrate the structure preservation, as well as the accuracy of the proposed port-Hamiltonian operator inference framework, through numerical experiments on a linear mass–spring-damper problem and a nonlinear Toda lattice problem.
97 MATHEMATICS AND COMPUTING↗
Strong Proton‐Phonon Coupling Drives Fast Ion Transport in Perovskites
Conduction of protons in solids is a cooperative process propelled by phonons, with molecular details obscured by the irregular movements in the thermal bath. It is shown that substitution with Y forms an imaginary phonon mode, instrumental for the function as proton conductor and effectively lowering the activation barrier for proton transport. To untangle the interplay in the exemplary proton conductor BaSn 0.9 Y 0.1 O 3 , its crystallographic structure is determined with high resolution neutron diffractometry and its phonon density of states with density functional theory calculations, experimentally validated by element specific nuclear resonant vibration spectroscopy. Based on phonon analysis, a quantitative transport model is present, which predicts the activation energy and performance by the ratio of ionic radii. Rather than individual vibrational modes, it is the oxygen sub-lattice which exerts its momentum on the protons. The extent of this momentum transfer is governed by the ratio of ionic radii. This model extends the transition state theory by the phonon-phonon interaction and complements the previously proposed idea that lattice dynamics is decisive for proton transport and specifies which properties of the material exactly define the vibration properties.
activation barriers↗
Weak phonon coupling to nematic quantum critical mode in BaFe 2 (As 1−𝑥 P 𝑥 ) 2
Here, in this work, we investigate the softening of the in-plane transverse acoustic phonon driven by electronic nematicity in BaFe 2 (As 1−𝑥 P 𝑥 ) 2 using inelastic x-ray scattering, with a focus on the optimally doped sample (𝑥 = 0.31) sample—a system exhibiting signatures of a putative nematic quantum critical point and minimal disorder among iron pnictides. We observe only a modest softening of the phonon frequency and no evidence of critical damping, suggesting that the nematic quantum critical fluctuations couple only weakly to the lattice from our quantum critical model. Given the close proximity of the structural and magnetic transition temperatures in the underdoped sample—which implies that spin-nematic fluctuations couple strongly to the lattice—we conjecture that the quantum critical nematic fluctuations are predominantly orbital in origin.
Wu, S. [University of California, Berkeley, CA (Un↗
Accessing the gluon momentum fraction of nucleons through the gradient flow
We calculate the gluon momentum fraction of the nucleon using lattice QCD, with a nonperturbative renormalization technique based on the gradient flow. The gluon momentum fraction is determined on a single Wilson-clover ensemble using 𝑁 𝑓 =2 +1 flavors with pion mass 358 MeV and lattice spacing 0.094 fm. We employ the variational method to reduce excited-state contamination and apply the distillation framework to ensure a large operator basis. To reduce systematic uncertainties, we apply Bayesian model averaging to all fit procedures. We apply matching coefficients to the flow-time dependent lattice results to recover the gluon momentum fraction in the $\overline{MS}$-scheme at 2 GeV. Our final result is ⟨𝑥⟩ 𝑔 (𝜇 =2 GeV) =0.482(35), where we quote only statistical uncertainties.
Lattice QCD↗
Inherent Anharmonicity of Harmonic Solids
Atomic vibrations, in the form of phonons, are foundational in describing the thermal behavior of materials. The possible frequencies of phonons in materials are governed by the complex bonding between atoms, which is physically represented by a spring-mass model that can account for interactions (spring forces) between the atoms (masses). The lowest-order, harmonic , approximation only considers linear forces between atoms and is thought incapable of explaining phenomena like thermal expansion and thermal conductivity, which are attributed to nonlinear, anharmonic , interactions. Here, we show that the kinetic energy of atoms in a solid produces a pressure much like the kinetic energy of atoms in a gas does. This vibrational or phonon pressure naturally increases with temperature, as it does in a gas and therefore results in a thermal expansion. Because thermal expansion thermodynamically defines a Grüneisen parameter γ , which is a typical metric of anharmonicity, we show that even a harmonic solid will necessarily have some anharmonicity. A consequence of this phonon pressure model is a harmonic estimation of the Grüneisen parameter as γ ≈ 3 / 2 3 − 4 x 2 / 1 + 2 x 2 , where x = v t / v l is the ratio of the transverse and longitudinal speeds of sound. We demonstrate the immediate utility of this model by developing a high-throughput harmonic estimate of lattice thermal conductivity that is comparable to other state-of-the-art estimations. By linking harmonic and anharmonic properties explicitly, this study provokes new ideas about the fundamental nature of anharmonicity, while also providing a basis for new material engineering design metrics.
42 ENGINEERING↗
Local structure, bonding, and asymmetry of $\mathrm{((NH_2)_2CH)PbBr_3}$, $\mathrm{CsPbBr_3}$, and $\mathrm{(CH_3MH_3)PbBr_3}$
Here we report local structure measurements for CsPbBr 3 and ((NH 2 ) 2 CH)PbBr 3 (FAPbBr 3 ) and compare them with recent results for (CH 3 NH 3 )PbBr 3 (MAPbBr 3 ). The Pb-Br bonding is similar for all three systems; the effective spring constants, κ, are comparable (ranging from 1.20 to 1.95 eV/Å 2 ), but small in magnitude indicating very soft materials. However, there are also important differences between the three systems. Static disorder is very small for CsPbBr3 but increases somewhat with the size of the organic molecular ions MA + and FA + . At room temperature, dynamic disorder dominates in all compounds. The thermal disorder of the Pb-Br pair distribution function (PDF), i.e., the Debye-Waller factor σ 2 follows a correlated Debye or Einstein model up to 300 K in CsPbBr 3 (orthorhombic phase), but for FAPbBr 3 and MAPbBr 3 , there is a break in the σ 2 (T) curve at the orthorhombic-tetragonal transition (o-t) near 150 K, indicating a small change in the spring constant κ. κ increases for MAPbBr 3 but decreases in FAPbBr 3 at this transition. These changes are attributed to changes in the H-bonding between Br - and MA + or FA + at this transition, as a result of librations or rotations of these molecular cations. In addition, the Pb-Br PDF becomes asymmetric at a relatively low temperature for FAPbBr 3 and MAPbBr 3 , while this effect is significantly smaller for CsPbBr 3 . Finally, we address the question of a model to explain the asymmetric PDF. Two main models are discussed in the literature, an anharmonic pair potential and a split-pair distribution, possibly driven by the presence of a lone pair on the Pb ion. We show that the fourth cumulant C 4 can differentiate between these two models and other possible models. Experimentally C 4 is positive at 250 K and above, for all three systems and that is inconsistent with a split-peak model, for which C 4 is negative for splittings larger than 0.12 Å.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗