Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solid dynamics”

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 469 records · Page 26

Dynamic characteristics of a variable-mass flexible missile: Dynamics of a two-stage variable-mass flexible rocket

The dynamic characteristics of two-stage slender elastic body were investigated. The first stage, containing a solid-fuel rocket, possesses variable mass while the second stage, envisioned as a flexible case, contains packaged instruments of constant mass. The mathematical formulation was in terms of vector equations of motion transformed by a variational principle into sets of scalar differential equations in terms of generalized coordinates. Solutions to the complete equations were obtained numerically by means of finite difference techniques. The problem has been programmed in the FORTRAN 4 language and solved on an IBM 360/50 computer. Results for limited cases are presented showing the nature of the solutions.

Meirovitch, L.↗

Three isoparametric solid elements for NASTRAN

Linear, quadratic, and cubic isoparametric hexahedral solid elements have been added to the element library of NASTRAN. These elements are available for static, dynamic, buckling, and heat-transfer analyses. Because the isoparametric element matrices are generated by direct numerical integration over the volume of the element, variations in material properties, temperatures, and stresses within the elements are represented in the computations. In order to compare the accuracy of the new elements, three similar models of a slender cantilever were developed, one for each element. All elements performed well. As expected, however, the linear element model yielded excellent results only when shear behavior predominated. In contrast, the results obtained from the quadratic and cubic element models were excellent in both shear and bending.

Johnson, S. E.↗

The spin-temperature theory of dynamic nuclear polarization and nuclear spin-lattice relaxation

A detailed derivation of the equations governing dynamic nuclear polarization (DNP) and nuclear spin lattice relaxation by use of the spin temperature theory has been carried to second order in a perturbation expansion of the density matrix. Nuclear spin diffusion in the rapid diffusion limit and the effects of the coupling of the electron dipole-dipole reservoir (EDDR) with the nuclear spins are incorporated. The complete expression for the dynamic nuclear polarization has been derived and then examined in detail for the limit of well resolved solid effect transitions. Exactly at the solid effect transition peaks, the conventional solid-effect DNP results are obtained, but with EDDR effects on the nuclear relaxation and DNP leakage factor included. Explicit EDDR contributions to DNP are discussed, and a new DNP effect is predicted.

Byvik, C. E.↗

Dynamics of volcanic plumes on Io

Ballistic and aerodynamic models are proposed to explain the dynamics of the volcanic plumes on Io. A ballistic model, which assumes that eruptive gas and solid-particle flows become completely and permanently decoupled over a small distance as compared with the maximum particle height, is shown to be unable to account for the observed shape of Plume 3. An aerodynamic model which assumes the entrainment and transportation of solid particles by the gas flow into the central portion of the top of the plume, with ballistic particle trajectories only after passage through the first part of their descent, is then developed. The model accounts for plume shape by the formation of a weak decompression shock near the edge of the plume stalk, which constrains horizontal expansion until a sufficient altitude is reached. Calculated plume dimensions are then compared with observations of Ionian Plumes 1 and 3.

Cook, A. F.↗

Hypersonic, stratified gas flows past an obstacle - Direct simulation Monte Carlo calculations

Monte Carlo particle dynamics are presently used to study the interaction of a stratified gas that impinges hypersonically on a solid body, whose size is comparable to the gas scale height. Attention is given to the results of these computations for two representative cases, where the obstacle adopted is a short, uniform cylindrical pipe whose upstream end is fully open, facing directly into the flow, while the downstream end is covered by a flat circular endplate containing an 'orifice' at its center. The cases are those of a completely closed and a fully open orifice; together, these bound the full range of obstructed hypersonic, stratified gas flows of interest.

Roberts, William W., Jr.↗

Space shuttle propulsion estimation development verification, volume 1

The results of the Propulsion Estimation Development Verification are summarized. A computer program developed under a previous contract (NAS8-35324) was modified to include improved models for the Solid Rocket Booster (SRB) internal ballistics, the Space Shuttle Main Engine (SSME) power coefficient model, the vehicle dynamics using quaternions, and an improved Kalman filter algorithm based on the U-D factorized algorithm. As additional output, the estimated propulsion performances, for each device are computed with the associated 1-sigma bounds. The outputs of the estimation program are provided in graphical plots. An additional effort was expended to examine the use of the estimation approach to evaluate single engine test data. In addition to the propulsion estimation program PFILTER, a program was developed to produce a best estimate of trajectory (BET). The program LFILTER, also uses the U-D factorized algorithm form of the Kalman filter as in the propulsion estimation program PFILTER. The necessary definitions and equations explaining the Kalman filtering approach for the PFILTER program, the models used for this application for dynamics and measurements, program description, and program operation are presented.

