Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuum modeling”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

RXTE Observations of the 1A 1118-61 in an Outburst, and the Discovery of a Cyclotron Line

We present the analysis of RXTE monitoring data obtained during the January 2009 outburst of the hard X-ray transient IA 1118-61. Using these observations the broadband (3.5-120 keV) spectrum of the source was measured for the first time ever. We have found that the broadband continuum spectrum of the source is similar to other accreting pulsars and is well described by several conventionally used phenomenological models. We have discovered that regardless of the applied continuum model, a prominent broad absorption feature at approx. 55 keV is observed. We interpret this feature as a Cyclotron Resonance Scattering Feature (CRSF). The observed CRSF energy is one of the highest known and corresponds to a magnetic field of B approx. 4.8 x 10(exp 12) G in the scattering region. Furthermore, our data suggests an iron emission line presence, which was not reported previously for lA 1118-61 as well. Timing properties of the source, including a strong spin-up, were found to be similar to those observed by CGRO/BATSE during the previous outburst, however the broadband capabilities of RXTE reveal a more complicated energy dependency of the pulse-profile.

Doroshenko, V.↗
NASA NTRS1982

Molecular dynamics simulation of propagating cracks

Steady state crack propagation is investigated numerically using a model consisting of 236 free atoms in two (010) planes of bcc alpha iron. The continuum region is modeled using the finite element method with 175 nodes and 288 elements. The model shows clear (010) plane fracture to the edge of the discrete region at moderate loads. Analysis of the results obtained indicates that models of this type can provide realistic simulation of steady state crack propagation.

Mullins, M.↗
NASA NTRS2018

Determination of Ballistic Limit for IM7/8552 Using LS-DYNA MAT 261

The goal of the NASA ACC High Energy Dynamic Impact Project is to determine the state of the art of dynamic fracture simulations for high velocity impact for composite fuselage shielding applications. Using a building block approach, several computational models considered under NASA ACC are being validated against test data, starting at unconfigured panels and progressing to configured panels under combined out-of-plane and in-plane loading due to ballistic impact. The computational models being evaluated in this project include MAT 162, MAT 213, MAT 261, SPG, and Peridynamics. In this paper, the simulation results using LS-Dyna Material MAT 261 are presented. In particular, a series of blind predictions for unconfigured panels were performed to determine the ballistic limit or V50 velocity. MAT 261 employs failure approach that is generally physically-based using fracture toughness criteria. The overall material model relies on typical ply-level stiffness properties, similar to MAT 162 and other composite continuum damage material models. The fracture toughness values are based on standard tests, and thus are not subject to extensive calibration. This approach is more efficient than performing extensive optimization studies for calibration of parameters. Also, this approach of relying on physical properties reduces the uncertainty of results, as questions concerning the quality and extent of the calibration studies is no longer relevant. However, it was found that carefully controlled coupon-level tests are needed to accurately obtain the required fracture toughness values. Additionally, it should be noted that there is one significant parameter in MAT 261 that does appear to require calibration, and that is the overall failure strain. This is the strain at which the element is deleted, and is not the same as the strain at which damage begins to accumulate. This failure strain is termed EFS (Effective Failure Strain), and is the maximum effective stain for element failure. Simulations have shown that this value will significantly affect impact response and failure. The paper presents the effect of this element failure strain parameter, along with possible uncertainties in fracture toughness values. With an adjusted appropriate value for EFS, it is seen that simulation results compare well with impact test data for predicted penetration velocity.

Byar, Alan↗
NASA NTRSDate not supplied

Understanding Variability in Mechanical Behavior of Adhesively Bonded Sandwich Joints due to Stochastic Material Damage

