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

Computational Materials Techniques for Thermal Protection Solutions: Materials and Process Design

Integrated computational materials techniques that span the atomistic and continuum scales have the potential to aid the design and manufacturing of thermal protection materials. Two cases demonstrating the practical application of these methods are discussed. Case one examines the selection of a high temperature coating for carbon/carbon, with the target application being a solar thermal propulsion heat exchanger. The performance of various refractory metal and metal-carbide coatings is characterized considering extreme thermal (3500 K) and chemical (hydrogen flows) conditions. The recession rate, hydrogen leakage, and likelihood of mechanical failure are characterized and provide directions for further experimental investigation. Case two examines the process optimization of a heat shield material composed of a woven silica fiber preform and cyanate ester resin. Frequently, internal voids were found to be present in this composite after the resin infusion and curing stages of manufacturing. Using the manufacturing conditions, computations indicate that both water adsorption and resin cure shrinkage are contributing factors to void formation. The results suggest an alternative process configuration for curing that would mitigate voids.

Haskins, Justin↗

Computational Materials Techniques for Thermal Protection Solutions: Materials and Process Design

Integrated computational materials techniques that span the atomistic and continuum scales have the potential to aid the design and manufacturing of thermal protection materials. Two cases demonstrating the practical application of these methods are discussed. Case one examines the selection of a high temperature coating for carbon/carbon, with the target application being a solar thermal propulsion heat exchanger. The performance of various refractory metal and metal-carbide coatings is characterized considering extreme thermal (3500 degrees Kelvin) and chemical (hydrogen flows) conditions. The recession rate, hydrogen leakage, and likelihood of mechanical failure are characterized and provide directions for further experimental investigation. Case two examines the process optimization of a heat shield material composed of a woven silica fiber preform and cyanate ester resin. Frequently, internal voids were found to be present in this composite after the resin infusion and curing stages of manufacturing. Using the manufacturing conditions, computations indicate that both water adsorption and resin cure shrinkage are contributing factors to void formation. The results suggest an alternative process configuration for curing that would mitigate voids.

Thermal Protection↗

Simulating Underexpanded Plumes in Martian and Lunar Environments

While the Loci/Chem Computational fluid dynamics (CFD) application has been validated and used extensively for simulating launch environments in atmospheric conditions, use of this tool for simulating supersonic plumes at much lower ambient pressure requires further validation. Simulations of underexpanded plumes have been performed using Loci/Chem at Martian pressure to predict several metrics for laminar, turbulent, and impinging plumes. The CFD results are compared with an experimental data set for low Reynolds number plumes at these conditions in order to evaluate the current capability of the Loci/Chem tool for these types of environments. The CFD validation results to date show reasonable agreement with the experimental data across all the metrics of interest for the configurations considered. The Loci/Chem-Boltzmann tool is a hybrid continuum/rarefied flow solver which extends modeling capabilities to very low-pressure environments such as those on the Moon. The Loci/Chem-Boltzmann solver uses a gradient-based continuum breakdown criterion to restrict solution of the computationally expensive Boltzmann equation to only a subset of the domain, while using the Navier-Stokes equations elsewhere. This CFD application is under active development, and current efforts toward establishing production ready capability for evaluating Lunar plume surface interactions are well under way. Initial simulations of an Apollo Lunar excursion module indicate that several regions of the flow require solving the Boltzmann equations due to extreme rarefaction. Early simulations are promising, indicating reasonable overall computational time for a full 3D human scale lander simulation.

Plume Surface Interaction↗

Simulating Underexpanded Jets in Martian and Lunar Environments

