Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model deformation”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Evidence of Quadrupole and Octupole Deformations in Zr 96 + Zr 96 and Ru 96 + Ru 96 Collisions at Ultrarelativistic Energies

In the hydrodynamic model description of heavy-ion collisions, the elliptic flow v 2 and triangular flow v 3 are sensitive to the quadrupole deformation β 2 and octupole deformation β 3 of the colliding nuclei. The relations between v n and β n have recently been clarified and were found to follow a simple parametric form. The STAR Collaboration has just published precision v n data from isobaric Ru 96 + Ru 96 and Zr 96 + Zr 96 collisions, where they observe large differences in central collisions v 2 , Ru > v 2 , Zr and v 3 , Ru < v 3 , Zr . Using a transport model simulation, we show that these orderings are a natural consequence of $\beta_{2,Ru} \gg \beta_{2,Zr}$ and $\beta_{3,Ru} \ll \beta_{3,Zr}$ We reproduce the centrality dependence of the v 2 ratio qualitatively and v 3 ratio quantitatively and extract values of β 2 and β 3 that are consistent with those measured at low-energy nuclear structure experiments. STAR data provide the first direct evidence of strong octupole correlations in the ground state of Zr 96 in heavy-ion collisions. Our analysis demonstrates that flow measurements in high-energy, heavy-ion collisions, especially using isobaric systems, are a new precision tool to study nuclear structure physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Cavitation in a soft porous material

Abstract We study the collapse and expansion of a cavitation bubble in a deformable porous medium. We develop a continuum-scale model that couples compressible fluid flow in the pore network with the elastic response of a solid skeleton. Under the assumption of spherical symmetry, our model can be reduced to an ordinary differential equation that extends the Rayleigh–Plesset equation to bubbles in soft porous media. The extended Rayleigh–Plesset equation reveals that finite-size effects lead to the breakdown of the universal scaling relation between bubble radius and time that holds in the infinite-size limit. Our data indicate that the deformability of the porous medium slows down the collapse and expansion processes, a result with important consequences for wide-ranging phenomena, from drug delivery to spore dispersion.

Leng, Yu↗
DOE OSTI2020

Structure of Nd 155 and Gd 163 from Cf 252 spontaneous fission

Background: A puzzle has arisen recently caused by the apparent shift in maximum deformation from the expected 66 Dy isotopic chain to the Nd 60 isotopic chain in the 82 155 Nd and 163 Gd, useful for constraining parameters in models that seek to answer the six proton shift in maximum deformation. Method: Data from the spontaneous fission of 252 Cf were taken by the Gammasphere detector array at Lawrence Berkeley National Laboratory to observe the excited states of 155 Nd and 163 Gd. Results: The structure of 163 Gd has been expanded with the addition of two new levels and three new γ rays, which are found to be consistent with previously published calculations and the structure of 165 Dy. In 155 Nd, nine new levels and 12 new γ rays are observed. The spins and parities of the previously known levels in 155 Nd have been reassigned from a ν3/2 - [521] ground state configuration to a ν5/2 + [642] isomeric configuration by comparison of these newly observed levels with levels in 153 Nd and 155 Sm. Conclusion: Further experimentation is required to determine the energy of the newly reassigned ν5/2+[642] level in 155 Nd with respect to the suspected ν3/2-[521] ground state. Additionally, more experiments should be conducted to further determine the structure of neutron rich nuclei, rarely produced in the spontaneous fission of 252 Cf, such as 163 Gd.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Self-equilibrated backstresses induce compensation between hardening and softening: Micromechanical and microstructural features

For crystal plasticity formulation neglecting the strain gradients at the dislocation scale, backstresses are necessary to satisfy the mesoscopic stress equilibrium and strain compatibility requirements. Backstresses in a material point and those between neighboring points are treated separately. Physically backstresses can dynamically induce either hardening or softening effects to adjust the dislocation hardening of slip systems, whereas the sum of the stresses resulting in backstresses in the representative volume element (RVE) is zero. Backstress hardening suppresses certain slip activities, while backstress softening triggers excess dislocations. A two-phase brass is studied. Backstresses make the hard phase deform smaller and the soft phase deform more than those without employing backstresses in a model. In the former, more homogeneous lattice strain and higher yield stress are obtained. Backstress has limited effects on the evolution of the critical resolved shear stress and crystallographic texture. Furthermore, it is a major cause of the banded structures featuring sharp crystallographic contrasts which are not parallel to the slip planes. Tests indicate that backstress induces a serrated behavior at the early stages of deformation and induces the geometrically necessary bands (GNBs), which are parallel to the macroscopic axes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Modeling active nematics via the nematic locking principle

