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 109 records · Page 6

Dynamics and stability of a tethered centrifuge in low earth orbit

The three-dimensional attitude dynamics of a spaceborne tethered centrifuge for artificial gravity experiments in low earth orbit is analyzed using two different methods. First, the tethered centrifuge is modeled as a dumbbell with a straight viscoelastic tether, point tip-masses, and sophisticated environmental models such as nonspherical gravity, thermal perturbations, and a dynamic atmospheric model. The motion of the centrifuge during spin-up, de-spin, and steady-rotation is then simulated. Second, a continuum model of the tether is developed for analyzing the stability of lateral tether oscillations. Results indicate that the maximum fluctuation about the 1-g radial acceleration level is less than 0.001 g; the time required for spin-up and de-spin is less than one orbit; and lateral oscillations are stable for any practical values of the system parameters.

Quadrelli, B. M.↗

Differential continuum damage mechanics models for creep and fatigue of unidirectional metal matrix composites

Three multiaxial isothermal continuum damage mechanics models for creep, fatigue, and creep/fatigue interaction of a unidirectional metal matrix composite volume element are presented, only one of which will be discussed in depth. Each model is phenomenological and stress based, with varying degrees of complexity to accurately predict the initiation and propagation of intergranular and transgranular defects over a wide range of loading conditions. The development of these models is founded on the definition of an initially transversely isotropic fatigue limit surface, static fracture surface, normalized stress amplitude function and isochronous creep damage failure surface, from which both fatigue and creep damage evolutionary laws can be obtained. The anisotropy of each model is defined through physically meaningful invariants reflecting the local stress and material orientation. All three transversely isotropic models have been shown, when taken to their isotropic limit, to directly simplify to previously developed and validated creep and fatigue continuum damage theories. Results of a nondimensional parametric study illustrate (1) the flexibility of the present formulation when attempting to characterize a large class of composite materials, and (2) its ability to predict anticipated qualitative trends in the fatigue behavior of unidirectional metal matrix composites. Additionally, the potential for the inclusion of various micromechanical effects (e.g., fiber/matrix bond strength, fiber volume fraction, etc.), into the phenomenological anisotropic parameters is noted, as well as a detailed discussion regarding the necessary exploratory and characterization experiments needed to utilize the featured damage theories.

Arnold, S. M.↗

Theoretical quasar emission-line ratios. VII - Energy-balance models for finite hydrogen slabs

The present energy balance calculations for finite, isobaric, hydrogen-slab quasar emission line clouds incorporate probabilistic radiative transfer (RT) in all lines and bound-free continua of a five-level continuum model hydrogen atom. Attention is given to the line ratios, line formation regions, level populations and model applicability results obtained. H lines and a variety of other considerations suggest the possibility of emission line cloud densities in excess of 10 to the 10th/cu cm. Lyman-beta/Lyman-alpha line ratios that are in agreement with observed values are obtained by the models. The observed Lyman/Balmer ratios can be achieved with clouds whose column depths are about 10 to the 22nd/sq cm.

Hubbard, E. N.↗

Finite Element Modeling of Plastic Deformation During Spin Forming of Aluminum 6061-O

Spin and flow forming are metal deformation techniques in which a disk or tube of material is radially thinned and axially lengthened over a rotating mandrel. A finite-element continuum model (FEM) based on commercial software (DEFORM®) was developed to determine the stress/strain distribution and damage accumulation in a flow-formed, near-net-shape part. Computational validation and experimental verification of the model was leveraged to assess the effectiveness of computer simulations. Quantitative comparisons between experimental and computational results show promise for applying a FEM to the flow forming process. Experimental failure locations correlated well with computed damage gradients, and there was good agreement between the measured and predicted roller forces during forming.

Elizabeth Urig↗

Application of advanced computational procedures for modeling solar-wind interactions with Venus: Theory and computer code

