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 181 records · Page 10

Simulating the Effect of Space Vehicle Environments on Directional Solidification of a Binary Alloy

Space microgravity missions are designed to provide a microgravity environment for scientific experiments, but these missions cannot provide a perfect environment, due to vibrations caused by crew activity, on-board experiments, support systems (pumps, fans, etc.), periodic orbital maneuvers, and water dumps. Therefore, it is necessary to predict the impact of these vibrations on space experiments, prior to performing them. Simulations were conducted to study the effect of the vibrations on the directional solidification of a dendritic alloy. Finite element ca!cu!attie?ls were dme with a simd2titcr based on a continuum model of dendritic solidification, using the Fractional Step Method (FSM). The FSM splits the solution of the momentum equation into two steps: the viscous intermediate step, which does not enforce continuity; and the inviscid projection step, which calculates the pressure and enforces continuity. The FSM provides significant computational benefits for predicting flows in a directionally solidified alloy, compared to other methods presently employed, because of the efficiency gains in the uncoupled solution of velocity and pressure. finite differences, arises when the interdendritic liquid reaches the eutectic temperature and concentration. When a node reaches eutectic temperature, it is assumed that the solidification of the eutectic liquid continues at constant temperature until all the eutectic is solidified. With this approach, solidification is not achieved continuously across an element; rather, the element is not considered solidified until the eutectic isotherm overtakes the top nodes. For microgravity simulations, where the convection is driven by shrinkage, it introduces large variations in the fluid velocity. When the eutectic isotherm reaches a node, all the eutectic must be solidified in a short period, causing an abrupt increase in velocity. To overcome this difficulty, we employed a scheme to numerically predict a more accurate value for the rate of change of fraction of liquid as the liquid in an element solidifies. The new method enables us to contrast results of simulations in which the alloy is subjected to no gravity or a steady-state acceleration versus simulations when the alloy is subjected to vibration disturbances; therefore, the effect of vibration disturbances can be assessed more accurately. To assess the impact of these vibration-perturbations, transient accelerometer data from a space shuttle mission are used as inputs for the simulation model. These on-orbit acceleration data were obtained from the Microgravity Science Division at Glenn Research Center (GRC- MSD) and are applied to the buoyancy term of the momentum equation in a simulation of a Pb-5.8 wt. % Sb alloy that solidifies in a thermal gradient of 4000 K/m and a translation velocity of 3 p d s . Figure 2 shows the vertical velocity of a node that begins in the all-liquid region and subsequently solidifies; the vibrations are applied at 5000 seconds in this simulation. An important difficulty, common to all solidification models based on finite elements or 2 The magnitudes of the velocity oscillations that are vibration-induced are very small and acceptable. The biggest concern is whether the concentration of the liquid near the dendrite tips is distorted because of the vibration-induced perturbations. Results for this case show no concentration oscillations present in the all-liquid region.

Westra, D. G.↗

Predicting the Influence of Nano-Scale Material Structure on the In-Plane Buckling of Orthotropic Plates

A multi-scale analysis of the structural stability of a carbon nanotube-polymer composite material is developed. The influence of intrinsic molecular structure, such as nanotube length, volume fraction, orientation and chemical functionalization, is investigated by assessing the relative change in critical, in-plane buckling loads. The analysis method relies on elastic properties predicted using the hierarchical, constitutive equations developed from the equivalent-continuum modeling technique applied to the buckling analysis of an orthotropic plate. The results indicate that for the specific composite materials considered in this study, a composite with randomly orientated carbon nanotubes consistently provides the highest values of critical buckling load and that for low volume fraction composites, the non-functionalized nanotube material provides an increase in critical buckling stability with respect to the functionalized system.

Gates, Thomas S.↗

The Effect of Chemical Functionalization on Mechanical Properties of Nanotube/Polymer Composites

The effects of the chemical functionalization of a carbon nanotube embedded in a nanotube/polyethylene composite on the bulk elastic properties are presented. Constitutive equations are established for both functionalized and non-functionalized nanotube composites systems by using an equivalent-continuum modeling technique. The elastic properties of both composites systems are predicted for various nanotube lengths, volume fractions, and orientations. The results indicate that for the specific composite material considered in this study, most of the elastic stiffness constants of the functionalized composite are either less than or equal to those of the non-functionalized composite.

