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 163 records · Page 9

Pressure measurements in a low-density nozzle plume for code verification

Measurements of Pitot pressure were made in the exit plane and plume of a low-density, nitrogen nozzle flow. Two numerical computer codes were used to analyze the flow, including one based on continuum theory using the explicit MacCormack method, and the other on kinetic theory using the method of direct-simulation Monte Carlo (DSMC). The continuum analysis was carried to the nozzle exit plane and the results were compared to the measurements. The DSMC analysis was extended into the plume of the nozzle flow and the results were compared with measurements at the exit plane and axial stations 12, 24 and 36 mm into the near-field plume. Two experimental apparatus were used that differed in design and gave slightly different profiles of pressure measurements. The DSMC method compared well with the measurements from each apparatus at all axial stations and provided a more accurate prediction of the flow than the continuum method, verifying the validity of DSMC for such calculations.

Penko, Paul F.↗

Static Aerodynamics of the Mars Exploration Rover Entry Capsule

The static aerodynamics for the Mars Exploration Rover (MER) aeroshell are presented. This aerodynamic database was an integral part of the end-to-end simulation used in pre- entry analysis for determining the MER entry design requirements for development of the MER entry system, as well as targeting the MER landing sites. The database was constructed using the same approach used for Mars Pathfinder (MPF). However, the MER aerodynamic database is of much higher fidelity and tailored to the MER entry trajectories. This set of data includes direct simulation Monte Carlo calculations covering the transitional regime of the entry trajectory and computational fluid dynamics calculations describing the aerodynamics in the hypersonic and supersonic continuum regimes. An overview of the methodology used to generate the data is given along with comparisons to important features in the MPF aerodynamics and related heritage data. The MER and MPF comparison indicates that trajectory specific data is required to properly model the flight characteristics of a.blunt entry capsule at Mars.

Schoenenberger, Mark↗

Decameter storm radiation. II

The height, size, directivity, and lifetime of decametric storm sources are investigated. The height of the continuum sources is determined from the rotation rate, and the gradient of electron density for the regions of storm origin is computed by assuming that the radiation originates at the plasma frequency. The mean angular size of the sources is found to be quite large (increasing with decreasing frequency), the storm continuum is strongly directive toward the disk center, and east-west asymmetry is observed at decameter wavelengths. Two distinct classes of type III bursts are observed during storms: 'off-fringe' (displaced in position from the continuum source) and 'on-fringe' (coinciding in position with the continuum). A model of the storm region is proposed in which the 'on-fringe' bursts originate in regions of diverging field lines above closed magnetic loops and the 'off-fringe' bursts are excited by energetic electron streams having access to open field lines at the base of the loops.

Gergely, T. E.↗

Ultraviolet and Visible Emission Mechanisms in Astrophysics

The project involved the study of ultraviolet (UV) and visible emission mechanisms in astrophysical and atmospheric environments. In many situations, the emission is a direct consequence of a charge transferring collision of an ion with a neutral with capture of an electron to an excited state of the product ion. The process is also important in establishing the ionization and thermal balance of an astrophysical plasma. As little of the necessary collision data are available, the main thrust of the project was the calculation of total and state-selective charge transfer cross sections and rate coefficients for a very large number of collision systems. The data was computed using modern explicit techniques including the molecular-orbital close-coupling (MOCC), classical trajectory Monte Carlo (CTMC), and continuum distorted wave (CDW) methods. Estimates were also made in some instances using the multichannel Landau-Zener (MCLZ) and classical over-the-barrier (COB) models. Much of the data which has been computed has been formatted for inclusion in a charge transfer database on the World Wide Web (cfadc.phy.ornl.gov/astro/ps/data/). A considerable amount of data has been generated during the lifetime of the grant. Some of it has not been analyzed, but it will be as soon as possible, the data placed on our website, and papers ultimately written.

Stancil, Phillip C.↗

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↗