As part of the Game Changing Development (GCD) Program, funded by NASA’s Space Technology Mission Directorate (STMD), the development of simulation capability for the prediction of extra-terrestrial Plume Surface Interaction (PSI) environments has been undertaken by the Fluid Dynamics Branch at NASA/MSFC. The GCD PSI Project, planned to be accomplished over a four year period, contains a Predictive Simulation Capability (PSC) Element focused on creating simulation capability for the reliable and accurate prediction of PSI in Martian (~650 Pa) and Lunar (vacuum) ambient environments. In addition to the PSC Element, the PSI Project also contains a companion Ground Testing Element for development of focused datasets for validation of predictive capability as well as a Flight-focused Instrumentation Element. This paper describes the activities and accomplishments in the past year for the Prediction of Plume Flow in low pressure environments component of the PSI Project. While the Loci/Chem Computational fluid dynamics (CFD) application has been validated and used extensively for simulating launch environments in atmospheric conditions, use of this tool for simulating supersonic plumes at Mars-like ambient pressure requires further validation. CFD simulations of underexpanded jets have been performed using Loci/Chem for Mars-like conditions to predict several metrics for steady laminar, turbulent, and impinging plumes. The CFD results are compared with an experimental data set for low Reynolds number plumes at these conditions in order to evaluate the current capability of the Loci/Chem tool for these types of environments. The CFD validation results to date show reasonable agreement with the experimental data across all of the metrics of interest for the configurations considered. The Loci/Chem-Boltzmann CFD application is a hybrid continuum/rarefied flow solver which extends modeling capabilities to very low pressure environments such as those on the Moon. The Loci/Chem-Boltzmann solver uses a gradient-based continuum breakdown criterion to restrict solution of the computationally expensive Boltzmann equation to only a subset of the domain, while using the Navier-Stokes equations elsewhere. This CFD application is under active development, and current efforts toward establishing production ready capability for evaluating Lunar plume surface interactions are well under way. Initial simulations of an Apollo LEM indicate that several regions of the flow require solving the Boltzmann equations due to extreme rarefaction. Early simulations are promising, indicating reasonable overall computational time for a full 3D human scale lander simulation.

Thomas Shurtz↗

Software for continuum modeling of controls-structures interactions

It is clear that computer software is needed to assist in the generation of the equations of motion for complex, flexible spacecraft. Daniel Poelaert of ESTEC has developed the software DISTEL with which he has modeled the structural dynamics for different satellites. He is interested in expanding the capabilities of DISTEL to include structural damping and control systems. Unfortunately, the software has not been released. The author has developed similar software, PDEMOD, which has been used to model the Spacecraft control Laboratory Experiment (SCOLE), the Solar Array Flight Experiment (SAFE), the Mini-MAST truss, and the LACE satellite. PDEMOD has been used also for optimal parameter estimation and integrated control-structures design. PDEMOD is also being extended to include structural damping and control systems which are imbedded into the same equations for the structural dynamics. This paper will address the formulation of the equations for the structural dynamics of spacecraft structures which are constructed of a 3-dimensional arrangement of rigid bodies and flexible beam elements. Control system dynamics are imbedded into the same equations so that model order reduction approximations are not necessary. The input data consists of the physical data of the elements and the topological information describing how the elements are connected. PDEMOD accomplishes the following: (1) automatically assembles the equations of motion for the entire structural model; (2) calculates the modal frequencies; (3) calculates the mode shapes; (4) generates perspective views of the mode shapes; and (5) forms selected transfer functions. The software PDEMOD continues to be developed to provide additional features to assist in analyzing and synthesizing control and structural systems for flexible spacecraft.

Taylor, Lawrence W., Jr.↗

Reentry heating and flow fields for Space Shuttle Orbiter at high angle of attack.

Rarefied flow fields and heating rates are calculated by a Monte Carlo direct simulation technique for several shuttle orbiter configurations at angles of attack between 40 and 60 deg. Continuum inviscid flow fields are determined by 2-D unsteady and 3-D steady finite-difference/artificial viscosity techniques used in combination with a 2-D shock layer analysis method. Continuum viscous flow fields and heating rates are computed by a boundary-layer integral matrix technique and an aerodynamic heating method.

Rochelle, W. C.↗

Further results on cooling mechanisms for soft X-ray flares

