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 379 records · Page 21

On equivalence of discrete-discrete and continuum-discrete design sensitivity analysis

Developments in design sensitivity analysis (DSA) method have been made using two fundamentally different approaches as shown. In the first approach, a discretized structural finite element model is used to carry out DSA. There are three different methods in the discrete DSA approach: finite difference, semi-analytical, and analytical methods. The finite difference method is a popular one due to its simplicity, but a serious shortcoming of the method is the uncertainty in the choice of a perturbation step size of design variables. In the semi-analytical method, the derivatives of stiffness matrix is computed by finite differences, whereas in the analytical method, the derivatives are obtained analytically. For the shape design variable, computation of analytical derivative of stiffness matrix is quite costly. Because of this, the semi-analytical method is a popular choice in discrete shape DSA approach. However, recently, Barthelemy and Haftka presented that the semi-analytical method can have serious accuracy problems for shape design variables in structures modeled by beam, plate, truss, frame, and solid elements. They found that accuracy problems occur even for a simple cantilever beam. In the second approach, a continuum model of the structure is used to carry out DSA.

Choi, Kyung K.↗

Multiscale Modeling of Damage Processes in fcc Aluminum: From Atoms to Grains

Molecular dynamics (MD) methods are opening new opportunities for simulating the fundamental processes of material behavior at the atomistic level. However, current analysis is limited to small domains and increasing the size of the MD domain quickly presents intractable computational demands. A preferred approach to surmount this computational limitation has been to combine continuum mechanics-based modeling procedures, such as the finite element method (FEM), with MD analyses thereby reducing the region of atomic scale refinement. Such multiscale modeling strategies can be divided into two broad classifications: concurrent multiscale methods that directly incorporate an atomistic domain within a continuum domain and sequential multiscale methods that extract an averaged response from the atomistic simulation for later use as a constitutive model in a continuum analysis.

Glaessgen, E. H.↗

A finite-difference program for stresses in anisotropic, layered plates in bending

The interlaminar stresses induced in a layered laminate that is bent into a cylindrical surface are studied. The laminate is modeled as a continuum, and the resulting elasticity equations are solved using the finite difference method. The report sets forth the mathematical framework, presents some preliminary results, and provides a listing and explanation of the computer program. Significant among the results are apparent symmetry relationships that will reduce the numerical size of certain problems and an interlaminar stress behavior having a sharp rise at the free edges.

Salamon, N. J.↗

The Mg II h and k lines. II - Comparison with synthesized profiles and Ca II K

Measured high-dispersion center and limb profiles of the solar Mg II h and k resonance lines are compared with synthetic spectra computed with a partial-redistribution formalism and based on several upper-photosphere and lower-chromosphere temperature distributions. Profiles of the analogously formed Ca II K resonance line are also synthesized for the same atmospheric models. The spectrum-synthesis approach is outlined, and the collisional and fixed radiative rates appropriate to the adopted model atoms and solar atmosphere are discussed. It is found that the HSRA and VAL models predict systematically lower intensities in the h, k, and K inner wings than observed and that models with a somewhat higher minimum temperature (about 4450 K) can reproduce the measured inner wings and limb darkening. A 'Ca II' solar model with a minimum temperature of 4450 K is proposed as an alternative to the class of models based on continuum observations.

Ayres, T. R.↗

Simulation of solidification in a Bridgman cell

Bridgman-type crystal growth techniques are attractive methods for producing homogeneous, high-quality infrared detector and junction device materials. However, crystal imperfections and interface shapes still must be controlled through modification of the temperature and concentration gradients created during solidification. The objective of this investigation was to study the temperature fields generated by various cell and heatpipe configurations and operating conditions. Continuum's numerical model of the temperature, species concentrations, and velocity fields was used to describe the thermal characteristics of Bridgman cell operation.

Dakhoul, Y. M.↗

Subcritical Fracture Propagation in Rocks: an Examination Using the Methods of Fracture Mechanics and Non-Destructive Testing

An experimental investigation of tensile rock fracture is presented with an emphasis on characterizing time dependent crack growth using the methods of fracture mechanics. Subcritical fracture experiments were performed in moist air on glass and five different rock types at crack velocities using the double torsion technique. The experimental results suggest that subcritical fracture resistance in polycrystals is dominated by microstructural effects. Evidence for gross violations of the assumptions of linear elastic fracture mechanics and double torsion theory was found in the tests on rocks. In an effort to obtain a better understanding of the physical breakdown processes associated with rock fracture, a series of nondestructive evaluation tests were performed during subcritical fracture experiments on glass and granite. Comparison of the observed process zone shape with that expected on the basis of a critical normal principal tensile stress criterion shows that the zone is much more elongated in the crack propagation direction than predicted by the continuum based microcracking model alone.