Active nematic systems consist of rod-like internally driven subunits that interact with one another to form large-scale coherent flows. They are important examples of far-from-equilibrium fluids, which exhibit a wealth of nonlinear behavior. This includes active turbulence, in which topological defects in the nematic order braid around one another in a chaotic fashion. One of the most studied examples of active nematics consists of a dense two-dimensional layer of microtubules, crosslinked by kinesin molecular motors that inject extensile deformations into the fluid. Though numerous theoretical studies have modeled microtubule-based active nematics, no consensus has emerged on how to fully and quantitatively capture the features of the experimental system. Here, to better understand the theoretical foundations for modeling this system, we propose a fundamental principle we call the nematic locking principle—individual microtubules cannot rotate without all neighboring microtubules also rotating. Physically, this is justified by the high density of the microtubules, their elongated nature, and their corresponding steric interactions. We assert that the nematic locking principle holds throughout the majority of the material, but breaks down in the neighborhood of topological defects and other regions of low density. We derive the most general nematic transport equation consistent with this principle and also derive the most general term that violates it, introducing fracturing into the material. We then examine the standard Beris–Edwards approach, commonly used to model this system, and show that it violates the nematic locking principle throughout the majority of the material due to fracturing. We then propose a modification to the Beris–Edwards model that enforces nematic locking nearly everywhere. This modification shuts off fracturing except in regions where the order parameter (a proxy for density) is reduced. In these regions fracturing is turned on. The resulting simulations in turn show strong nematic locking throughout the bulk of the material, with narrow bands of fracturing, consistent with experimental observation. One additional advantage of enforcing nematic locking is that nontrivial stationary state solutions, common in Beris–Edwards simulations but not seen in experiments, are eliminated.