A mesoscale progressive damage analysis was applied to the analysis of recently tested adhesively bonded sandwich panels, in which ten nominally identical tension loaded panels failed in three different failure modes, with peakloads ranging from 13.0 –16.1 kips. A parametric finite element model for progressive damage analysis of an adhesively bonded joint (ABJ) considering matrix and fiber damage, and delamination cracks in unidirectional and fabric composite sections of the joint, as well as core crushing, and adhesive fracture was developed. The model was developed using Abaqus/Explicit and utilized continuum damage mechanics material models for intralaminar damage in the facesheet and doubler plies, cohesive elements for delamination and adhesive debonding, and discretized cells considering plasticity to model crushing of the aluminum honeycomb core. Analysis predictions obtained with a pristine model agreed well with the average experimental peak load, strain in the joint, and post-test damage states. Then, various methods to incorporate stochastic material variability into ABJ model were evaluated. These methods included simulation of pre-existing matrix damage, pre-existing adhesive porosity, and uniform modification of adhesive properties. These material and manufacturing defects changed the damage mode predictions in some cases. The analysis of stochastic strength due to the considered defects revealed unintuitive damage mode interaction in the ABJ, explaining some structure-property relationships observed experimentally.

Richard Larson↗
NASA NTRS2006

Molecular-dynamics Simulation-based Cohesive Zone Representation of Intergranular Fracture Processes in Aluminum

A traction-displacement relationship that may be embedded into a cohesive zone model for microscale problems of intergranular fracture is extracted from atomistic molecular-dynamics simulations. A molecular-dynamics model for crack propagation under steady-state conditions is developed to analyze intergranular fracture along a flat 99 [1 1 0] symmetric tilt grain boundary in aluminum. Under hydrostatic tensile load, the simulation reveals asymmetric crack propagation in the two opposite directions along the grain boundary. In one direction, the crack propagates in a brittle manner by cleavage with very little or no dislocation emission, and in the other direction, the propagation is ductile through the mechanism of deformation twinning. This behavior is consistent with the Rice criterion for cleavage vs. dislocation blunting transition at the crack tip. The preference for twinning to dislocation slip is in agreement with the predictions of the Tadmor and Hai criterion. A comparison with finite element calculations shows that while the stress field around the brittle crack tip follows the expected elastic solution for the given boundary conditions of the model, the stress field around the twinning crack tip has a strong plastic contribution. Through the definition of a Cohesive-Zone-Volume-Element an atomistic analog to a continuum cohesive zone model element - the results from the molecular-dynamics simulation are recast to obtain an average continuum traction-displacement relationship to represent cohesive zone interaction along a characteristic length of the grain boundary interface for the cases of ductile and brittle decohesion. Keywords: Crack-tip plasticity; Cohesive zone model; Grain boundary decohesion; Intergranular fracture; Molecular-dynamics simulation

Yamakov, Vesselin I.↗
NASA NTRS1989

A study of material damping in large space structures

A constitutive model was developed for predicting damping as a function of damage in continuous fiber reinforced laminated composites. The damage model is a continuum formulation, and uses internal state variables to quantify damage and its subsequent effect on material response. The model is sensitive to the stacking sequence of the laminate. Given appropriate baseline data from unidirectional material, and damping as a function of damage in one crossply laminate, damage can be predicted as a function of damage in other crossply laminates. Agreement between theory and experiment was quite good. A micromechanics model was also developed for examining the influence of damage on damping. This model explicitly includes crack surfaces. The model provides reasonable predictions of bending stiffness as a function of damage. Damping predictions are not in agreement with the experiment. This is thought to be a result of dissipation mechanisms such as friction, which are not presently included in the analysis.

Highsmith, A. L.↗
NASA NTRS1994

Millimeter continuum measurements of circumstellar dust around very young low-mass stars

We investigate the question of disk formation during the protostar phase. We build on the results of Keene and Masson (1990) whose analysis of L1551 showed the millimeter continuum emission comes from both an unresolved circumstellar component, i.e., a disk and an extended cloud core. We model the dust continuum emission from the cloud core and show how it is important at 1.3 mm but negligible at 2.7 mm. Combining new 2.7 mm Owens Valley Interferometer data of IRAS-Dense cores with data from the literature we conclude that massive disks are also seen toward a number of other sources. However, 1.3 mm data from the IRAM 30 m telescope for a larger sample shows that massive disks are relatively rare, occurring around perhaps 5% of young embedded stars. This implies that either massive disks occur briefly during the embedded phase or that relatively few young stars form massive disks. At 1.3 mm the median flux of IRAS-Dense cores is nearly the same as T Tauri stars in the sample of Beckwith et al. (1990). We conclude that the typical disk mass during the embedded phase is nearly the same or less than the typical disk mass during the T Tauri phase.

