SEARCH · Engineering Papers

Results for “finite density”

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 451 records · Page 25

All sources NASA NTRS DOE OSTI

NASA NTRS2011

Numerically Exact Computer Simulations of Light Scattering by Densely Packed, Random Particulate Media

Direct computer simulations of electromagnetic scattering by discrete random media have become an active area of research. In this progress review, we summarize and analyze our main results obtained by means of numerically exact computer solutions of the macroscopic Maxwell equations. We consider finite scattering volumes with size parameters in the range, composed of varying numbers of randomly distributed particles with different refractive indices. The main objective of our analysis is to examine whether all backscattering effects predicted by the low-density theory of coherent backscattering (CB) also take place in the case of densely packed media. Based on our extensive numerical data we arrive at the following conclusions: (i) all backscattering effects predicted by the asymptotic theory of CB can also take place in the case of densely packed media; (ii) in the case of very large particle packing density, scattering characteristics of discrete random media can exhibit behavior not predicted by the low-density theories of CB and radiative transfer; (iii) increasing the absorptivity of the constituent particles can either enhance or suppress typical manifestations of CB depending on the particle packing density and the real part of the refractive index. Our numerical data strongly suggest that spectacular backscattering effects identified in laboratory experiments and observed for a class of high-albedo Solar System objects are caused by CB.

Dlugach, Janna M.↗

NASA NTRS1989

Graphical postprocessing for 3-D mesh quality evaluation

An important objective of three-dimensional graphical finite-element postprocessing is to indicate to the engineer the accuracy of analysis results. The inclusion of mesh quality sensors permits a subjective evaluation of the adequacy of a single analysis being interpreted. For graphical approaches, both strain-energy-density gradients and discontinuities of unsmoothed responses and their gradients have proved to be effective sensors. Interactive graphical tools which can display discontinuity information effectively are described; these are essentially different from the ordinary methods used for the viewing of smoothed results.

Panthaki, M. J.↗

NASA NTRS1990

Simulation of leading-edge vortex flows

An implicit upwind-relaxation finite-difference algorithm solving the incompressible Navier-Stokes equations is employed to simulate low-speed, three-dimensional, laminar, leading-edge vortex flows over three round-edged low-aspect-ratio wings. The effects of grid density, angle of attack, Reynolds number, and wing planform on the flowfield structures and integral values are studied. Computed results are presented and compared with experimental data.

Hsu, C.-H.↗

NASA NTRS1991

The free compressible viscous vortex

The present study investigates the effects of compressibility on free (unsteady) viscous heat-conducting vortices. Analytical solutions are found in the limit of large but finite Reynolds number and small but finite Mach number. It is shown that the spreading of the vortex causes a radial flow. This flow is given by the solution of an ordinary differential equation, which gives the dependence of the radial velocity on the tangential velocity, density, and temperature profiles of the vortex. Estimates of the radial velocity found by solving this equation are found to be in good agreement with numerical solutions of the full equations. The equations for the viscous evolution are expanded in powers of Mach number to obtain detailed analytical solutions. It is shown that swirling axisymmetric compressible flows generate negative radial velocities far from the vortex core owing to viscous effects, regardless of the initial distributions of vorticity, density, and entropy.

Colonius, Tim↗

NASA NTRS1991

Fragmentation in stellar collapse

A three dimensional Cartesian code was developed to model molecular cloud collapse and protostellar formation. The basic Eulerian scheme employs spatially second-order accurate, finite difference methods to advance the fluid variables in time. The hydrodynamic equations governing the collapse include the effects of self-gravity, rotation, and radiative transfer. Supplementary equations include Poisson's equation which relates the gravitational potential to the density distribution and a mean intensity equation derived with the Eddington approximation for radiative transfer. The gas pressures and internal energies are determined from equations of state derived for a molecular cloud composed of hydrogen, helium, and heavier elements.

Myhill, E. M.↗

NASA NTRS1993

Small rocket flowfield diagnostic chambers