Computational procedures are developed and applied to the prediction of solar wind interaction with nonmagnetic terrestrial planet atmospheres, with particular emphasis to Venus. The theoretical method is based on a single fluid, steady, dissipationless, magnetohydrodynamic continuum model, and is appropriate for the calculation of axisymmetric, supersonic, super-Alfvenic solar wind flow past terrestrial planets. The procedures, which consist of finite difference codes to determine the gasdynamic properties and a variety of special purpose codes to determine the frozen magnetic field, streamlines, contours, plots, etc. of the flow, are organized into one computational program. Theoretical results based upon these procedures are reported for a wide variety of solar wind conditions and ionopause obstacle shapes. Plasma and magnetic field comparisons in the ionosheath are also provided with actual spacecraft data obtained by the Pioneer Venus Orbiter.

Stahara, S. S.↗

The Long-stable Hard State of XTE J1752-223 and the Disk Truncation Dilemma

The degree to which the thin accretion disks of black hole X-ray binaries are truncated during hard spectral states remains a contentious open question in black hole astrophysics. During its singular observed outburst in 2009–2010, the black hole X-ray binary XTE J1752−223 spent ∼1 month in a long-stable hard spectral state at a luminosity of ∼0.02–0.1 L Edd . It was observed with 56 RXTE pointings during this period, with simultaneous Swift-XRT daily coverage during the first 10 days of the RXTE observations. While reflection modeling has been extensively explored in the analysis of these data, there is disagreement surrounding the geometry of the accretion disk and corona implied by the reflection features. We reexamine the combined, high signal-to-noise, simultaneous Swift and RXTE observations, and perform extensive reflection modeling with the latest relxill suite of reflection models, including newer high disk density models. We show that reflection modeling requires that the disk be within ∼5 R ISCO during the hard spectral state, while weaker constraints from the thermal disk emission imply higher truncation (Rin = 6–80 R ISCO ). We also explore more complex coronal continuum models, allowing for two Comptonization components instead of one, and show that the reflection features still require only a mildly truncated disk. Finally we present a full comparison of our results to previous constraints found from analyses of the same data set.

Accretion↗

Calculation of solar wind flows about terrestrial planets

A computational model was developed for the determination of the plasma and magnetic field properties of the global interaction of the solar wind with terrestrial planetary magneto/ionospheres. The theoretical method is based on an established single fluid, steady, dissipationless, magnetohydrodynamic continuum model, and is appropriate for the calculation of supersonic, super Alfvenic solar wind flow past terrestrial planets. A summary is provided of the important research results.

Stahara, S. S.↗

A Theoretical Investigation of the Plausibility of Reactions Between Ammonia and Carbonyl Species (Formaldehyde, Acetaldehyde, and Acetone) in Interstellar Ice Analogs at Ultracold Temperatures

We have reexamined the reaction between formaldehyde and ammonia, which was previously studied by us and other workers in modestly sized cluster calculations. Larger model systems with up to 12H2O were employed, and reactions of two more carbonyl species, acetaldehyde and acetone, were also carried out. Calculations were performed at the B3LYP/6-31+G** level with bulk solvent effects treated with a polarizable continuum model; limited MP2/6-31+G** calculations were also performed. We found that while the barrier for the concerted proton relay mechanism described in previous work remains modest, it is still prohibitively high for the reaction to occur under the ultracold conditions that prevail in dense interstellar clouds. However, a new pathway emerged in more realistic clusters that involves at least one barrierless step for two of the carbonyl species considered here: ammonia reacts with formaldehyde and acetaldehyde to form a partial charge transfer species in small clusters (4H2O) and a protonated hydroxyamino intermediate species in large clusters (9H2O, 12H2O); modest barriers that decrease sharply with cluster size are found for the analogous processes for the acetone-NH3 reaction. Furthermore, if a second ammonia replaces one of the water molecules in calculations in the 9H2O clusters, deprotonation can occur to yield the same neutral hydroxyamino species that is formed via the original concerted proton relay mechanism. In at least one position, deprotonation is barrierless when zero-point energy is included. In addition to describing the structures and energetics of the reactions between formaldehyde, acetaldehyde, and acetone with ammonia, we report spectroscopic predictions of the observable vibrational features that are expected to be present in ice mixtures of different composition.