Terebey, S.↗
NASA NTRS1993

Rarefied-flow Shuttle aerodynamics model

A rarefied-flow shuttle aerodynamic model spanning the hypersonic continuum to the free molecule-flow regime was formulated. The model development has evolved from the High Resolution Accelerometer Package (HiRAP) experiment conducted on the Orbiter since 1983. The complete model is described in detail. The model includes normal and axial hypersonic continuum coefficient equations as functions of angle-of-attack, body flap deflection, and elevon deflection. Normal and axial free molecule flow coefficient equations as a function of angle-of-attack are presented, along with flight derived rarefied-flow transition bridging formulae. Comparisons are made with data from the Operational Aerodynamic Design Data Book (OADDB), applicable wind-tunnel data, and recent flight data from STS-35 and STS-40. The flight-derived model aerodynamic force coefficient ratio is in good agreement with the wind-tunnel data and predicts the flight measured force coefficient ratios on STS-35 and STS-40. The model is not, however, in good agreement with the OADDB. But, the current OADDB does not predict the flight data force coefficient ratios of either STS-35 or STS-40 as accurately as the flight-derived model. Also, the OADDB differs with the wind-tunnel force coefficient ratio data.

Blanchard, Robert C.↗
NASA NTRSDate not supplied

Higher-Order Approximations for Stabilizing Zero-Energy Modes in Peridynamics Crystal Plasticity Models with Large Horizon Interactions

The non-ordinary state-based peridynamics theory combines non-local dynamic techniques with a desirable correspondence material principle, allowing for the use of continuum mechanics constitutive models. Such an approach presents a unique capability for solving problems involving discontinuities (e.g., strain localization, fracture, and fragmentation). However, the correspondence-based peridynamics models often suffer from zero-energy mode instabilities in numerical implementation, primarily due to the approximations of the non-local deformation gradient tensor. This paper focuses on a computational scheme for eliminating the zero-energy mode oscillations using a choice of influence functions that improve the truncation error in a higher-order Taylor series expansion of the deformation gradient. The novelty here is a tensor-based derivation of the linear constraint equations, which can be used to systematically identify the particle interaction weight functions for various user-specified horizon radii. In this paper, the proposed higher-order stabilization scheme is demonstrated for multi-dimensional examples involving polycrystalline and composite microstructures, along with comparisons against conventional finite element methods. The proposed stabilization scheme is shown to be highly effective in suppressing the spurious zero-energy mode oscillations in all numerical examples while enabling efficient simulations of strain localizations across material interfaces.

Non-Ordinary State-Based Peridynamics↗
NASA NTRS2010

Investigation of a Macromechanical Approach to Analyzing Triaxially-Braided Polymer Composites

A macro level finite element-based model has been developed to simulate the mechanical and impact response of triaxially-braided polymer matrix composites. In the analytical model, the triaxial braid architecture is simulated by using four parallel shell elements, each of which is modeled as a laminated composite. The commercial transient dynamic finite element code LS-DYNA is used to conduct the simulations, and a continuum damage mechanics model internal to LS-DYNA is used as the material constitutive model. The material stiffness and strength values required for the constitutive model are determined based on coupon level tests on the braided composite. Simulations of quasi-static coupon tests of a representative braided composite are conducted. Varying the strength values that are input to the material model is found to have a significant influence on the effective material response predicted by the finite element analysis, sometimes in ways that at first glance appear non-intuitive. A parametric study involving the input strength parameters provides guidance on how the analysis model can be improved.

