Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite Volume”

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 1,477 records · Page 82

Effects of Microstructural Variability on Thermo-Mechanical Properties of a Woven Ceramic Matrix Composite

The objectives of this paper include identifying important architectural parameters that describe the SiC/SiC five-harness satin weave composite and characterizing the statistical distributions and correlations of those parameters from photomicrographs of various cross sections. In addition, realistic artificial cross sections of a 2D representative volume element (RVE) are generated reflecting the variability found in the photomicrographs, which are used to determine the effects of architectural variability on the thermo-mechanical properties. Lastly, preliminary information is obtained on the sensitivity of thermo-mechanical properties to architectural variations. Finite element analysis is used in combination with a response surface and it is shown that the present method is effective in determining the effects of architectural variability on thermo-mechanical properties.

Goldsmith, Marlana B.↗

Modeling Pore Closure in the Hot Isostatic Pressing of Additively Manufactured Inconel 718 Samples

Additive manufacturing provides opportunity for a new world of possibilities for metallic component design and fabrication; however, as-built part quality is notoriously variable, anisotropic, and degraded by porous defects. To improve quality, parts are commonly post processed using a variety of techniques including hot isostatic pressing (HIP). During HIP, high temperature and pressure is employed to improve microstructure and reduce pore volume by inducing visco-plastic deformation in the metallic substrate. Validated computational capabilities can provide predictions to guide the selection of processing condition, expected pore reduction, and the effects of varying entrapped gas in the pores. This discussion describes a micro-scale finite element structural model developed to predict pore closure due to HIP in as-built additively manufactured Inconel 718 samples produced by the laser powder bed fusion process. The model accounts for plastic deformation, creep, and internal pressure for a single pore. Pores are characterized in physical samples before and after HIP, and model predictions are compared to this measured data.

defect model prediction↗

Interplanetary Program to Optimize Simulated Trajectories (IPOST). Volume 1: User's guide

IPOST is intended to support many analysis phases, from early interplanetary feasibility studies through spacecraft development and operations. The IPOST output provides information for sizing and understanding mission impacts related to propulsion, guidance, communications, sensor/actuators, payload, and other dynamic and geometric environments. IPOST models three degree of freedom trajectory events, such as launch/ascent, orbital coast, propulsive maneuvering (impulsive and finite burn), gravity assist, and atmospheric entry. Trajectory propagation is performed using a choice of Cowell, Encke, Multiconic, Onestep, or Conic methods. The user identifies a desired sequence fo trajectory events, and selects which parameters are independent (controls) and dependent (targets), as well as other constraints and the coat function. Targeting and optimization is performed using the Stanford NPSOL algorithm. IPOST structure allows sub-problems within a master optimization problem to aid in the general constrained parameter optimization solution. An alternate optimization method uses implicit simulation and collocation techniques.

Hong, P. E.↗

Heat transfer, thermal control, and heat pipes

This volume provides information on recent progress in spacecraft thermal control and the supporting disciplines of conduction, thermal radiation, and heat pipe theory and application. Four problem areas are considered: conduction heat transfer, radiation heat transfer, thermal control, and heat pipes. The topics covered include finite-element methodology for transient conduction/forced-convection thermal analysis; effects of surface finish on thermal contact resistance between different materials; mathematical models for wide-band nongray gas radiation in spherical and cylindrical geometries; thermal design, analysis and testing of the Shuttle remote manipulator arm; porous heat pipe; and transient behavior of liquid trap heat-pipe thermal diodes. Also discussed is the thermal design concept for a high-resolution UV spectrometer.

Olstad, W. B.↗

Preliminary Soot Computations Based on a Model Aircraft Combustor With OpenNCC

In a preliminary effort to extend our current capabilities used in the design & testing of aero-engine combustors to particulate emissions, we incorporated several existing 2- equation soot models into the Open National Combustion Code (OpenNCC). The paper presents results from a validation study based on a well-characterized, swirl-stabilized, turbulent sooting flame generated by a model aero-engine combustor developed at DLR, Germany. The calculations are based on a Reynolds-Averaged Navier Stokes (RANS) simulation with standard k-e turbulence model & a finite-rate chemical kinetic mechanism. A separate transport equation is solved for all individual species involved in ethylene/air combustion. The predicted velocity and temperature are mostly in good agreement with the measured data and reasonable agreement was found for soot. While the calculated SVF (Soot Volume Fraction) behavior is similar to some other predicted results reported in literature, the measured soot lasts farther downstream than the calculated soot indicating a need for further improvement in the soot models employed.