Development of a Global 3-D magnetohydrodynamic computational model for solar wind-cometary and planetary studies

The development of a global, 3-D magnetohydrodynamic computational model is described. This model is used to quantitatively describe the detailed continuum field and plasma interaction process of the solar wind with cometary and planetary magneto/ionopause shapes. The present solar wind/terrestrial planet interaction (Level 1) model (which is based on an axisymmetric gas dynamic plus frozen field approximation to the full MHD equations) is extended to a full 3-D gas dynamic (Level 2) approximation. A mass loading capability in the Level 1 interaction model is implemented.

Stahara, S. S.↗

The Development of Computational Techniques for the Nonlinear Analysis of Composite Structures at High Temperature

The objective of this research project was to develop robust and efficient computational tools for deformation and life analysis techniques. This included: development of viscoplastic continuum models, associated numerical integration techniques, suitable micromechanics models/methods, and coupled deformation and damage algorithms. The work was involved with models/methods applicable to composite structures on both microscale and macroscale levels. The microscale analysis utilized Aboudi's generalized method of cells (GMS) micromechanical model as well as finite element models of similar composite architectures, i.e., square and hexagonal fiber pack models.

Keith, Theo G., Jr.↗

A flight experiment to measure rarefied-flow aerodynamics

A flight experiment to measure rarefied-flow aerodynamics of a blunt lifting body is being developed by NASA. This experiment, called the Rarefied-Flow Aerodynamic Measurement Experiment (RAME), is part of the Aeroassist Flight Experiment (AFE) mission, which is a Pathfinder design tool for aeroassisted orbital transfer vehicles. The RAME will use flight measurements from accelerometers, rate gyros, and pressure transducers, combined with knowledge of AFE in-flight mass properties and trajectory, to infer aerodynamic forces and moments in the rarefied-flow environment, including transition into the hypersonic continuum regime. Preflight estimates of the aerodynamic measurements are based upon environment models, existing computer simulations, and ground test results. Planned maneuvers at several altitudes will provide a first-time opportunity to examine gas-surface accommondation effects on aerodynamic coefficients in an environment of changing atmospheric composition. A description is given of the RAME equipment design.

Blanchard, Robert C.↗

Boundary formulations for three-dimensional continuum structural shape sensitivity analysis

The direct, singular, boundary element analysis formulation is shown to provide a basis for a computationally efficient and accurate shape design sensitivity analysis approach for the structural response of three-dimensional solid objects. The theoretical formulation for surface displacement and traction component sensitivities, and all components of the stress tensor is presented along with a formulation for the recovery of displacement and stress components in the interior of the object under consideration. Discussion of computational issues related to the overall efficiency of these formulations is given, along with numerical results to demonstrate the accuracy and efficiency of this approach.

Kane, J. H.↗

Scattered light in the echelle modes of the Goddard high-resolution spectrograph aboard the Hubble Space Telescope. I - Analysis of prelaunch calibration data

A study of scattered light in the echelle mode of the Goddard high-resolution spectrograph (GHRS), one of the first-generation science instruments aboard the Hubble Space Telescope, is reported. The data analyzed are from the GHRS ground-based test program and include observations with the 500 channel Digicon detectors through the small science aperture of UV emission-line sources, continuum sources, and continuum sources viewed through molecular absorption cells. The two largest contributors to scattered light are the echelle and cross-disperser, the effects of which have been quantified. A computer simulation using the derived scattering functions demonstrates that the combined effects of these sources of scattered light adequately explain the observed level of interorder light when viewing UV continuum sources. Observing strategies and data reduction algorithms for measuring and removing the effects of the scattered light are discussed.

Cardelli, Jason A.↗

Updates on the Predictive Materials Modeling Software Tools

