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 1,513 records · Page 84

Systematics of E 2 strength in the s d shell with the valence-space in-medium similarity renormalization group

Recent developments in ab initio nuclear theory demonstrate promising results in medium- to heavy-mass nuclei. A particular challenge for many of the many-body methodologies, however, is an accurate treatment of the electric-quadrupole, E2, strength associated with collectivity. The valence-space in-medium similarity renormalization group (VS-IMSRG) is a particularly powerful method for accessing medium- and high-mass nuclei but has been found to underpredict E2 strengths. The purpose of this work is to evaluate the isospin dependence of this underprediction. We perform a systematic comparison of VS-IMSRG calculations with available literature. We make use of isoscalar and isovector contributions to the E2 matrix elements to assess isoscalar and isovector contributions to the missing strength. It is found that the E2 strength is consistent throughout T z =|12|, T z =|1|, T z =|32|, and T z =2 pairs within the sd shell. Furthermore, no isovector contribution to the deficiency is identified. A comparison with toy-models and coupled-cluster calculations is used to discuss potential origins of the missing strength, which arises from missing many-particle, many-hole excitations out of the model space. The absence of any significant isovector contribution to the missing E2 strength indicates that the E2 strength discrepancy, and therefore any correction, is largely independent of the isospin of the nuclei in question.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Deep anistropic shell program for tire analysis

A finite element program was constructed to model the mechanical response of a tire, treated as a deep anisotropic shell, to specified static loads. The program is based on a Sanders Budiansky type shell theory with the effects of transverse shear deformation and bending-extensional coupling included. A displacement formulation is used together with a total Lagrangian description of the deformation. Sixteen-node quadrilateral elements with bicubic shape functions are employed. The Noor basis reduction technique and various type of symmetry considerations serve to improve the computational efficiency.

Andersen, C. M.↗

Mirror nucleon removal reactions in p-shell nuclei

Nucleon removal reactions have been shown to be an effective tool for studying the single particle structure of nuclei. This work continues efforts to experimentally probe and benchmark the reaction and structure models used to calculate the removal reaction cross sections when using microscopic nuclear structure inputs. Here, three different single nucleon removal reactions were performed, from p -shell nuclei with masses A = 7, 9, and 10. The residual nuclei from the reactions were detected in coincidence with γ rays to determine partial cross sections to individual final states. The eikonal direct-reaction model is combined with overlap functions and residual nucleus densities from microscopic, variational Monte Carlo calculations to provide consistent nuclear structure input to the partial cross section calculations. Comparisons of measured and calculated cross sections, including for mirror reactions, are presented. The analysis of the partial cross sections leading to the ground states shows a similar behavior to the one observed from analyses of inclusive cross sections using shell model nuclear structure input: the theoretical description of the removal process is in better agreement with the data when removing weakly bound nucleons, than when removing well-bound ones. The two mirror reaction pairs presented here show consistent results between the respective members of the pairs. The results obtained for the population of the excited states, however, show a systematically different trend that appears connected to the structure part of the calculation. Additional cases are needed to better understand the respective roles of structure and dynamical effects in the deviations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The Rings Around the Egg Nebula

We present an eccentric binary model for the formation of the proto-planetary nebula CRL 2688 (the Egg Nebula) that exhibits multiple concentric shells. Given the apparent regularity of the structure in the Egg Nebula, we postulate that the shells are caused by the periodic passages of a companion star. Such an orbital period would have to lie in the range of 100-500 yr, the apparent time that corresponds to the spacing between the rings. We assume, in this model, that an asymptotic giant branch (AGB) star, which is the origin of the matter within the planetary nebula, loses mass in a spherically symmetric wind. We further suppose that the AGB star has an extended atmosphere (out to approximately 10 stellar radii) in which the outflow speed is less than the escape speed; still farther out, grains form and radiation pressure accelerates the grains along with the trapped gas to the escape speed. Once escape speed has been attained, the presence of a companion star will not significantly affect the trajectories of the matter leaving in the wind and the mass loss will be approximately spherically symmetric. On the other hand, if the companion star is sufficiently close that the Roche lobe of the AGB star moves inside the extended atmosphere, then the slowly moving material will be forced to flow approximately along the critical potential surface (i.e., the Roche lobe) until it flows into the potential lobe of the companion star. Therefore, in our model, the shells are caused by periodic cessations of the isotropic wind rather than by any periodic enhancement in the mass-loss process. We carry out detailed binary evolution calculations within the context of this scenario, taking into account the nuclear evolution and stellar wind losses of the giant as well as the effects of mass loss and mass transfer on the evolution of the eccentric binary orbit. From the initial binary parameters that we find are required to produce a multiple concentric shell nebula and the known properties of primordial binaries, we conclude that approximately 0.3% of all planetaries should go through a phase with multiple concentric shells.

