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.

SourceResearch connectionFrom yearThrough year

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

At least 1,405 records · Page 78

All sources NASA NTRS DOE OSTI

NASA NTRS1985

Computation of transonic vortex flows past delta wings Integral equation approach

The steady full-potential equation is written in the form of Poisson's equation, and the solution of the velocity field is expressed in terms of an integral equation. The solution consists of a surface integral of vorticity distribution on the wing and its free-vortex sheets and a volume integral of source distribution within a volume around the wing and its free-vortex sheets. The solution is obtained through successive iteration cycles. The source distribution is computed by using a mixed finite-difference scheme of the Murman-Cole type. The method is applied to delta wings. Numerical examples show that a conical shock is captured on the suction side of the wing. It is attached to the lower surface of the leading-edge vortex but does not necessarily reach to the wing surface.

Kandil, O. A.↗

NASA NTRS1996

Numerical Investigation of Slag Behavior for RSRM

It is known that the flow field of the redesigned solid rocket motor (RSRM) is very complicated due to the complex characteristics of turbulent multi-phase flow, chemical reaction, particle combustion, evaporation, breakup and agglomeration etc. It requires multi-phase calculations, chemical reaction simulation, and particle combustion, evaporation, and breakup models to obtain a better understanding of thermophysics for the RSRM design using numerical methods. Also, the slag buildup due to the molten particles is another factor affecting the performance of the RSRM. To achieve this goal, the volume of fluid (VOF) method is used to capture the free surface motion so as to simulate the accumulation of the molten particles (slag) of the RSRM. A finite rate chemistry model is used to simulate the chemical reaction effects. For multi-phase calculations, the Hermsen combustion model is used for the aluminum particle combustion analysis and the Taylor Analogy Breakup (TAB) model is used for the particle breakup analysis. An interphase mas-exchange model introduced by Spalding is used for the evaporation calculation. The particle trajectories are calculated using a one-step implicit method for several groups of particle sizes by which the drag forces and heat fluxes are then coupled with the gas phase equations. The preliminary results predicted a reasonable physical simulation of the particle effects using a simple two dimensional solid rocket motor configuration. It shows that the AL/AL2O3 particle sizes are reduced due to the combustion, evaporation, and breakup. The flow field is disturbed by the particles. Mach number distributions in the nozzle are deformed due to the effect of particle concentrations away from the center line.

Liaw, P.↗

NASA NTRS1996

Progressive Damage Analysis of Laminated Composite (PDALC)-A Computational Model Implemented in the NASA COMET Finite Element Code

A method for analysis of progressive failure in the Computational Structural Mechanics Testbed is presented in this report. The relationship employed in this analysis describes the matrix crack damage and fiber fracture via kinematics-based volume-averaged variables. Damage accumulation during monotonic and cyclic loads is predicted by damage evolution laws for tensile load conditions. The implementation of this damage model required the development of two testbed processors. While this report concentrates on the theory and usage of these processors, a complete list of all testbed processors and inputs that are required for this analysis are included. Sample calculations for laminates subjected to monotonic and cyclic loads were performed to illustrate the damage accumulation, stress redistribution, and changes to the global response that occur during the load history. Residual strength predictions made with this information compared favorably with experimental measurements.

Lo, David C.↗

NASA NTRS1998

Progressive Damage Analysis of Laminated Composite (PDALC) (A Computational Model Implemented in the NASA COMET Finite Element Code)

A method for analysis of progressive failure in the Computational Structural Mechanics Testbed is presented in this report. The relationship employed in this analysis describes the matrix crack damage and fiber fracture via kinematics-based volume-averaged damage variables. Damage accumulation during monotonic and cyclic loads is predicted by damage evolution laws for tensile load conditions. The implementation of this damage model required the development of two testbed processors. While this report concentrates on the theory and usage of these processors, a complete listing of all testbed processors and inputs that are required for this analysis are included. Sample calculations for laminates subjected to monotonic and cyclic loads were performed to illustrate the damage accumulation, stress redistribution, and changes to the global response that occurs during the loading history. Residual strength predictions made with this information compared favorably with experimental measurements.

Coats, Timothy W.↗

NASA NTRS2007

High-LET Patterns of DSBs in DNA Loops, the HPRT Gene and Phosphorylation Foci

