Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SHELL MODELS”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Nucleon-pair approximation for nuclei from spherical to deformed regions

Here we model low-lying states of atomic nuclei in the nucleon-pair approximation of the shell model, using three approaches to select collective nucleon pairs: the generalized seniority scheme, the conjugate gradient method, and the Hartree-Fock approach. We find the collective pairs obtained from the generalized seniority scheme provide a good description for nearly spherical nuclei, and those from the conjugate gradient method or the Hartree-Fock approach work well for transitional and deformed nuclei. Our NPA calculations using collective pairs with angular momenta 0, 2, and 4 (denoted by S D G pairs) reproduce the nuclear shape evolution in the N = 26 isotones, Ca 46 , Ti 48 , Cr 50 , and Fe 52 , and yield good agreement with full configuration-interaction calculations of low-lying states in medium-heavy transitional and deformed nuclei: Ti 44 – 48 , Cr 48 , Cr 50 , Fe 52 , Zn 60 – 64 , Ge 64 , 66 , Mo 84 , and Xe 108 – 112 . Finally, using the S D G I -pair approximation we describe low-lying states of Ba 112 , 114 , cases difficult to reach by conventional configuration-interaction methods.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Union of rotational and vibrational modes in generator-coordinate-type calculations, with application to neutrinoless double-β decay

Good many-body methods for medium and heavy nuclei are important. Here we combine rotational motion with vibrational modes in a novel generator-coordinate method (GCM): starting from a mean-field solution (Hartree-Fock-Bogoliubov), we create nonorthogonal reference states from low-lying quasiparticle Tamm-Dancoff modes, and then project onto states of good angular momentum and particle number. Here, the results we benchmark against full shell model calculations. Even with just a few such modes we find improvement over standard GCM calculations in excitation spectra. We also find significant improvement in 0νββ nuclear matrix elements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Persistence of collectivity in the low-lying states of 30,31 Na inside the 𝑁 = 20 island of inversion

Near the 𝛽-stability line, nuclei with magic numbers are expected to exhibit spherical ground states. However, mass measurements of neutron-rich nuclei near the 𝑁 = 20 island of inversion have indicated an excess in the binding energies. These results, together with subsequent spectroscopic studies and theoretical works, point to the structural evolution from normal configurations to intruder-deformed configurations characterized by particle-hole excitations across the shell gap. Despite evidence of intruder-dominant ground states for the 30,31 Na isotopes in the vicinity of 𝑁 = 20, a question remains as to whether the deformation persists into the ground-state bands beyond the first excited states. Here, in this work, we report on a heavy-ion inelastic scattering measurement performed with GRETINA, the TRIPLEX device, and the S800 spectrograph to study the low-lying excitations in neutron-rich 30,31 Na . The observed gamma rays for the (4 + ) → (3 + ), (3 + ) → 2$^+_{g.s.}$, and (4 + ) → 2$^+_{g.s.}$ decays in 30 Na and for the (7/2 + ) → (5/2 + ) and (5/2 + )→3/2$^+_{g.s.}$ decays in 31 Na were used to extract the reduced transition probabilities 𝐵⁡(𝐸⁢2↑). We report the first measurement of the 𝐵⁡(𝐸⁢2↑) between the 3/2$^+_{g.s.}$ and the (7/2 + ) states in 31 Na . The results from this study are compared to shell-model calculations, confirming the persistence of large collectivity in the low-lying excited states consistent with the formation of the well-deformed ground-state bands in 30,31 Na .

collective models↗
DOE OSTI2020

Shell Element Material Model Verification Problems for DYNA3D: Part II

A suite of supplementary shell element material model verification tests were developed for the explicit finite element program DYNA3D, in continuation of the work performed in LLNL-TR-792469 (Shell Element Material Model Verification Problems for DYNA3D). The testing procedure developed in the preceding report is extended, and used to verify all remaining untested shell model features and inputs. A collection of 78 feature-specific verification tests are proposed, some of which are applicable to multiple material models, though many are specialized to a particular model. The collective suite of tests cover all 27 currently available material models in DYNA3D. In the course of developing the proposed test suite, 24 separate shell model-related issues were identified and have subsequently been resolved.

42 ENGINEERING↗
DOE OSTI2025

Puzzling 𝐵⁡(𝐸⁢2;0 + → 2 + ) strength in the proton dripline nucleus 36 Ca