Raju, Manthena↗

Vibrational Power Flow Analysis of Rods and Beams

A new method to model vibrational power flow and predict the resulting energy density levels in uniform rods and beams is investigated. This method models the flow of vibrational power in a manner analogous to the flow of thermal power in a heat conduction problem. The classical displacement solutions for harmonically excited, hysteretically damped rods and beams are used to derive expressions for the vibrational power flow and energy density in the rod and beam. Under certain conditions, the power flow in these two structural elements will be shown to be proportional to the energy density gradient. Using the relationship between power flow and energy density, an energy balance on differential control volumes in the rod and beam leads to a Poisson's equation which models the energy density distribution in the rod and beam. Coupling the energy density and power flow solutions for rods and beams is also discussed. It is shown that the resonant behavior of finite structures complicates the coupling of solutions, especially when the excitations are single frequency inputs. Two coupling formulations are discussed, the first based on the receptance method, and the second on the travelling wave approach used in Statistical Energy Analysis. The receptance method is the more computationally intensive but is capable of analyzing single frequency excitation cases. The traveling wave approach gives a good approximation of the frequency average of energy density and power flow in coupled systems, and thus, is an efficient technique for use with broadband frequency excitation.

Wohlever, James Christopher↗

A Hybrid Transfer Function Procedure for Broadband Auralization Within Small Flight Vehicle Interiors

Auralization of the sound fields within flight vehicle interiors is of interest to aircraft designers considering crew and passenger ride quality. Computer models can be used to predict the filtering effects that the fuselage structure and interior volume have on the exterior sound field in order to simulate the resulting cabin interior noise. However, this can become expensive when a wide frequency range of analysis suitable for passenger auralization is desired. This is particularly the case when using deterministic modeling methods such as finite element analysis (FEA) wherein prohibitively high levels of detail and finite element discretizations are required to resolve the wavelengths at high frequencies. On the contrary, Statistical Energy Analysis (SEA), while not capturing important modal behaviors at low frequencies, is better suited for higher frequency analyses where structural and interior dynamics exhibit high modal overlap and spatially averaged physical quantities suffice. For this reason, a hybrid procedure is considered here that utilizes FE modal solutions at lower frequencies and SEA results at higher frequencies. A transfer function dataset containing the exterior to interior vibroacoustic filtering effect is calculated accordingly and subsequently stored to be used for a variety of exterior dynamic load cases. At low frequencies, the transfer function dataset detail is related to the level of FE model refinement whereas at high frequencies it functions as an equalizer array with a level of detail corresponding to the SEA subsystem partitioning scheme. The above procedure is described herein and demonstrated on a six passenger flight vehicle with eight propulsors in a hover flight condition.

auralization↗

A Hybrid Transfer Function Procedure for Broadband Auralization Within Small Flight Vehicle Interiors

Auralization of the sound fields within flight vehicle interiors is of interest to aircraft designers considering crew and passenger ride quality. Computer models can be used to predict the filtering effects that the fuselage structure and interior volume have on the exterior sound field in order to simulate the resulting cabin interior noise. However, this can become expensive when a wide frequency range of analysis suitable for passenger auralization is desired. This is particularly the case when using deterministic modeling methods such as finite element analysis (FEA) wherein prohibitively high levels of detail and finite element discretizations are required to resolve the wavelengths at high frequencies. On the contrary, Statistical Energy Analysis (SEA), while not capturing important modal behaviors at low frequencies, is better suited for higher frequency analyses where structural and interior dynamics exhibit high modal overlap and spatially averaged physical quantities suffice. For this reason, a hybrid procedure is considered here that utilizes FE modal solutions at lower frequencies and SEA results at higher frequencies. A transfer function dataset containing the exterior to interior vibroacoustic filtering effect is calculated accordingly and subsequently stored to be used for a variety of exterior dynamic load cases. At low frequencies, the transfer function dataset detail is related to the level of FE model refinement whereas at high frequencies it functions as an equalizer array with a level of detail corresponding to the SEA subsystem partitioning scheme. The above procedure is described herein and demonstrated on a six passenger flight vehicle with eight propulsors in a hover flight condition.

auralization↗

