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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

NASA NTRS1973

A highly stable explicit integration technique for computational continuum mechanics

A user-oriented subroutine package is built around a highly stable explicit integration algorithm for solution of large order systems of ordinary differential equations, as result from discretization of initial-boundary value problems in continuum mechanics. Fast and accurate solutions, for problems in laminar and turbulent, two and three-dimensional viscous flow fields and multi-dimensional transient heat transfer, are presented using this algorithm, as embodied within a general purpose finite element computer program.

Baker, A. J.↗

NASA NTRS1972

COMOC: Thermal analysis variant. User's manual

The Thermal Analysis Variant of the COMOC (computational continuum mechanics) computer system solves problems involving transient heat conduction and convection in stationary continua spanning arbitrarily irregular two-dimensional and axisymmetric solution domains. COMOC is based upon a finite element solution algorithm for the energy equation, and solves for the transient nodal temperature distribution using a highly stable and automatic explicit integration procedure. COMOC is extensively user-oriented, requires minimal input, and no a priori knowledge concerning the stability character of the differential equation system. It can readily output computed data in user-specified format fields, that geometrically resemble the solution domain discretization (for rapid engineering evaluation). Complete information is provided for applying COMOC to a specific problem.

Bauer, A. M.↗

NASA NTRS1989

Rarefied gas dynamics

The current state of those aspects of rarefied gas dynamics research that appear to be most important to research planned over the next decade is evaluated. These aspects encompass assessments of computational rarefied-gas dynamics (CRGD) that will allows their use as surrogates for experiments, the development of hybris-flowfield computational techniques matching continuum computations with particle computations, and the validation of CRGD through the results of experimental studies of Knudsen layers in simple flows. The design of surfaces for the achievement of stable, low-momentum and thermal accommodation coefficients will be a major priority, together with theorization and experimentation on evaporation and condensation effects close to surfaces.

Muntz, E. P.↗

NASA NTRS2012

Multiscale Modeling of Ultra High Temperature Ceramics (UHTC) ZrB2 and HfB2: Application to Lattice Thermal Conductivity

We are developing a multiscale framework in computational modeling for the ultra high temperature ceramics (UHTC) ZrB2 and HfB2. These materials are characterized by high melting point, good strength, and reasonable oxidation resistance. They are candidate materials for a number of applications in extreme environments including sharp leading edges of hypersonic aircraft. In particular, we used a combination of ab initio methods, atomistic simulations and continuum computations to obtain insights into fundamental properties of these materials. Ab initio methods were used to compute basic structural, mechanical and thermal properties. From these results, a database was constructed to fit a Tersoff style interatomic potential suitable for atomistic simulations. These potentials were used to evaluate the lattice thermal conductivity of single crystals and the thermal resistance of simple grain boundaries. Finite element method (FEM) computations using atomistic results as inputs were performed with meshes constructed on SEM images thereby modeling the realistic microstructure. These continuum computations showed the reduction in thermal conductivity due to the grain boundary network.

Lawson, John W.↗

NASA NTRS2012

Multiscale Modeling of UHTC: Thermal Conductivity

Lawson, John W.↗

NASA NTRS1987

The 2D Rydberg series in Al I

High quality ab initio electonic structure calculations were performed on the 2D Rydberg series in Al I. The configuration 3s3p2(2D) is shown to contribute substantially to the lowest four 2D Rydberg states. The same configuration also contributes substantially to a 2D state embedded in the ionization continuum. Computed oscillator strengths for the first six members of the 2D Rydberg transitions are given: these should be of substantially high accuracy than currently available values.

Taylor, Peter R.↗

NASA NTRS1988

The 2D Rydberg series in Al I

High-quality ab initio electronic structure calculations have been performed on the 2D Rydberg series in Al I. The configuration 3s3p2(2D) is shown to contribute substantially to the lowest four 2D Rydberg states. The same configuration also contributes substantially to a 2D state embedded in the ionization continuum. Computed oscillator strengths for the first six members of the 2D Rydberg transitions are given: these should be of substantially higher accuracy than currently available values.

Taylor, Peter R.↗

NASA NTRS1978

The dynamics and control of large flexible space structures. Part B: Development of continuum model and computer simulation