Harpaz, Amos↗

Modeling and analysis of the space shuttle nose-gear tire with semianalytic finite elements

A computational procedure is presented for the geometrically nonlinear analysis of aircraft tires. The Space Shuttle Orbiter nose gear tire was modeled by using a two-dimensional laminated anisotropic shell theory with the effects of variation in material and geometric parameters included. The four key elements of the procedure are: (1) semianalytic finite elements in which the shell variables are represented by Fourier series in the circumferential direction and piecewise polynominals in the meridional direction; (2) a mixed formulation with the fundamental unknowns consisting of strain parameters, stress-resultant parameters, and generalized displacements; (3) multilevel operator splitting to effect successive simplifications, and to uncouple the equations associated with different Fourier harmonics; and (4) multilevel iterative procedures and reduction techniques to generate the response of the shell. Numerical results of the Space Shuttle Orbiter nose gear tire model are compared with experimental measurements of the tire subjected to inflation loading.

Kim, Kyun O.↗

Effect of longitudinal stringers on sound transmission into a thin cylindrical shell

In the context of the transmission of airborne noise into an aircraft fuselage, a mathematical model is presented for the transmission of airborne noise into a stiffened cylindrical shell. The stiffeners are longitudinal stringers and are modeled as discrete structural elements. The numerical cases examined were typical of a narrow-bodied jet transport fuselage. The stringers appeared to raise the cylinder transmission loss in the mass-controlled region, although they produced dips at the stringer resonances. The ring-frequency dip in transmission loss, which is characteristic of monocoque shells, was found to still be present. There appeared to be a small increase in transmission loss as the number of stringers was increased. The effect of high damping was also investigated.

Koval, L. R.↗

INTERNAL STRUCTURE OF ASTEROIDS HAVING SURFACE SHEDDING DUE TO ROTATIONAL INSTABILITY

Surface shedding of an asteroid is a failure mode where surface materials fly off due to strong centrifugal forces beyond the critical spin period, while the internal structure does not deform significantly. This paper proposes a possible structure of an asteroid interior that leads to surface shedding due to rapid rotation rates. A rubble pile asteroid is modeled as a spheroid composed of a surface shell and a concentric internal core, the entire assembly called the test body. The test body is assumed to be uniformly rotating around a constant rotation axis. We also assume that while the bulk density and the friction angle are constant, the cohesion of the surface shell is different from that of the internal core. First, developing an analytical model based on limit analysis, we provide the upper and lower bounds for the actual surface shedding condition. Second, we use a Soft-sphere Discrete Element Method (SSDEM) to study dynamical deformation of the test body due to a quasi-static spin-up. In this paper we show the consistency of both approaches. Additionally, the SSDEM simulations show that the initial failure always occurs locally and not globally. In addition, as the core becomes larger, the size of lofted components becomes smaller. These results imply that if there is a strong core in a progenitor body, surface shedding is the most likely failure mode.

79 ASTRONOMY AND ASTROPHYSICS↗

Effects of inclination angle on the spectra of X-ray binaries

Iron emission and absorption features from nine X-ray binaries are compared with a model that predicts characteristics of iron K-shell features corresponding to different viewing angles for the accretion disk. In this model, emission and absorption features arise from an accretion disk corona whose structure is computed under the assumptions of ionization, thermal, and hydrostatic balance of gas illuminated by the central continuum source. Hence the equivalent widths of the emission and absorption features are predicted without the introduction of any free parameters to control the physical conditions in the disk corona. Our data/model comparison suggests an explanation for the equivalent widths of the iron K-alpha emission and the depth of the absorption edge in terms of (1) the angle of inclination between the observer's line of sight and the plane of the accretion disk and (2) the shape of the incident source spectrum.

Vrtilek, S. D.↗

Finite Element Based Model of Fiber-optic Strain Sensor Mounted on a Semi-monocoque Structure

Fiber-optic strain sensing (FOSS) provides an opportunity for additional awareness and can enable other mitigations such as flutter suppression or gust load alleviation but requires an accurate computational model as part of a design process that integrates structures and controls. A modeling approach for fiber-optic strain sensors has been developed using a shell element formulation in a detailed finite element model. The new FOSS modeling approach is designed to better scale to more complex realistic structures and fiber layouts than previous approaches. The modeling method was applied to the X-56A vehicle and results were compared to experimental static ground-test data. Additionally, the experimental sensor noise of these fiber optic sensors was assessed for the correlations in time and correlations along the fiber. The noise does show significant correlations both in time and along the fibers. These correlations of the noise should also be considered in modeling and assessing algorithms that will use FOSS.