Chen, Lina↗

Development of a computational model for predicting solar wind flows past nonmagnetic terrestrial planets

A computational model for the determination of the detailed plasma and magnetic field properties of the global interaction of the solar wind with nonmagnetic terrestrial planetary obstacles is described. The theoretical method is based on an established single fluid, steady, dissipationless, magnetohydrodynamic continuum model, and is appropriate for the calculation of supersonic, super-Alfvenic solar wind flow past terrestrial ionospheres.

Stahara, S. S.↗

Creating a Simple Single Computational Approach to Modeling Rarefied and Continuum Flow About Aerospace Vehicles

We proposed to create a single computational code incorporating methods that can model both rarefied and continuum flow to enable the efficient simulation of flow about space craft and high altitude hypersonic aerospace vehicles. The code was to use a single grid structure that permits a smooth transition between the continuum and rarefied portions of the flow. Developing an appropriate computational boundary between the two regions represented a major challenge. The primary approach chosen involves coupling a four-speed Lattice Boltzmann model for the continuum flow with the DSMC method in the rarefied regime. We also explored the possibility of using a standard finite difference Navier Stokes solver for the continuum flow. With the resulting code we will ultimately investigate three-dimensional plume impingement effects, a subject of critical importance to NASA and related to the work of Drs. Forrest Lumpkin, Steve Fitzgerald and Jay Le Beau at Johnson Space Center. Below is a brief background on the project and a summary of the results as of the end of the grant.

Goldstein, David B.↗

Application of a global solar wind/planetary obstacle interaction computational model: Earth, Venus, Mars, Jupiter and Saturn studies

This work relates to a series of collaborative investigations involving the application of a computational model for the determination of the detailed plasma and magnetic field properties associated with the global interaction of the solar wind with various planetary obstacles throughout the solar system. The theoretical method is based on an established single fluid, steady, dissipationless, magnetohydrodynamic continuum model, and is appropriate for the calculation of supersonic, super-Alfvenic solar wind flow past planetary obstacles. The investigations undertaken relate to studies of various solar wind interaction phenomena with Venus, Earth, Jupiter, and Saturn.

Stahara, S. S.↗

Solar wind flow past Venus - Theory and comparisons

Advanced computational procedures are applied to an improved model of solar wind flow past Venus to calculate the locations of the ionopause and bow wave and the properties of the flowing ionosheath plasma in the intervening region. The theoretical method is based on a single-fluid, steady, dissipationless, magneto-hydrodynamic continuum model and is appropriate for the calculation of axisymmetric supersonic, super-Alfvenic solar wind flow past a nonmagnetic planet possessing a sufficiently dense ionosphere to stand off the flowing plasma above the subsolar point and elsewhere. Determination of time histories of plasma and magnetic field properties along an arbitrary spacecraft trajectory and provision for an arbitrary oncoming direction of the interplanetary solar wind have been incorporated in the model. An outline is provided of the underlying theory and computational procedures, and sample comparisons of the results are presented with observations from the Pioneer Venus orbiter.

Spreiter, J. R.↗

Multi-Scale Modeling of Cross-Linked Nanotube Materials

The effect of cross-linking single-walled carbon nanotubes on the Young's modulus of a nanotube-reinforced composite is modeled with a multi-scale method. The Young's modulus is predicted as a function of nanotube volume fraction and cross-link density. In this method, the constitutive properties of molecular representative volume elements are determined using molecular dynamics simulation and equivalent-continuum modeling. The Young's modulus is subsequently calculated for cross-linked nanotubes in a matrix which consists of the unreacted cross-linking agent. Two different cross-linking agents are used in this study, one that is short and rigid (Molecule A), and one that is long and flexible (Molecule B). Direct comparisons between the predicted elastic constants are made for the models in which the nanotubes are either covalently bonded or not chemically bonded to the cross-linking agent. At a nanotube volume fraction of 10%, the Young's modulus of Material A is not affected by nanotube crosslinking, while the Young's modulus of Material B is reduced by 64% when the nanotubes are cross-linked relative to the non-cross-linked material with the same matrix.