Rogers, Robert M.↗

Investigation of Conjugate Heat Transfer in Turbine Blades and Vanes

We report on work carried out to develop a 3-D coupled Finite Volume/BEM-based temperature forward/flux back (TFFB) coupling algorithm to solve the conjugate heat transfer (CHT) which arises naturally in analysis of systems exposed to a convective environment. Here, heat conduction within a structure is coupled to heat transfer to the external fluid which is convecting heat into or out of the solid structure. There are two basic approaches to solving coupled fluid structural systems. The first is a direct coupling where the solution of the different fields is solved simultaneously in one large set of equations. The second approach is a loose coupling strategy where each set of field equations is solved to provide boundary conditions for the other. The equations are solved in turn until an iterated convergence criterion is met at the fluid-solid interface. The loose coupling strategy is particularly attractive when coupling auxiliary field equations to computational fluid dynamics codes. We adopt the latter method in which the BEM is used to solve heat conduction inside a structure which is exposed to a convective field which in turn is resolved by solving the NASA Glenn compressible Navier-Stokes finite volume code Glenn-HT. The BEM code features constant and bi-linear discontinuous elements and an ILU-preconditioned GMRES iterative solver for the resulting non-symmetric algebraic set arising in the conduction solution. Interface of flux and temperature is enforced at the solid/fluid interface, and a radial-basis function scheme is used to interpolated information between the CFD and BEM surface grids. Additionally, relaxation is implemented in passing the fluxes from the conduction solution to the fluid solution. Results from a simple test example are reported.

Kassab, A. J.↗

Convective Heat Transfer in the Reusable Solid Rocket Motor of the Space Transportation System

This simulation involved a two-dimensional axisymmetric model of a full motor initial grain of the Reusable Solid Rocket Motor (RSRM) of the Space Transportation System (STS). It was conducted with CFD (computational fluid dynamics) commercial code FLUENT. This analysis was performed to: a) maintain continuity with most related previous analyses, b) serve as a non-vectored baseline for any three-dimensional vectored nozzles, c) provide a relatively simple application and checkout for various CFD solution schemes, grid sensitivity studies, turbulence modeling and heat transfer, and d) calculate nozzle convective heat transfer coefficients. The accuracy of the present results and the selection of the numerical schemes and turbulence models were based on matching the rocket ballistic predictions of mass flow rate, head end pressure, vacuum thrust and specific impulse, and measured chamber pressure drop. Matching these ballistic predictions was found to be good. This study was limited to convective heat transfer and the results compared favorably with existing theory. On the other hand, qualitative comparison with backed-out data of the ratio of the convective heat transfer coefficient to the specific heat at constant pressure was made in a relative manner. This backed-out data was devised to match nozzle erosion that was a result of heat transfer (convective, radiative and conductive), chemical (transpirating), and mechanical (shear and particle impingement forces) effects combined.

Ahmad, Rashid A.↗

Reduction of Wake-Stator Interaction Noise Using Passive Porosity

The present study was conducted to assess the potential of Passive Porosity Technology as a mechanism to reduce interaction noise in turbomachinery by reducing the fluctuating forces acting on the vane surfaces. To do so, a typical fan stator airfoil was subjected to the effects of a transversely moving wake; time histories of the primitive aerodynamic variables, obtained from Computational Fluid Dynamics (CFD) solutions, were then input into an acoustic prediction code. This procedure was performed on the solid airfoil to obtain a baseline, and on a series of porous configurations in order to isolate those that yield maximum noise reductions without compromising the aerodynamic performance of the stator. It was found that communication between regions of high pressure differential - made possible by the use of passive porosity - is necessary to significantly alter the noise radiation pattern of the stator airfoil. In general, noise reductions were obtained for those configurations incorporating passive porosity in the region between x/c is approximately 0.15 on the suction side of the airfoil and x/c is approximately 0.20 on the pressure side. Reductions in overall radiated noise of approximately 1.0 dB were obtained. The noise benefit increased to about 2.5 dB when the effects of loading noise alone were considered.

Tinetti, Ana F.↗