Viscel: A general purpose computer program for analysis of linear viscoelastic structures, volume 2

The VISCEL program is a general purpose computer program developed for equilibrium analysis of linear viscoelastic structures. The program is written in FORTRAN 5 language to operate on the Univac 1108 computer under the EXEC 8 operating system. The program, an extension of the linear equilibrium problem solver ELAS, is an updated and extended version of its earlier form written for the IBM 7094 computer. Finite element matrix displacement approach coupled with the synchronized material property concept, utilizing incremental time steps, was adopted for the solution presented. The step-by-step procedure involves solution of recursive equations in the time domain, which takes into account the memory of material properties. Incremental and accumulative displacements and stresses are obtained at the end of each time step. In order to minimize the extent of computations resulting from accumulative effects of material memory, the program provides an option which enables the employment of constant time steps in the logarithmic scale. Program documentation is presented.

Gupta, K. K.↗

Modeling of Fillets in Thin-Walled Structures for Dynamic Analysis

Relatively new developments in manufacturing methods have made it possible to produce machined parts with a wall thickness of less than 0.010". While the parts are being machined. Fillets of relatively large radius to thickness ratios are created by the milling tool. This paper discusses an accurate new technique for finite-element modeling of fillets for dynamic and stiffness analysis using efficient plate elements rather than more computationally intensive, high density solid meshes, A simple filleted cantilever beam of 0.040" thickness, 1.6" length, 1" depth, and 0.250" radius in the corner was modeled using both solid elements and several different plate element geometries. The finite element results were then compared with static and modal testing of a machined sample. The highest-scoring plate element technique uses a bridge of elements that is constructed through the tangent point of the Fillet radius such that the thickness of the element matches tile volume of the fillet. This model produces errors of less than 4.6 percent for static loading and less than 7.4 percent for modal analysis. This simple plate element technique will prove critical for efficient, timely, and accurate dynamic analysis of complex thin-walled structures.

Seugling, Richard M.↗

The Effect of Fiber Strength Stochastics and Local Fiber Volume Fraction on Multiscale Progressive Failure of Composites

Continuous fiber unidirectional polymer matrix composites (PMCs) can exhibit significant local variations in fiber volume fraction as a result of processing conditions that can lead to further local differences in material properties and failure behavior. In this work, the coupled effects of both local variations in fiber volume fraction and the empirically-based statistical distribution of fiber strengths on the predicted longitudinal modulus and local tensile strength of a unidirectional AS4 carbon fiber/ Hercules 3502 epoxy composite were investigated using the special purpose NASA Micromechanics Analysis Code with Generalized Method of Cells (MAC/GMC); local effective composite properties were obtained by homogenizing the material behavior over repeating units cells (RUCs). The predicted effective longitudinal modulus was relatively insensitive to small (~8%) variations in local fiber volume fraction. The composite tensile strength, however, was highly dependent on the local distribution in fiber strengths. The RUC-averaged constitutive response can be used to characterize lower length scale material behavior within a multiscale analysis framework that couples the NASA code FEAMAC and the ABAQUS finite element solver. Such an approach can be effectively used to analyze the progressive failure of PMC structures whose failure initiates at the RUC level. Consideration of the effect of local variations in constituent properties and morphologies on progressive failure of PMCs is a central aspect of the application of Integrated Computational Materials Engineering (ICME) principles for composite materials.

multiscale analysis↗

Design and Development of NEA Scout Solar Sail Deployer Mechanism

The 6U (approximately10cm x 20cm x 30cm) cubesat Near Earth Asteroid (NEA) Scout, projected for launch in September 2018 aboard the maiden voyage of the Space Launch System (SLS), will utilize a solar sail as its main method of propulsion throughout its approximately 3 year mission to a near earth asteroid. Due to the extreme volume constraints levied onto the mission, an acutely compact solar sail deployment mechanism has been designed to meet the volume and mass constraints, as well as provide enough propulsive solar sail area and quality in order to achieve mission success. The design of such a compact system required the development of approximately half a dozen prototypes in order to identify unforeseen problems and advance solutions. Though finite element analysis was performed during this process in an attempt to quantify forces present within the mechanism during deployment, both the boom and the sail materials do not lend themselves to achieving high-confidence results. This paper focuses on the obstacles of developing a solar sail deployment mechanism for such an application and the lessons learned from a thorough development process. The lessons presented here will have significant applications beyond the NEA Scout mission, such as the development of other deployable boom mechanisms and uses for gossamer-thin films in space.