Goldberg, Robert K.↗
NASA NTRS1990

An isolated, well-defined infrared cirrus cloud

High spatial resolution H I measurements of a high-latitude FIR cirrus cloud (l = 9 deg, b = 51 deg) were obtained with the Arecibo 310 m telescope. The cloud is well isolated spatially and has a velocity well-separated from other H I clouds, providing an ideal laboratory for studying diffuse interstellar clouds. The Iv(100)/N(H) and Iv(60)/Iv(100) ratios are consistent with the Draine and Anderson model which includes a distribution of grain sizes down to 3 A. The small grains undergo transient heating and affect the 60 micron FIR continuum emission. Current models, however, underpredict Iv(12)/Iv(100) and Iv(25)/Iv(100) by more than a factor of 2. The effects small grains have on the thermal equilibrium of diffuse interstellar clouds are discussed.

Herter, T.↗
NASA NTRS1993

Simultaneous ROSAT, Ginga, VLA, IUE, and optical observations of the bright quasar H1821+643

Simultaneous optical, UV, and X-ray spectroscopy of the low-redshift high-luminosity QSO H1821+643 is presented in order to examine the connection between the optical UV bump and the soft X-ray excess. The correlation of the short-term variability in the optical/UV and soft X-ray band is examined, and the applicability of the accretion disk models to the simultaneous continuum energy distribution is tested. By fitting the continuum energy distribution with the thin accretion disk model, it was found that bare disk models cannot easily account for the emission in both the UV and the steep excess in the soft X-ray band. The intrinsic X-ray spectrum of H1821+643 exhibited a typical power-law slope above about 0.3 keV, as detected by ROSAT and Ginga. The first detection of H1821+643 in the radio band is reported. This radio-quiet QSO was found to have a substantial radio flux at 20 cm and a typical steep spectral slope.

Kolman, Michiel↗
NASA NTRS1995

Water vapor absorption in the atmospheric window at 239 GHz

Absolute absorption rates of pure water vapor and mixtures of water vapor and nitrogen have been measured in the atmospheric window at 239 GHz. The dependence on pressure as well as temperature has been obtained. The experimental data are compared with several theoretical or empirical models, and satisfactory agreement is obtained with the models involving a continuum; in the case of pure water vapor, the continuum contribution based upon recent theoretical developments gives good results. The temperature dependence is stronger than that proposed in a commonly used atmospheric transmission model.

Bauer, A.↗
NASA NTRS2013

Vertical Settling and Radial Segregation of Large Dust Grains in the Circumstellar Disk of the Butterfly Star

Circumstellar disks are considered as the environment for the formation of planets. The growth of dust grains in these disks is the first step in the core accretion-gas capture planet formation scenario. Indicators as well as evidence of disk evolution can be traced in spatially resolved images and the spectral energy distribution (SED) of these objects. Aims. We aim at developing a model for the dust phase of the edge-on orientated circumstellar disk of the Butterfly Star which allows one to fit observed multi-wavelength images as well as the SED simultaneously. Our model is based on spatially resolved high angular resolution observations at 1.3 mm, 894 μm, 2.07 μm, 1.87 μm, 1.60 μm and 1.13 μm as well as an extensively covered SED ranging from 12 μm to 2.7 mm, including a detailed spectrum obtained with the Spitzer Space Telescope in the range from 12 μm to 38 μm. A parameter study based on a grid search method, involving the detailed analysis of every parameter, was performed to constrain the disk parameters and find the best-fit model for the independent observations. The individual observations were modeled simultaneously, using our continuum radiative transfer code MC3D. We derived a model that is capable of reproducing all of the observations of the disk at the same time. We find quantitative evidence for grain growth up to ~100 μm sized particles, vertical settling of larger dust grains toward the disk midplane, and radial segregation of the latter toward the central star. Our results are consistent with current theoretical models for the evolution of circumstellar disks and the early stages of planet formation.

Radiative transfer↗
NASA NTRSDate not supplied