Swanson, P. L.↗

Transport properties of droplet clusters in gravity-free fields

Clusters of liquid droplets are suspended in an atmosphere of saturated vapor and are subjected to an external force field. This system can be modeled as a continuum whose macroscopic properties may be determined by applying the generalized theory of Taylor dispersion.

Brenner, Howard↗

Hypersonic rarefied-flow aerodynamics inferred from Shuttle Orbiter acceleration measurements

Data obtained from multiple flights of sensitive accelerometers on the Space Shuttle Orbiter during reentry have been used to develop an improved aerodynamic model for the Orbiter normal- and axial-force coefficients in hypersonic rarefied flow. The lack of simultaneous atmospheric density measurements was overcome in part by using the ratio of normal-to-axial acceleration, in which density cancels, as a constraint. Differences between the preflight model and the flight-acceleration-derived model in the continuum regime are attributed primarily to real gas effects. New insights are gained into the variation of the force coefficients in the transition between the continuum regime and free molecule flow.

Blanchard, R. C.↗

H-beta photometry of V1010 Ophiuchi

An analysis is presented of previously unpublished intermediate and narrow band H-beta light curves of the short period eclipsing binary V1010 Oph. The temperature of the primary is determined directly by comparing the UV continuum with Kurucz model atmospheres and solve the broadband lightcurve using the Wilson-Devinney program. The present solution has the primary filling its Roche lobe while the secondary is detached from its lobe. An asymmetry in the light curve in which first quadrature (0.25 P) is brighter than the second quadrature (0.75 P) may be the result of local heating of the inner surface of the secondary by mass transfer from the primary. The present semidetached solution indicates that the system is probably evolved.

Corcoran, M. F.↗

Inelastic Deformation of Metal Matrix Composites: Plasticity and Damage Mechanisms - Part 1

The deformation mechanisms of a Ti 15-3/SCS6 (SiC fiber) metal matrix composite (MMC) were investigated using a combination of mechanical measurements and microstructural analysis. The objectives were to evaluate the contributions of plasticity and damage to the overall inelastic response, and to confirm the mechanisms by rigorous microstructural evaluations. The results of room temperature experiments performed on 0 degree and 90 degree systems primarily are reported in this report. Results of experiments performed on other laminate systems and at high temperatures will be provided in a forthcoming report. Inelastic deformation of the 0 degree MMC (fibers parallel to load direction) was dominated by the plasticity of the matrix. In contrast, inelastic deformations of the 90 degree composite (fibers perpendicular to loading direction) occurred by both damage and plasticity. The predictions of a continuum elastic plastic model were compared with experimental data. The model was adequate for predicting the 0 degree response; however, it was inadequate for predicting the 90 degree response largely because it neglected damage. The importance of validating constitutive models using a combination of mechanical measurements and microstructural analysis is pointed out. The deformation mechanisms, and the likely sequence of events associated with the inelastic deformation of MMCs, are indicated in this paper.

Majumdar, B. S.↗

An explicit solution to the optimal LQG problem for flexible structures with collocated rate sensors

We present a class of compensators in explicit form (not requiring numerical computer calculations) for stabilizing flexible structures with collocated rate sensors. They are based on the explicit solution, valid for both Continuum and FEM Models, of the LQG problem for minimizing mean square rate. They are robust with respect to system stability (will not destabilize modes even with mismatch of parameters), can be instrumented in state space form suitable for digital controllers, and can be specified directly from the structure modes and mode 'signature' (displacement vectors at sensor locations). Some simulation results are presented for the NASA LaRC Phase-Zero Evolutionary Model - a modal Trust model with 86 modes - showing damping ratios attainable as a function of compensator design parameters and complexity.

Balakrishnan, A. V.↗

Similarity and scale in catchment storm response

Until recently, very little progress had been made in understanding the relationship between small-scale variability of topography, soil, and rainfalls and the storm response seen at the catchment scale. The work reviewed here represents the first attempt at a systematic theoretical framework for such understanding in the context of surface runoff generation by different processes. The parameterization of hydrological processes over a range of scales is examined, and the concept of the 'representative elementary area' (REA) is introduced. The REA is a fundamental scale for catchment modeling at which continuum assumptions can be applied for the spatially variable controls and parameters, and spatial patterns no longer have to be considered explicitly. The investigation of scale leads into the concept of hydrologic similarity in which the effects of the environmental controls on runoff generation and flood frequency response be investigated independently of catchment scale. The paper reviews the authors' initial results and hopefully will motivate others to also investigate the issues of hydrologic scale and similarity.

Wood, Eric F.↗