Updates on NASA‘s efforts to build a Predictive Material Modeling (PMM) framework from the micro-scale to the macro-scale are presented in this abstract. The PMM effort is part of the Entry Systems Modeling (ESM) project under NASA’s Game Changing Development (GCD) program. To reduce the need for extensive testing and accelerate the design cycle process, ESM is developing simulation and modeling tools that enable the characterization of the properties of thermal protection materials and their response to extremely hot plasma. The Porous Microstructure Analysis (PuMA) software has been developed to compute effective material properties and perform material response simulations on digitized microstructures of porous media. PuMA is able to import three-dimensional digital images obtained from X-ray microtomography or to generate artificial microstructures that mimic real materials. PuMA also provides a module for interactive 3D visualizations. Version 3, which was recently released as open-source, includes modules to compute simple morphological properties such as porosity, volume fractions, pore diameter, and specific surface area. Additional capabilities include the determination of effective thermal and electrical conductivity (both radiative and solid conduction - including the ability to simulate local anisotropy for the latter); effective diffusivity and tortuosity from the continuum to the rarefied regime; techniques to determine the local material orientation, as well as mechanical properties (elasticity coefficients), and permeability. Computed properties are then used to inform a macro-scale material response model, such as those implemented in the Porous material Analysis Toolbox based on OpenFOAM (PATO) software developed within ESM. The computational model in PATO is a generic heat and mass transfer model for porous reactive materials containing several solid phases and a single gas phase. The detailed chemical interactions occurring between the solid phases and the gas phase are modeled at the pore scale, assuming Local Thermal Equilibrium. Recent efforts include the development of a mechanical erosion model as well as a unified model allowing an intrinsic coupling between fluid and material. Comparison to flight data (Mars Science Laboratory [MSL] Entry Descent and Landing Instrument [MEDLI] and Mars 2020 MEDLI2) is critical in order to validate these computational tools. Examples of ablative material response using the code will be presented, including 3D simulations of the full-scale heatshield of the MSL capsule. The simulations demonstrated the ability of the modern material response code, PATO, to handle the material response of geometrically complex and large domains through the use of massively parallel computations.

material modeling↗

Statistical equilibrium calculations for silicon in early-type model stellar atmospheres

Line profiles of 36 multiplets of silicon (Si) II, III, and IV were computed for a grid of model atmospheres covering the range from 15,000 to 35,000 K in effective temperature and 2.5 to 4.5 in log (gravity). The computations involved simultaneous solution of the steady-state statistical equilibrium equations for the populations and of the equation of radiative transfer in the lines. The variables were linearized, and successive corrections were computed until a minimal accuracy of 1/1000 in the line intensities was reached. The common assumption of local thermodynamic equilibrium (LTE) was dropped. The model atmospheres used also were computed by non-LTE methods. Some effects that were incorporated into the calculations were the depression of the continuum by free electrons, hydrogen and ionized helium line blocking, and auto-ionization and dielectronic recombination, which later were found to be insignificant. Use of radiation damping and detailed electron (quadratic Stark) damping constants had small but significant effects on the strong resonance lines of Si III and IV. For weak and intermediate-strength lines, large differences with respect to LTE computations, the results of which are also presented, were found in line shapes and strengths. For the strong lines the differences are generally small, except for the models at the hot, low-gravity extreme of our range. These computations should be useful in the interpretation of the spectra of stars in the spectral range B0-B5, luminosity classes III, IV, and V.

Kamp, L. W.↗

Modeling Deformable Linear Objects for Autonomous Robotic Outfitting of Lunar Surface Systems