Jeffrey Ouellette↗

Finite Element Based Model of Fiber-optic Strain Sensor Mounted on a Semi-monocoque Structure

Fiber-optic strain sensing (FOSS) provides an opportunity for additional awareness and can enable other mitigations such as flutter suppression or gust load alleviation but requires an accurate computational model as part of a design process that integrates structures and controls. A modeling approach for fiber-optic strain sensors has been developed using a shell element formulation in a detailed finite element model. The new FOSS modeling approach is designed to better scale to more complex realistic structures and fiber layouts than previous approaches. The modeling method was applied to the X-56A vehicle and results were compared to experimental static ground-test data. Additionally, the experimental sensor noise of these fiber optic sensors was assessed for the correlations in time and correlations along the fiber. The noise does show significant correlations both in time and along the fibers. These correlations of the noise should also be considered in modeling and assessing algorithms that will use FOSS.

Jeffrey Ouellette↗

Improved Subcell Model for the Prediction of Braided Composite Response

In this work, the modeling of triaxially braided composites was explored through a semi-analytical discretization. Four unique subcells, each approximated by a "mosaic" stacking of unidirectional composite plies, were modeled through the use of layered-shell elements within the explicit finite element code LS-DYNA. Two subcell discretizations were investigated: a model explicitly capturing pure matrix regions, and a novel model which absorbed pure matrix pockets into neighboring tow plies. The in-plane stiffness properties of both models, computed using bottom-up micromechanics, correlated well to experimental data. The absorbed matrix model, however, was found to best capture out-of- plane flexural properties by comparing numerical simulations of the out-of-plane displacements from single-ply tension tests to experimental full field data. This strong correlation of out-of-plane characteristics supports the current modeling approach as a viable candidate for future work involving impact simulations.

Cater, Christopher R.↗

On the Origin of Fast Radio Bursts (FRBs)

We derive stringent constraints on the persistent source associated with FRB 121102: size 0.3<R{sub 17.5} =(R/10{sup 17.5} cm)<3, age <10{sup 2.5} year, energy E≈10{sup 49}(ε{sub e}/0.2GeV){sup 3} erg, characteristic electron energy 0.1⩽ε{sub e}/1 GeV⩽0.5; the radiating plasma is confined by a cold plasma of mass M{sub c}<10{sup −1.5}R{sub 17.5}{sup 4} M{sub ⊙}; these properties are inconsistent with typical “magnetar wind nebulae” model predictions. The fact that ε{sub e}∼m{sub p}c{sup 2} suggests that the hot plasma was created by the ejection of a mildly relativistic, M≈E/c{sup 2}≈10{sup −5} M{sub ⊙} shell, which propagated into an extended ambient medium or collided with a pre-ejected shell. Independent of the persistent source model, we suggest a physical mechanism for the generation of fast radio bursts (FRBs): the ejection from an underlying compact object, R{sub s}=10{sup 6}R{sub s,6} cm, of highly relativistic shells with energy E{sub s}=10{sup 41}E{sub 41} erg and Lorentz factor γ{sub s}=10{sup 3}E{sub 41}{sup 1/8}R{sub s,6}{sup −3/8}, into a surrounding e − p plasma with density n∼10{sup −1} cm{sup −3} (consistent with that inferred for the persistent source). For E {sub s} similar to observed FRB energies, plasma conditions appropriate for strong synchrotron maser emission at ν{sub coh.}≈0.5E{sub 41}{sup 1/4}R{sub s,6}{sup −3/4} GHz are formed. A significant fraction of the deposited energy is converted to an FRB with duration R{sub s}/c, accompanied by ∼10 MeV gamma-rays carrying less energy than the FRB. The inferred energy and mass associated with the source suggest some type of a “weak stellar explosion,” where a neutron star is formed with relatively low mass and energy ejection. However, the current upper limit on R does not allow one to rule out M{sub c}∼1 M{sub ⊙}, or the ejection of a larger mass well before the ejection of the confining shell.

79 ASTRONOMY AND ASTROPHYSICS↗

Prediction of Redox Potentials for U, Np, Pu, and Am in Aqueous Solution

