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,657 records · Page 92

A comprehensive alpha-heating model for inertial confinement fusion

In this paper, a comprehensive model is developed to study alpha-heating in inertially confined plasmas. It describes the time evolution of a central low-density hot spot confined by a compressible shell, heated by fusion alphas, and cooled by radiation and thermal losses. The model includes the deceleration, stagnation, and burn phases of inertial confinement fusion implosions, and is valid for sub-ignited targets with ≤10× amplification of the fusion yield from alpha-heating. The results of radiation-hydrodynamic simulations are used to derive realistic initial conditions and dimensionless parameters for the model. It is found that most of the alpha energy (~90%) produced before bang time is deposited within the hot spot mass, while a small fraction (~10%) drives mass ablation off the inner shell surface and its energy is recycled back into the hot spot. Of the bremsstrahlung radiation emission, ~40% is deposited in the hot spot, ~40% is recycled back in the hot spot by ablation off the shell, and ~20% leaves the hot spot. We show here that the hot spot, shocked shell, and outer shell trajectories from this analytical model are in good agreement with simulations. Finally, a detailed discussion of the effect of alpha-heating on the hydrodynamics is also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electron temperature measurements inside the ablating plasma of gas-filled hohlraums at the National Ignition Facility

The first measurement of the electron temperature (T{sub e}) inside a National Ignition Facility hohlraum is obtained using temporally resolved K-shell X-ray spectroscopy of a mid-Z tracer dot. Both isoelectronic- and interstage-line ratios are used to calculate the local T{sub e} via the collisional–radiative atomic physics code SCRAM [Hansen et al., High Energy Density Phys 3, 109 (2007)]. The trajectory of the mid-Z dot as it is ablated from the capsule surface and moves toward the laser entrance hole (LEH) is measured using side-on x-ray imaging, characterizing the plasma flow of the ablating capsule. Data show that the measured dot location is farther away from the LEH in comparison to the radiation-hydrodynamics simulation prediction using HYDRA [Marinak et al., Phys. Plasmas 3, 2070 (1996)]. To account for this discrepancy, the predicted simulation T{sub e} is evaluated at the measured dot trajectory. The peak T{sub e}, measured to be 4.2 keV ± 0.2 keV, is ∼0.5 keV hotter than the simulation prediction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

User-Defined Material Model for Progressive Failure Analysis

An overview of different types of composite material system architectures and a brief review of progressive failure material modeling methods used for structural analysis including failure initiation and material degradation are presented. Different failure initiation criteria and material degradation models are described that define progressive failure formulations. These progressive failure formulations are implemented in a user-defined material model (or UMAT) for use with the ABAQUS/Standard1 nonlinear finite element analysis tool. The failure initiation criteria include the maximum stress criteria, maximum strain criteria, the Tsai-Wu failure polynomial, and the Hashin criteria. The material degradation model is based on the ply-discounting approach where the local material constitutive coefficients are degraded. Applications and extensions of the progressive failure analysis material model address two-dimensional plate and shell finite elements and three-dimensional solid finite elements. Implementation details and use of the UMAT subroutine are described in the present paper. Parametric studies for composite structures are discussed to illustrate the features of the progressive failure modeling methods that have been implemented.

Knight, Norman F. Jr.↗

The 3D inelastic analysis methods for hot section components

A two-year program to develop advanced 3D inelastic structural stress analysis methods and solution strategies for more accurate and cost effective analysis of combustors, turbine blades, and vanes is described. The approach was to develop a matrix of formulation elements and constitutive models. Three constitutive models were developed in conjunction with optimized iterating techniques, accelerators, and convergence criteria within a framework of dynamic time incrementing. Three formulation models were developed: an eight-noded midsurface shell element; a nine-noded midsurface shell element; and a twenty-noded isoparametric solid element. A separate computer program has been developed for each combination of constitutive model-formulation model. Each program provides a functional stand alone capability for performing cyclic nonlinear structural analysis. In addition, the analysis capabilities incorporated into each program can be abstracted in subroutine form for incorporation into other codes or to form new combinations.

Mcknight, R. L.↗

Theoretical Development of an Orthotropic Elasto-Plastic Generalized Composite Material Model

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites is becoming critical as these materials are gaining increased usage in the aerospace and automotive industries. While there are several composite material models currently available within LSDYNA (Livermore Software Technology Corporation), there are several features that have been identified that could improve the predictive capability of a composite model. To address these needs, a combined plasticity and damage model suitable for use with both solid and shell elements is being developed and is being implemented into LS-DYNA as MAT_213. A key feature of the improved material model is the use of tabulated stress-strain data in a variety of coordinate directions to fully define the stress-strain response of the material. To date, the model development efforts have focused on creating the plasticity portion of the model. The Tsai-Wu composite failure model has been generalized and extended to a strain-hardening based orthotropic yield function with a nonassociative flow rule. The coefficients of the yield function, and the stresses to be used in both the yield function and the flow rule, are computed based on the input stress-strain curves using the effective plastic strain as the tracking variable. The coefficients in the flow rule are computed based on the obtained stress-strain data. The developed material model is suitable for implementation within LS-DYNA for use in analyzing the nonlinear response of polymer composites.