The equations of motion of an arbitrary flexible body in orbit were derived. The model includes the effects of gravity with all its higher harmonics. As a specific example, the motion of a long, slender, uniform beam in circular orbit was modelled. The example considers both the inplane and three dimensional motion of the beam in orbit. In the case of planar motion with only flexible vibrations, the pitch motion is not influenced by the elastic motion of the beam. For large values of the square of the ratio of the structural modal frequency to the orbital angular rate the elastic motion was decoupled from the pitch motion. However, for small values of the ratio and small amplitude pitch motion, the elastic motion was governed by a Hill's 3 term equation. Numerical simulation of the equation indicates the possibilities of instability for very low values of the square of the ratio of the modal frequency to the orbit angular rate. Also numerical simulations of the first order nonlinear equations of motion for a long flexible beam in orbit were performed. The effect of varying the initial conditions and the number of modes was demonstrated.

Bainum, P. M.↗

NASA NTRS1993

Analysis of the anisotropic viscoplastic-damage response of composite laminates - Continuum basis and computational algorithms

The mathematical structure underlying the rate equations of a recently-developed constitutive model for the coupled viscoplastic-damage response of anisotropic composites is critically examined. In this regard, a number of tensor projection operators have been identified, and their properties were exploited to enable the development of a general computational framework for their numerical implementation using the Euler fully-implicit integration method. In particular, this facilitated (i) the derivation of explicit expressions of the (consistent) material tangent stiffnesses that are valid for both three-dimensional as well as subspace (e.g. plane stress) formulations, (ii) the implications of the symmetry or unsymmetry properties of these tangent operators from a thermodynamic standpoint, and (iii) the development of an effective time-step control strategy to ensure accuracy and convergence of the solution. In addition, the special limiting case of inviscid elastoplasticity is treated. The results of several numerical simulations are given to demonstrate the effectiveness of the schemes developed.

Saleeb, A. F.↗

NASA NTRS1997

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

NASA NTRS2021

An Aerodynamics Database for Aerocapture of a Conceptual Morphable Entry System at Mars

An aerodynamics database for aerocapture of a conceptual morphable entry system at Mars was developed. The vehicle is envisioned to employ deflections of its geometry to modulate angle of attack and sideslip angle for trajectory control. Modified Newtonian methods were used to compute hypersonic continuum aerodynamics. Direct simulation Monte Carlo calculations were carried out to compute data in the rarefied transitional regime. Free-molecular theory was used to compute values in the free-molecular regime. A bridging function method was then used to bridge data from all three flow regimes together and produce a database that can be used for future trajectory simulations to assess the vehicle’s performance. The vehicle was found to have the control authority in hypersonic continuum to effect an 18 deg change in the trim angle of attack using a 10 deg deflection of the pitching control surfaces. Similarly, a 6 deg change in the trim sideslip angle could be achieved with a 5 deg deflection of the yawing control surfaces. Different configurations of the vehicle were found to be statically stable in pitch and yaw in the continuum regime, but unstable in the rarefied and free-molecular regimes.

aerocapture↗

NASA NTRS2017

Multiscale Fiber Kinking: Computational Micromechanics and a Mesoscale Continuum Damage Mechanics Models

In this work, the fiber kinking phenomenon, which is known as the failure mechanism that takes place when a fiber reinforced polymer is loaded under longitudinal compression, is studied. A computational micromechanics model is employed to interrogate the assumptions of a recently developed mesoscale continuum damage mechanics (CDM) model for fiber kinking based on the deformation gradient decomposition (DGD) and the LaRC04 failure criteria.

Herraez, Miguel↗

NASA NTRS2018

Prediction of Thermal Protection System Material Permeability and Hydraulic Tortuosity Factor Using Direct Simulation Monte Carlo

Carbon preforms used in Thermal Protection System (TPS) materials are 80 to 90% porous, allowing for boundary layer and pyrolysis gases to flow through the porous regions. The bulk material properties such as permeability and hydraulic tortuosity factor affect the transport of the boundary layer gases. The use of Direct Simulation Monte Carlo along with the Klinkenberg permeability formulation allows us to compute the continuum permeability and Knudsen correction factor for flow in the transition regime. In this work, we have computed the permeability for two types of carbon preforms, namely, Morgan Felt and FiberForm, and assessed the effect of orientation on the permeability. Since both the materials are anisotropic, the permeability was found to depend on orientation, wherein, the materials are more permeable in the in-plane orientation than the through-thickness orientation. The through-thickness orientation was also more tortuous compared to the in-plane material orientation. Compared to Morgan Felt, FiberForm is less permeable, in both, through thickness and in-plane directions.

Jambunathan, Revathi↗

NASA NTRS1987

Computational fluid dynamics near the continuum limit

Two dimensional Navier-Stokes equations for a perfect gas are solved for hypersonic flow over a flat plate at Reynolds numbers ranging from continuum to the free molecule flow. In the transition flow regime, new slip boundary conditions are introduced, which reduce to the well-known slip conditions of Maxwell at small Knudsen numbers and yield the correct shear stress and heat transfer in the limiting case of free molecule flow. Comparison of the computed results with the existing experimental data and Monte Carlo calculations indicates that the continuum Navier-Stokes equations give surprisingly realistic results throughout the transition flow regime, when the new slip boundary conditions are used.

