Search NASASearch

SEARCH · Search NASA

Results for “Computational solid mechanics”

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.

386 records · Page 8

Boundary Layer Analysis of Shock Tube Flows

Shock tubes offer a controlled environment to reproduce kinetic and radiative phenomena characteristic of atmospheric entry flows under ground-test conditions. The boundary layer developing behind the incident shock wave determines the available test time, influences particle residence times important for similarity scaling, and can affect radiative energy transport. In this study, we couple a quasi-1D space marcher with the compressible boundary-layer equations to numerically compute the post-shock flow in a shock frame of reference for various test gas mixtures representative of different planetary atmospheres. The solvers are individually verified against CFD simulations and analytical correlations available in the literature. Coupled solutions are computed for a finite-rate chemistry in the boundary layer and a non-catalytic isothermal wall. Results yield refined estimates of the maximum separation distance, as well as insights into concentration profiles of relevant species within the boundary layer.

Andrea Fagnani

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation

Advances in Tactical Rocket Propulsion

An AGARD Colloquium on 'Progress in Tactical Rocket Propulsion' was held at La Jolla, California, on April 22 and 23, 1965. The unclassified papers, or unclassified versions of classified papers, presented at this meeting form the contents of the present AGARD conference proceedings, which has been identified by a new title that reflects more properly the scope of the remaining technical papers. The articles contained in this volume represent significant evaluations of the status of solid-propellant, liquid-propellant and hybrid-propellant rocket development. Together with the reviews, they should provide a useful reference source for research workers and students alike.

S S Penner

Future Aerospace Technology in the Service of the Alliance: Sustained Hypersonic Flight - Volume 3

Sustained Hypersonic Flight capability is envisioned for several future military and civilian applications, such as long range immediate reaction reconnaissance, high speed interception of air targets, long range precision strike against hardened or time critical targets, and access to space. The inherent reduction in time-to-target and low vulnerability will permit new operational tactics. The symposium outlined mission needs and operational scenarios for hypersonic vehicles. Addressed were: Technological issues and challenges in external hypersonic aerodynamics and design, propulsion and engine/airframe integration, military ramjet applications, overall system design including structures and materials development needs, and test facilities. Ground based test facilities cannot provide full simulation at Mach numbers above 5. Therefore, the use of Computational Fluid Dynamics design tools is essential. For routine use of the computational tools in analysis, design and optimization, it is necessary to reduce the time for the entire computational process by two orders of magnitude. Hypersonic boundary layer transition remains a critical design issue because of the important impact on engine drag and on heating, which can affect the choice of materials and thermal protection systems. In addition, reliable prediction of scramjet net thrust is an absolute must in resolving hypersonic air breathing vehicle design issues. Due to current facility and computational shortfalls, the development of future hypersonic flight systems requires research flight tests in the technology areas of boundary layer transition and air-breathing propulsion engine performance. For sustained hypersonic flight beyond Mach 6, the supersonic combustion ramjet (scramjet) engine is the only choice for the near future. Only this air-breathing concept offers a significant promise of large reductions in required propellant fractions, increased payload fractions, and reduced size vehicles, together with a foreseeable technological feasibility. Airframe/engine integration, combustor design and thermal management are the predominant engineering tasks. Fuels, hydrogen or hydrocarbon, must be matched to the operational needs of military or civil use. Experience in existing ramjet propelled missiles capable of speeds up to Mach 4 can support the development effort. The potential mission and cost benefits of sustained hypersonic flight to both military and civil applications are tremendous. From the budget point of view, the possibility of sharing development costs between military and civil programs offers a specific advantage.

Combustion chambers

The Evolution of the Moon and the Terrestrial Planets

The thermal evolutions of the Moon, Mars, Venus, and Mercury are calculated theoretically starting from cosmochemical condensation models. An assortment of geological, geochemical, and geophysical data are used to constrain both the present-day temperatures and the thermal histories of the planets' interiors. Such data imply that the planets were heated during or shortly after formation and that all the terrestrial planets started their differentiations early in their history. Initial temperatures and core formation play the most important roles in the early differentiation. The size of the planet is the primary factor in determining its present-day thermal state. The Moon, smallest in size, is characterized as a differentiated body with a crust, a thick solid mantle, and an interior region which may be partially molten. It is presently cooling rapidly and is relatively inactive tectonically. Mercury, which probably has a large core, may have a 500-km-thick solid lithosphere and a partially molten core, if it is assumed that some heat sources exist in the core. If this is not the case, the planet's interior temperatures are everywhere below the melting curve for iron. The thermal evolution is dominated by the core separation and the high conductivity of iron that makes up the bulk of Mercury. Mars, intermediate in size, is assumed to have differentiated an Fe-FeS core. While the formation of an early crust is evident, large-scale melting and differentiation of the mantle silicates has occurred at least up until 1 b.y. ago. Present-day temperature profiles indicate moderate tectonic activity at the present time. Venus is characterized as a planet not unlike the Earth in many respects. Core formation has occurred probably during the first billion years after the formation. The present-day temperatures indicate a partially molten upper mantle overlain by a 100-km-thick lithosphere and a molten Fe-Ni core. We can expect that today Venus may have tectonic processes similar to the Earth's.