Instrumented and optically-accessible rocket chambers are being developed to be used for diagnostics of small rocket (less than 440 N thrust level) flowfields. These chambers are being tested to gather local fluid dynamic and thermodynamic flowfield data over a range of test conditions. This flowfield database is being used to better understand mixing and heat transfer phenomena in small rockets, influence the numerical modeling of small rocket flowfields, and characterize small rocket components. The diagnostic chamber designs include: a chamber design for gathering wall temperature profiles to be used as boundary conditions in a finite element heat flux model; a chamber design for gathering inner wall temperature and static pressure profiles; and optically-accessible chamber designs, to be used with a suite of laser-based diagnostics for gathering local species concentration, temperature, density, and velocity profiles. These chambers were run with gaseous hydrogen/gaseous oxygen (GH2/GO2) propellants, while subsequent versions will be run on liquid oxygen/hydrocarbon (LOX/HC) propellants. The purpose, design, and initial test results of these small rocket flowfield diagnostic chambers are summarized.

Morren, Sybil↗

NASA NTRS1997

The Reduction of Ducted Fan Engine Noise Via A Boundary Integral Equation Method

The development of a Boundary Integral Equation Method (BIEM) for the prediction of ducted fan engine noise is discussed. The method is motivated by the need for an efficient and versatile computational tool to assist in parametric noise reduction studies. In this research, the work in reference 1 was extended to include passive noise control treatment on the duct interior. The BEM considers the scattering of incident sound generated by spinning point thrust dipoles in a uniform flow field by a thin cylindrical duct. The acoustic field is written as a superposition of spinning modes. Modal coefficients of acoustic pressure are calculated term by term. The BEM theoretical framework is based on Helmholtz potential theory. A boundary value problem is converted to a boundary integral equation formulation with unknown single and double layer densities on the duct wall. After solving for the unknown densities, the acoustic field is easily calculated. The main feature of the BIEM is the ability to compute any portion of the sound field without the need to compute the entire field. Other noise prediction methods such as CFD and Finite Element methods lack this property. Additional BIEM attributes include versatility, ease of use, rapid noise predictions, coupling of propagation and radiation both forward and aft, implementable on midrange personal computers, and valid over a wide range of frequencies.

Tweed, J.↗

NASA NTRS2002

Numerical Modeling of Cavitating Venturi: A Flow Control Element of Propulsion System

In a propulsion system, the propellant flow and mixture ratio could be controlled either by variable area flow control valves or by passive flow control elements such as cavitating venturies. Cavitating venturies maintain constant propellant flowrate for fixed inlet conditions (pressure and temperature) and wide range of outlet pressures, thereby maintain constant, engine thrust and mixture ratio. The flowrate through the venturi reaches a constant value and becomes independent of outlet pressure when the pressure at throat becomes equal to vapor pressure. In order to develop a numerical model of propulsion system, it is necessary to model cavitating venturies in propellant feed systems. This paper presents a finite volume model of flow network of a cavitating venturi. The venturi was discretized into a number of control volumes and mass, momentum and energy conservation equations in each control volume are simultaneously solved to calculate one-dimensional pressure, density, and flowrate and temperature distribution. The numerical model predicts cavitations at the throat when outlet pressure was gradually reduced. Once cavitation starts, with further reduction of downstream pressure, no change in flowrate is found. The numerical predictions have been compared with test data and empirical equation based on Bernoulli's equation.

Majumdar, Alok↗

NASA NTRS2011

Structural Analysis of Pressurized Small Diameter Lines in a Random Vibration Environment