This paper presents structural models of deformable linear objects (DLOs). DLOs are a subclass of deformable objects that encompasses common outfitting elements such as cables and ropes. Models are validated through hardware experiments, and integration in a robotic autonomy architecture for space environments is discussed. A persistent human presence on the lunar surface is one of the next major milestones in space exploration. This requires the development of robust extraplanetary construction technologies including structures and materials modeling and robotic systems. Previous robotic construction technology development has primarily focused on structural assembly, with significantly less focus on robotically performed outfitting tasks to instantiate subsystems providing power, data, life support, etc. These tasks involve manipulation of highly flexible elements, which are difficult to model, such as cable harnesses, ropes, and hoses. Robotic manipulation of DLOs, especially cable harnesses, is an active area of research as cable harnesses are essential for providing power and data to space assets. DLO models that can be used for robot manipulator trajectory generation are necessary for autonomous operation of lunar infrastructure. There are many proposed methods for modeling DLOs, and they primarily fall into three types: 1) discrete model-based, 2) continuum model-based, and 3) Neural Network-based. These types each have pros and cons, and the tradeoff between model accuracy and computational speed informs which type should be used. An understanding of this trade-off is imperative for real-time control of autonomous systems. High computational requirements reduce the speed of the model, making real-time control difficult, while accuracy is critical to preventing collisions. Discrete models, such as a mass-spring multibody representation, require relatively few calculations, and accuracy is directly tied to the step size of the discretization. Continuum models, such as a B-spline representation or a Cosserat rod model (a mix of continuous and discrete), are more informed of the structural properties of the cable and are much more accurate than a rigid body mass-spring model, but at significant computational cost. A Neural Network approach can provide an online solution with very few computational steps, but properly generating training data can be difficult and validation for an in-space application is not trivial. This paper explores the trade-off between different modeling approaches and compares accuracy and computational speed/complexity of the three types mentioned above. Model accuracy is evaluated using a cable in a static configuration. True cable shape is obtained using a depth camera for RGB images and point-cloud segmentation. The purpose of this experiment is to evaluate the trade-offs of different approaches to the DLO modeling problem. Understanding the tradeoffs between different cable modeling techniques paves the way for developing robotic control and planning architectures necessary for real-time manipulation of DLOs for lunar infrastructure outfitting. Real-time control is required for robotic systems to be able to actively manipulate a cable in a harsh environment where model and sensor errors compound, and environmental conditions can cause significant disturbances. Cable routing must be performed in areas with high density of objects/obstacles: through truss structures, near solar panels or mirror arrays, next to bundles of electrical equipment. Understanding the best way to plan and manipulate a cable without disrupting the environment or damaging the cable is imperative to robotic outfitting operations on the lunar surface.

Amy M Quartaro↗

Two-dimensional hybrid continuum/particle approach for rarefied flows

A hybrid numerical technique previously developed for one-dimensional rarefied gas flows is generalized to two dimensions. The method is based on the fact that the flowfield that develops near a body in a rarefied gas typically contains local regions of continuum, transitional and free-molecular flow. By utilizing the solution technique most appropriate for each region and coupling the techniques in an interface region where both are applicable, more computationally efficient solutions can be obtained, or equivalently, more complex flowfields can be analyzed. The present method combines finite difference solution of the Navier-Stokes equations in the continuum regions, with Direct Simulation Monte Carlo in the more rarefied regions. The two schemes are coupled interactively via a general conservative flux boundary condition. The method is tested by application to the model problem of pressure-driven rarefied flow through a slit. Results show the hybrid scheme offers a speedup of a factor of nearly two for the nominal conditions considered, due to the decrease in the size of the Direct Simulation domain.

Wadsworth, Dean C.↗

Numerical Study of a Continuum Sonic Jet Interacting with a Rarefield Flow

The results of a numerical study with flow and boundary conditions based on an experiment of a continuum sonic jet interacting with rarefied flow about a sharp leading edge flat plate at zero incidence are presented. Comparisons are made between computational fluid dynamics (CFD) and direct simulation Monte Carlo (DSMC) solutions which provide an assessment of applying each technique to the flow conditions of the experiment. An analysis of the CFD results revealed a correlation between the interaction interface of the jet continuum breakdown surface and a non-dimensional parameter derived from jet and free stream flow conditions. Using the breakdown surface from the correlation, the continuum jet was uncoupled from the interaction, thus allowing an uncoupled CFD-DSMC solution to be obtained. Also, a nearest neighbor collision algorithm, similar to the subcell technique, was implemented in the DSMC solution technique. The comparison between CFD and DSMC results shows good qualitative agreement in the interaction region and good quantitative agreement elsewhere.