Polymer Matrix Composites↗

Theoretical Development of an Orthotropic Elasto-Plastic Generalized Composite Material Model

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites is becoming critical as these materials are gaining increased usage in the aerospace and automotive industries. While there are several composite material models currently available within LS-DYNA (Registered), there are several features that have been identified that could improve the predictive capability of a composite model. To address these needs, a combined plasticity and damage model suitable for use with both solid and shell elements is being developed and is being implemented into LS-DYNA as MAT_213. A key feature of the improved material model is the use of tabulated stress-strain data in a variety of coordinate directions to fully define the stress-strain response of the material. To date, the model development efforts have focused on creating the plasticity portion of the model. The Tsai-Wu composite failure model has been generalized and extended to a strain-hardening based orthotropic material model with a non-associative flow rule. The coefficients of the yield function, and the stresses to be used in both the yield function and the flow rule, are computed based on the input stress-strain curves using the effective plastic strain as the tracking variable. The coefficients in the flow rule are computed based on the obtained stress-strain data. The developed material model is suitable for implementation within LS-DYNA for use in analyzing the nonlinear response of polymer composites.

Plasticity↗

Fossil Shell in 3C 84 as TeV γ-Ray Emitter and Cosmic-Ray Accelerator

We explore physical properties of the shocked external medium (i.e., a shell) in 3C 84 associated with the recurrent radio lobe born around 1960. In the previous work of Ito et al., we investigated a dynamical and radiative evolution of such a shell after the central engine stops the jet launching, and we found that a fossil-shell emission overwhelms that of the rapidly fading radio lobe. We apply this model to 3C 84 and find the following. (1) The fossil shell made of shocked diffuse ambient matter with the number density of 0.3 cm{sup −3} radiates bright inverse-Compton (IC) emission with the seed photons of the radio emission from the central compact region and the IC emission is above the sensitivity threshold of the Cherenkov Telescope Array (CTA). (2) When the fossil shell is produced in a geometrically thick ionized plasma with the number density of 10{sup 3} cm{sup −3}, the field strength in the shell may reach about 17 mG in the presence of magnetic fields amplification and the radio emission becomes comparable to the sensitivity of deep imaging VLBI observations. The possible production of ultra high-energy cosmic-rays (UHECRs) in the dense shocked plasma is also argued.

79 ASTRONOMY AND ASTROPHYSICS↗

Glass shell manufacturing in space

Highly-uniform, hollow glass spheres (shells), which are used for inertial confinement fusion targets, were formed from metal-organic gel powder feedstock in a vertical furnace. As a result of the rapid pyrolysis caused by the furnace, the gel is transformed to a shell in five distinct stages: (a) surface closure of the porous gel; (b) generation of a closed-cell foam structure in the gel; (c) spheridization of the gel and further expansion of the foam; (d) coalescence of the closed-cell foam to a single-void shell; and (e) fining of the glass shell. The heat transfer from the furnace to the falling gel particle was modeled to determine the effective heating rate of the gel. The model predicts the temperature history for a particle as a function of mass, dimensions, specific heat, and absorptance as well as furnace temperature profile and thermal conductivity of the furnace gas. A model was developed that predicts the gravity-induced degradation of shell concentricity in falling molten shells as a function of shell characteristics and time.

Downs, R. L.↗

Quantification of Processing Effects on Filament Wound Pressure Vessels

A computational simulation procedure is described which is designed specifically for the modeling and analysis of filament wound pressure vessels. Cylindrical vessels with spherical or elliptical end caps can be generated automatically. End caps other than spherical or elliptical may be modeled by varying circular sections along the x-axis according to the C C! end cap shape. The finite element model generated is composed of plate type quadrilateral shell elements on the entire vessel surface. This computational procedure can also be sued to generate grid, connectivity and material cards (bulk data) for component parts of a larger model. These bulk data are assigned to a user designated file for finite element structural/stress analysis of composite pressure vessels. The procedure accommodates filament would pressure vessels of all types of shells-of-revolution. It has provisions to readily evaluate initial stresses due to pretension in the winding filaments and residual stresses due to cure temperature.

Aiello, Robert A.↗

Quantification of Processing Effects on Filament Wound Pressure Vessels