We present new results obtained with our model based on the track structure and chromatin geometry that predicts the DSB spatial and genomic distributions in a cell nucleus with the full genome represented. The model generates stochastic patterns of DSBs in the physical space of the nucleus filled with the realistic configuration of human chromosomes. The model was re-used to find the distribution of DSBs in a physical volume corresponding to a visible phosphorylation focus believed to be associated with a DSB. The data shows whether there must more than one DSB per foci due to finite size of the visible focus, even if a single DSB is radiochemically responsible for the phosphorylation of DNA in its vicinity. The same model can predict patterns of closely located DSBs in a given gene, or in a DNA loop, one of the large-scale chromatin structures. We demonstrated for the example of the HPRT gene, how different sorts of radiation lead to proximity effect in DSB locations, which is important for modeling gene deletions. The spectrum of intron deletions and total gene deletions was simulated for the HPRT gene. The same proximity effect of DSBs in a loop can hinder DSB restitutions, as parts of the loop between DSBs is deleted with a higher likelihood. The distributions of DSBs and deletions of DNA in a loop are presented.

Ponomarev, Artem L.↗

NASA NTRS1985

Study for the optimization of a transport aircraft wing for maximum fuel efficiency. Volume 1: Methodology, criteria, aeroelastic model definition and results

Work performed in the design of a transport aircraft wing for maximum fuel efficiency is documented with emphasis on design criteria, design methodology, and three design configurations. The design database includes complete finite element model description, sizing data, geometry data, loads data, and inertial data. A design process which satisfies the economics and practical aspects of a real design is illustrated. The cooperative study relationship between the contractor and NASA during the course of the contract is also discussed.

Radovcich, N. A.↗

NASA NTRS1976

Numerical simulation of low gravity draining

A boundary value problem was solved numerically for a liquid that is assumed to be inviscid and incompressible, having a motion that is irrotational and axisymmetric, and having a constant (5 degrees) solid-liquid contact angle. The avoidance of excessive mesh distortion, encountered with strictly Lagrangian or Eulerian kinematics, was achieved by introducing an auxiliary kinematic velocity field along the free surface in order to vary the trajectories used in integrating the ordinary differential equations simulating the moving boundary. The computation of the velocity potential was based upon a nonuniform triangular mesh which was automatically revised to varying depths to accommodate the motion of the free surface. These methods permitted calculation of draining induced axisymmetric slosh through the many (or fractional) finite amplitude oscillations that can occur depending upon the balance of draining, gravitational, and surface tension forces. Velocity fields, evolution of the free surface with time, and liquid residual volumes were computed for three and one half decades of Weber number and for two Bond numbers, tank fill levels, and drain radii. Comparisons with experimental data are very satisfactory.

Bizzell, G. D.↗

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 NTRS1984

Development and characterization of orthotropic-birefringent materials

Daniel, I. M.↗

NASA NTRS1989

Preliminary design of a large tetrahedral truss/hexagonal heatshield panel aerobrake

An aerobrake structural concept is introduced which consists of two primary components: (1) a lightweight erectable tetrahedral support truss; and (2) sandwich hexagonal heatshield panels which, when attached to the truss, form a continuous impermeable aerobraking surface. Generic finite element models and a general analysis procedure to design tetrahedral truss/hexagonal heatshield panel aerobrakes is developed, and values of the aerobrake design parameters which minimize mass and packaging volume for a 120-foot-diameter aerobrake are determined. Sensitivity of the aerobrake design to variations in design parameters is also assessed. The results show that a 120-foot-diameter aerobrake is viable using the concept presented (i.e., the aerobrake mass is less than or equal to 15 percent of the payload spacecraft mass). Minimizing the aerobrake mass (by increasing the number of rings in the support truss) however, leads to aerobrakes with the highest part count.

Dorsey, John T.↗

NASA NTRS1998

Effects of regenerator geometry on pulse tube refrigerator performance