Glass, Christoper E.↗

Extension of a Kinetic Approach to Chemical Reactions to Electronic Energy Levels and Reactions Involving Charged Species With Application to DSMC Simulations

The ability to compute rarefied, ionized hypersonic flows is becoming more important as missions such as Earth reentry, landing high mass payloads on Mars, and the exploration of the outer planets and their satellites are being considered. Recently introduced molecular-level chemistry models that predict equilibrium and nonequilibrium reaction rates using only kinetic theory and fundamental molecular properties are extended in the current work to include electronic energy level transitions and reactions involving charged particles. These extensions are shown to agree favorably with reported transition and reaction rates from the literature for nearequilibrium conditions. Also, the extensions are applied to the second flight of the Project FIRE flight experiment at 1634 seconds with a Knudsen number of 0.001 at an altitude of 76.4 km. In order to accomplish this, NASA's direct simulation Monte Carlo code DAC was rewritten to include the ability to simulate charge-neutral ionized flows, take advantage of the recently introduced chemistry model, and to include the extensions presented in this work. The 1634 second data point was chosen for comparisons to be made in order to include a CFD solution. The Knudsen number at this point in time is such that the DSMC simulations are still tractable and the CFD computations are at the edge of what is considered valid because, although near-transitional, the flow is still considered to be continuum. It is shown that the inclusion of electronic energy levels in the DSMC simulation is necessary for flows of this nature and is required for comparison to the CFD solution. The flow field solutions are also post-processed by the nonequilibrium radiation code HARA to compute the radiative portion of the heating and is then compared to the total heating measured in flight.

Liechty, Derek S.↗

Extension of a Kinetic Approach to Chemical Reactions to Electronic Energy Levels and Reactions Involving Charged Species with Application to DSMC Simulations

The ability to compute rarefied, ionized hypersonic flows is becoming more important as missions such as Earth reentry, landing high mass payloads on Mars, and the exploration of the outer planets and their satellites are being considered. Recently introduced molecular-level chemistry models that predict equilibrium and nonequilibrium reaction rates using only kinetic theory and fundamental molecular properties are extended in the current work to include electronic energy level transitions and reactions involving charged particles. These extensions are shown to agree favorably with reported transition and reaction rates from the literature for near-equilibrium conditions. Also, the extensions are applied to the second flight of the Project FIRE flight experiment at 1634 seconds with a Knudsen number of 0.001 at an altitude of 76.4 km. In order to accomplish this, NASA's direct simulation Monte Carlo code DAC was rewritten to include the ability to simulate charge-neutral ionized flows, take advantage of the recently introduced chemistry model, and to include the extensions presented in this work. The 1634 second data point was chosen for comparisons to be made in order to include a CFD solution. The Knudsen number at this point in time is such that the DSMC simulations are still tractable and the CFD computations are at the edge of what is considered valid because, although near-transitional, the flow is still considered to be continuum. It is shown that the inclusion of electronic energy levels in the DSMC simulation is necessary for flows of this nature and is required for comparison to the CFD solution. The flow field solutions are also post-processed by the nonequilibrium radiation code HARA to compute the radiative portion.

Liechty, Derek S.↗

Apparent solar temperature enhancement due to large-amplitude waves

The effect of slow-mode acoustic-type MHD waves propagating outward in the solar atmosphere on the temperature structures predicted by empirical models is investigated analytically. A model is constructed, and numerical results are presented for wing intensities, line profiles, temperature enhancements, waves with higher energy flux, temperature depression, and the Si continuum. The flux in the MgII and CaII UV resonance lines is found to be increased relative to that in the IR continuum, leading to model temperatures which depend systematically on which observations are used in the computation. It is suggested that mechanical heating may take place in smaller regions such as flux tubes rather than uniformly over the surface.

Kalkofen, W.↗