The pressurization and propellant feed lines for the Ares 1 Upper Stage Reaction and Roll Control Systems (ReCS and RoCS) were required to be in a high g-load random vibration flight environment. The lines connected the system components and were filled with both liquid hydrazine and gaseous helium. They are considered small and varied between one fourth to one inch in diameter. The random vibration of the lines was considered to be base excitation through the mating components and mounting hardware. It was found that reducing the amount of support structure for the lines added flexibility to the system and improved the line stresses from random vibration, but caused higher stresses from the static g-loads. The locations and number of brackets were optimized by analyzing the mode shapes of the lines causing high stresses. The use of brackets that only constrain motion in the direction of concern further reduced the stresses in the lines. Finite element analysis was used to perform the analysis. The lines were pre-stressed by temperature and internal pressure with fluid and insulation included as non-structural mass. Base excitation was added to the model using Power Spectral Density (PSD) data for the expected flight loads. The random vibration and static g-load cases were combined to obtain the total stress in the lines. This approach advances the state of the art in line analysis by using FEA to predict the stresses in the lines and to optimize the entire system based on the expected flight environment. Adding flexibility to lines has been used in piping system for temperature loads, but in flight environments flexibility has been limited for the static stresses. Adding flexibility to the system in a flight environment by reducing brackets has the benefit of reducing stresses and weight

Davis, Mark↗

NASA NTRS2019

Effect of Phenolic Matrix Microcracking on the Structural Response of a 3-D Woven Thermal Protection System

The effect of microcracking in the phenolic matrix of a three-dimensional woven thermal protection system (TPS) and the resulting material stiffness reduction was studied via a comparison of finite element analysis results from a linear analysis and an iterative linear analysis. A TPS is necessary to protect space vehicles from the aerodynamic heating of planetary entry. The Heatshield for Extreme Entry Environment Technology (HEEET) project has developed a TPS for use in high heat-flux and pressure missions. The material is a dual-layer continuous dry weave, which is then infiltrated with a low-density phenolic resin matrix to form a rigid ablator. The phenolic resin matrix does not have structural load transfer requirements, and testing has shown that the phenolic resin can fully satisfy thermal requirements when the matrix contains microcracks. Due to high stresses in the through-the-thickness direction of the material, phenolic microcracks may form in the matrix material, which would result in a reduction of stiffness. An exploratory study was conducted to determine if reduction in material stiffness would change the load paths and/or decrease the structural margins. A comparison was performed between a linear finite element analysis that did not take into account phenolic microcracking and an iterative linear finite element analysis that accounted for propagation of phenolic microcracking. Four subcases using varying assumptions were analyzed and the results indicate that the assumed strength at which the phenolic microcracking propagates was the critical parameter for determining the extent of microcracking in the phenolic matrix. Phenolic microcracking does not have an adverse effect on the structural response of the test article and is not a critical failure.

Langston, Sarah L.↗

NASA NTRS1976

A three-dimensional, time-dependent model of Mobile Bay

A three-dimensional, time-variant mathematical model for momentum and mass transport in estuaries was developed and its solution implemented on a digital computer. The mathematical model is based on state and conservation equations applied to turbulent flow of a two-component, incompressible fluid having a free surface. Thus, bouyancy effects caused by density differences between the fresh and salt water, inertia from thare river and tidal currents, and differences in hydrostatic head are taken into account. The conservation equations, which are partial differential equations, are solved numerically by an explicit, one-step finite difference scheme and the solutions displayed numerically and graphically. To test the validity of the model, a specific estuary for which scaled model and experimental field data are available, Mobile Bay, was simulated. Comparisons of velocity, salinity and water level data show that the model is valid and a viable means of simulating the hydrodynamics and mass transport in non-idealized estuaries.

Pitts, F. H.↗

NASA NTRS1981

Development of orthotropic birefringent materials for photoelastic stress analysis

Materials were selected and fabrication procedures developed for orthotropic birefringent materials. An epoxy resin (Maraset 658/558 system) was selected as the matrix material. Fibers obtained from style 3733 glass cloth and type 1062 glass roving were used as reinforcement. Two different fabrication procedures were used. In the first one, layers of unidirectional fibers removed from the glass cloth were stacked, impregnated with resin, bagged and cured in the autoclave at an elevated temperature. In the second procedure, the glass roving was drywound over metal frames, impregnated with resin and cured at room temperature under pressure and vacuum in an autoclave. Unidirectional, angle-ply and quasi-isotropic laminates of two thicknesses and with embedded flaws were fabricated. The matrix and the unidirectional glass/epoxy material were fully characterized. The density, fiber volume ratio, mechanical, and optical properties were determined. The fiber volume ratio was over 0.50. Birefringent properties were in good agreement with predictions based on a stress proportioning concept and also, with one exception, with properties predicted by a finite element analysis.