Frankland, S. J. V.↗

Multi-Scale Modeling of a Graphite-Epoxy-Nanotube System

A multi-scale method is utilized to determine some of the constitutive properties of a three component graphite-epoxy-nanotube system. This system is of interest because carbon nanotubes have been proposed as stiffening and toughening agents in the interlaminar regions of carbon fiber/epoxy laminates. The multi-scale method uses molecular dynamics simulation and equivalent-continuum modeling to compute three of the elastic constants of the graphite-epoxy-nanotube system: C11, C22, and C33. The 1-direction is along the nanotube axis, and the graphene sheets lie in the 1-2 plane. It was found that the C11 is only 4% larger than the C22. The nanotube therefore does have a small, but positive effect on the constitutive properties in the interlaminar region.

Frankland, S. J. V.↗

Mathematical modeling of spinning elastic bodies for modal analysis.

The problem of modal analysis of an elastic appendage on a rotating base is examined to establish the relative advantages of various mathematical models of elastic structures and to extract general inferences concerning the magnitude and character of the influence of spin on the natural frequencies and mode shapes of rotating structures. In realization of the first objective, it is concluded that except for a small class of very special cases the elastic continuum model is devoid of useful results, while for constant nominal spin rate the distributed-mass finite-element model is quite generally tractable, since in the latter case the governing equations are always linear, constant-coefficient, ordinary differential equations. Although with both of these alternatives the details of the formulation generally obscure the essence of the problem and permit very little engineering insight to be gained without extensive computation, this difficulty is not encountered when dealing with simple concentrated mass models.

Likins, P. W.↗

Particle kinetic simulation of high altitude hypervelocity flight

In this grant period, the focus has been on enhancement and application of the direct simulation Monte Carlo (DSMC) particle method for computing hypersonic flows of re-entry vehicles. Enhancement efforts dealt with modeling gas-gas interactions for thermal non-equilibrium relaxation processes and gas-surface interactions for prediction of vehicle surface temperatures. Both are important for application to problems of engineering interest. The code was employed in a parametric study to improve future applications, and in simulations of aeropass maneuvers in support of the Magellan mission. Detailed comparisons between continuum models for internal energy relaxation and DSMC models reveals that several discrepancies exist. These include definitions of relaxation parameters and the methodologies for implementing them in DSMC codes. These issues were clarified and all differences were rectified in a paper (Appendix A) submitted to Physics of Fluids A, featuring several key figures in the DSMC community as co-authors and B. Haas as first author. This material will be presented at the Fluid Dynamics meeting of the American Physical Society on November 21, 1993. The aerodynamics of space vehicles in highly rarefied flows are very sensitive to the vehicle surface temperatures. Rather than require prescribed temperature estimates for spacecraft as is typically done in DSMC methods, a new technique was developed which couples the dynamic surface heat transfer characteristics into the DSMC flow simulation code to compute surface temperatures directly. This model, when applied to thin planar bodies such as solar panels, was described in AIAA Paper No. 93-2765 (Appendix B) and was presented at the Thermophysics Conference in July 1993. The paper has been submitted to the Journal of Thermophysics and Heat Transfer. Application of the DSMC method to problems of practical interest requires a trade off between solution accuracy and computational expense and limitations. A parametric study was performed and reported in AIAA Paper No. 93-2806 (Appendix C) which assessed the accuracy penalties associated with simulations of varying grid resolution and flow domain size. The paper was also presented at the Thermophysics Conference and will be submitted to the journal shortly. Finally, the DSMC code was employed to assess the pitch, yaw, and roll aerodynamics of the Magellan spacecraft during entry into the Venus atmosphere at off-design attitudes. This work was in support of the Magellan aerobraking maneuver of May 25-Aug. 3, 1993. Furthermore, analysis of the roll characteristics of the configuration with canted solar panels was performed in support of the proposed 'Windmill' experiment. Results were reported in AIAA Paper No. 93-3676 (Appendix D) presented at the Atmospheric Flight Mechanics Conference in August 1993, and were submitted to Journal of Spacecraft and Rockets.