The redox properties of the actinides in aqueous solution are important for fuel production/reprocessing and understanding the environmental impact of nuclear waste. The redox potentials for U, Np, Pu, and Am in oxidation states from 0 up to VII (as appropriate) in aqueous solutions have been predicted at the density functional theory level with the B3LYP functional, Stuttgart small core pseudopotential basis sets for the actinides, and explicit (30H 2 O molecules)/implicit treatment of the aqueous solvent using the self-consistent reaction field COSMO and SMD approaches for the implicit solvation. The predictions of the structural parameters of clusters incorporating first and second solvation shells are consistent with the available experimental data., Our results are typically within 0.2 V of the available experimental data using two explicit solvation shells with an implicit solvent model. The use of the PW91 functional substantially improved the prediction of the Pu(VI/V) redox couple. The redox couples for An(VI/IV) and An(V/IV) which involve the addition of protons and removal of the actinyl oxygens led to slightly larger differences from experiment. Here, the An(IV/0) and An(III/0) couples were reliably predicted with our approach. Predictions of the unknown An(II/I) redox potentials were negative, consistent with expectations, and predictions for unknown An(VII/VI), An(III/II), and An(II/0) redox couples improve prior estimates.

Actinides↗

Small Angle Neutron Scattering Study of Complex Coacervate Micelles and Hydrogels Formed from Ionic Diblock and Triblock Copolymers

A complex coacervate is a fluid phase that results from the electrostatic interactions between two oppositely charged macromolecules. The nature of the coacervate core structure of hydrogels and micelles formed from complexation between pairs of diblock or triblock copolymers containing oppositely charged end-blocks as a function of polymer and salt concentration was investigated. Both ABA triblock copolymers of poly[(allyl glycidyl ether)-b-(ethylene oxide)-b-(allyl glycidyl ether)] and analogous poly[(allyl glycidyl ether)-b-(ethylene oxide)] diblock copolymers, which were synthesized to be nearly one-half of the symmetrical triblock copolymers, were studied. The poly(allyl glycidyl ether) blocks were functionalized with either guanidinium or sulfonate groups via postpolymerization modification. Mixing of oppositely charged block copolymers resulted in the formation of nanometer-scale coacervate domains. Small angle neutron scattering (SANS) experiments were used to investigate the size and spacing of the coacervate domains. The SANS patterns were fit using a previously vetted, detailed model consisting of polydisperse core–shell micelles with a randomly distributed sphere or body-centered cubic (BCC) structure factor. For increasing polymer concentration, the size of the coacervate domains remained constant while the spatial extent of the poly(ethylene oxide) (PEO) corona decreased. However, increasing salt concentration resulted in a decrease in both the coacervate domain size and the corona size due to a combination of the electrostatic interactions being screened and the shrinkage of the neutral PEO blocks. Additionally, for the triblock copolymers that formed BCC ordered domains, the water content in the coacervate domains was calculated to increase from approximately 16.8% to 27.5% as the polymer concentration decreased from 20 to 15 wt %.

36 MATERIALS SCIENCE↗

Activation of ultrathin SrTiO 3 with subsurface SrRuO 3 for the oxygen evolution reaction

Here, water electrolysis occurs via the oxygen/hydrogen evolution reactions. Achieving sufficiently high oxygen evolution reaction (OER) activity while maintaining stability of the active catalyst surface is a primary challenge of designing OER catalysts. Often high electrocatalytic activity is accompanied by structural instability. This is the case for the metallic perovskite oxide SrRuO 3 (SRO), which exhibits rapid dissolution under OER conditions, both as thin-films and as nanoparticles. On the other hand, large band-gap perovskite oxides such as SrTiO 3 (STO) are inactive for OER in the dark but are stable. In this work, we demonstrate that burying as little as one unit cell of SRO beneath an ultrathin STO capping layer activates the otherwise inert electrocatalyst and gives excellent stability. Using density functional theory (DFT) calculations, we find that such a chemical modification of STO by sub-surface SRO introduces new 4d electronic states including Ru states within the STO band gap, raising the energy level of the electrons, changing the electronic hybridization, and facilitating an easy transfer to the adsorbed intermediates. We validate this hypothesis experimentally using atomically precise heteroepitaxial deposition. We find that a single-unit-cell layer of SRO is sufficient to activate the topmost STO layer towards OER; burying SRO underneath two unit cells of STO protects the inherently unstable SRO against corrosion during OER. Generally, our layered heterostructures are model systems of oxide core–shell structures, where an unstable catalyst (core) is protected against degradation and electronically activates the otherwise inactive shell. As demonstrated in this work, the growth of ultrathin heterostructures establishes a rigorous platform for screening core/shell nanoparticle combinations/thicknesses, and elucidates sub-surface activation mechanism for achieving stable and active oxide electrocatalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Comparative stereodynamics in molecule-atom and molecule-molecule rotational energy transfer: NO(A{sup 2}Σ{sup +}) + He and D{sub 2}