Recent measurements of the 𝐸⁢2 transition rate from the ground state to the first 2 + excited state of the proton dripline nucleus 36 Ca show an unusual pattern when compared to its isotopic neighbor 38 Ca: despite having a higher 𝐸 𝑥 ⁡(2$^{+}_{1}$) excitation energy, the 𝐵⁡(𝐸⁢2;0$^{+}_{1}$ → 2$^{+}_{1}$) rate in 36 Ca is larger. The question that naturally arises is to what extent this observation can be attributed to the unbound character of the 2$^{+}_{1}$ state. To understand the influence of the continuum space on the low-energy properties of 36 Ca, we carried out Gamow shell model calculations that can account for the continuum coupling effects associated with the occupation of unbound 𝑓⁡𝑝 shells. We found that in the threshold 2 + state, 36 Ca is spatially diffused, which impacts the observed 𝐵⁡(𝐸⁢2) trend.

20 ≤ A ≤ 38↗
DOE OSTI2022

α-cluster structure of 18 Ne

In this work, we investigate α-clustering in 18 Ne and compare it with what is known about clustering in the mirror nucleus 18 O. The excitation function for resonant elastic scattering of 14 O + α was measured in inverse kinematics using the TexAT active-target detector system, covering excitation energies from 8 to 17 MeV. The data were analyzed using a multi-channel R-matrix approach. Detailed spectroscopic information—including excitation energies, spin and parity assignments, and partial α and total widths—was extracted from the R-matrix analysis and compared with theoretical models and previously available data. Correspondence between levels in 18 O and 18 Ne is established. An extensive shell-model analysis of the 18 O– 18 Ne mirror pair shows very good agreement between theory and experiment, especially in the distribution of clustering strength. The comparison reveals that certain states, particularly at higher excitation energies, exhibit significantly stronger α-cluster character than predicted. This suggests that these states have a collective structure aligned with the corresponding α-reaction channel.

6 ≤ A ≤ 19↗
DOE OSTI2017

SS-HORSE Method for Analysis of Resonances: Charged-Particle Scattering

The method for analyzing resonance states based on the harmoni©oscillator representation of scattering equations (HORSE) and on analytic properties of partial-wave scattering amplitudes which has been proposed earlier, is generalized to the case of charged-particle scattering. The method in question is tested by applying it to the model problem of pα scattering and can be used to study resonance states on the basis of microscopic calculations within various versions of the shell model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Quantum magic and multipartite entanglement in the structure of nuclei

Motivated by the Gottesman-Knill theorem, we present a detailed study of the quantum complexity of p -shell and s d -shell nuclei. Valence-space nuclear shell-model wave functions generated by the code are mapped to qubit registers using the Jordan-Wigner mapping (12 qubits for the p shell and 24 qubits for the s d shell), from which measures of the many-body entanglement ( n -tangles) and magic (nonstabilizerness) are determined. While exact evaluations of these measures are possible for nuclei with a modest number of active nucleons, Monte Carlo simulations are required for the more complex nuclei. The broadly applicable Pauli-string I ̂ Z ̂ exact (PSIZe) Markov chain Monte Carlo (MCMC) technique is introduced to accelerate the evaluation of measures of magic in deformed nuclei (with hierarchical wave functions), by factors of ≈ 8 for some nuclei. Significant multinucleon entanglement is found in the s d shell, dominated by proton-neutron configurations, along with significant measures of magic. This is evident not only for the deformed states, but also for nuclei on the path to instability via regions of shape coexistence and level inversion. These results indicate that quantum-computing resources will accelerate precision simulations of such nuclei and beyond. Published by the American Physical Society 2025

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

SS-HORSE method for studying resonances

A new method for analyzing resonance states based on the Harmonic-Oscillator Representation of Scattering Equations (HORSE) formalism and analytic properties of partial-wave scattering amplitudes is proposed. The method is tested by applying it to the model problem of neutral-particle scattering and can be used to study resonance states on the basis of microscopic calculations performed within various versions of the shell model.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1978

Finite element dynamic analysis on CDC STAR-100 computer

Computational algorithms are presented for the finite element dynamic analysis of structures on the CDC STAR-100 computer. The spatial behavior is described using higher-order finite elements. The temporal behavior is approximated by using either the central difference explicit scheme or Newmark's implicit scheme. In each case the analysis is broken up into a number of basic macro-operations. Discussion is focused on the organization of the computation and the mode of storage of different arrays to take advantage of the STAR pipeline capability. The potential of the proposed algorithms is discussed and CPU times are given for performing the different macro-operations for a shell modeled by higher order composite shallow shell elements having 80 degrees of freedom.

Noor, A. K.↗
NASA NTRS2012

Application of Out-of-Plane Warping to Control Rotor Blade Twist