M Nafi Toksöz

Parallel Processing and Expert Systems

Whether it be monitoring the thermal subsystem of Space Station Freedom, or controlling the navigation of the autonomous rover on Mars, NASA missions in the 90's cannot enjoy an increased level of autonomy without the efficient use of expert systems. Merely increasing the computational speed of uniprocessors may not be able to guarantee that real time demands are met for large expert systems. Speed-up via parallel processing must be pursued alongside the optimization of sequential implementations. Prototypes of parallel expert systems have been built at universities and industrial labs in the U.S. and Japan. The state-of-the-art research in progress related to parallel execution of expert systems was surveyed. The survey is divided into three major sections: (1) multiprocessors for parallel expert systems; (2) parallel languages for symbolic computations; and (3) measurements of parallelism of expert system. Results to date indicate that the parallelism achieved for these systems is small. In order to obtain greater speed-ups, data parallelism and application parallelism must be exploited.

Jerry C Yan

Mechanics of Preloaded Bolt Tensile Loading With Focus on Load Introduction Factor

The bolt tensile and joint separation loads are directly influenced by the locations at which the external loads enter the clamped members of a preloaded bolted joint (PBJ) and the associated load-paths through the joint. This physical load introduction mechanism affecting the bolt tensile loading is typically represented in the bolt tensile load equation, in part, by a load introduction factor (LIF), which was shown by H.M. Lee of Marshall Spaceflight Center to be a natural product of the bolt tensile load equation using a linear spring stiffness model. This LIF, being a function of load-path stiffness, has subsequently been denoted as the stiffness-based LIF (SBLIF), providing a framework to calculate the LIF using whatever load-path stiffness approximations are appropriate. Expanding upon the work of Lee, it is shown that the SBLIF and the joint stiffness factor are functions of the stiffnesses of the same load-paths and regions within a PBJ, and thus they should not be treated as independent variables. Mathematical expressions for the SBLIF are presented. Comparisons are shown between the analytically calculated SBLIF, the analytically calculated geometric LIF (GLIF), which is a simple clamped-member thickness ratio, the experimentally derived LIF, and the LIF determined by finite element analysis (FEA). Using experiment and FEA as a benchmark, the SBLIF, using traditional load-path stiffness approximations, enables a more accurate prediction of bolt tensile loading than the GLIF, although it can be unconservative near joint separation. The GLIF generally attributes more of the externally applied tensile load to the bolt than does the SBLIF, potentially resulting in heavier and/or more costly bolted joints. Mathematical relationships between the SBLIF and the GLIF are developed. Supplemental material is provided in the appendixes where the historical practice of using the joint compressive stiffness in place of the joint tensile stiffness is evaluated. The appendixes include step-by-step examples demonstrating the calculation of the SBLIF using traditional stiffness approximations and conclude with the development of the joint diagram in terms of the SBLIF, culminating into formulas for the key features of a joint diagram, which is useful for programming.

Load Path

Softening the Gap between Wöhler and Paris – New Approaches for Fatigue Analysis –

Fatigue analysis tools can vary across industries. For example, automotive engineers often use the Wöhler (S-N) approach to design for safe-life, while aerospace engineers prioritize damage tolerance and inspection intervals, relying instead on crack growth models such as Paris’ law. Although both approaches may deal with the control of cracks in similar materials, their analysis tools and material characterizations are fundamentally distinct. This divide mirrors the classic split between stress-based strength analysis and linear elastic fracture mechanics. However, modern nonlinear models that incorporate material softening, such as cohesive laws, blur this boundary and capture fracture behaviors across scales. This presentation describes the CF23 fatigue model, which uses cohesive softening to link S-N crack initiation with crack propagation rates. CF23 spans the full fatigue spectrum, from initial propagation transients to steady-state growth and threshold conditions, offering a unified framework that bridges Wöhler and Paris-based methodologies. Example applications include fatigue crack propagation transients in adhesive interfaces and skin/stiffener separation.