Analysis of data from the soft X-ray proportional counter spectrometer on OSO-4 has been continued along the lines of earlier work. It was noticed that the time profiles of X-ray counts from the OSO-4 instrument's 1 to 3 Angstrom detectors for two events showed a very rapid initial decay after maximum, followed by a much more slowly falling section. At least in the case of the larger of the two events, the change of slope seems quite definitely to be discontinuous. Temperatures and emission measures for both flares have been derived, by fitting an emission function to eight-channel count histograms which form the output of the 1 to 3 Angstrom detectors, analyzed into eight energy intervals by pulse-height analysis. The computer program is described. It employs temperature and continuum emission measures in the Culhane-Acton formula and the flux in the 6.7 KeV iron-line feature as free parameters, adjusting them by small amounts in successive iterations until the original historgram is approximately reproduced. A chi-squared is used to examine the agreement between histograms and terminates the iteration accordingly. Account is taken of the proportional counters' energy resolution in the fitting process.

Craig, I. J. D.↗

Design sensitivity analysis and optimization tool (DSO) for sizing design applications

The DSO tool, a structural design software system that provides the designer with a graphics-based menu-driven design environment to perform easy design optimization for general applications, is presented. Three design stages, preprocessing, design sensitivity analysis, and postprocessing, are implemented in the DSO to allow the designer to carry out the design process systematically. A framework, including data base, user interface, foundation class, and remote module, has been designed and implemented to facilitate software development for the DSO. A number of dedicated commercial software/packages have been integrated in the DSO to support the design procedures. Instead of parameterizing an FEM, design parameters are defined on a geometric model associated with physical quantities, and the continuum design sensitivity analysis theory is implemented to compute design sensitivity coefficients using postprocessing data from the analysis codes. A tracked vehicle road wheel is given as a sizing design application to demonstrate the DSO's easy and convenient design optimization process.

Chang, Kuang-Hua↗

Cosmic Metal Production and the Contribution of QSO Absorption Systems to the Ionizing Background

The recent discovery by Cowie et al. (1995) and Tytler et al. (1995) of metals in the Ly alpha clouds shows that the intergalactic medium (IGM) at high redshift is contaminated by the products of stars and suggests that ionizing photons from massive star formation may be a significant contributor to the UV background radiation at early epochs. We assess the validity of the stellar photoionization hypothesis. Based on recent computations of metal yields and 0-star Lyman continuum (Lyc) fluxes, we find that 0.2 percent of the rest-mass energy of the metals produced is radiated as Lyc. By modeling the transfer of ionizing radiation through the IGM and the rate of chemical enrichment, we demonstrate that the background intensity of photons at 1 ryd that accompanies the production of metals in the Ly alpha forest clouds may be significant, approaching 0.5 x 10(exp -21) ergs cm squared s(-1) Hz(-1) sr(-1) at z approximately equals 3 if the Lyc escape fraction is greater than of equal to 0.25. Together with quasars, massive stars could then, in principle, provide the hydrogen and helium Lyc photons required to ionize the universe at high redshifts. We propose that observations of the He2 Gunn-Peterson effect and of the metal ionization states of the Ly alpha forest and Lyman-limit absorbers should show the signature of a stellar spectrum. We also note that the stellar photoionization model fails if a large fraction of the UV radiation emitted from stars cannot escape into the IGM, as suggested by the recent Hopkins Ultraviolet Telescope observations by Leitherer et al. (1995) of low-redshift starburst galaxies, or if most of the metals observed at z is approximately 3 were produced at much earlier epochs.

Madau, Piero↗

Estimating Drag and Heating Coefficients for Hollow Reentry Objects in Transitional Flow Using DSMC

In NASA’s Object Reentry Survival Analysis Tool (ORSAT), aerodynamic drag and aerothermal heating coefficients are computed for each of the free-molecular, continuum, and transitional flow regimes using analytical and semi-analytical methods. These methods are typically limited to convex, blunt objects (such as spheres) and are applied to other objects such as boxes and cylinders using multiplicative “shape factors” to account for the different behavior. Previous literature has analyzed the aerodynamic and aerothermodynamic properties of flow around sharp-edged objects like boxes and cylinders in transitional flow, though only those objects with solid external boundaries. However, many reentry objects we have encountered in real spacecraft have been hollow (i.e., with the potential to allow flow through them). We present here preliminary results from analyses performed using the NASA Direct Simulation Monte Carlo (DSMC) Analysis Code (DAC) on hollow cylinders and boxes (with varying wall thickness-diameter ratio).

