Search NASA⌕ Search

SEARCH · Search NASA

Results for “shell model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17

Radiative transfer in circumstellar dust shells

An iterative numerical procedure is developed to calculate the radiation field as a function of position, direction, and wavelength within a spherically symmetric circumstellar dust shell. The dust shell is assumed to consist of gray, isotropically scattering dust particles in thermal equilibrium with the radiation field and to be characterized by seven parameters: radius of the central star, inner radius of the shell, outer radius of the shell, total optical depth of the shell, an index which specifies the density distribution, albedo of the dust particles, and temperature of the central star. When the shell is close to the star, the angular distribution of the stellar component of the radiation field is determined by direct integration along lines of sight to the stellar surface, and the back-warming effect of the shell is included in the calculations to determine the net stellar flux. Using this procedure, the temperature distributions and radiation fields within several model dust shells are determined and used to calculate for each shell the spectral-energy and spatial-energy distributions of the radiation emitted at 11 wavelengths from 0.4 to 20 microns. For shells of moderate optical depth, it is found that the radiation field is quite anisotropic, especially at shorter wavelengths, and that the shape of the spectral-energy distribution is strongly dependent on all seven model parameters. The Eddington factor as a function of wavelength and optical depth is calculated for each model.

Harvel, C. A.↗

Skin-Stiffener Debond Prediction Based on Computational Fracture Analysis

Interlaminar fracture mechanics has proven useful for characterizing the onset of delaminations in composites and has been used with limited success primarily to investigate onset in fracture toughness specimens and laboratory size coupon type specimens. Future acceptance of the methodology by industry and certification authorities however, requires the successful demonstration of the methodology on structural level. For this purpose a panel was selected that is reinforced with stringers. Shear loading causes the panel to buckle and the resulting out-of-plane deformations initiate skin/stringer separation at the location of an embedded defect. For finite element analysis, the panel and surrounding load fixture were modeled with shell elements. A small section of the stringer foot and the panel in the vicinity of the embedded defect were modeled with a local 3D solid model. Across the width of the stringer foot the mixed-mode strain energy release rates were calculated using the virtual crack closure technique. A failure index was calculated by correlating the results with the mixed-mode failure criterion of the graphite/epoxy material. For small applied loads the failure index is well below one across the entire width. With increasing load the failure index approaches one first near the edge of the stringer foot from which delamination is expected to grow. With increasing delamination lengths the buckling pattern of the panel changes and the failure index increases which suggests that rapid delamination growth from the initial defect is to be expected.

Krueger, Ronald↗

Analysis of Composite Panel-Stiffener Debonding Using a Shell/3D Modeling Technique

Interlaminar fracture mechanics has proven useful for characterizing the onset of delaminations in composites and has been used with limited success primarily to investigate onset in fracture toughness specimens and laboratory size coupon type specimens. Future acceptance of the methodology by industry and certification authorities however, requires the successful demonstration of the methodology on structural level. For this purpose a panel was selected that was reinforced with stringers. Shear loading cases the panel to buckle and the resulting out-of-plane deformations initiate skin/stringer separation at the location of an embedded defect. For finite element analysis, the panel and surrounding load fixture were modeled with shell element. A small section of the stringer foot and the panel in the vicinity of the embedded defect were modeled with a local 3D solid model. A failure index was calculated by correlating computed mixed-mode failure criterion of the graphite/epoxy material.

Krueger, Ronald↗

Delayed detonation models for normal and subluminous type Ia sueprnovae: Absolute brightness, light curves, and molecule formation