This paper gives results of the cooling performance of a double-inlet pulse tube refrigerator using various regenerators. The same pulse tube was used for all the experiments and measured 4.76 mm in diameter and 46.2 mm in length. A commercial linear compressor with a swept volume of 4 cm3 was used in these experiments. The operating conditions were held constant at a mean pressure of 2.0 MPa and a frequency of 54 Hz. Using finite difference software called REGEN3.1, developed at NIST, and recent experiment results, we optimized a series of regenerators based on dimensions, materials and screen packing. The values used for calculating the thermal conduction through stacked screens by REGEN3.1 were based on recent experimental results from NIST. The regenerator tubes were designed using 316 stainless steel and titanium materials. The regenerator matrices investigated were 400-mesh and 500-mesh stainless steel screen. The valve settings for both orifices were adjusted to minimize the no-load temperature for all regenerators. A cooling capacity curve from 0 to 3 W was also determined. The performance of the pulse tube refrigerator using the different regenerators is discussed. The experimental results from the various regenerators are evaluated and compared with their corresponding numerically calculated coefficient of performance (COP) and regenerator design as determined by REGEN3.1.

NASA Discipline Life Sciences Technologies↗

NASA NTRS2019

Nanopore-Templated Silver Nanoparticle Arrays Photopolymerized in Zero-Mode Wavelengths

In situ fabrication of nanostructures within a solid-polymer electrolyte confined to subwavelength-diameter nanoapertures is a promising approach for producing nanomaterials for nanophotonic and chemical sensing applications. The solid-polymerelectrolyte can be patterned by lithographic photopolymerization of poly(ethyleneglycol) diacrylate (PEGDA)-based silver cation (Ag+)-containing polyelectrolyte. Here, we present a new method for fabricating nanopore-templated Ag nanoparticle (AgNP)arrays by in situ photopolymerization using a zero-mode waveguide (ZMW) array to simultaneously template embedded AgNPs and control the spatial distribution of the optical field used for photopolymerization. The approach starts with an array of nanopores fabricated by sequential layer-by-layer deposition and focused ion beam milling. These structures have an optically transparent bottom, allowing access of the optical radiation to the attoliter-volume ZMW region to photopolymerize a PEGDA monomer solution containing AgNPs and Ag+. The electric field intensity distribution is calculated for various ZMW optical cladding layer thicknesses using finite-element simulations, closely following the light-blocking efficiency of the optical cladding layer. The fidelity of the polyelectrolyte nanopillar pattern was optimized with respect to experimental conditions, including the presence or absence of Ag+ and AgNPs and the concentrations of PEGDA and Ag+. The self-templated approach for photo patterning high-resolution photolabile polyelectrolyte nanostructures directly within a ZMW array could lead to a new class of metamaterials formed by embedding metal nanoparticles within a dielectric in a well-defined spatial array.

Zero-mode waveguide↗

NASA NTRS1977

Thermal expansion of an epoxy-glass microsphere composite

The thermal expansion of a composite of epoxy (diglycidyl ether of bisphenol A) and solid glass microspheres was investigated. The microspheres had surfaces which were either untreated or treated with a silicone release agent, an epoxy coupling agent, or a general purpose silane coupling agent. Both room temperature (about 300 K) and elevated temperature (about 475 K) cures were used for the epoxy. Two microsphere size ranges were used, about 50 microns, which is applicable in filled moldings, and about 125 microns, which is applicable as bond line spacers. The thermal expansion of the composites was measured from 300 to 350 K or from 300 to 500 K, depending on the epoxy cure temperature. Measurements were made on composites containing up to .6 volume fraction microspheres. Two predictive models, which required only the values of thermal expansion of the polymer and glass and their specific gravities, were tested against the experimental data. A finite element analysis was made of the thermal strain of a composite cell containing a single microsphere surrounded by a finite-thickness interface.

Price, H. L.↗

NASA NTRS1990

Metal matrix composites microfracture: Computational simulation

Fiber/matrix fracture and fiber-matrix interface debonding in a metal matrix composite (MMC) are computationally simulated. These simulations are part of a research activity to develop computational methods for microfracture, microfracture propagation and fracture toughness of the metal matrix composites. The three-dimensional finite element model used in the simulation consists of a group of nine unidirectional fibers in three by three unit cell array of SiC/Ti15 metal matrix composite with a fiber volume ration of 0.35. This computational procedure is used to predict the fracture process and establish the hierarchy of fracture modes based on strain energy release rate. It is also used to predict stress redistribution to surrounding matrix-fibers due to initial and progressive fracture of fiber/matrix and due to debonding of fiber-matrix interface. Microfracture results for various loading cases such as longitudinal, transverse, shear and bending are presented and discussed. Step-by-step procedures are outlined to evaluate composite microfracture for a given composite system.

