Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computing continuum”

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

Comprehensive Micromechanics-Analysis Code - Version 4.0

Version 4.0 of the Micromechanics Analysis Code With Generalized Method of Cells (MAC/GMC) has been developed as an improved means of computational simulation of advanced composite materials. The previous version of MAC/GMC was described in "Comprehensive Micromechanics-Analysis Code" (LEW-16870), NASA Tech Briefs, Vol. 24, No. 6 (June 2000), page 38. To recapitulate: MAC/GMC is a computer program that predicts the elastic and inelastic thermomechanical responses of continuous and discontinuous composite materials with arbitrary internal microstructures and reinforcement shapes. The predictive capability of MAC/GMC rests on a model known as the generalized method of cells (GMC) - a continuum-based model of micromechanics that provides closed-form expressions for the macroscopic response of a composite material in terms of the properties, sizes, shapes, and responses of the individual constituents or phases that make up the material. Enhancements in version 4.0 include a capability for modeling thermomechanically and electromagnetically coupled ("smart") materials; a more-accurate (high-fidelity) version of the GMC; a capability to simulate discontinuous plies within a laminate; additional constitutive models of materials; expanded yield-surface-analysis capabilities; and expanded failure-analysis and life-prediction capabilities on both the microscopic and macroscopic scales.

Arnold, S. M.↗

Some experience with arc-heater simulation of outer planet entry radiation

An electric arc heater was operated at 800 amperes and 100,000 pa (1 atm) with hydrogen, helium, and two mixtures of hydrogen and helium. A VUV-scanning monochromator was used to record the spectra from an end view while a second spectrometer was used to determine the plasma temperature using hydrogen continuum radiation at 562 nm. Except for pure helium, the plasma temperature was found to be too low to produce significant helium radiation, and the measured spectra were primarily the hydrogen spectra with the highest intensity in the pure hydrogen case. A radiation computer code was used to compute the spectra for comparison to the measurements and to extend the study to simulation of outer planet entry radiation. Conductive cooling prevented ablation of phenolic carbon material samples mounted inside the arc heater during a cursory attempt to produce radiation absorption by ablation gases.

Wells, W. L.↗

Formation of the O I resonance triplet and intercombination doublet in the solar chromosphere

Spectrum synthesis calculations are presented for the O I resonance triplet at 1304 A and the intercombination doublet at 1358 A for the solar atmosphere and several variants, allowing for triplet fluorescence by the Ly-beta emission of H I. Profiles, synthesized from a seven-level plus continuum O I atom are compared to observations taken with the high-resolution spectrometer on OSO 8. It is found that the O I triplet emission is dominated by the Ly-beta fluorescence and that the agreement between observations and profiles computed with current chromospheric models is much improved over earlier studies.

Skelton, D. L.↗

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.↗

Turbomachinery

The discipline research in turbomachinery, which is directed toward building the tools needed to understand such a complex flow phenomenon, is based on the fact that flow in turbomachinery is fundamentally unsteady or time dependent. Success in building a reliable inventory of analytic and experimental tools will depend on how the time and time-averages are treated, as well as on who the space and space-averages are treated. The raw tools at disposal (both experimentally and computational) are truly powerful and their numbers are growing at a staggering pace. As a result of this power, a case can be made that a situation exists where information is outstripping understanding. The challenge is to develop a set of computational and experimental tools which genuinely increase understanding of the fluid flow and heat transfer in a turbomachine. Viewgraphs outline a philosophy based on working on a stairstep hierarchy of mathematical and experimental complexity to build a system of tools, which enable one to aggressively design the turbomachinery of the next century. Examples of the types of computational and experimental tools under current development at Lewis, with progress to date, are examined. The examples include work in both the time-resolved and time-averaged domains. Finally, an attempt is made to identify the proper place for Lewis in this continuum of research.

Simoneau, Robert J.↗

Accuracy of the Burnett equations for hypersonic real gas flows