Haas, Brian L.↗

The New BeppoSAX Observation of the Brightest X-Ray Quasar at Redshift

This grant was to support the reduction and analysis of our approved SAX observation of the high redshift (z=3.2) blazar PKS 2126-158. This is the brightest quasar at z greater than 3 and has been intensively studied in X-ray, since the first Einstein detection. In 1994 Elvis et al., discovered a strong low energy cutoff in this object, which could imply either quasar frame photoelectric absorption by a column of 0.8-2.7 x 1e22 cm-2 cold gas, or a lower column of cold gas at z=0. Subsequent ASCA observations of this object, could not definitely address this issue, nor could establish whether the curvature of the low energy portion of the spectrum was due to pure photoelectric absorption (considerably exceeding the Galactic value along the line of sight) or to an intrinsic continuum curvature. We proposed to observe PKS 2126-158 with BeppoSAX, to try to solve this puzzle (thanks to the broadband of BeppoSAX: 0.1-250 keV). PKS 2126 was observed by BeppoSAX on May 1999, with a MECS exposure of 100 ks. We have reduced and analyzed the BeppoSAX data, and compared them with a Chandra ACIS observation of the same object, taken only 6 months apart (Nov. 1999). We have recently finished to write a paper on the BeppoSAX data only, that concentrate on the properties of the X-ray absorber, which is highly requested by our SAX data, independently on the continuum model adopted. The paper (P.I.F. Fiore) will be submitted to APJ in the next few days. A second paper on the combined BeppoSAX and Chandra data, and based on the broad band spectral energy distribution of this quasar, is currently in preparation. Our main results, on the X-ray absorber, are: (a) the presence of an X-ray absorber is confirmed, indipendently on the continuum adopted (simple power law, or curved continuum); (b) if the absorber is not significantly ionized, then the BeppoSAX data do prefer a low redshift absorber; (c) if the gas is ionized, then it can be located in the quasar environment, but its metal abundances must be lower than 0.2 times solar (because of the absence of a strong FeK absorption edges that would follow at approx. 2 keV, where the MECS response is the highest).

Nicastro, Fabrizio↗

An Assessment of Five Modeling Approaches for Thermo-Mechanical Stress Analysis of Laminated Composite Panels

A study is made of the effects of variation in the lamination and geometric parameters, and boundary conditions of multi-layered composite panels on the accuracy of the detailed response characteristics obtained by five different modeling approaches. The modeling approaches considered include four two-dimensional models, each with five parameters to characterize the deformation in the thickness direction, and a predictor-corrector approach with twelve displacement parameters. The two-dimensional models are first-order shear deformation theory, third-order theory; a theory based on trigonometric variation of the transverse shear stresses through the thickness, and a discrete layer theory. The combination of the following four key elements distinguishes the present study from previous studies reported in the literature: (1) the standard of comparison is taken to be the solutions obtained by using three-dimensional continuum models for each of the individual layers; (2) both mechanical and thermal loadings are considered; (3) boundary conditions other than simply supported edges are considered; and (4) quantities compared include detailed through-the-thickness distributions of transverse shear and transverse normal stresses. Based on the numerical studies conducted, the predictor-corrector approach appears to be the most effective technique for obtaining accurate transverse stresses, and for thermal loading, none of the two-dimensional models is adequate for calculating transverse normal stresses, even when used in conjunction with three-dimensional equilibrium equations.

Noor, A. K.↗