Atmospheric Trace Gas Abundances and Stable Isotope Ratios via IR-LIF

We propose to develop new technologies with support provided by PIDDP that will enable the in situ measurements of abundances and stable isotope ratios in important radiatively and biogenically active gases such as carbon dioxide, carbon monoxide, water, methane, nitrous oxide, and hydrogen sulfide to very high precision (0.1 per mil or better for the isotopic ratios, for example). Such measurements, impossible at present, could provide pivotal new constraints on the global (bio)geochemical budgets of these critical species, and could also be used to examine the dynamics of atmospheric transport on Mars, Titan, and other solar system bodies. We believe the combination of solid state light sources with imaging of the IR laser induced fluorescence (IR-LIF) via newly available detector arrays will make such in situ measurements possible for the first time. Even under ambient terrestrial conditions, the LIF yield from vibrational excitation of species such as water and carbon dioxide should produce emission measures well in excess of ten billion photons/sec from samples volumes of order 1 c.c. These count rates can, in principle, yield detection limits into the sub-ppt range that are required for the in situ isotopic study of atmospheric trace gases. While promising, such technologies are relatively immature, but developing rapidly, and there are a great many uncertainties regarding their applicability to in situ IR-LIF planetary studies. We therefore feel PIDDP support will be critical to developing these new tools, and propose a three-year program to combine microchip near-IR lasers with low background detection axes and state-of-the-art HgCdTe detectors developed for astronomical spectroscopy to investigate the sensitivity of IR-LIF under realistic planetary conditions, to optimize the optical pumping and filtering schemes for important species, and to apply the spectrometer to the non-destructive measurement of stable isotopes in a variety of test samples. These studies form the necessary precursors to the development of compact, lightweight stable isotope/trace gas sensors for future planetary missions.

Blake, Geoffrey A.↗

A Change of Inertia-Supporting the Thrust Vector Control of the Space Launch System

The Space Launch System (SLS) is America's next launch vehicle. To utilize the vehicle more economically, heritage hardware from the Space Transportation System (STS) will be used when possible. The Solid Rocket Booster (SRB) actuators could possibly be used in the core stage of the SLS. The dynamic characteristics of the SRB actuator will need to be tested on an Inertia Load Stand (ILS) that has been converted to Space Shuttle Main Engine (SSME). The inertia on the pendulum of the ILS will need to be changed to match the SSME inertia. In this testing environment an SRB actuator can be tested with the equivalent resistence of an SSME.

Dziubanek, Adam J.↗

Material Response Modeling of Ablative Thermal Protection Systems using PATO

Developing thermal protection systems (TPS) for future space vehicles involves an extensive design and test cycle. A key component of the cycle is determining margins for a safe design that rely on confidence in the performance of the TPS. Predicting the complicated multiphysics phenomena that occur during atmospheric entry requires high-fidelity modeling tools. To this end, the Porous-material Analysis Toolbox based on OpenFOAM (PATO) has been developed. PATO is an open-source software for computing material response of reactive porous materials submitted to high-temperature environments. Recent applications of PATO include computation of the full heatshield 3D material response from the Mars Science Laboratory atmospheric entry and ablation during arc jet testing. Current efforts are underway to loosely couple PATO with other discipline specialized codes including hypersonic computational fluid dynamics (CFD) to assess the effects of pyrolysis-gas blowing into the boundary layer and computational solid mechanics to address modeling of mechanical erosion. Surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation are also being added.

Thermal Protection Systems↗

Material Response Modeling of Ablative Thermal Protection Systems using PATO

Developing thermal protection systems (TPS) for future space vehicles involves an extensive design and test cycle. A key component of the cycle is determining margins for a safe design that rely on confidence in the performance of the TPS. Predicting the complicated multiphysics phenomena that occur during atmospheric entry requires high-fidelity modeling tools. To this end, the Porous-material Analysis Toolbox based on OpenFOAM (PATO) has been developed. PATO is an open-source software for computing material response of reactive porous materials submitted to high-temperature environments. Recent applications of PATO include computation of the full heatshield 3D material response from the Mars Science Laboratory atmospheric entry and ablation during arc jet testing. Current efforts are underway to loosely couple PATO with other discipline specialized codes including hypersonic computational fluid dynamics (CFD) to assess the effects of pyrolysis-gas blowing into the boundary layer and computational solid mechanics to address modeling of mechanical erosion. Surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation are also being added.