Odegard, G. M.↗

Viewpoint 9--molecular structure of aqueous interfaces

In this review we summarize recent progress in our understanding of the structure of aqueous interfaces emerging from molecular level computer simulations. It is emphasized that the presence of the interface induces specific structural effects which, in turn, influence a wide variety of phenomena occurring near the phase boundaries. At the liquid-vapor interface, the most probable orientations of a water molecule is such that its dipole moment lies parallel to the interface, one O-H bond points toward the vapor and the other O-H bond is directed toward the liquid. The orientational distributions are broad and slightly asymmetric, resulting in an excess dipole moment pointing toward the liquid. These structural preferences persist at interfaces between water and nonpolar liquids, indicating that the interactions between the two liquids in contact are weak. It was found that liquid-liquid interfaces are locally sharp but broadened by capillary waves. One consequence of anisotropic orientations of interfacial water molecules is asymmetric interactions, with respect to the sign of the charge, of ions with the water surface. It was found that even very close to the surface ions retain their hydration shells. New features of aqueous interfaces have been revealed in studies of water-membrane and water-monolayer systems. In particular, water molecules are strongly oriented by the polar head groups of the amphiphilic phase, and they penetrate the hydrophilic head-group region, but not the hydrophobic core. At infinite dilution near interfaces, amphiphilic molecules exhibit behavior different from that in the gas phase or in bulk water. This result sheds new light on the nature of hydrophobic effect in the interfacial regions. The presence of interfaces was also shown to affect both equilibrium and dynamic components of rates of chemical reactions. Applications of continuum models to interfacial problems have been, so far, unsuccessful. This, again, underscores the importance of molecular-level information about interfaces.

Review, Tutorial↗

The Dynamics of Coupled Oscillator Phase Control

Arrays of coupled oscillators have been proposed as means of realizing high power rf sources via coherent spatial power combining. In such applications, a uniform phase distribution over the aperture is usually desired. However, it has been shown that by detuning some of the oscillators away from the oscillation frequency of the ensemble of oscillators, one may achieve other useful aperture phase distributions. Of particular interest among those achievable are linear phase distributions because these result in steering of the output rf beam away from the broadside direction. The theory describing the behavior of such arrays of coupled oscillators is quite complicated since the phenomena involved are inherently nonlinear. However, a simplified theory has been developed which facilitates intuitive understanding. This simplified theory is based on a "continuum model" in which the aperture phase is represented by a continuous function of the aperture coordinates. A challenging aspect of the development of this theory is the derivation of appropriate boundary conditions at the edges or ends of the array.

Pogorzelski, R. J.↗

Buckling of Carbon Nanotube-Reinforced Polymer Laminated Composite Materials Subjected to Axial Compression and Shear Loadings

A multi-scale method to predict the stiffness and stability properties of carbon nanotube-reinforced laminates has been developed. This method is used in the prediction of the buckling behavior of laminated carbon nanotube-polyethylene composites formed by stacking layers of carbon nanotube-reinforced polymer with the nanotube alignment axes of each layer oriented in different directions. Linking of intrinsic, nanoscale-material definitions to finite scale-structural properties is achieved via a hierarchical approach in which the elastic properties of the reinforced layers are predicted by an equivalent continuum modeling technique. Solutions for infinitely long symmetrically laminated nanotube-reinforced laminates with simply-supported or clamped edges subjected to axial compression and shear loadings are presented. The study focuses on the influence of nanotube volume fraction, length, orientation, and functionalization on finite-scale laminate response. Results indicate that for the selected laminate configurations considered in this study, angle-ply laminates composed of aligned, non-functionalized carbon nanotube-reinforced lamina exhibit the greatest buckling resistance with 1% nanotube volume fraction of 450 nm uniformly-distributed carbon nanotubes. In addition, hybrid laminates were considered by varying either the volume fraction or nanotube length through-the-thickness of a quasi-isotropic laminate. The ratio of buckling load-to-nanotube weight percent for the hybrid laminates considered indicate the potential for increasing the buckling efficiency of nanotube-reinforced laminates by optimizing nanotube size and proportion with respect to laminate configuration.

Riddick, J. C.↗

Two Suzaku Observations of the Cyclotron Line Source 4U 1907+09