We compute optical and infrared light curves of the pulsating class of delayed detonation models for Type Ia supernovae (SN Ia's) using an elaborate treatment of the Local Thermodynamic Equilbrium (LTE) radiation transport, equation of state and ionization balance, expansion opacity including the cooling by CO, Co(+), and SiO, and a Monte Carlo gamma-ray deposition scheme. The models have an amount of Ni-56 in the range from approximately or equal to 0.1 solar mass up to 0.7 solar mass depending on the density at which the transition from a deflagration to a detonation occurs. Models with a large nickel production give light curves comparable to those of typical Type Ia supernovae. Subluminous supernovae can be explained by models with a low nickel production. Multiband light curves are presented in comparison with the normally bright event SN 1992bc and the subluminous events Sn 1991bg and SN 1992bo to establish the principle that the delayed detonation paradigm in Chandrasekhar mass models may give a common explosion mechanism accounting for both normal and subluminous SN Ia's. Secondary IR-maxima are formed in the models of normal SN Ia's as a photospheric effect if the photospheric radius continues to increase well after maximum light. Secondary maxima appear later and stronger in models with moderate expansion velocities and with radioactive material closer to the surface. Model light curves for subluminous SN Ia's tend to show only one 'late' IR-maximum. In some delayed detonation models shell-like envelopes form, which consist of unburned carbon and oxygen. The formation of molecules in these envelopes is addressed. If the model retains a C/O-envelope and is subluminous, strong vibration bands of CO may appear, typically several weeks past maximum light. CO should be very weak or absent in normal Sn Ia's.

Hoflich, P.↗

Well-posedness of a model for structural acoustic coupling in a cavity enclosed by a thin cylindrical shell

A fully coupled mathematical model describing the interactions between a vibrating thin cylindrical shell and enclosed acoustic field is presented. Because the model will ultimately be used in control applications involving piezoceramic actuators, the loads and material contributions resulting from piezoceramic patches bonded to the shell are included in the discussion. Theoretical and computational issues lead to the consideration of a weak form of the modeling set of partial differential equations (PDE's) and through the use of a semigroup formulation, well-posedness results for the system model are obtained.

Banks, H. T.↗

Panel-Stiffener Debonding and Analysis Using a Shell/3D Modeling Technique

A shear loaded, stringer reinforced composite panel is analyzed to evaluate the fidelity of computational fracture mechanics analyses of complex structures. Shear loading causes the panel to buckle. The resulting out-of-plane deformations initiate skin/stringer separation at the location of an embedded defect. The panel and surrounding load fixture were modeled with shell elements. A small section of the stringer foot, web and noodle as well as the panel skin near the delamination front were modeled with a local 3D solid model. Across the width of the stringer foot, the mixed-mode strain energy release rates were calculated using the virtual crack closure technique. A failure index was calculated by correlating the results with a mixed-mode failure criterion of the graphite/epoxy material. The objective was to study the effect of the fidelity of the local 3D finite element model on the computed mixed-mode strain energy release rates and the failure index.

Krueger, Ronald↗

Panel Stiffener Debonding Analysis using a Shell/3D Modeling Technique

A shear loaded, stringer reinforced composite panel is analyzed to evaluate the fidelity of computational fracture mechanics analyses of complex structures. Shear loading causes the panel to buckle. The resulting out -of-plane deformations initiate skin/stringer separation at the location of an embedded defect. The panel and surrounding load fixture were modeled with shell elements. A small section of the stringer foot, web and noodle as well as the panel skin near the delamination front were modeled with a local 3D solid model. Across the width of the stringer fo to, the mixed-mode strain energy release rates were calculated using the virtual crack closure technique. A failure index was calculated by correlating the results with a mixed-mode failure criterion of the graphite/epoxy material. The objective was to study the effect of the fidelity of the local 3D finite element model on the computed mixed-mode strain energy release rates and the failure index.

Krueger, Ronald↗

Assessment of computational models for multilayered anisotropic plates

The accuracy of the static and vibrational responses of multilayer anisotropic plates, predicted by eight modeling approaches based on two-dimensional shear-deformation theories, was investigated as a function of the lamination and geometric parameters of the shells. The models used included two-dimensional global approximation models, two-dimensional discrete-layer models, and predictor-corrector models. It is shown that the accuracy of the stress and vibrational responses predicted by the different models is strongly dependent on the significance of the transverse shear deformations.

Noor, Ahmed K.↗

The chromosphere of VV cephei and the distribution of circumstellar dust around red giants and supergiants

The work on this project has followed two separate paths of inquiry. The first project was entitled 'the Chromosphere of VV Cephei.' The examination of the archival spectra revealed significant changes in the spectra. Therefore, we obtained additional observing time with IUE to monitor the system during the summer of 1991. Short-term changes continue to be seen in both the overall spectrum and individual line profiles. Work continues on this object. The second project was entitled 'the Distribution of Circumstellar Dust around Red Giants and Supergiants.' A number of cool evolved stars are surrounded by dust shells of sufficient angular size as to appear extended in the IRAS survey data. The aim of this project has been to convolve the predictions of the flux distribution from model dust shells with the IRAS beam profiles in order to reproduce the observed IRAS scans. At the time of the last status report, the cross-scan profiles of the IRAS detectors had just been added to the modeling procedure. For scans in which the star passed near the detector center, there was no significant variation in predicted scan profile for different detectors. Scans in which the detector did not pass over the bright central star had been anticipated to be particularly useful in determining the dust distribution; however, significant differences in the predicted scan profiles were seen for different detector profiles. For this reason, and due to the cross-talk effects discussed in the previous status report, further work on the scans not including a central star has been postponed in favor of further analysis of scans passing over the central star.

Bauer, Wendy Hagen↗

Probing the limits of statistical neutron capture for the r process: Experimental constraints on 141 Cs nuclear level densities

The r-process abundance peaks, particularly near mass number A ∼ 130, reflect underlying nuclear structure effects such as closed neutron shells, yet modeling the nucleosynthesis in this region remains hindered by uncertain neutron-capture rates. These rates are especially sensitive to nuclear level densities (NLDs) and γ-ray strength functions of neutron-rich nuclei, where experimental data are scarce. We present the first experimental constraint on the NLD of 141 Cs using the β-Oslo method, extending sensitivity to the neutron-rich regime near the N = 82 closed shell. Our data allow for critical calibration of microscopic NLD models and reveal that 141 Cs lies near the limit of statistical model applicability. Using this experimental input, we evaluate radiative neutron-capture rates across neighboring isotones using both Hauser–Feshbach (HF) and High Fidelity Resonance (HFR) models. Our results show order-of-magnitude rate increases for nuclei along the N = 86 line, signaling a transition to resonance-dominated capture in this region. These findings underscore the importance of constraining NLDs to improve r-process reaction network predictions, particularly in environments where the validity of statistical models breaks down.

Nuclear level density↗

On sound transmission into an orthotropic shell

A mathematical model is presented for the transmission of airborne noise through the walls of an orthotropic cylindrical shell. Parameters were varied to see how orthotropicity affected noise transmission. When compared to that for an isotropic shell, the cylinder transmission loss was found to be quite sensitive to the ratio of circumferential to axial modulus of elasticity. A modulus ratio greater than unity appears to enhance transmission loss in the mass-controlled region, while a ratio less than unity degrades it. Below the ring frequency, the trends appear to be reversed. The cylinder transmission loss appears to be relatively insensitive to changes in the shear modulus

Koval, L. R.↗

Simplified projection on total spin zero for state preparation on quantum computers

Here, we introduce a simple algorithm for projecting on J = 0 states of a many-body system by performing a series of rotations to remove states with angular momentum projections greater than zero. Existing methods rely on unitary evolution with the two-body operator J 2 , which when expressed in the computational basis contains many complicated Pauli strings requiring Trotterization and leading to very deep quantum circuits. Our approach performs the necessary projections using the one-body operators J x and J z . By leveraging the method of Cartan decomposition, the unitary transformations that perform the projection can be parametrized as a product of a small number of two-qubit rotations, with angles determined by an efficient classical optimization. Given the reduced complexity in terms of gates, this approach can be used to prepare approximate ground states of even-even nuclei by projecting onto the J = 0 component of deformed Hartree-Fock states. We estimate the resource requirements in terms of the universal gate set {H,S, CNOT ,T} and briefly discuss a variant of the algorithm that projects onto J = 1/2 states of a system with an odd number of fermions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Dielectric coated wire antennas

An electrically thin dielectric insulating shell on an antenna composed of electrically thin circular cylindrical wires is examined. A moment method solution is obtained, and the insulating shell is modeled by equivalent volume polarization currents. These polarization currents are related in a simple manner to the surface charge density on the wire antenna. In this way the insulating shell causes no new unknowns to be introduced, and the size of the impedance matrix is the same as for the uninsulated wires. The insulation is accounted for entirely through a modification of the symmetric impedance matrix. This modification influences the current distribution, impedance, efficiency, field patterns, and scattering properties. The theory is compared with measurement for dielectric coated antennas in air.

Richmond, J. H.↗

Nonlinear analysis of damaged stiffened fuselage shells subjected to combined loads

The results of an analytical study of the nonlinear response of stiffened fuselage shells with long cracks are presented. The shells are modeled with a hierarchical modeling strategy that accounts for global and local response phenomena accurately. Results are presented for internal pressure and mechanical bending loads. The effects of crack location and orientation on shell response are described. The effects of mechanical fasteners on the response of a lap joint and the effects of elastic and elastic-plastic material properties on the buckling response of tension-loaded flat panels with cracks are also addressed.

Starnes, James H., Jr.↗

Ensemble Simulations of the 2012 July 12 Coronal Mass Ejection with the Constant-turn Flux Rope Model

Flux-rope-based magnetohydrodynamic modeling of coronal mass ejections (CMEs) is a promising tool for prediction of the CME arrival time and magnetic field at Earth. In this work, we introduce a constant-turn flux rope model and use it to simulate the 2012 July 12 16:48 CME in the inner heliosphere. We constrain the initial parameters of this CME using the graduated cylindrical shell (GCS) model and the reconnected flux in posteruption arcades. We correctly reproduce all the magnetic field components of the CME at Earth, with an arrival time error of approximately 1 hr. We further estimate the average subjective uncertainties in the GCS fittings by comparing the GCS parameters of 56 CMEs reported in multiple studies and catalogs. We determined that the GCS estimates of the CME latitude, longitude, tilt, and speed have average uncertainties of 5°.74, 11°.23, 24°.71, and 11.4%, respectively. Using these, we have created 77 ensemble members for the 2012 July 12 CME. We found that 55% of our ensemble members correctly reproduce the sign of the magnetic field components at Earth. We also determined that the uncertainties in GCS fitting can widen the CME arrival time prediction window to about 12 hr for the 2012 July 12 CME. On investigating the forecast accuracy introduced by the uncertainties in individual GCS parameters, we conclude that the half-angle and aspect ratio have little impact on the predicted magnetic field of the 2012 July 12 CME, whereas the uncertainties in longitude and tilt can introduce relatively large spread in the magnetic field predicted at Earth.

Talwinder Singh↗

Experiments on stiffened conical shell structures using cast epoxy models

Description of a casting technique for fabricating high-quality plastic structural models, and review of results regarding the use of such specimens to parametrically study the effect of base ring stiffness on the critical buckling pressure of a ring-stiffened conical shell. The fabrication technique involves machining a metal mold to the desired configuration and vacuum-drawing the plastic material into the mold. A room-temperature curing translucent thermoset epoxy was the casting material selected. A shell of revolution computer program which employs a nonlinear axisymmetric prebuckling strain field to obtain a bifurcation buckling solution was used to guide the selection of congifurations tested. The shell experimentally exhibited asymmetric collapse behavior, and the ultimate load was considerably higher than the analytical bifurcation prediction. The asymmetric buckling mode shape, however, initially appeared at a pressure near the analysis bifurcation solution.

Williams, J. G.↗

A Melting Layer Model for Passive/Active Microwave Remote Sensing Applications: Model Formulation and Comparison with Observations - Part 1

In this study, a 1-D steady-state microphysical model which describes the vertical distribution of melting precipitation particles is developed. The model is driven by the ice-phase precipitation distributions just above the freezing level at applicable gridpoints of "parent" 3-D cloud-resolving model (CRM) simulations. It extends these simulations by providing the number density and meltwater fraction of each particle in finely separated size categories through the melting layer. The depth of the modeled melting layer is primarily determined by the initial material density of the ice-phase precipitation. The radiative properties of melting precipitation at microwave frequencies are calculated based upon different methods for describing the dielectric properties of mixed phase particles. Particle absorption and scattering efficiencies at the Tropical Rainfall Measuring Mission Microwave Imager frequencies (10.65 to 85.5 GHz) are enhanced greatly for relatively small (approx. 0.1) meltwater fractions. The relatively large number of partially-melted particles just below the freezing level in stratiform regions leads to significant microwave absorption, well-exceeding the absorption by rain at the base of the melting layer. Calculated precipitation backscatter efficiencies at the Precipitation Radar frequency (13.8 GHz) increase in proportion to the particle meltwater fraction, leading to a "bright-band" of enhanced radar reflectivities in agreement with previous studies. The radiative properties of the melting layer are determined by the choice of dielectric models and the initial water contents and material densities of the "seeding" ice-phase precipitation particles. Simulated melting layer profiles based upon snow described by the Fabry-Szyrmer core-shell dielectric model and graupel described by the Maxwell-Garnett water matrix dielectric model lead to reasonable agreement with radar-derived melting layer optical depth distributions. Moreover, control profiles that do not contain mixed-phase precipitation particles yield optical depths that are systematically lower than those observed. Therefore, the use of the melting layer model to extend 3-D CRM simulations appears justified, at least until more realistic spectral methods for describing melting precipitation in high-resolution, 3-D CRM's are implemented.

Olson, William S.↗

Parameter identification of material constants in a composite shell structure

One of the basic requirements in engineering analysis is the development of a mathematical model describing the system. Frequently comparisons with test data are used as a measurement of the adequacy of the model. An attempt is typically made to update or improve the model to provide a test verified analysis tool. System identification provides a systematic procedure for accomplishing this task. The terms system identification, parameter estimation, and model correlation all refer to techniques that use test information to update or verify mathematical models. The goal of system identification is to improve the correlation of model predictions with measured test data, and produce accurate, predictive models. For nonmetallic structures the modeling task is often difficult due to uncertainties in the elastic constants. A finite element model of the shell was created, which included uncertain orthotropic elastic constants. A modal survey test was then performed on the shell. The resulting modal data, along with the finite element model of the shell, were used in a Bayes estimation algorithm. This permitted the use of covariance matrices to weight the confidence in the initial parameter values as well as confidence in the measured test data. The estimation procedure also employed the concept of successive linearization to obtain an approximate solution to the original nonlinear estimation problem.

Martinez, David R.↗