Thermal Protection Systems↗

Comparison of In-situ Near Infrared Melt Pool Imagery to Optical Microscopy Measurements

Additive manufacturing (AM) is a rapidly growing technology. An area of major importance is the integrity and repeatability of AM parts. The goal is to reduce obstacles to certify AM built parts to allow for use in critical aerospace applications. In-situ nondestructive evaluation sensors can be used for build assessment and can potentially play a key role in certifying AM parts. For example, melt pool features are understood to have a strong correlation to microstructural defects and the use of a near infrared (NIR) camera can be used to record the melt pool, cooling areas, and temperature gradients during the build. This work explores the use of a low cost NIR camera to obtain single line track imagery of the Ti-6Al-4V melt pools for various processing parameters. The NIR camera is radiometrically calibrated and configured in-line with the laser source to obtain high resolution imagery of the melt pool shape and dynamics. The challenge to measure melt pool shapes is to identify the transition points between the metal solid to liquid phase. Factors for melt pool measurements such as thermal camera pixel resolution, surface emissivity, and blurring due to the laser beam movement are discussed. Lastly, the melt pool imagery are compared to optical microscopy measurements for validation.

additive manufacturing↗

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO). PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Material Response↗

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO) [1,2,3]. PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) [4] CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) [5] program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield [6]. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem [7]. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Thermal Protection Systems↗

A Coupled Fluid-Structure Interaction Analysis of Solid Rocket Motor with Flexible Inhibitors

A capability to couple NASA production CFD code, Loci/CHEM, with CFDRC's structural finite element code, CoBi, has been developed. This paper summarizes the efforts in applying the installed coupling software to demonstrate/investigate fluid-structure interaction (FSI) between pressure wave and flexible inhibitor inside reusable solid rocket motor (RSRM). First a unified governing equation for both fluid and structure is presented, then an Eulerian-Lagrangian framework is described to satisfy the interfacial continuity requirements. The features of fluid solver, Loci/CHEM and structural solver, CoBi, are discussed before the coupling methodology of the solvers is described. The simulation uses production level CFD LES turbulence model with a grid resolution of 80 million cells. The flexible inhibitor is modeled with full 3D shell elements. Verifications against analytical solutions of structural model under steady uniform pressure condition and under dynamic condition of modal analysis show excellent agreements in terms of displacement distribution and eigen modal frequencies. The preliminary coupled result shows that due to acoustic coupling, the dynamics of one of the more flexible inhibitors shift from its first modal frequency to the first acoustic frequency of the solid rocket motor.

Yang, H. Q.↗

Coupled Fluid-Structure Interaction Analysis of Solid Rocket Motor with Flexible Inhibitors

Flexible inhibitors are generally used in solid rocket motors (SRMs) as a means to control the burning of propellant. Vortices generated by the flow of propellant around the flexible inhibitors have been identified as a driving source of instabilities that can lead to thrust oscillations in launch vehicles. Potential coupling between the SRM thrust oscillations and structural vibration modes is an important risk factor in launch vehicle design. As a means to predict and better understand these phenomena, a multidisciplinary simulation capability that couples the NASA production CFD code, Loci/CHEM, with CFDRC's structural finite element code, CoBi, has been developed. This capability is crucial to the development of NASA's new space launch system (SLS). This paper summarizes the efforts in applying the coupled software to demonstrate and investigate fluid-structure interaction (FSI) phenomena between pressure waves and flexible inhibitors inside reusable solid rocket motors (RSRMs). The features of the fluid and structural solvers are described in detail, and the coupling methodology and interfacial continuity requirements are then presented in a general Eulerian-Lagrangian framework. The simulations presented herein utilize production level CFD with hybrid RANS/LES turbulence modeling and grid resolution in excess of 80 million cells. The fluid domain in the SRM is discretized using a general mixed polyhedral unstructured mesh, while full 3D shell elements are utilized in the structural domain for the flexible inhibitors. Verifications against analytical solutions for a structural model under a steady uniform pressure condition and under dynamic modal analysis show excellent agreement in terms of displacement distribution and eigenmode frequencies. The preliminary coupled results indicate that due to acoustic coupling, the dynamics of one of the more flexible inhibitors shift from its first modal frequency to the first acoustic frequency of the solid rocket motor. This insight could have profound implications for SRM and flexible inhibitor designs for current and future launch vehicles including SLS.

Yang, H. Q.↗