Suzaku observed the accreting pulsar 4U 1907+09 in May 2006 and in April 2007. We present and compare the light curves and broad band spectra of both observations. The approximately 10 mCrab source is strongly variable and shows a rich phenomenology. While the 2006 observation is dominated by typical, non-absorption-related dipping behavior, the 2007 observation shows equally typical flaring behavior in both the XIS and PIN light curves. Above 2 keV the average combined XIS and PIN broad band spectrum can be well described by standard pulsar continuum models, modified by the fundamental cyclotron resonant scattering feature at approximately 19 keV. We model the cyclotron line using the standard phenomenological Gaussian model as well as a new physically motivated Green's functions based approach (Schoenherr et al., 2007, A&A 472, 353). Below 2 keV the XIS spectra allow for the first detection of a soft excess in this source, in addition to providing the best determination of the parameters of the weak iron line at 6.4 keV (EW approximately 70 eV) to date.

Pottschmidt, Katja↗

Broad-Band Spectroscopy of Hercules X-1 with Suzaku

Hercules X-1 was observed with Suzaku in the main-on state from 2005 to 2010. The 0.4- 100 keV wide-band spectra obtained in four observations showed a broad hump around 4-9 keV in addition to narrow Fe lines at 6.4 and 6.7 keV. The hump was seen in all the four observations regardless of the selection of the continuum models. Thus it is considered a stable and intrinsic spectral feature in Her X-1. The broad hump lacked a sharp structure like an absorption edge. Thus it was represented by two different spectral models: an ionized partial covering or an additional broad line at 6.5 keV. The former required a persistently existing ionized absorber, whose origin was unclear. In the latter case, the Gaussian fitting of the 6.5-keV line needs a large width of sigma = 1.0-1.5 keV and a large equivalent width of 400-900 eV. If the broad line originates from Fe fluorescence of accreting matter, its large width may be explained by the Doppler broadening in the accretion flow. However, the large equivalent width may be inconsistent with a simple accretion geometry.

X-rays: stars↗

Searching for Decaying Dark Matter in Deep XMM-Newton Observation of the Draco Dwarf Spheroidal

We present results of a search for the 3.5 keV emission line in our recent very long (approx. 1.4 Ms) XMM-Newton observation of the Draco dwarf spheroidal galaxy. The astrophysical X-ray emission from such dark matter-dominated galaxies is faint, thus they provide a test for the dark matter origin of the 3.5 keV line previously detected in other massive, but X-ray bright objects, such as galaxies and galaxy clusters. We do not detect a statistically significant emission line from Draco; this constrains the lifetime of a decaying dark matter particle to tau >(7-9) × 10(exp 27) s at 95% CL (combining all three XMM-Newton cameras; the interval corresponds to the uncertainty of the dark matter column density in the direction of Draco). The PN camera, which has the highest sensitivity of the three, does show a positive spectral residual (above the carefully modeled continuum) at E = 3.54 +/- 0.06 keV with a 2.3(sigma) significance. The two MOS cameras show less-significant or no positive deviations, consistently within 1(sigma) with PN. Our Draco limit on tau is consistent with previous detections in the stacked galaxy clusters, M31 and the Galactic Centre within their 1 − 2(sigma) uncertainties, but is inconsistent with the high signal from the core of the Perseus cluster (which has itself been inconsistent with the rest of the detections). We conclude that this Draco observation does not exclude the dark matter interpretation of the 3.5 keV line in those objects.

dark matter↗

Disk-Wind Connection During the Heartbeats of GRS 1915+105

Disk and wind signatures are seen in the soft state of Galactic black holes, while the jet is seen in the hard state. Here we study the disk-wind connection in the Rho class of variability in GRS 1915+105 using a joint NuSTAR-Chandra observation. The source shows 50 s limit cycle oscillations. By including new information provided by the reflection spectrum and using phase-resolved spectroscopy, we find that the change in the inner disk inferred from the blackbody emission is not matched by reflection measurements. The latter is almost constant, independent of the continuum model. The two radii are comparable only if the disk temperature color correction factor changes, an effect that could be due to the changing opacity of the disk caused by changes in metal abundances. The disk inclination is similar to that inferred from the jet axis, and oscillates by approx.10 deg. The simultaneous Chandra data show the presence of two wind components with velocities between 500 and 5000 km s(exp. −1), and possibly two more with velocities reaching 20,000 km s(exp. −1) (approx. 0.06 c). The column densities are approx. 5 × 10(exp. 22) cm(exp. −2). An upper limit to the wind response time of 2 s is measured, implying a launch radius of less than 6 × 10(exp. 10) cm. The changes in wind velocity and absorbed flux require the geometry of the wind to change during the oscillations, constraining the wind to be launched from a distance of 290-1300 r (sub g) from the black hole. Both data sets support fundamental model predictions in which a bulge originates in the inner disk and moves outward as the instability progresses.