Mitchell, Kevin A. [Univ. of California, Merced, C↗
DOE OSTI2023

Gapped lineon and fracton models on graphs

We introduce a $\mathbb{Z}_N$ stabilizer code that can be defined on any spatial lattice of the form $\Gamma\times C_{L_z}$, where $\Gamma$ is a general graph. We also present the low-energy limit of this stabilizer code as a Euclidean lattice action, which we refer to as the anisotropic $\mathbb{Z}_N$ Laplacian model. It is gapped, robust (i.e., stable under small deformations), and has lineons. Its ground state degeneracy (GSD) is expressed in terms of a "mod $N$-reduction" of the Jacobian group of the graph $\Gamma$. In the special case when space is an $L\times L\times L_z$ cubic lattice, the logarithm of the GSD depends on $L$ in an erratic way and grows no faster than $O(L)$. Furthermore, we also discuss another gapped model, the $\mathbb{Z}_N$ Laplacian model, which can be defined on any graph. It has fractons and a similarly strange GSD.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

A Modified Network Alteration Model to Predict Quasi-Static Behavior of the Cross-Linked Polymers During Hydrolytic Aging

Abstract The network alteration theory is utilized to capture the effects of hydrolytic aging on the mechanical integrity of the polymer matrix. Hydrolytic degradation generally considered as a result of two simultaneous chemical phenomena; (i) reduction of the cross-links, and (ii) active chain loss. Assuming the alteration of the chain density is identical to the change of the number average molecular weight along the aging trajectory, the strain energy of polymer matrix can be rewritten as a function of deformation, deformation history, storage time, and aging temperature. The proposed model includes only five physically inspired material parameters. Thus, while it is computationally efficient, it shows good agreement with own experimental data, which performed on various range of accelerated aging temperatures and times. With respect to its computational efficiency, simplicity, accuracy, and interpret-ability, the model is a good choice for advanced implementations in FE programs.

Bahrololoumi, Amir↗
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 OSTI2020

Crystallographic orientation-dependent strain hardening in a precipitation-strengthened Al-Cu alloy

While the strengthening of Al-Cu alloys due to precipitation has been extensively studied, the effect of crystallographic orientation of the matrix and precipitates, as well as precipitate morphology, on the strain hardening behavior is not well understood. In this paper, we investigate this effect with in situ neutron diffraction during deformation of an Al-Cu alloy (206) after multiple aging treatments. Precipitate-dislocation interactions were found to change from precipitate shearing for microstructures predominantly containing GPI and θ" precipitates to Orowan looping for microstructures with primarily θ' and θ precipitates. Notably, significant anisotropy in strain hardening behavior was observed when θ' precipitates were present, which was attributed to crystallographic orientation dependent load transfer from the Al matrix to the θ' precipitates. The anisotropic load transfer is hypothesized to be caused by the extent of rotation of high aspect-ratio θ' precipitates, owing to dislocations looping around them during plastic deformation of the matrix. Predictions from an analytical model describing the anisotropic magnitude of load transfer from precipitate rotation agree well with experimental results, successfully validating the precipitate rotation hypothesis and explaining the anisotropic strain hardening behavior. This model allows for the prediction of stresses separately in the precipitate and matrix phases as a function of crystallographic orientation, only given the bulk mechanical properties.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Can classical DEM simultaneously capture compressibility and flowability of milled biomass?

Accurate prediction of the rheological behavior of biomass is essential for the design and operation of hoppers, feeders, and storage systems in biorefineries. This study examines whether the classical, coarse-grained discrete element method (DEM) formulation can simultaneously reproduce the compressibility and flowability of milled herbaceous biomass, using Miscanthus × giganteus as a representative material. The model represents particles as rigid spheres interacting through Hertz-Mindlin elastic-frictional contacts augmented with an area-dependent cohesion term. Laboratory cyclic compression and wedge-shaped hopper discharge experiments were used as calibration benchmarks. Although the model can independently reproduce each behavior by appropriately tuning particle Young's modulus E and cohesion energy density k, an extensive parametric investigation comprising more than 600 simulations reveals that the optimal parameter regions for compression and hopper flow are distinct and non-overlapping in (E, k) space. Surrogate surface analysis further shows that the corresponding objective-function valleys exhibit similar trends but are approximately parallel and spatially offset, precluding a unified calibration within the explored domain. Sensitivity analysis indicates that compressibility is governed predominantly by stiffness and cohesion, whereas the slope of the mass flow rate-opening relation in hopper discharge is primarily controlled by tangential friction. Extensions incorporating particle size distribution and clumped-sphere representations do not eliminate the incompatibility. These results systematically reveal, for the first time, the structural limitation of simplified DEM formulations in representing biomass rheological behavior, underscoring the necessity for models incorporating additional physical mechanisms, such as particle deformability or enhanced interlocking, to achieve unified predictive capability for biomass handling behavior.

09 BIOMASS FUELS↗
DOE OSTI2025

Probing Cretaceous-Paleogene crustal thickness in southern Tibet using quartz-zircon chronobarometry

Knowledge of the crustal thickness history of southern Tibet during the India-Asia collision is key to understanding what deformation mechanisms accommodated northward propagating crustal shortening. Thermoisotopic models suggest that a relatively thin margin of southern Asia persisted from ca. 200 to 45 Ma whereas trace element paleodepth proxies are interpreted to indicate a thinning phase from ca. 100 to 65 Ma when the crust reached ∼30 km-thickness. As neither of these methods directly measures crustal thickness, resolution of this conflict awaits development of a method that can. Here, in this study, we place bounds on the thickening history of southern Tibet using coupled quartz-zircon thermobarometry and U-Pb geochronology of granitoid plutons in the Gangdese batholith. We find that crustal thicknesses were as high (or higher) as 60–75 km at 65 Ma, or 10–15 Ma prior to the onset of hard continental collision. Magmatic inflation in the lower crust likely contributed to crustal thickening and thermal weakening of the lower crust, suggesting that pure shear in the lower crust was likely the primary accommodation mechanism for N-S shortening during collision. These data are in partial agreement with results of thermoisotopic models of crustal thickness but contradict thickness histories derived from empirical trace element proxies.

Chronobarometry↗
DOE OSTI2023

Microphysical Modeling of Fault Slip and Stability Transition in Hydrothermal Conditions

Abstract Field and laboratory observations indicate that the frictional behaviors of faults depend on hydrothermal conditions. We extend the microphysical Chen‐Niemeijer‐Spiers (CNS) model to hydrothermal conditions by using the observed temperature variation of indentation hardness to infer the temperature dependence of a microphysical parameter . This parameter is assumed constant in previous versions of the CNS model. A simple spring‐slider system is used to simulate the fault system and investigate the steady‐state frictional behaviors of wet granite gouges. Our numerical results quantitatively reproduce experimental data showing the frictional‐plastic transition. The results also describe the transition from velocity‐strengthening at low temperatures (<160°C), to velocity‐weakening at intermediate temperatures (160°C–370°C), then back to velocity‐strengthening at high temperatures (>370°C). In our extended CNS model, these results suggest that the dominant shear deformation mechanism does transition from frictional granular flow to fully plastic creep with increasing temperature.

58 GEOSCIENCES↗
DOE OSTI2025

Pomeranchuk instability from electronic correlations in CsTi 3 Bi 5 kagome metal

Electronic nematicity, the spontaneous breaking of rotational symmetry, has emerged as a key instability in correlated quantum systems. CsTi 3 Bi 5 , a kagome metal of the AV 3 Sb 5 (A = K, Rb, Cs) family, hosts rich unconventional electronic phases, yet the origin of its nematicity remains unsettled. Here, we combine polarization-dependent angle-resolved photoemission spectroscopy with functional renormalization group calculations on a fully interacting ab initio model. We reveal an orbital-selective nematic deformation in the low-energy band structure and identify a finite angular momentum (d-wave) Pomeranchuk instability driven by electronic correlations in specific orbital channels and detuning from Van Hove singularities. Our results establish a direct link between orbital selectivity and symmetry-breaking instabilities in CsTi 3 Bi 5 , providing a microscopic framework for nematic order in kagome systems.

Bigi, Chiara [Synchrotron SOLEIL, Saint-Aubin (Fra↗
DOE OSTI2024

Modeling of Sagging for 3D Printed Layers During the Curing Process

Here, a printed layer of silica-enforced poly(dimethylsiloxane)-co-(diphenylsiloxane) is modeled as a two-phase system consisting of air and polymer with an interface set up during the printing process. The structural geometry changes mostly due to the action of surface tension, while all material properties are strongly temperature dependent. Polymer flow is described using equations of the extended Herschel–Bulkley model, with parameters strongly dependent on temperature and degree of curing. Parameters of the model are determined using flow sweep measurements and separate experiments with vertical structure sagging at different temperatures. The curing process is modeled using dependencies between the curing rate, degree of curing, and temperature obtained in studies by differential scanning calorimetry. The developed model is used for simulations of printed structure deformations with different initial and boundary conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Shape evolution of neutron-rich 106,108,110 Mo isotopes in the triaxial degree of freedom

Background: Neutron-rich nuclei with mass number between 100 and 110 attract much attention, since several kinds of shapes, such as spherical, prolate, oblate, and triaxial shapes, are predicted. In particular, for neutron-rich Mo isotopes, different models predict different magnitudes and rigidity of triaxial deformation. Previous interpretations of experimental results based solely on low-lying $2^+_2$ states are insufficient to distinguish between the rigid triaxial shape, $\gamma$ vibration, or $\gamma$-soft rotor. Purpose: The shape evolution of 106 Mo, 108 Mo, and 110 Mo is investigated through their $2^+_1$-state lifetimes, decay-branching ratios of the $2^+_2$ state, and energies of the low-lying collective excited states with $K^π = 0^+, 2^+$, and $4^+$. Method: $\beta$-delayed $\gamma$-ray spectroscopy was employed for neutron-rich Nb and Zr isotopes produced at the RIKEN RI Beam Factory to populate excited states in 106 Mo, 108 Mo, and 110 Mo. The EUroball-RIKEN Cluster Array was used for high-resolution $\gamma$-ray detection and lifetimes of the $2^+_1$ states were determined using the UK fast-timing array of LaBr 3 (Ce) detectors. Results: New $\gamma$-ray transitions and levels are reported, including newly assigned $0^+_2$ states in 108,110 Mo. Quadrupole deformations were obtained for 106,108,110 Mo from their $2^+_1$ energies and lifetimes. The $\beta$-delayed neutron-emission probabilities of 108 Nb and 110 Nb were determined by examining the $\gamma$ rays of their respective daughter decays. Conclusions: In this work, the even-odd energy staggering in the $2^+_2$ band was compared with typical patterns of the $\gamma$-vibrational band, rigid triaxial rotor, and $\gamma$-soft rotor. The very small even-odd staggering of 106 Mo, 108 Mo, and 110 Mo favors a $\gamma$-vibrational band assignment. The kinematic moment of inertia for the $2^+_2$ band showed a trend similar to the ground-state band, which is as expected for the $\gamma$-vibrational band. Beyond-mean-field calculations employing the constrained Hartree-Fock-Bogoliubov and local quasiparticle-random-phase approximation method using the SLy5 + T interaction reproduced the ground and $2^+_2$ bands in 106 Mo and 108 Mo. The collective wave functions are consistent with the interpretation of the $2^+_2$ band as the $\gamma$-vibrational band of the prolate shape. However, the staggering pattern observed in 110 Mo differs from the one suggested in the calculations which predict a $\gamma$-soft rotor. There was no experimental indication of the oblate shape or the $\gamma$-soft rotor predicted in these Mo isotopes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Updated universal relations for tidal deformabilities of neutron stars from phenomenological equations of state

Equation of state (EOS) insensitive relations, so-called universal relations, between the neutron star (NS) compactness, its multipolar tidal deformability coefficients, and between the tidal parameters for binary systems are essential to break degeneracies in gravitational wave data analysis. Here, we validate and recalibrate these universal relations using a large set of almost 2 million phenomenological EOSs that are consistent with current observations. In doing so, we extend universal relations to a larger region of the EOS parameter space, most notably to softer EOSs and larger compactnesses. We show that waveform models that neglect higher-than-leading-order tidal deformations of the NSs accumulate as much as 3.5 radians of dephasing from 20Hz to merger. Further, we also perform a full Bayesian parameter estimation of the GW170817 data, and we compare the NS radius constraints produced using universal relations from the literature and the updated fits we propose here. We find that the new fits yield a NS radius that is smaller by about 500 meters. This difference is less than the statistical uncertainty on the radius at the signal-to-noise-ratio of GW170817, but it is significantly larger than the precision anticipated for next-generation detectors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

(2 + δ)-dimensional theory of the electromechanics of lipid membranes: Electrostatics

The coupling of electric fields to the mechanics of lipid membranes gives rise to intriguing electromechanical behavior, as, for example, evidenced by the deformation of lipid vesicles in external electric fields. Electromechanical effects are relevant for many biological processes, such as the propagation of action potentials in axons and the activation of mechanically gated ion channels. Currently, a theoretical framework describing the electromechanical behavior of arbitrarily curved and deforming lipid membranes does not exist. Purely mechanical models commonly treat lipid membranes as two-dimensional surfaces, ignoring their finite thickness. While holding analytical and numerical merit, this approach cannot describe the coupling of lipid membranes to electric fields and is thus unsuitable for electromechanical models. In a sequence of articles, we derive an effective surface theory of the electromechanics of lipid membranes, called the (2 + δ)-dimensional theory, which has the advantages of surface descriptions while accounting for finite thickness effects. The present article proposes a generic dimension reduction procedure relying on low-order spectral expansions. This procedure is applied to the electrostatics of lipid membranes to obtain the (2 + δ)-dimensional theory that captures potential differences across and electric fields within lipid membranes. Finally, this model is tested on different geometries relevant for lipid membranes, showing good agreement with the corresponding three-dimensional electrostatics theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