New Developments in the CompDam Fiber Kinking Model for the Interaction of Kinking and Splitting Cracks with Application to Open Hole Compression Specimens

This paper describes two new additions to the NASA fiber kinking model: 1) A failure index for fiber kinking is developed to facilitate interaction between kinking and matrix cracking and 2) A criterion for fiber rupture is introduced to account for fiber breakage before or after fiber kinking. The model formulation is described with the aid of graphical Considered constructions to elucidate the behavior. Verification studies are conducted with single element models to show the behavior of the model in light of the new developments and to evaluate the effect of numerical parameters of the model. The model is used to analyze open hole compression specimens with soft, quasi-isotropic, and hard layups. Analysis results are shown for a conventional continuum damage mechanics model as well as the proposed model to highlight the differences between the two. Comparison with test data highlights the capabilities and remaining gaps in the quantitative prediction of damage evolution in open hole compression. Overall, the results show that the proposed model is capable of representing the kinematics of fiber kinking, and specifically the role of fiber kinking in open hole compression.

Andrew Cole Bergan↗
NASA NTRS1984

Apparent solar temperature enhancement due to large-amplitude waves

The effect of slow-mode acoustic-type MHD waves propagating outward in the solar atmosphere on the temperature structures predicted by empirical models is investigated analytically. A model is constructed, and numerical results are presented for wing intensities, line profiles, temperature enhancements, waves with higher energy flux, temperature depression, and the Si continuum. The flux in the MgII and CaII UV resonance lines is found to be increased relative to that in the IR continuum, leading to model temperatures which depend systematically on which observations are used in the computation. It is suggested that mechanical heating may take place in smaller regions such as flux tubes rather than uniformly over the surface.

Kalkofen, W.↗
NASA NTRSDate not supplied

New Developments in the CompDam Fiber Kinking Model for the Interaction

This paper describes two new additions to the CompDam fiber kinking model: 1) a failure index for fiber kinking to facilitate interaction between kinking and matrix cracking and 2) a strain-based criterion for fiber rupture to account for fiber breakage before or after fiber kinking. The model formulation is described with the aid of the graphical Considère construction, which is used to determine the conditions for fiber kinking. Verification studies are conducted with single element models to show the behavior of the model in light of the new developments and to evaluate the effect of numerical parameters of the model. The model is used to analyze open hole compression (OHC) specimens with soft, quasi-isotropic, and hard layups. Analysis results are shown for a conventional continuum damage mechanics model using a trilinear constitutive law for fiber damage and for the proposed model to highlight the differences between the two. Comparison with test data highlights the capabilities and remaining gaps in the quantitative prediction of the damage and failure response. While both models can predict many features of the OHC response, accurate prediction of kinking onset in the OHC specimens remains elusive. The results suggest that it is necessary to predict the combined effects of compression and shear loading as well as the structural effects (e.g., ply thickness, strain gradient effect) accurately to match the experimental observations of fiber kinking onset in the three OHC laminates.

Andrew C. Bergan↗
NASA NTRS2002

Constitutive Modeling of Nanotube-Reinforced Polymer Composites

In this study, a technique is presented for developing constitutive models for polymer composite systems reinforced with single-walled carbon nanotubes (SWNT). Because the polymer molecules are on the same size scale as the nanotubes, the interaction at the polymer/nanotube interface is highly dependent on the local molecular structure and bonding. At these small length scales, the lattice structures of the nanotube and polymer chains cannot be considered continuous, and the bulk mechanical properties can no longer be determined through traditional micromechanical approaches that are formulated by using continuum mechanics. It is proposed herein that the nanotube, the local polymer near the nanotube, and the nanotube/polymer interface can be modeled as an effective continuum fiber using an equivalent-continuum modeling method. The effective fiber serves as a means for incorporating micromechanical analyses for the prediction of bulk mechanical properties of SWNT/polymer composites with various nanotube lengths, concentrations, and orientations. As an example, the proposed approach is used for the constitutive modeling of two SWNT/polyimide composite systems.

Odegard, G. M.↗