Marichalar, Jeremiah J.↗

Developing a Pyrolysis Gas Thermal Blocking Model for Reentry Demise

In NASA’s Object Reentry Survival Analysis Tool (ORSAT), aerodynamic drag and aerothermal heating coefficients are computed for each of the free-molecular, continuum, and transitional flow regimes using analytical and semi-analytical methods. These heating coefficients were derived for typical metallic materials that melt and do not have a strong gas-phase contribution to the flow in the boundary layer. Modern satellites typically feature fiber-reinforced polymer (FRP) components, such as solar array booms, facesheets of sandwich panels, or overwraps for composite-overwrapped pressure vessels (COPV). These FRP materials do not behave the same as metals in the reentry environment, but instead will pyrolyze and develop significant volumes of gas into the boundary layer. Accurately predicting the reentry demise of FRP components is critical to assessing the reentry casualty risk for modern spacecraft. Research in recent years has shown that this demisability can depend heavily on how the expulsion of gaseous pyrolysis products through the outer surface of the material affects the heat flux at the surface. The ODPO has been developing a reduced-order model of the effect of pyrolysis gas blowing on the heat flux based on correlations between a blowing factor and a non-dimensional heat flux to be incorporated in the upcoming version 7.3 of the Object Reentry Survivability Analysis Tool (ORSAT). This presentation discusses the progress of this development project and the challenges remaining for generalizing the model across families of FRP materials.

Benton Greene↗

Equivalent-Continuum Modeling With Application to Carbon Nanotubes

A method has been proposed for developing structure-property relationships of nano-structured materials. This method serves as a link between computational chemistry and solid mechanics by substituting discrete molecular structures with equivalent-continuum models. It has been shown that this substitution may be accomplished by equating the vibrational potential energy of a nano-structured material with the strain energy of representative truss and continuum models. As important examples with direct application to the development and characterization of single-walled carbon nanotubes and the design of nanotube-based devices, the modeling technique has been applied to determine the effective-continuum geometry and bending rigidity of a graphene sheet. A representative volume element of the chemical structure of graphene has been substituted with equivalent-truss and equivalent continuum models. As a result, an effective thickness of the continuum model has been determined. This effective thickness has been shown to be significantly larger than the interatomic spacing of graphite. The effective thickness has been shown to be significantly larger than the inter-planar spacing of graphite. The effective bending rigidity of the equivalent-continuum model of a graphene sheet was determined by equating the vibrational potential energy of the molecular model of a graphene sheet subjected to cylindrical bending with the strain energy of an equivalent continuum plate subjected to cylindrical bending.

Odegard, Gregory M.↗

Fast Particle Methods for Multiscale Phenomena Simulations

We are developing particle methods oriented at improving computational modeling capabilities of multiscale physical phenomena in : (i) high Reynolds number unsteady vortical flows, (ii) particle laden and interfacial flows, (iii)molecular dynamics studies of nanoscale droplets and studies of the structure, functions, and evolution of the earliest living cell. The unifying computational approach involves particle methods implemented in parallel computer architectures. The inherent adaptivity, robustness and efficiency of particle methods makes them a multidisciplinary computational tool capable of bridging the gap of micro-scale and continuum flow simulations. Using efficient tree data structures, multipole expansion algorithms, and improved particle-grid interpolation, particle methods allow for simulations using millions of computational elements, making possible the resolution of a wide range of length and time scales of these important physical phenomena.The current challenges in these simulations are in : [i] the proper formulation of particle methods in the molecular and continuous level for the discretization of the governing equations [ii] the resolution of the wide range of time and length scales governing the phenomena under investigation. [iii] the minimization of numerical artifacts that may interfere with the physics of the systems under consideration. [iv] the parallelization of processes such as tree traversal and grid-particle interpolations We are conducting simulations using vortex methods, molecular dynamics and smooth particle hydrodynamics, exploiting their unifying concepts such as : the solution of the N-body problem in parallel computers, highly accurate particle-particle and grid-particle interpolations, parallel FFT's and the formulation of processes such as diffusion in the context of particle methods. This approach enables us to transcend among seemingly unrelated areas of research.