Burnett solutions for the structure of normal shock waves for Mach numbers ranging from 1.1 to 50 have been obtained by using the time-dependent equation set in conjunction with modern computational fluid dynamic techniques. The current work seeks to provide corrected and improved high-resolution direct simulation Monte Carlo (DSMC) results for hypersonic shocks in monatomic gases. These DSMC results will serve as a reliable test base for assessing the accuracy of various continuum equations more advanced than Navier-Stokes, such as the Burnett equations. The main conclusion of Fisko and Chapman (1988) that the Burnett equations are more accurate than the Navier-Stokes equations for all gases and at all Mach numbers in matching DSMC shock wave results is reassesed by comparing the improved and corrected DSMC results with continuum Burnett results. Various Burnett results for shocks in a diatomic gas, nitrogen, which include improvements to the work of Lumpkin and Chapman (1989) are presented. It is concluded that the Burnett equations represent a significant improvement over Navier-Stokes equations for flows exhibiting translational nonequilibrium.

Lumpkin, Forrest E., III↗

Multilevel optimization using a continuum model for structures

The new concept of using a continuum model in the multilevel decomposition of structural optimization problems is presented. The practicality of this concept is demonstrated by application in a scheme for the optimization of beam-like space structures. This scheme is tested against traditional optimization procedures for savings in computational cost. Results from both optimization methods are presented for comparison.

Yates, K.↗

Probabilistic finite elements

In the Probabilistic Finite Element Method (PFEM), finite element methods have been efficiently combined with second-order perturbation techniques to provide an effective method for informing the designer of the range of response which is likely in a given problem. The designer must provide as input the statistical character of the input variables, such as yield strength, load magnitude, and Young's modulus, by specifying their mean values and their variances. The output then consists of the mean response and the variance in the response. Thus the designer is given a much broader picture of the predicted performance than with simply a single response curve. These methods are applicable to a wide class of problems, provided that the scale of randomness is not too large and the probabilistic density functions possess decaying tails. By incorporating the computational techniques we have developed in the past 3 years for efficiency, the probabilistic finite element methods are capable of handling large systems with many sources of uncertainties. Sample results for an elastic-plastic ten-bar structure and an elastic-plastic plane continuum with a circular hole subject to cyclic loadings with the yield stress on the random field are given.

Belytschko, Ted↗

The structure and spectrum of a colliding-cloud system and its possible relationship to QSOs

A collision between two gas clouds with initial densities of approximately 10 million per cu cm, velocities of about 1000 km/s, and radii of approximately 1 pc is investigated quantitatively by coupling a calculation of the radiation spectrum resulting from the anticipated shock fronts with a computation for the conversion of this high-energy radiation into optical emission in adjacent photoionized regions. The detailed structure of the colliding clouds is discussed, and the effects of an ambient magnetic field are considered. The combined emission-line spectrum is presented along with continuum emission estimates for thermal, synchrotron, and very-high-energy bremsstrahlung mechanisms. It is shown that significant continua can be produced over the range from 300 microns to 3 keV, including a blackbody contribution from a high-density neutral region between the shock fronts, free-free and free-bound radiation from the cooling zones directly behind the shocks, and free-free, free-bound, and two-photon radiation from the photoionized regions immediately ahead of and behind the cooling zones. The theoretical spectrum of the structure resulting from the collision is found to be similar in general and in some details to those observed for typical quasars.

Daltabuit, E.↗

Atomization simulations using an Eulerian-VOF-Lagrangian method

This paper summarizes the technical development and validation of a multiphase computational fluid dynamics (CFD) numerical method using the volume-of-fluid (VOF) model and a Lagrangian tracking model which can be employed to analyze general multiphase flow problems with free surface mechanism. The gas-liquid interface mass, momentum and energy conservations are modeled by continuum surface mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed flow regimes. The objectives of the present study are to develop and verify the fractional volume-of-fluid cell partitioning approach into a predictor-corrector algorithm and to demonstrate the effectiveness of the present innovative approach by simulating benchmark problems including the coaxial jet atomization.

Chen, Yen-Sen↗

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↗

Limb-darkening scans of Jupiter