accretion↗

3D Mapping of the Neutral X-Ray Absorption in the Local Interstellar Medium: The Gaia and XMM-Newton Synergy

We present a three-dimensional map of the hydrogen density distribution in the Galactic interstellar medium. The hydrogen-equivalent column densities were obtained from the Exploring the X-ray Transient and variable Sky project (EXTRAS) which provides equivalent NH values from X-ray spectral fits of observations within the XMM-Newton Data Release. EXTRAS include multiple fits for each source, allowing an accurate determination of the equivalent column densities, which depends on the continuum modelling of the spectra. A cross-correlation between the EXTRAS catalogue and the first Gaia Data Release was performed in order to obtain accurate parallax and distance measurements. We use a Bayesian method explained in Rezaei Kh. et al. (2017) in order to predict the most probable distribution of the density at any arbitrary point, even for lines of sight along which there are no initial observation. The resulting map shows small-scale density structures which could not have been modelled by using analytic density profiles. In this paper, we present a proof of concept of the kind of science possible with the synergy of these catalogues. However, given the systematic uncertainties connected to the source identification and to the dependence of NH on the spectral model, the present maps should be considered qualitatively at this point.

Gatuzz, Efrain↗

Broad-band Spectroscopy of Hercules X-1 with Suzaku

Hercules X-1 was observed with Suzaku in the main-on state from 2005 to 2010. The 0.4–100 keV wide-band spectra obtained in four observations showed a broad hump around 4–9 keV in addition to narrow Fe lines at 6.4 and 6.7 keV. The hump was seen in all the four observations regardless of the selection of the continuum models. Thus it is considered a stable and intrinsic spectral feature in Her X-1. The broad hump lacked a sharp structure like an absorption edge. Thus it was represented by two different spectral models: an ionized partial covering or an additional broad line at 6.5 keV. The former required a persistently existing ionized absorber, whose origin was unclear. In the latter case, the Gaussian fitting of the 6.5-keV line needs a large width of σ = 1.0–1.5 keV and a large equivalent width of 400–900 eV. If the broad line originates from Fe fluorescence of accreting matter, its large width may be explained by the Doppler broadening in the accretion flow. However, the large equivalent width may be inconsistent with a simple accretion geometry.

Asami, Fumi↗

Radiation Interaction with Particulates in Planetary Entry Flows

A hypersonic particle phase solver based on a continuum model is implemented in the NASA HyperSolve CFD code and used to predict particle phase conditions in a dusty Mars entry scenario. Black-body radiative emission from high temperature dust particles is predicted for a Mars relevant entry problem, and the emission integrated to an entry vehicle using a ray tracing approach to estimate surface heat flux. The surface heating due to this mechanism is found to be small relative to convective and radiation from the gas phase emission. Additionally, the effect of atmospheric haze particulates vaporizing in the Titan atmosphere is studied. The increase in carbon availability due to the composition of the haze particles is found to increase the gas radiative heating to the vehicle by up to 3% on the forebody of a representative Titan entry vehicle with approximately 0.15% haze in the freestream by mass.

Andrew Hinkle↗

Assessment of Radiation Scattering by Dust in Mars Entry

The impact of gas phase radiation scattering by dust particle presence in Mars entry scenarios is assessed using a suite of numerical tools. A hypersonic gas and dilute particle multiphase flow is simulated using a mixed continuum model and an efficient size distribution quadrature scheme. Radiative transport is simulated with a hybrid ray-tracing and P1 spherical harmonics model. The NASA HARA radiation code and a Lorentz-Mie scattering model are used to predict the gas-particle mixture radiative properties. Application of the numerical suite to a representative case based on the Mars 2020 entry shows that, even at conditions in excess of the strongest recorded dust storm, the augmentation in surface radiative heat flux due to scattering alone is lower than 1%, and therefore negligible. Predictions with even higher dust loading conditions confirm that even highly uncertain dust mass profiles are unlikely to be a concern for the scattering augmentation mechanism.