Koumoutsakos, P.↗

Equivalent-Continuum Modeling of Nano-Structured Materials

A method has been developed for modeling structure-property relationships of nano-structured materials. This method serves as a link between computational chemistry and solid mechanics by substituting discrete molecular structures with an equivalent-continuum model. It has been shown that this substitution may be accomplished by equating the vibrational potential energy of a nano-structured material with the strain energy of representative truss and continuum models. As an important example with direct application to the development and characterization of single-walled carbon nanotubes, the model has been applied to determine the effective continuum geometry of a graphene sheet. A representative volume element of the equivalent-continuum model has been developed with an effective thickness. This effective thickness has been shown to be similar to, but slightly smaller than, the interatomic spacing of graphite.

Odegard, Gregory M.↗

Multiscale Analysis of Structurally-Graded Microstructures Using Molecular Dynamics, Discrete Dislocation Dynamics and Continuum Crystal Plasticity

A multiscale modeling methodology is developed for structurally-graded material microstructures. Molecular dynamic (MD) simulations are performed at the nanoscale to determine fundamental failure mechanisms and quantify material constitutive parameters. These parameters are used to calibrate material processes at the mesoscale using discrete dislocation dynamics (DD). Different grain boundary interactions with dislocations are analyzed using DD to predict grain-size dependent stress-strain behavior. These relationships are mapped into crystal plasticity (CP) parameters to develop a computationally efficient finite element-based DD/CP model for continuum-level simulations and complete the multiscale analysis by predicting the behavior of macroscopic physical specimens. The present analysis is focused on simulating the behavior of a graded microstructure in which grain sizes are on the order of nanometers in the exterior region and transition to larger, multi-micron size in the interior domain. This microstructural configuration has been shown to offer improved mechanical properties over homogeneous coarse-grained materials by increasing yield stress while maintaining ductility. Various mesoscopic polycrystal models of structurally-graded microstructures are generated, analyzed and used as a benchmark for comparison between multiscale DD/CP model and DD predictions. A final series of simulations utilize the DD/CP analysis method exclusively to study macroscopic models that cannot be analyzed by MD or DD methods alone due to the model size.

Saether, Erik↗

High-voltage spark atomic emission detector for gas chromatography

A dc-powered, double-gap, miniature nanosecond spark source for emission spectrochemical analysis of gas chromatographic effluents is described. The spark is formed between two thoriated tungsten electrodes by the discharge of a coaxial capacitor. The spark detector is coupled to the gas chromatograph by a heated transfer line. The gas chromatographic effluent is introduced into the heated spark chamber where atomization and excitation of the effluent occurs upon breakdown of the analytical gap. A microcomputer-controlled data acquisition system allows the implementation of time-resolution techniques to distinguish between the analyte emission and the background continuum produced by the spark discharge. Multiple sparks are computer averaged to improve the signal-to-noise ratio. The application of the spark detector for element-selective detection of metals and nonmetals is reported.

Calkin, C. L.↗

Explicit Hilbert-space representations of atomic and molecular photoabsorption spectra - Computational studies of Stieltjes-Tchebycheff functions

Computational methods are reported for construction of discrete and continuum Schroedinger states in atoms and molecules employing explicit Hilbert space procedures familiar from bound state studies. As theoretical development, the Schroedinger problem of interest is described, the Cauchy-Lanczos bases and orthonormal polynomials used in constructing L-squared Stieltjes-Tchebycheff (ST) approximations to the discrete and continuum states are defined, and certain properties of these functions are indicated. Advantages and limitations of the ST approach to spectral studies relative to more conventional calculations are discussed, and aspects of the approach in single-channel approximations to larger molecules are described. Procedures are indicated for construction of photoejection anisotropies and for performing coupled-channel calculations employing the ST formalism. Finally, explicit descriptive intercomparisons are made of the nature and diagnostic value of ST functions with more conventional scattering functions.

Hermann, M. R.↗