cohesive elements

Battery Pack Shape Optimization using Transient Heat Conduction Coupled with Cell-Discharge Analysis

Battery electric systems exhibit significant time-dependence, especially when evaluated in the context of an aircraft mission profile with continually changing power demands. Additionally, when evaluating battery-powered aircraft concepts, it is important to accurately compute the temperature of the batteries and properly characterize the thermal response of the system. The temperature of the batteries has a significant impact on cell performance, in addition to safety considerations of maintaining battery temperatures below their operating limit. Because of these considerations, battery models for preliminary design and optimization of aircraft should include the capability to accurately compute the temperature distribution within the battery pack. Furthermore, battery pack designs should be as light-weight as possible to maximize the pack energy density, while also considering battery temperature limits. Here, we demonstrate a simultaneous trajectory and shape optimization of a battery pack concept, using a transient heat transfer finite element model coupled with a time-varying cell-discharge battery model to provide this capability. Including the transient finite element problem in the loop enables accurate temperatures that can be passed back to the cell discharge model, while the cell discharge model can supply the finite element model with time-varying heat boundary conditions to the finite element problem, further benefiting the fidelity of the thermal response of the batteries. We first demonstrate the coupling capability between the battery cell-discharge model and the transient finite-element heat transfer through an optimization which computes the optimal current profile for the battery pack while ensuring the battery temperatures remain below their operational limit. We then build on this optimization by adding shape optimization to the problem, which allows us to consider a composite objective function which also minimizes the mass of the battery pack, while also producing an optimal current discharge profile.

Optimization

Battery simulation program

Computer program calculates spacecraft battery energy at specific times dictated by input sequence of recharge, discharge, and no activity phases. Recharge phases are assumed during times not specified, unless batteries are fully charged. Warnings are printed when energy falls below specified level. Program assumes two identical batteries.

Abreu, A.

Battery-Charge-State Model

Charge-state model for lead/acid batteries proposed as part of effort to make equivalent of fuel gage for battery-powered vehicles. Models based on equations that approximate observable characteristics of battery electrochemistry. Uses linear equations, easier to simulate on computer, and gives smooth transitions between charge, discharge, and recuperation.

Vivian, H. C.

Battery Pack Shape Optimization using Transient Heat Conduction Coupled with Cell-Discharge Analysis

Battery electric systems exhibit significant time-dependence, especially when evaluated in the context of an aircraft mission profile with continually changing power demands. Additionally, when evaluating battery-powered aircraft concepts, it is important to accurately compute the temperature of the batteries and properly characterize the thermal response of the system. The temperature of the batteries has a significant impact on cell performance, in addition to safety considerations of maintaining battery temperatures below their operating limit. Because of these considerations, battery models for preliminary design and optimization of aircraft should include the capability to accurately compute the temperature distribution within the battery pack. Furthermore, battery pack designs should be as light-weight as possible to maximize the pack energy density, while also considering battery temperature limits. Here, we demonstrate a simultaneous trajectory and shape optimization of a battery pack concept, using a transient heat transfer finite element model coupled with a time-varying cell-discharge battery model to provide this capability. The transient finite-element analysis is done using TACS, and the cell-discharge battery model uses OpenMDAO and dymos. Including the transient finite element problem in the loop enables accurate temperatures that can be passed back to the cell discharge model, while the cell discharge model can supply the finite element model with time-varying heat boundary conditions, further benefiting the fidelity of the thermal response of the batteries. We first demonstrate the coupling capability between the battery cell-discharge model and the transient finite-element heat transfer through an optimization which computes the optimal current profile for the battery pack while ensuring the battery temperatures remain below their operational limit. We then build on this optimization by adding shape optimization to the problem, which allows us to consider a composite objective function which also minimizes the mass of the battery pack, while also producing an optimal current discharge profile.

Optimization

Trajectory Engineering with Modular Patched Conics for Entry Systems and TPS (TEMPEST)