A computational simulation procedure is described which is designed specifically for the modeling and analysis of filament wound pressure vessels. Cylindrical vessels with spherical or elliptical end caps can be generated automatically. End caps other than spherical or elliptical may be modeled by varying circular sections along the x-axis according to the end cap shape. The finite element model generated is composed of plate type quadrilateral shell elements on the entire vessel surface. This computational procedure can also be used to generate grid, connectivity and material cards (bulk data) for component parts of a larger model. These bulk data are assigned to a user designated file for finite element structural/stress analysis of composite pressure vessels. The procedure accommodates filament wound pressure vessels of all types of shells-of -revolution. It has provisions to readily evaluate initial stresses due to pretension in the winding filaments and residual stresses due to cure temperature.

Aiello, Robert A.↗

A Numerical Study of the Liquid Motion in Titan’s Subsurface Ocean

An ocean filled with liquid water lies beneath the icy surface of several Jovian and Saturnian moons. In such an ocean, the currents are driven by various phenomena such as the tidal forcing, the deformation of the ice shell lying at its top, the temperature gradient resulting from the surface and bottom heat fluxes…The flow induced by the first two forcings can be modelled by means of a 2D depth-averaged model, while the third one generates horizontal and vertical density variations whose effects can only be captured by a 3D baroclinic model. We study the tides of Titan’s subsurface ocean and the impact of the ice shell on the liquid motion by means of the Second-generation Louvain-la-Neuve Ice–ocean Model, SLIM (https://www.climate.be/slim). The impact of the ice shell lying at the top of the ocean is modelled by a surface friction term and surface pressure terms. The latter are a function of the difference between the ocean elevation and the vertical displacement of the shell and the time derivative of this difference. Because of Titan’s appreciable obliquity (0.306°), the tidal motion expected (and found) is similar to the Europa tidal scenario described by Tyler (2008): the surface elevation consists of two bulges rotating around Titan and the associated depth-averaged velocity field consists of two gyres, separated by an area of high speed flow, whose centre follows a sinusoidal path centred on the equator. The ice shell damps the surface motion, thus slowing down the flow, without significantly modifying the spatial patterns of these fields. The depth of the ocean and the mechanical characteristics of the ice shell being poorly constrained, a sensitivity analysis is conducted. The depth-averaged flow slows down when the depth is increased and a lag appears in the tidal phase but the tidal range remains similar. The ice shell mechanical characteristics influences both the elevation and depth-averaged velocity fields in terms of magnitude but does not modify the spatial patterns of these fields. The influence of the surface heat flux is studied by means of the 3D baroclinic version of SLIM. The heat flux derived from Titan’s topography by Kvorka et al. (2018) is used as surface boundary condition for the temperature equation while a uniform bottom heat flux is implemented. Its value is computed assuming that the heat budget of the ocean is at equilibrium. These boundary conditions cause density variations, which impact the hydrodynamics of the ocean. While the flow velocity induced by these variations is two orders of magnitude smaller than the tidal flow, its orientation is time-independent, hence impacting the orientation of the velocity field. Although the variations of ocean surface elevation and speed with respect to the shell mechanical properties can be larger than those induced by the surface heat flux, taking into account the latter results in large variations of the velocity field global patterns, which was not observed when modifying the shell mechanical properties. Future studies should therefore focus on modelling the surface and bottom heat fluxes while uncertainties about the mechanical characteristics of the shell can be tolerated.

Titan↗

Plates and shells containing a surface crack under general loading conditions

Various through and part-through crack problems in plates and shells are considered. The line-spring model of Rice and Levy is generalized to the skew-symmetric case to solve surface crack problems involving mixed-mode, coplanar crack growth. Compliance functions are introduced which are valid for crack depth to thickness ratios at least up to .95. This includes expressions for tension and bending as well as expressions for in-plane shear, out-of-plane shear, and twisting. Transverse shear deformation is taken into account in the plate and shell theories and this effect is shown to be important in comparing stress intensity factors obtained from the plate theory with three-dimensional solutions. Stress intensity factors for cylinders obtained by the line-spring model also compare well with three-dimensional solution. By using the line-spring approach, stress intensity factors can be obtained for the through crack and for part-through crack of any crack front shape, without recalculation integrals that take up the bulk of the computer time. Therefore, parameter studies involving crack length, crack depth, shell type, and shell curvature are made in some detail. The results will be useful in brittle fracture and in fatigue crack propagation studies. All problems considered are of the mixed boundary value type and are reducted to strongly singular integral equations which make use of the finite-part integrals of Hadamard. The equations are solved numerically in a manner that is very efficient.

Joseph, Paul F.↗

Copper fine-structure K-shell electron impact ionization cross sections for fast-electron diagnostic in laser-solid experiments

The K-shell electron impact ionization (EII) cross section, along with the K-shell fluorescence yield, is one of the key atomic parameters for fast-electron diagnostic in laser-solid experiments through the K-shell emission cross section. In addition, copper is a material that has been often used in those experiments because it has a maximum total K-shell emission yield. Furthermore, in a campaign dedicated to the modeling of the K lines of astrophysical interest (Palmeri et al., 2012), the K-shell fluorescence yields for the K-vacancy fine-structure atomic levels of all the copper isonuclear ions have been calculated. In this study, the K-shell EII cross sections connecting the ground and the metastable levels of the parent copper ions to the daughter ions K-vacancy levels considered in Palmeri et al. (2012) have been determined. The relativistic distorted-wave (DW) approximation implemented in the FAC atomic code has been used for the incident electron kinetic energies up to 10 times the K-shell threshold energies. Moreover, the resulting DW cross sections have been extrapolated at higher energies using the asymptotic form proposed by Davies et al. (2013)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Spatiotemporal Formation and Growth Kinetics of Polyelectrolyte Complex Micelles with Millisecond Resolution

We have directly observed the in situ self-assembly kinetics of polyelectrolyte complex (PEC) micelles by synchrotron time-resolved small-angle X-ray scattering, equipped with a stopped-flow device that provides millisecond temporal resolution. A synthesized neutral-charged diblock polycation and homopolyanion that we have previously investigated as a model charge-matched, core-shell micelle system were selected for this work. The initial micellization of the oppositely charged polyelectrolytes was completed within the dead time of mixing of 100 ms, followed by micelle growth and equilibration up to several seconds. By combining the structural evolution of the radius of gyration (R-g) with complementary molecular dynamics simulations, we show how the self-assemblies evolve incrementally in size over time through a two-step kinetic process: first, oppositely charged polyelectrolyte chains pair to form nascent aggregates that immediately assemble into spherical micelles, and second, these PEC micelles grow into larger micellar entities. Furthermore, this work has determined one possible kinetic pathway for the initial formation of PEC micelles, which provides useful physical insights for increasing fundamental understanding self-assembly dynamics, driven by polyelectrolyte complexation that occurs on ultrafast time scales.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A low-cost approach to electronic excitation energies based on the driven similarity renormalization group

Here, we propose an economical state-specific approach to evaluate electronic excitation energies based on the driven similarity renormalization group truncated to second order (DSRG-PT2). Starting from a closed-shell Hartree–Fock wave function, a model space is constructed that includes all single or single and double excitations within a given set of active orbitals. The resulting VCIS-DSRG-PT2 and VCISD-DSRG-PT2 methods are introduced and benchmarked on a set of 28 organic molecules [M. Schreiber et al. , J. Chem. Phys. 128, 134110 (2008)]. Taking CC3 results as reference values, mean absolute deviations of 0.32 and 0.22 eV are observed for VCIS-DSRG-PT2 and VCISD-DSRG-PT2 excitation energies, respectively. Overall, VCIS-DSRG-PT2 yields results with accuracy comparable to those from time-dependent density functional theory using the B3LYP functional, while VCISD-DSRG-PT2 gives excitation energies comparable to those from equation-of-motion coupled cluster with singles and doubles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hyperspectral interference tomography of nacre

Significance We invented an optical technique—hyperspectral interference tomography—that rapidly and nondestructively extracts nanoscale structural information across large samples of nacre (mother-of-pearl) and other layered materials by combining multiangle and polarization-resolved hyperspectral imaging with optical-interference modeling. We investigated nacre in mollusk shells from two different species, red abalone and rainbow abalone, and discovered a previously unknown relationship between the age of the mollusk and the thickness of aragonite tablets in nacre. Hyperspectral interference tomography will have applications in climate science, since nacre tablet thickness in fossil shells is a proxy for ancient seawater temperature, and in bioinspired mechanics, because the layered structure of nacre inspires engineered materials with exceptional strength and toughness.

47 OTHER INSTRUMENTATION↗

Multi-Fidelity Bayesian Optimization with Gaussian Processes for Double Shell Inertial Confinement Fusion Target Design

Reliable, secure access to energy is a major focus for national security efforts. One potential route to such energy is through fusion reactions in inertial confinement fusion (ICF) experiments. Such experiments are carried out at facilities such as the National Ignition Facility (NIF) in Livermore, California, where high powered lasers are used to compress a DT fuel-containing target to the necessary high temperature, high pressure conditions. These experiments are limited in number, which creates a heavy dependence on high fidelity predictive physics simulations and analysis performed “pre shot,” or before the experiment occurs. Many of these simulations in higher dimensions (2D and 3D) are computationally expensive, so finding optimal simulation-based designs presents its own challenges. In this work, we present our multi-fidelity Bayesian optimization with Gaussian processes (GPs) for ICF double shell targets, where a 1D surrogate model is used to help find a 2D surrogate model, enabling us to find optimal targets in the higher fidelity (2D), while saving computational cost.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