Sobey, Alexander R.↗

Rotational cars application to simultaneous and multiple-point temperature and concentration determination in a turbulent flow

Coherent Anti-stokes Raman Scattering (CARS) from the pure rotational Raman lines of N2 is employed to measure the instantaneous rotational temperature of N2 gas at room temperature and below with good spatial resolution. A broad-bandwidth dye laser is used to obtain the entire rotational spectrum from a signal laser pulse; the CARS signal is then dispersed by a spectrograph and recorded on an optical multichannel analyzer. A best-fit temperature is found in several seconds with the aid of a computer for each experimental spectrum by a least squares comparison with calculated spectra. The model used to calculate the theoretical spectra incorporates the temperature and pressure dependence of the pressure-broadened rotational Raman lines, includes the nonresonant background susceptibility, and assumes that the pump laser has a finite linewidth. Temperatures are fit to experimental spectra recorded over the temperature range of 135 to 296K, and over the pressure range of 0.13 to 15.3 atm. In addition to the spatially resolved single point work, we have used multipoint CARS to obtain information from many spatially resolved volume elements along a cylindrical line (0.1 x 0.1 x 2.0 mm). We also obtained qualitative information on the instantaneous species concentration and temperature at 20 spatially resolved volume elements (0.1 x 0.1 x 0.1 mm) along a line.

Snow, J. B.↗

Molecular-dynamics Simulation-based Cohesive Zone Representation of Intergranular Fracture Processes in Aluminum

A traction-displacement relationship that may be embedded into a cohesive zone model for microscale problems of intergranular fracture is extracted from atomistic molecular-dynamics simulations. A molecular-dynamics model for crack propagation under steady-state conditions is developed to analyze intergranular fracture along a flat 99 [1 1 0] symmetric tilt grain boundary in aluminum. Under hydrostatic tensile load, the simulation reveals asymmetric crack propagation in the two opposite directions along the grain boundary. In one direction, the crack propagates in a brittle manner by cleavage with very little or no dislocation emission, and in the other direction, the propagation is ductile through the mechanism of deformation twinning. This behavior is consistent with the Rice criterion for cleavage vs. dislocation blunting transition at the crack tip. The preference for twinning to dislocation slip is in agreement with the predictions of the Tadmor and Hai criterion. A comparison with finite element calculations shows that while the stress field around the brittle crack tip follows the expected elastic solution for the given boundary conditions of the model, the stress field around the twinning crack tip has a strong plastic contribution. Through the definition of a Cohesive-Zone-Volume-Element an atomistic analog to a continuum cohesive zone model element - the results from the molecular-dynamics simulation are recast to obtain an average continuum traction-displacement relationship to represent cohesive zone interaction along a characteristic length of the grain boundary interface for the cases of ductile and brittle decohesion. Keywords: Crack-tip plasticity; Cohesive zone model; Grain boundary decohesion; Intergranular fracture; Molecular-dynamics simulation

Yamakov, Vesselin I.↗

Scattering and Diffraction of Electromagnetic Radiation: An Effective Probe to Material Structure