The goal of this ongoing study is to develop and demonstrate the feasibility of a blade actuation system to dynamically change the twist, and/or the camber, of an airfoil section and, consequently, alter the in-flight aerodynamic loading on the blade for efficient flight control. The required analytical and finite element tools are under development to enable an accurate and comprehensive aeroelastic assessment of the current Full-Blade Warping and 3D Warping Actuated Trailing Edge Flap concepts. The feasibility of the current concepts for swashplateless rotors and higher harmonic blade control is also being investigated. In particular, the aim is to complete the following objectives, some of which have been completed (as noted below) and others that are currently ongoing: i) Develop a Vlasov finite element model and validate against the ABAQUS shell models (completed). ii) Implement the 3D warping actuation concept within the comprehensive analysis code DYMORE. iii) Perform preliminary aeroelastic simulations of blades using DYMORE with 3D warping actuation: a) Investigate the blade behavior under 1 per/rev actuation. Determine whether sufficient twist can be generated and sustained to achieve primary blade control. b) Investigate the behavior of a trailing edge flap configuration under higher harmonic excitations. Determine how much twist can be obtained at the harmonics 2-5 per/rev. iv) Determine actuator specifications such as the power required, load and displacements, and identify the stress and strain distributions in the actuated blades. In general, the completion of Item ii) above will give an additional research capability in rotorcraft dynamics analyses, i.e., the capability to calculate the rotor blade twist due to warping, something that is not currently available in any of the existing comprehensive rotorcraft analyses.

VanWeddingen, Yannick↗
DOE OSTI2017

Development of a Full Ice-cream Cone Model for Halo Coronal Mass Ejections

It is essential to determine three-dimensional parameters (e.g., radial speed, angular width, and source location) of coronal mass ejections (CMEs) for the space weather forecast. In this study, we investigate which cone type represents a halo CME morphology using 29 CMEs (12 Solar and Heliospheric Observatory (SOHO) /Large Angle and Spectrometric Coronagraph (LASCO) halo CMEs and 17 Solar Terrestrial Relations Observatory ( STEREO )/Sun–Earth Connection Coronal and Heliospheric Investigation COR2 halo CMEs) from 2010 December to 2011 June. These CMEs are identified as halo CMEs by one spacecraft ( SOHO or one of STEREO A and B ) and limb ones by the other spacecraft (One of STEREO A and B or SOHO ). From cone shape parameters of these CMEs, such as their front curvature, we find that the CME observational structures are much closer to a full ice-cream cone type than a shallow ice-cream cone type. Thus, we develop a full ice-cream cone model based on a new methodology that the full ice-cream cone consists of many flat cones with different heights and angular widths to estimate the three-dimensional parameters of the halo CMEs. This model is constructed by carrying out the following steps: (1) construct a cone for a given height and angular width, (2) project the cone onto the sky plane, (3) select points comprising the outer boundary, and (4) minimize the difference between the estimated projection speeds with the observed ones. By applying this model to 12 SOHO /LASCO halo CMEs, we find that 3D parameters from our method are similar to those from other stereoscopic methods (i.e., a triangulation method and a Graduated Cylindrical Shell model).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Proton and neutron contributions to the quadrupole transition strengths in 39 Ca and 39 K studied by lifetime measurements of mirror transitions

The E2 transition matrix elements of isobaric multiplets are expected to follow a linear trend as a function of isospin projection. However, measurements of the 2 + → 0 + transitions in the A = 38 triplet of Ca, K, and Ar show a deviation from this trend with an enhanced transition strength in 38 Ca with respect to its mirror 38 Ar. We have studied analogue 11/2 – → 7/2 – E2 transitions in 39 Ca and its mirror partner 39 K to determine if this enhancement persists in neighboring Ca isotopes. Recoil-distance lifetime measurements of 39 Ca and 39 K were performed utilizing a 42 Sc secondary beam, the TRIPLEX plunger, the GRETINA array, and the S800 spectrograph. Our data provide a lifetime measurement of the (11/2 – ) state in 39 Ca as well as an improved lifetime result for the (9/2 – ) state, while the 39 K data are used to validate the present analysis. Furthermore, a comparison of the present data to shell-model calculations suggests an enhanced transition strength in 39 Ca, pointing to both proton and neutron contributions to core excitations across the Z = N = 20 shell gaps in close proximity to 40 Ca.

39 ≤ A ≤ 58↗
DOE OSTI2018

Modeling and simulation of CO 2 capture using semipermeable elastic microcapsules