Representative center-limb photometric measurements of several areas on Jupiter are presented in the form of limb-darkening scans taken through narrow-band interference filters which isolate continuum regions near the 3-0 and 4-0 hydrogen quadrupole bands (8200 and 6300 A, respectively). Results of scans measured through filters which isolate two regions dominated by methane absorption (6190 and 7250 A) are also discussed. Limb-darkening curves at all four wavelengths are computed for semiinfinite homogeneous isotropically scattering atmospheres with different values of the single-scattering albedo of the cloud particles. It is found that albedo values between 0.985 and 0.989 match the data over most of the northern component of the South Equatorial Belt. The single-scattering albedo values are used to show that the equatorial regions of Jupiter which emit at 5 microns may be relatively clear while the tropical regions may be relatively cloudy.

Pilcher, C. B.↗

Decoupled predictions of radiative heating in air using a particle simulation method

A particle simulation technique with decoupled radiation is used to estimate the radiative emission along the stagnation streamline and the radiative heating at the stagnation point of a blunt-nosed vehicle entering the earth's atmosphere at hypersonic speed. The direct simulation Monte Carlo (DSMC) method is used to compute the fluid mechanics of the weakly ionized flow. The radiative emission is computed using the NEQAIR computer code with the radiation decoupled from the flowfield solution. It is concluded that the new models make it possible to decrease the predicted total radiative heating at the stagnation point of the vehicle by a factor of 15. A comparison of the DSMC approach with a continuum flow model shows that the total predicted radiative heating estimates agree within a factor of 2.

Boyd, Iain D.↗

Numerical modeling of spray combustion with an advanced VOF method

This paper summarizes the technical development and validation of a multiphase computational fluid dynamics (CFD) numerical method using the volume-of-fluid (VOF) model and a Lagrangian tracking model which can be employed to analyze general multiphase flow problems with free surface mechanism. The gas-liquid interface mass, momentum and energy conservation relationships are modeled by continuum surface mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed flow regimes. The objectives of the present study are to develop and verify the fractional volume-of-fluid cell partitioning approach into a predictor-corrector algorithm and to demonstrate the effectiveness of the present approach by simulating benchmark problems including laminar impinging jets, shear coaxial jet atomization and shear coaxial spray combustion flows.

Chen, Yen-Sen↗

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.↗

Evaluation of out-of-core computer programs for the solution of symmetric banded linear equations

FORTRAN coded out-of-core equation solvers that solve using direct methods symmetric banded systems of simultaneous algebraic equations. Banded, frontal and column (skyline) solvers were studied as well as solvers that can partition the working area and thus could fit into any available core. Comparison timings are presented for several typical two dimensional and three dimensional continuum type grids of elements with and without midside nodes. Extensive conclusions are also given.

Dunham, R. S.↗

Application of a vectorized particle simulation in high-speed near-continuum flow

An efficient particle simulation technique, developed for use on vector architecture based supercomputers for studying hypersonic rarefied gas flows is employed to simulate the complex wake generated by Mach six flow over a 10 deg half-angle wedge for freestream Reynolds numbers of 1780 and 3560. Data obtained are compared against higher Reynolds number experimental results. Simulations utilized as many as 10 to the 5th computational cells and 10 to the 7th simulated particles having power-law interaction potentials. A code performance of 1.8-2.4 microsec of Cray-2 CPU time to process a single particle per timestep is achieved. Diffuse adiabatic and isothermal wedge surface models are used in this investigation. Although the wedge geometry is two-dimensional, the simulation incorporates a width-wise direction, resulting in a three-dimensional computation.

Woronowicz, Michael S.↗

DSMC analysis in a heterogeneous parallel computing environment

A methodology for implementing parallel DSMC codes in a heterogeneous computing environment is described. The methodology involves the use of a common message-passing software library together with recently developed software that handles the actual interprocessor communications in a standard manner across a variety of computing platforms. Benchmark tests using a simple DSMC model problem were performed on an Intel iPSC/860, a Cray-YMP and a group of Sun workstations. The approach was found to give speedups that scaled linearly with problem size on all the computing platforms tested. This methodology was then incorporated into a production-type DSMC code to allow the simulation of problems that would not otherwise have been practical. The application of this production code to simulations of hypersonic shear flows and shock-lip interactions under near-continuum conditions is described. Synchronous and asynchronous models for implementing parallelism into DSMC simulations are also described and both models are shown to produce the same steady-state result.

Wilmoth, R. G.↗