Evidence for elemental variation in the ejecta of the Tycho supernova remnant

We present an X-ray study of the Tycho supernova remnant utilizing archival data from the high-resolution imagers (HRIs) on Einstein and ROSAT, the low-energy imaging telescopes (LEITs) on EXOSAT, and spectral data from the Broad Band X-ray Telescope (BBXRT). We have made use of the differing HRI bandpasses to construct images of Tycho in two spectral bands, o.7-1.8 keV and 1.8-4.5 keV. We find that the two images differ, with the harder image showing enhanced emission along much of the south, west, and north periphery. There appears to be enhanced soft emission in the interior and in one particular knot of emission in the southeast. Besides continuum (which we model here as thermal bremsstrahlung emission), we believe the hard image shows primarily the distribution of high-ionization Si and S K-shell lines which lie in the 1.8-2.6 keV band, while the softer image has contributions from Si as well as Fe XVII to Fe XXIV L-shell lines in the 0.7-1.4 keV band. Guided by the results of nonequilibrium ionization modeling of the BBXRT spectral data, we interpret the observed contrast in hard and soft X-ray emission in terms of variations in abundance, ionization timescale, and temperature. The most likely explanation for the spectral differences are spatial variations of the relative abundances of Si, S, and Fe.

Vancura, O.↗

XTE Proposal #20102--"SS 433's High Energy Spectrum"

We observed the jet-producing compact binary system SS 433 with RXTE during three multiwavelength campaigns, the first in conjunction with ASCA observations, the second simultaneous with a VLA-VLBA-MERLIN campaign, and the third associated with a Nobeyama millimeter-band campaign. All these campaigns included optical observations. Occurring at different jet precession and binary phases, the observations also monitored the system during a radio flare. The data provide SS 433's X-ray spectrum over more than an energy decade, and track the spectral variations as the X-ray source was partially eclipsed. The continuum can be modeled as a power law with an exponential cutoff, which can be detected to approximately 50 keV. Strong line emission is evident in the 5-10 keV range which can be modeled as a broad line whose energy is precession independent and a narrow line whose energy does vary with jet precession phase; this line model is clearly an over simplification since the PCA does not have sufficient energy resolution to detect the lines ASCA observed. The eclipses are deeper at high energy and at jet precession phases when the jets are more inclined towards and away from us. A large radio flare occurred between two sets of X-ray monitoring observations; an X-ray observation at the peak of the flare found a softer spectrum with a flux approximately 1/3 that of the quiescent level.

Band, David L.↗

An Intercomparison of Radiation Codes for Retrieving Upper Tropospheric Humidity in the 6.3-micron Band: A Report from the 1st GVaP Workshop

An intercomparison of radiation codes used in retrieving upper tropospheric humidity (UTH) from observations in the v2 (6.3 microns) water vapor absorption band was performed. This intercomparison is one part of a coordinated effort within the GEWEX Water Vapor Project (GVaP) to assess our ability to monitor the distribution and variations of upper tropospheric moisture from space-borne sensors. A total of 23 different codes, ranging from detailed line-by-line (LBL) models, to coarser resolution narrow-band (NB) models, to highly-parameterized single-band (SB) models participated in the study. Forward calculations were performed using a carefully selected set of temperature and moisture profiles chosen to be representative of a wide range of atmospheric conditions. The LBL model calculations exhibited the greatest consistency with each other, typically agreeing to within 0.5 K in terms of the equivalent blackbody brightness temperature (T(sub b)). The majority of NB and SB models agreed to within +/- 1 K of the LBL models, although a few older models exhibited systematic T(sub b) biases in excess of 2 K. A discussion of the discrepancies between various models, their association with differences in model physics (e.g. continuum absorption), and their implications for UTH retrieval and radiance assimilation is presented.

Soden, B.↗

Concept for Determining the Life of Ceramic Matrix Composites Using Nondestructive Characterization Techniques

Damping measurements have been taken on ceramic matrix composite (CMC) turbopump blisks in the as fabricated, post proof testing, and post turbopump testing conditions. These results indicate that damping is able to quantify fatigue of the CMC blisk. This gives hope for the potential of determining the actual and residual life of CMC materials using a combination of nondestructive techniques. If successful, then this new paradigm for life prediction of CMCs could revolutionize the approach for designing and servicing CMC components, thereby significantly reducing costs for design, development, health monitoring, and maintenance of CMC components and systems. The Nondestructive Characterization (NDC) life prediction approach would complement life prediction using micromechanics and continuum finite element models. This paper reports on the initial concept of NDC life prediction and how changes in damping and ultrasonic elastic modulus data have established the concept as a possibility.

Effinger, Michael↗