Mital, Subodh K.↗

NASA NTRS1990

Metal matrix composites microfracture - Computational simulation

Mital, S. K.↗

NASA NTRS1972

VISCEL: A general-purpose computer program for analysis of linear viscoelastic structures (user's manual), volume 1

This program, an extension of the linear equilibrium problem solver ELAS, is an updated and extended version of its earlier form (written in FORTRAN 2 for the IBM 7094 computer). A synchronized material property concept utilizing incremental time steps and the finite element matrix displacement approach has been adopted for the current analysis. A special option enables employment of constant time steps in the logarithmic scale, thereby reducing computational efforts resulting from accumulative material memory effects. A wide variety of structures with elastic or viscoelastic material properties can be analyzed by VISCEL. The program is written in FORTRAN 5 language for the Univac 1108 computer operating under the EXEC 8 system. Dynamic storage allocation is automatically effected by the program, and the user may request up to 195K core memory in a 260K Univac 1108/EXEC 8 machine. The physical program VISCEL, consisting of about 7200 instructions, has four distinct links (segments), and the compiled program occupies a maximum of about 11700 words decimal of core storage.

Gupta, K. K.↗

NASA NTRS1993

Analysis of microstrip dipoles and slots transversely coupled to a microstrip line using the FDTD method

Printed dipole elements and their complement, linear slots, are elementary radiators that have found use in low-profile antenna arrays. Low-profile antenna arrays, in addition to their small size and low weight characteristics, offer the potential advantage of low-cost, high-volume production with easy integration with active integrated circuit components. The design of such arrays requires that the radiation and impedance characteristics of the radiating elements be known. The FDTD (Finite-Difference Time-Domain) method is a general, straight-forward implementation of Maxwell's equations and offers a relatively simple way of analyzing both printed dipole and slot elements. Investigated in this work is the application of the FDTD method to the analysis of printed dipole and slot elements transversely coupled to an infinite transmission line in a multilayered configuration. Such dipole and slot elements may be used in dipole and slot series-fed-type linear arrays, where element offsets and interelement line lengths are used to obtain the desired amplitude distribution and beam direction, respectively. The design of such arrays is achieved using transmission line theory with equivalent circuit models for the radiating elements. In an equivalent circuit model, the dipole represents a shunt impedance to the transmission line, where the impedance is a function of dipole offset, length, and width. Similarly, the slot represents a series impedance to the transmission line. The FDTD method is applied to single dipole and slot elements transversely coupled to an infinite microstrip line using a fixed rectangular grid with Mur's second order absorbing boundary conditions. Frequency-dependent circuit and scattering parameters are obtained by saving desired time-domain quantities and using the Fourier transform. A Gaussian pulse excitation is applied to the microstrip transmission line, where the resulting reflected signal due to the presence of the radiating element is used to determine the equivalent element impedance.

Tulintseff, A. N.↗

NASA NTRS1992

The effects of uneven fiber spacing on thermal residual stresses in a unidirectional SCS-6Ti-15-3 laminate

High residual stresses develop in SCS-6/Ti-15-3 composites during cooldown from the fabrication temperature; these residual stresses can effect the mechanical and physical properties of the composite. Discrete fiber-matrix finite element models were used to study the residual stresses due to the temperature change during the fabrication process, including the effects of uneven fiber spacing, the free surface, and increased fiber volume fractions. To accurately model the effects of the free surface, it is only necessary to model one fiber through the thickness. Below the first ply, the analysis predicts stress distributions that are identical to the infinite array predictions. For uneven fiber space less then 0.042 mm in an interior ply, the maximum hoop stress was predicted to occur between fibers within a ply and to increase as the fiber spacing decreased. The maximum hoop stress correlated well with the observed radial cracking between fibers. For the case of touching fibers, the analysis predicted tensile radial stresses at the fiber-matrix debonding during the fabrication cooldown. Identical trends were predicted for uneven fiber spacing in surface plies with slightly greater values of maximum stresses. The analysis predicted matrix yielding to occur upon cooldown when the edge-to-edge fiber spacing was less than or equal to 0.022 mm. The stress distributions predicted for increasing fiber volume fractions were similar to those predicted for decreasing the fiber spacing for two adjacent fibers within a ply.

Bigelow, C. A.↗