We present a crossed molecular beam scattering study, using velocity-map ion-imaging detection, of state-to-state rotational energy transfer for NO(A{sup 2}Σ{sup +}) in collisions with the kinematically identical colliders He and D{sub 2}. We report differential cross sections and angle-resolved rotational angular momentum polarization moments for transfer of NO(A, v = 0, N = 0, j = 0.5) to NO(A, v = 0, N′ = 3, 5-12) in collisions with He and D{sub 2} at respective average collision energies of 670 cm{sup −1} and 663 cm{sup −1}. Quantum scattering calculations on a literature ab initio potential energy surface for NO(A)-He [J. Kłos et al., J. Chem. Phys. 129, 244303 (2008)] yield near-quantitative agreement with the experimental differential scattering cross sections and good agreement with the rotational polarization moments. This confirms that the Kłos et al. potential is accurate within the experimental collisional energy range. Comparison of the experimental results for NO(A) + D{sub 2} and He collisions provides information on the hitherto unknown NO(A)-D{sub 2} potential energy surface. The similarities in the measured scattering dynamics of NO(A) imply that the general form of the NO(A)-D{sub 2} potential must be similar to that calculated for NO(A)-He. A consistent trend for the rotational rainbow maximum in the differential cross sections for NO(A) + D{sub 2} to peak at more forward angles than those for NO(A) + He is consistent with the NO(A)-D{sub 2} potential being more anisotropic with respect to NO(A) orientation. No evidence is found in the experimental measurements for coincident rotational excitation of the D{sub 2}, consistent with the potential having low anisotropy with respect to D{sub 2}. The NO(A) + He polarization moments deviate systematically from the predictions of a hard-shell, kinematic-apse scattering model, with larger deviations as N′ increases, which we attribute to the shallow gradient of the anisotropic repulsive NO(A)-He potential energy surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Measurement of the structure function of the nearly free neutron using spectator tagging in inelastic H 2 ( e , e ' p s ) X scattering with CLAS

In this study, much less is known about neutron structure than that of the proton due to the absence of free neutron targets. Neutron information is usually extracted from data on nuclear targets such as deuterium, requiring corrections for nuclear binding and nucleon off-shell effects. These corrections are model dependent and have significant uncertainties, especially for large values of the Bjorken scaling variable x. As a consequence, the same data can lead to different conclusions, for example, about the behavior of the d quark distribution in the proton at large x.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Phase transitions in core-collapse supernova matter at sub-saturation densities

Phase transitions in hot, dense matter in the collapsing cores of massive stars have an important impact on the core-collapse supernova mechanism as they absorb heat, disrupt homology, and so weaken the developing shock. We perform a three-dimensional, finite temperature Skyrme-Hartree-Fock (SHF) study of inhomogeneous nuclear matter to determine the critical density and temperature for the phase transition between the pasta phase and homogeneous matter and its properties. We employ four different parametrizations of the Skyrme nuclear energy-density functional, SkM*, SLy4, NRAPR, and SQMC700, which span a range of saturation-density symmetry energy behaviors constrained by a variety of nuclear experimental probes. For each of these interactions we calculate free energy, pressure, entropy, and chemical potentials in the range of particle number densities where the nuclear pasta phases are expected to exist, 0.02–0. 12 fm -3 , temperatures 2–8 MeV, and a proton fraction of 0.3. We find unambiguous evidence for a first-order phase transition to uniform matter, unsoftened by the presence of the pasta phases. No conclusive signs of a first-order phase transition between the pasta phases is observed, and it is argued that the thermodynamic quantities vary continuously right up to the first-order phase transition to uniform matter. We compare our results with thermodynamic spinodals calculated using the same Skyrme parametrizations, finding that the effect of short-range Coulomb correlations and quantum shell effects included in our model leads to the pasta phases existing at densities up to 0. 01 fm -3 above the spinodal boundaries, thus increasing the transition density to uniform matter by the same amount. The transition density is otherwise shown to be insensitive to the symmetry energy at saturation density within the range constrained by the concordance of a variety of experimental constraints, and can be taken to be a well determined quantity.

79 ASTRONOMY AND ASTROPHYSICS↗