Daniel, I. M.↗

NASA NTRS1982

Application of two-point difference schemes to the conservative Euler equations for one-dimensional flows

An implicit finite difference method is presented for obtaining steady state solutions to the time dependent, conservative Euler equations for flows containing shocks. The method used the two-point differencing approach of Keller with dissipation added at supersonic points via the retarded density concept. Application of the method to the one-dimensional nozzle flow equations for various combinations of subsonic and supersonic boundary conditions shows the method to be very efficient. Residuals are typically reduced to machine zero in approximately 35 time steps for 50 mesh points. It is shown that the scheme offers certain advantages over the more widely used three-point schemes, especially in regard to application of boundary conditions.

Wornom, S. F.↗

NASA NTRS1984

Development and characterization of orthotropic-birefringent materials

Daniel, I. M.↗

NASA NTRS1984

Dynamics of the Venus ionosphere - A two-dimensional model study

A two-dimensional model of the ionosphere of Venus which simulates ionospheric dynamics by self-consistently solving the plasma equations of motion, including the inertial term, in finite difference form has been constructed. The model, which is applied over the solar zenith angle range extending from 60 to 140 deg and the altitude range 100 to 480 km, simulates the measured horizontal velocity field quite satisfactorily. The ion density field is somewhat overestimated on the dayside, because of the choice model neutral atmosphere, and underestimated on the nightside, because of setting the ionopause height at too low an altitude. It is concluded that solar photoionization on the dayside and ion recombination on the nightside are the processes mainly responsible for accelerating the plasma to the observed velocities. The plasma flow appears to be sufficient to maintain the nightside ionosphere at or near the observed median level of ion densities.

Whitten, R. C.↗

NASA NTRS1986

Interactive graphics and analysis accuracy

An important objective of graphical finite element postprocessing is the facility to indicate to the engineer the accuracy of analysis results. The inclusion of mesh quality sensors permits a subjective evaluation of the adequacy of a single analysis being interpreted. For graphical approaches, both strain energy density gradients and discontinuities of unsmoothed responses and their gradients have proved to be effective sensors. Graphical tools which can display discontinuity information effectively are described; these are essentially different from the ordinary methods used for the viewing of smoothed results.

Abel, J. F.↗

NASA NTRS1993

On the nonlinear interfacial instability of rotating core-annular flow

The interfacial stability of rotating core-annular flows is investigated. The linear and nonlinear effects are considered for the case when the annular region is very thin. Both asymptotic and numerical methods are used to solve the flow in the core and film regions which are coupled by a difference in viscosity and density. The long-term behavior of the fluid-fluid interface is determined by deriving its nonlinear evolution in the form of a modified Kuramoto-Sivashinsky equation. We obtain a generalization of this equation to three dimensions. The flows considered are applicable to a wide array of physical problems where liquid films are used to lubricate higher or lower viscosity core fluids, for which a concentric arrangement is desired. Linearized solutions show that the effects of density and viscosity stratification are crucial to the stability of the interface. Rotation generally destabilizes non-axisymmetric disturbances to the interface, whereas the centripetal forces tend to stabilize flows in which the film contains the heavier fluid. Nonlinear affects allow finite amplitude helically travelling waves to exist when the fluids have different viscosities.

Coward, Aidrian V.↗

NASA NTRS1993

Noniterative implicit method for tracking particles in mixed Lagrangian-Eulerian formulations

The existing implicit methods for the current initial value problems (IVPs) concerning particle-laden flows are complicated and iterative in nature. This paper presents a noniterative implicit method which can be used with pressure-based as well as with density-based algorithms. The method is illustrated by analyzing a dilute dispersion of noninteracting solid particles in an isothermal flow in a passage bounded by one straight wall and one wavy wall, in which all particles are spherical and have a finite velociy relative to the continuum phase at the inflow boundary.

Shih, T. I.-P.↗