Brief Presenter Biography (35 word limit): Bohdan Wesely is an Aerospace Engineer in the Entry Systems and Technology Division at Ames. He has worked on a variety of projects for NASA including integrated TPS (thermal protection system) flight hardware deliveries and testing services for commercial partners. Introduction: TEMPEST is a new trajectory analysis framework that is designed to fill the gap between dedicated flight mechanics tools and aerothermal and TPS sizing tools. The project started as an SJSU master’s thesis and has since evolved into a general conceptual design tool capable of studying a wide variety of entry problems. Development is ongoing in the Entry Systems and Technology Division at NASA ARC. Why TEMPEST: Space missions involving entry into a planetary atmosphere involve a series of unique requirements across multiple disciplines. Whether it is traditional entry descent and landing (EDL), or aerocapture, the vehicle must navigate to its target landing location or orbit state, and the TPS must protect the payload during entry. The design process typically involves iterative handoffs between various flight mechanics, flow solver, and material response level tools. During the early conceptual phase, a wide variety of feasible trajectories are simulated in a Monte Carlo scenario which broadly satisfy the mission or landing requirements. Next, computational fluid dynamics (CFD), direct simulation Monte Carlo (DSMC), and other flow solver analyses are performed at various key trajectory points to generate an aero-database, heating and TPS design requirements also emerge at this stage. At this point, with updated aerodynamics from the various flow solvers, trajectories can be re-run, this in turn can change the required freestream conditions for the CFD tools, and as a project progresses, these analyses converge, and uncertainty is reduced. However, there is always a “hand-off” occurring between two inherently coupled phenomena. Analysis Description: One of the goals with TEMPEST is to use a variety of first principles estimation methods coupled with an atmosphere model to predict vehicle aerothermodynamics across the entire flight regime while propagating a 3 or 6 degree of freedom (DoF) trajectory. Aerodynamics methods include modified Newtonian, Maxwell and Cercignani- Lampis-Lord (CLL) for continuum, transitional, and free molecular flow regimes. Aerothermodynamics include boundary layer and reference enthalpy methods, and Mutation++ for non-equilibrium chemistry modeling. TEMPEST is also capable of stitching multiple trajectory segments together to study mission scenarios like multi-pass aerocapture and aero-gravity assists. Most of the program is implemented in MATLAB using modern system objects, it relies on several C++ shared libraries for supporting tools like Gmsh, the Global Reference Atmospheric Model (GRAM), and Mutation++. The various first principles aerothermal estimation methods are discretized across either a structured axisymmetric panel mesh or an unstructured tri-mesh generated from an open-source tool such as Gmsh, this allows solutions on the same mesh to be compared across tools such as CB-Aero. CFD Coupling. A physics-aware, gaussian process CFD anchoring scheme is proposed to adjust the various first principles methods as a CFD database is populated. One goal for this anchoring module is to inform the project where CFD should be run. Full knowledge of the entire trajectory, atmosphere, and aerothermodynamics allows for easier identification of high sensitivity areas and uncertainty quantification. While the first principles effects are well known and proven accurate in existing tools such as CB- Aero and Cart3D, a physics aware CFD anchoring scheme increases tool credibility across a project lifecycle. Material Response Modeling. Correct TPS sizing is critical for optimizing mass for science payloads and ensuring mission success. The process typically involves a thermal analysis along the trajectory with surface heating environments as a boundary condition. Several design constraints are maximum bondline temperature and maximum recession with various margining techniques. The material response tool FIAT, developed out of NASA Ames, is currently being integrated into the TEMPEST environment. TPS recession, shape change, mass loss, and mass property alteration are all factors that can perturb an entry trajectory. For missions like Mars 2020, recession was minimal and was safely handled separately as a post process. For missions such as Jupiter Galileo with a high TPS mass fraction or asteroid entries, recession plays a major role. The proposed fully coupled scheme is to use an epoch-based approach where the trajectory integration is halted after a recession threshold, the energy balance and FIAT are solved at each panel, the mesh, aerodynamics, and mass properties are updated, and the trajectory continues. Several computational tradeoffs have been made during the development of TEMPEST to limit the cost of a single trajectory and preserve its utility as a conceptual, rapid iteration tool. Conclusion: Development of TEMPEST is ongoing and the project is still in its infancy. This talk aims to showcase its unique capabilities to support future NASA entry systems missions.

Bohdan O Wesely

Impact of Inertia, Friction, and Backlash upon Force Control in Telemanipulation