Andrew Hinkle↗

Continuum Damage Mechanics Models for the Analysis of Progressive Failure in Open-Hole Tension Laminates

The performance of a state-of-the-art continuum damage mechanics model for interlaminar damage, coupled with a cohesive zone model for delamination is examined for failure prediction of quasi-isotropic open-hole tension laminates. Limitations of continuum representations of intra-ply damage and the effect of mesh orientation on the analysis predictions are discussed. It is shown that accurate prediction of matrix crack paths and stress redistribution after cracking requires a mesh aligned with the fiber orientation. Based on these results, an aligned mesh is proposed for analysis of the open-hole tension specimens consisting of different meshes within the individual plies, such that the element edges are aligned with the ply fiber direction. The modeling approach is assessed by comparison of analysis predictions to experimental data for specimen configurations in which failure is dominated by complex interactions between matrix cracks and delaminations. It is shown that the different failure mechanisms observed in the tests are well predicted. In addition, the modeling approach is demonstrated to predict proper trends in the effect of scaling on strength and failure mechanisms of quasi-isotropic open-hole tension laminates.

Song, Kyonchan↗

Applicability of the Continuum-Shell Theories to the Mechanics of Carbon Nanotubes

Validity of the assumptions relating the applicability of continuum shell theories to the global mechanical behavior of carbon nanotubes is examined. The present study focuses on providing a basis that can be used to qualitatively assess the appropriateness of continuum-shell models for nanotubes. To address the effect of nanotube structure on their deformation, all nanotube geometries are divided into four major classes that require distinct models. Criteria for the applicability of continuum models are presented. The key parameters that control the buckling strains and deformation modes of these classes of nanotubes are determined. In an analogy with continuum mechanics, mechanical laws of geometric similitude are presented. A parametric map is constructed for a variety of nanotube geometries as a guide for the applicability of different models. The continuum assumptions made in representing a nanotube as a homogeneous thin shell are analyzed to identify possible limitations of applying shell theories and using their bifurcation-buckling equations at the nano-scale.

Harik, V. M.↗

Multiscale Fiber Kinking: Computational Micromechanics and a Mesoscale Continuum Damage Mechanics Models

In this work, the fiber kinking phenomenon, which is known as the failure mechanism that takes place when a fiber reinforced polymer is loaded under longitudinal compression, is studied. A computational micromechanics model is employed to interrogate the assumptions of a recently developed mesoscale continuum damage mechanics (CDM) model for fiber kinking based on the deformation gradient decomposition (DGD) and the LaRC04 failure criteria.

Herraez, Miguel↗

Analysis of Open Hole Compression Specimens Using the CompDam Continuum Damage Mechanics Model

This report documents a validation study on the standard Open Hole Compression (OHC) laminate test specimen conducted as part of the NASA Advanced Composites Project (ACP). Tests were conducted on OHC specimens with hard, quasi-isotropic, and soft layups using digital image correlation and X-Ray computed tomography to capture the structural response and damage evolution. Progressive damage models were constructed for use with the CompDam continuum damage mechanics code following the best practices established during the ACP. Detailed interrogation of the analysis results and comparison with experimental measurements provide a basis for assessing the capability of the modeling approach for OHC. The structural response is found to be captured well, with strengths predicted within 3\% of the experimental values for hard and quasi-isotropic laminates. In the soft laminate, the model predicts failure to be more brittle than the nonlinear, ductile response that was measured. Damage states extracted from the models at the same load level as test measurements are overlaid to show directly the similarities and differences between test and analysis results. Studying the damage evolution predicted by the analysis reveals that the failure process is a competition between fiber damage and delamination/sub-laminate buckling, with fiber damage dominating the collapse in the hard laminate and sub-laminate buckling governing in the soft laminate. Finally, a series of parametric studies varying numerical solution parameters (mesh size, mass scaling) and physical properties (fiber direction compressive strength and toughness) reveal sensitivities and deficiencies of the model. To the authors' knowledge, this study is the first for OHC specimens to include detailed evaluation of damage mode interactions, direct overlay of predicted and measured damage states, and sensitivity of the predicted results to difficult-to-measure fiber direction material properties.

Continuum damage mechanics↗