Maccormack, Robert W.↗

NASA NTRS2019

Aerodynamics for the ADEPT SR-1 Flight Experiment

Adaptable, Deployable, Entry, and Placement Technology (ADEPT) is a combination of a heatshield and an aerodynamic decelerator for atmospheric entry applications. The ADEPT Sounding Rocket (SR)-1 mission was a suborbital flight experiment of an 0.7 m-diameter ADEPT to verify system-level performance and to characterize dynamic stability behavior. The aerodynamic database for ADEPT SR-1 was constructed from non-continuum and continuum flowfield computations, along with data from recent ADEPT ground testing and the IRVE-3 flight test vehicle. High-altitude (free-molecular and transitional regimes) data were generated using DSMC methods. Pre-flight predictions of continuum static aerodynamics coefficients were derived from Reynolds-Averaged Navier-Stokes solutions at conditions along a design trajectory, with comparisons to available ground test data of the nano-ADEPT geometry. Dynamic pitch damping characteristics were taken from functional forms developed for the IRVE-3 flight test vehicle through ballistic range testing. Comparison of pre-flight predictions to post-flight reconstruction of aerodynamic force and moment coefficients is presented.

Korzun, Ashley M.↗

NASA NTRS1980

Impact melt generation and transport

The results from the first two calculations in a series of continuum mechanics computer code calculations, investigating the effects of variations in impactor mass and velocity on the generation and transport of impact melt, are reported. In the present calculations, the impactor is modeled as a spherical iron projectile with a mass of one trillion grams, and the target as a gabbroic anorthosite (GA) half-space, where the cases calculated have impact velocities of 5 and 15.8 km/sec. Early-time ejection velocities are 1-2 km/sec in both cases. The first calculation results in 0.07 projectile masses of GA being partly or completely melted, with all the melted GA being ejected from the crater, and a maximum impact range for the ejected melted material of 30 km. The second calculation yields 10.4 projectile masses of melted GA, 50% of which is ejected from the crater to ranges of up to about 130 km. Peak shock pressure attenuation with depth is reported for both cases, and transient cavity dynamics are described and compared to that for surface and near-surface explosions.

Orphal, D. L.↗

NASA NTRS1981

Z-model analysis of impact cratering - An overview

The Maxwell Z-Model has been applied to two continuum mechanics computer calculations: (1) a laboratory-scale impact of an aluminum projectile into plasticene clay, and (2) a planetary-scale impact of an iron meteor into gabbroic anorthosite. The material flow in the cratering flow field may be well approximated by incompressible flow for most of the excavation stage of crater growth. The center of the flow field is located beneath, not at, the surface. Soon after energy partitioning is complete, Z can assume values less than 2.0 associated with the initial directedness of the projectile's momentum. The Z-Model parameters are time dependent during a significant portion of the crater growth time, and Z increases steadily with time from about 2.0 or slightly less at the beginning of the excavation stage to level off at values in the neighborhood of about 3.0 before the excavation stage is half-over.

Austin, M. G.↗

NASA NTRS2001

General Multimechanism Reversible-Irreversible Time-Dependent Constitutive Deformation Model Being Developed

Since most advanced material systems (for example metallic-, polymer-, and ceramic-based systems) being currently researched and evaluated are for high-temperature airframe and propulsion system applications, the required constitutive models must account for both reversible and irreversible time-dependent deformations. Furthermore, since an integral part of continuum-based computational methodologies (be they microscale- or macroscale-based) is an accurate and computationally efficient constitutive model to describe the deformation behavior of the materials of interest, extensive research efforts have been made over the years on the phenomenological representations of constitutive material behavior in the inelastic analysis of structures. From a more recent and comprehensive perspective, the NASA Glenn Research Center in conjunction with the University of Akron has emphasized concurrently addressing three important and related areas: that is, 1) Mathematical formulation; 2) Algorithmic developments for updating (integrating) the external (e.g., stress) and internal state variables; 3) Parameter estimation for characterizing the model. This concurrent perspective to constitutive modeling has enabled the overcoming of the two major obstacles to fully utilizing these sophisticated time-dependent (hereditary) constitutive models in practical engineering analysis. These obstacles are: 1) Lack of efficient and robust integration algorithms; 2) Difficulties associated with characterizing the large number of required material parameters, particularly when many of these parameters lack obvious or direct physical interpretations.

Saleeb, A. F.↗