The mechanical behavior of master controllers of telemanipulators has been a concern of both designers and implementors of telerobotic systems. In general, the literature recommends that telemanipulator systems be constructed that minimize inertia, friction, and backlash in an effort to improve telemanipulative performance. For the most part, these recommendations are founded upon theoretical analysis or simply intuition. Although these recommendations are not challenged on their merit, the material results are measured of building and fielding telemanipulators that possess less than ideal mechanical behaviors. Experiments are described in which forces in a mechanical system with human input are evaluated as a function of mechanical characteristics such as inertia, friction, and backlash. Results indicate that the ability of the human to maintain gripping forces was relatively unaffected by dynamic characteristics in the range studied, suggesting that telemanipulator design in this range should be based on task level force control requirements rather than human factors.

Neil A Duffie

In-Situ Scanning Electron Microscope Experiments for Microscale Mechanical Testing and Validated Modeling of Fiber Reinforced Thermoplastics

A novel, in-situ, scanning electron microscope (SEM) mechanical testing capability for materials at the microscale which provides experimental validation to a machine learning (ML) toolset for full-field validation of physics-based micromechanics models is being developed by researchers at NASA Glenn Research Center. These are enabling technologies for the integration of multiscale digital twins for materials into system level models which will result in the improved performance, material discovery, reduced production cost and time, rapid characterization, and prognostic structural health monitoring (SHM) for materials and structures for extreme environments in support of NASA space exploration missions. In order to bridge the material structure-to-system gap for digital twins, physics-based models must be experimentally validated at multiple length scales. Seminal microscale experiments, conducted at the Air Force Research Laboratory (AFRL), were limited to transverse compression of single-layer, unidirectional thermoset polymer matrix composite (PMC) micropillar specimens [1]. The early phases of the current project followed those initial results and setup to reproduce the compression testing of PMC material on the custom-built piezoelectric actuated micromechanical testing rig built by MicroTesting Solutions LLC. In this work, samples of thermoplastic PMC material were first machined into 3 mm cubes, and then further machining and final milling was done using a Focused Ion Beam (FIB). The initial experiment was done on a pillar roughly 20 µm x 20 µm x 40 µm tall. Additional pillars were milled with final sizes ranging from 20 µm x 20 µm x 40 µm tall to 40 µm x 40 µm x 65 µm tall. A speckle pattern for in-situ full-field measurements using Digital Image Correlation (DIC) was applied with platinum, which was coated on the surface, and then the FIB was used to mill away some of the coating to produce an irregular pattern of Pt on the pillar surface. The samples were loaded into the custom testing rig and placed into the SEM and loaded under compression until failure. Images were collected in the SEM during testing. Post-processing of the images was conducted using DIC to obtain full-field displacement and strain measurements elucidating the role of the matrix as well as fiber-fiber interaction at the microscale within the composite subjected to compression loading well into the non-linear regime of the material. Moreover, the evolution of fiber-matrix debonding and matrix cracking is observed in-situ at the microscale. This data, along with images segmented with a newly developed ML toolset [2], was used to create and validate physics-based micromechanics models. An image of the failed micropillar is shown in Figure 1. The techniques developed in the initial compression experiment was tailored to the validation needs of the models and expanded to include different sized samples as well as possibly tension and fatigue.

Laura Wilson

Experimental Investigations of Low-Energy (4 to 40 eV) Collisions of O - ( 2 P) Ions and O( 3 P) Atoms With Surfaces

Using a newly-developed, magnetically confined source, low-energy, ground state oxygen negative ions and neutral atoms are generated. The energy range is variable, and atom and neutrals have been generated at energies varying from 2 eV to 40 eV and higher. It was found that the interaction of these low-energy species with a solid magnesium fluoride target leads to optical emissions in the (at least) visible and infrared regions of the spectrum. Researchers describe y details of the photodetachment source, and present spectra of the neutral and ion glows in the wavelength range 250 to 850 nm (for O - ) and 600 to 850 nm (for O), and discuss the variability of the emissions for incident energies between 4 and 40 eV.

A Chutjian

Radiation Temperature and Extinction of Transient Gaseous Diffusion Flames in Microgravity

The absence of buoyancy-induced flows in micro-g and the resulting increase in the reactant residence time significantly alters the fundamentals of many combustion processes. Substantial differences between 1-g and micro-g flames have been reported in experiments on candle flames, flame spread over solids, droplet combustion and others. These differences are more basic than just in the visible flame shape. Longer residence times and higher concentration of combustion products in the flame zone create a thermochemical environment which changes the flame chemistry and the heat and mass transfer processes. Processes such as flame radiation (and its interaction with flame chemistry), that are often ignored under normal gravity, become very important and sometimes even controlling. This is particularly true for conditions at extinction of a micro-g diffusion flame.

Arvind Atreya