Very high rates of CO2 absorption have been demonstrated using semipermeable microcapsules filled with environmentally benign sodium carbonate solutions (Vericella et al., Nature communications 6, 2015). However, such capsules have also exhibited several traits, including water loss/uptake, elastic swelling and buckling, and precipitation of solids, which make them more complicated to engineer as a carbon capture material. To address this, a mechanistic model for mass transfer and chemical reaction has been developed, which accounts for these behaviors, and can be used as a predictive tool to explore capsule performance under different CO2 capture scenarios. The model uses a modified film theory to describe CO2 and H2O mass transfer rates, a concentration based description of vapor-liquid-solid equilibrium inside the capsule, and an elastic description of the capsule shell. Model predictions compare favorably to literature data for precipitating carbonates, as well as recent data for microcapsule CO2 absorption. The model was used to assess practical issues of capsule water transfer during a temperature swing CO2 absorption/desorption cycle. Model predictions indicate that capsules filled with high weight percent, precipitating sodium carbonate solution could be significantly or completely dehydrated during CO2 absorption. Subsequent rehydration of lean capsules at elevated temperature showed strong sensitivity to both gas humidity, and capsule shell elasticity.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1988

Part 1: Aspects of lithospheric evolution on Venus. Part 2: Thermal and collisional histories of chondrite parent bodies

The geological evolution of distinctly different kinds of solar system objects is addressed. Venus has been observed over the past decade by orbital radars on both American and Soviet spacecraft. These surface measurements provide clues to the structure and evolution of the lithosphere. The parent bodies of chondritic meteorites, thought to resemble asteroids, represent the other end of the size spectrum of terrestrial objects. Their early thermal and collisional histories may be constrained by the chemical and textural record preserved in meteorite samples. Impact craters on Venus have been observed by the Soviet Venera 15/16 spacecraft. A formalism is presented by which the size-frequency distribution of impact craters may be used to estimate upper bounds on the mean global rates of volcanic resurfacing and lithospheric recycling on that planet over the past several hundred million years. The impact crater density reported from Venera observations, if valid for the entire Venus surface, indicates a mean volcanic flux no greater than 2 cu km/y, corresponding to a maximum average rate of resurfacing of about 4 km/b.y. For the lowest estimated mean crater retention age of the surface of Venus imaged by Venera 15/16, the rate of lithospheric recycling on Venus does not exceed 1.5 sq km/y. Ordinary chondrite meteorites show textural and chemical patterns indicative of varying intensities of thermal metamorphism. The conventional onion-shell model, which envisions highly metamorphosed material in the core and less intensely heated rocks near the surface, predicts an inverse relation between peak temperature and cooking rate, but none has been observed. A metamorphosed-planetesimal model is devised to explain this discrepancy, whereby heating occurs in planetesimals a few kilometers in radius which then accrete to form 100-km-radius parent bodies. Cooling rates are then randomly controlled by burial depth. Thermal and collisional constraints are examined, and the model is found to be applicable only to highly insulating Al-26-rich planetesimals that remain closely aggregated upon accretion. An alternative model is presented, in which onion-shell parent bodies are collisionally fragmented during metamorphism and then gravitationally reassembled. If reassembly times are short, then cooling rates would be determined by burial depth in the reaccreted parent body. This model, unlike previous ones, can explain both coherent and incoherent cooling of Breccia clasts by collisions during or after metamorphism, respectively.

Grimm, Robert E.↗
DOE OSTI2023

Projection algorithm for state preparation on quantum computers

We present an efficient method to prepare states of a many-body system on quantum hardware, first isolating individual quantum numbers and then using time evolution to isolate the energy. Our method in its simplest form requires only one additional auxiliary qubit. The total time evolved for an accurate solution is proportional to the ratio of the spectrum range of the trial state to the gap to the lowest excited state, a substantial improvement over other projection algorithms, and the accuracy increases exponentially with the time evolved. Isolating the quantum numbers is efficient because of the known eigenvalues and increases the gap thus shortening the propagation time required. The success rate of the algorithm, or the probability of producing the desired state, is a simple function of measurement times and phases and is dominated by the square overlap of the original state to the desired state. We present examples from the nuclear shell model and the Heisenberg model. In conclusion, we compare this algorithm to previous algorithms for short evolution times and discuss potential further improvements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Evolution of low-lying M1 modes in germanium isotopes

Magnetic dipole strength functions are determined for the series of germanium isotopes from N=Z=32 to N = 48 on the basis of a large number of transition strengths calculated within the shell model. The evolution of the strength with increasing neutron number in the 1g 9/2 orbital is analyzed. A bimodal structure comprising an enhancement toward low transition energy and a resonance in the region of the scissors mode is identified. The low-energy enhancement is strongest near closed shells, in particular at the almost completely filled 1g 9/2 orbital, while the scissorslike resonance is most pronounced in the middle of the open shell, which correlates with the magnitude of the also deduced electric quadrupole transition strengths. The results are consistent with previous findings for the shorter series of iron isotopes and prove the occurrence and correlation of the two low-lying magnetic dipole modes as a global structural feature.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