Scattered electromagnetic waves from material bodies of different forms contain, in an intricate way, precise information on the intrinsic, geometrical and physical properties of the objects. Scattering theories, ever deepening, aim to provide dependable interpretation and prediction to the complicated interaction of electromagnetic radiation with matter. There are well-established multiple-scattering formulations based on classical electromagnetic theories. An example is the Generalized Multi-particle Mie-solution (GMM), which has recently been extended to a special version ̶ the GMM-PA approach, applicable to finite periodic arrays consisting of a huge number (e.g., >>106) of identical scattering centers [1]. The framework of the GMM-PA is nearly complete. When the size of the constituent unit scatterers becomes considerably small in comparison with incident wavelength, an appropriate array of such small element volumes may well be a satisfactory representation of a material entity having an arbitrary structure. X-ray diffraction is a powerful characterization tool used in a variety of scientific and technical fields, including material science. A diffraction pattern is nothing more than the spatial distribution of scattered intensity, determined by the distribution of scattering matter by way of its Fourier transform [1]. Since all linear dimensions entered into Maxwell's equations are normalized by wavelength, an analogy exists between optical and X-ray diffraction patterns. A large set of optical diffraction patterns experimentally obtained can be found in the literature [e.g., 2,3]. Theoretical results from the GMM-PA have been scrutinized using a large collection of publically accessible, experimentally obtained Fraunhofer diffraction patterns. As far as characteristic structures of the patterns are concerned, theoretical and experimental results are in uniform agreement; no exception has been found so far. Closely connected with the spatial distribution of scattered intensities are cross sections, such as for extinction, scattering, absorption, and radiation pressure, as a critical type of key quantity addressed in most theoretical and experimental studies of radiative scattering. Cross sections predicted from different scattering theories are supposed to be in general agreement. For objects of irregular shape, the GMM-PA solutions can be compared with the highly flexible Discrete Dipole Approximation (DDA) [4,5] when dividing a target to no more than ~106 unit cells. Also, there are different ways to calculate the cross sections in the GMM-PA, providing an additional means to examine the accuracy of the numerical solutions and to unveil potential issues concerning the theoretical formulations and numerical aspects. To solve multiple scattering by an assembly of material volumes through classical theories such as the GMM-PA, the radiative properties of the component scatterers, the complex refractive index in particular, must be provided as input parameters. When using a PA to characterize a material body, this involves the use of an adequate theoretical tool, an effective medium theory, to connect Maxwell's phenomenogical theory with the atomistic theory of matter. In the atomic theory, one regards matter as composed of interacting particles (atoms and molecules) embedded in the vacuum [6]. However, the radiative properties of atomic-scaled particles are known to be substantially different from bulk materials. Intensive research efforts in the fields of cluster science and nanoscience attempt to bridge the gap between bulk and atom and to understand the transition from classical to quantum physics. The GMM-PA calculations, which place virtually no restriction on the component-particle size, might help to gain certain insight into the transition.

Xu, Yu-Lin↗

Unified Aeroacoustics Analysis for High Speed Turboprop Aerodynamics and Noise: Development of Theory for Blade Loading, Wakes, and Noise - Volume 1

A unified theory for the aerodynamics and noise of advanced turboprops are presented. Aerodynamic topics include calculation of performance, blade load distribution, and non-uniform wake flow fields. Blade loading can be steady or unsteady due to fixed distortion, counter-rotating wakes, or blade vibration. The aerodynamic theory is based on the pressure potential method and is therefore basically linear. However, nonlinear effects associated with finite axial induction and blade vortex flow are included via approximate methods. Acoustic topics include radiation of noise caused by blade thickness, steady loading (including vortex lift), and unsteady loading. Shielding of the fuselage by its boundary layer and the wing are treated in separate analyses that are compatible but not integrated with the aeroacoustic theory for rotating blades.

Hanson, D. B.↗

Reliability-based failure analysis of brittle materials

The reliability of brittle materials under a generalized state of stress is analyzed using the Batdorf model. The model is modified to include the reduction in shear due to the effect of the compressive stress on the microscopic crack faces. The combined effect of both surface and volume flaws is included. Due to the nature of fracture of brittle materials under compressive loading, the component is modeled as a series system in order to establish bounds on the probability of failure. A computer program was written to determine the probability of failure employing data from a finite element analysis. The analysis showed that for tensile loading a single crack will be the cause of total failure but under compressive loading a series of microscopic cracks must join together to form a dominant crack.

Powers, Lynn M.↗

Interplanetary program to optimize simulated trajectories (IPOST). Volume 4: Sample cases

The Interplanetary Program to Optimize Simulated Trajectories (IPOST) is intended to support many analysis phases, from early interplanetary feasibility studies through spacecraft development and operations. The IPOST output provides information for sizing and understanding mission impacts related to propulsion, guidance, communications, sensor/actuators, payload, and other dynamic and geometric environments. IPOST models three degree of freedom trajectory events, such as launch/ascent, orbital coast, propulsive maneuvering (impulsive and finite burn), gravity assist, and atmospheric entry. Trajectory propagation is performed using a choice of Cowell, Encke, Multiconic, Onestep, or Conic methods. The user identifies a desired sequence of trajectory events, and selects which parameters are independent (controls) and dependent (targets), as well as other constraints and the cost function. Targeting and optimization are performed using the Standard NPSOL algorithm. The IPOST structure allows sub-problems within a master optimization problem to aid in the general constrained parameter optimization solution. An alternate optimization method uses implicit simulation and collocation techniques.

Hong, P. E.↗