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,369 records · Page 76

Attenuation behavior of solid dense random media at microwave frequencies

To better understand scattering from nontenuous dense random media such as sea ice and snow, attenuation measurements have been performed on two different types of random media with ka values ranging from 0.5 to 0.7, and 1.5 to 2.1. Experimental results are presented for wave propagation in plane slabs of finite thickness composed of a random distribution of identical, finite scatterers, and a random distribution of scatterers with narrow size distribution. The observed behavior is described in terms of attenuation versus volume fraction, and the behavior of attenuation versus frequency. Results presented are compared to the behavior reported by some earlier experiments where the medium properties are different.

Nadimi, S. A.↗

Structural Anomalies Detected in Ceramic Matrix Composites Using Combined Nondestructive Evaluation and Finite Element Analysis (NDE and FEA)

Most reverse engineering approaches involve imaging or digitizing an object and then creating a computerized reconstruction that can be integrated, in three dimensions, into a particular design environment. The rapid prototyping technique builds high-quality physical prototypes directly from computer-aided design files. This fundamental technique for interpreting and interacting with large data sets is being used here via Velocity2 (an integrated image-processing software, ref. 1) using computed tomography (CT) data to produce a prototype three-dimensional test specimen model for analyses. A study at the NASA Glenn Research Center proposes to use these capabilities to conduct a combined nondestructive evaluation (NDE) and finite element analysis (FEA) to screen pretest and posttest structural anomalies in structural components. A tensile specimen made of silicon nitrite (Si3N4) ceramic matrix composite was considered to evaluate structural durability and deformity. Ceramic matrix composites are being sought as candidate materials to replace nickel-base superalloys for turbine engine applications. They have the unique characteristics of being able to withstand higher operating temperatures and harsh combustion environments. In addition, their low densities relative to metals help reduce component mass (ref. 2). Detailed three-dimensional volume rendering of the tensile test specimen was successfully carried out with Velocity2 (ref. 1) using two-dimensional images that were generated via computed tomography. Subsequent, three-dimensional finite element analyses were performed, and the results obtained were compared with those predicted by NDE-based calculations and experimental tests. It was shown that Velocity2 software can be used to render a three-dimensional object from a series of CT scan images with a minimum level of complexity. The analytical results (ref. 3) show that the high-stress regions correlated well with the damage sites identified by the CT scans and the experimental data. Furthermore, modeling of the voids collected via NDE offered an analytical advantage that resulted in more accurate assessments of the material s structural strength. The top figure shows a CT scan image of the specimen test section illustrating various hidden structural entities in the material and an optical image of the test specimen considered in this study. The bottom figure represents the stress response predicted from the finite element analyses (ref .3 ) for a selected CT slice where it clearly illustrates the correspondence of the high stress risers due to voids in the material with those predicted by the NDE. This study is continuing, and efforts are concentrated on improving the modeling capabilities to imitate the structural anomalies as detected.

Abdul-Aziz, Ali↗

Thermal Pollution Mathematical Model: Verification of One-Dimensional Numerical Model at Lake Keowee - Volume 2

A one dimensional model for studying the thermal dynamics of cooling lakes was developed and verified. The model is essentially a set of partial differential equations which are solved by finite difference methods. The model includes the effects of variation of area with depth, surface heating due to solar radiation absorbed at the upper layer, and internal heating due to the transmission of solar radiation to the sub-surface layers. The exchange of mechanical energy between the lake and the atmosphere is included through the coupling of thermal diffusivity and wind speed. The effects of discharge and intake by power plants are also included. The numerical model was calibrated by applying it to Cayuga Lake. The model was then verified through a long term simulation using Lake Keowee data base. The comparison between measured and predicted vertical temperature profiles for the nine years is good. The physical limnology of Lake Keowee is presented through a set of graphical representations of the measured data base.

Lee, S. S.↗

Thermal-capillary model for Czochralski growth of semiconductor materials

The success of efficiently calculating the temperature field, crystal radius, melt mensicus, and melt/solid interface in the Czochralski crystal growth system by full finite-element solution of the government thermal-capillary model is demonstrated. The model predicts realistic response to changes in pull rate, melt volume, and the thermal field. The experimentally observed phenomena of interface flipping, bumping, and the difficulty maintaining steady-state growth as the melt depth decreases are explained by model results. These calculations will form the basis for the first quantitative picture of Cz crystal growth. The accurate depiction of the melt meniscus is important in calculating the crystal radius and solidification interface. The sensitivity of the results to the equilibrium growth angle place doubt on less sophisticated attempts to model the process without inclusion of the meniscus. Quantitative comparison with experiments should be possible once more representation of the radiation and view factors in the thermal system and the crucible are included. Extensions of the model in these directions are underway.

Derby, J. J.↗

Micromechanics of crenulated fibers in carbon/carbon composites

The influence of crenulated noncircular fibers on the micromechanical stress states due to a transverse strain and to a temperature change in carbon/carbon composites is examined using the finite element method. Stresses at the interface of both fully bonded and fully disbonded fibers having two crenulation amplitudes and with two fiber volume fractions are presented. In each case, these interface stresses are compared to stresses at the interface of circular fibers which have the same degree of disbond and fiber volume fraction and are under the same loading conditions. For the disbonded cases, deformed meshes showing locations of fiber/matrix contact are also included. In addition to the interface stress states, selected composite properties are also computed and compared in each case examined. Interest in studying noncircular fibers stems from a desire to increase the transverse properties of carbon/carbon by introducing a mechanical interlocking between the fiber and the matrix. Results presented here indicate that this interlocking does in fact occur. Evidence from the interface stress data suggests, however, that any possible advantage of this interlocking may be outweighed by the disadvantage of stress concentrations which arise at the interface due to the crenulated geometry of the fibers.

Carapella, E. E.↗

Development of Long-Lifetime Pulsed Gas Valves for Pulsed Electric Thrusters

It is advantageous for gas-fed pulsed electric thrusters to employ pulsed valves so propellant is only flowing to the device during operation. The propellant utilization of the thruster will be maximized when all the gas injected into the thruster is acted upon by the fields produced by the electrical pulse. Gas that is injected too early will diffuse away from the thruster before the electrical pulse can act to accelerate the propellant. Gas that is injected too late will miss being accelerated by the already-completed electrical pulse. As a consequence, the valve must open quickly and close equally quickly, only remaining open for a short duration. In addition, the valve must have only a small amount of volume between the sealing body and the thruster so the front and back ends of the pulse are as coincident as possible with the valve cycling, with very little latent propellant remaining in the feed lines after the valve is closed. For a real mission of interest, a pulsed thruster can be expected to pulse at least 10(exp 10) - 10(exp 11) times, setting the range for the number of times a valve must open and close. The valves described in this paper have been fabricated and tested for operation in an inductive pulsed plasma thruster (IPPT) for in-space propulsion. In general, an IPPT is an electrodeless space propulsion device where a capacitor is charged to an initial voltage and then discharged, producing a high-current pulse through a coil. The field produced by this pulse ionizes propellant, inductively driving current in a plasma located near the face of the coil. Once the plasma is formed, it can be accelerated and expelled at a high exhaust velocity by the electromagnetic Lorentz body force arising from the interaction of the induced plasma current and the magnetic field produced by the current in the coil. The valve characteristics needed for the IPPT application require a fast-acting valve capable of a minimum of 10(exp 10) valve actuation cycles. Since even 10(exp 9) cycles is well above anything demonstrated, this lower value was selected as the design point for the present work. The valve seal must remain leak-tight throughout operation, and the body must maintain a low internal leakage at relatively high operating temperatures. The full set of design requirements used for this program are summarized in Table 1. In this work, we describe two pulsed gas valves that have been fabricated to have long lifetime and demonstrate the characteristics listed above. The first is a miniaturized, conventional electromagnet-based valve while the second is a piezoelectric-based valve design. The conventional valve, shown in Fig. 1, is opened by use of a solenoid electromagnetic actuator. When current is applied to the solenoid coil, magnetic forces pull the plunger away from the valve seat, allowing fluid to flow through the valve. Removal of electrical current permits the spring and fluid pressure to seat the plunger, halting the flow of fluid. The valve body is fabricated from 304L corrosion resistant steel (CRES) and while the parts that form the magnetic circuit are fabricated from 430 CRES. This material does not have optimum magnetic properties, but its corrosion resistance permits incorporation into a design without requiring an additional plating process. A viton O-ring compound (Parker V0884-75), selected for its mechanical strength at elevated temperatures, was used for the valve seat seal. The design was based solely on the use of analytical sizing calculations, as opposed to a more rigorous finite element analysis. While this valve is small and relatively lightweight, it does not represent a design that is optimized for mass and/or a given volume envelope. The piezoelectric valve is a "puller" valve design. Applying a voltage to the piezo crystal causes it to elongate and pull a pintle off the seat, opening the valve. The valve seal consists of the pintle with an external, spherically-formed tip fabricated from Inconel 625 and the seat fabricated from Inconel 718 and possessing a spherical surface that permits the pintle to align with the valve body centerline after each cycle, maintaining the circumferential metal-on-metal seal. The contact stress on the valve pintle is expected in the worst case to be 96.5 MPa (14 ksi), which is well below the work hardening value of 413.7 MPa (60 ksi) for Inconel 615. Results presented in this work will compare the performance of each valve, both to each other and to the valve design requirements given in Table 1. In general, both valves meet all the performance envelope requirements, but the piezoelectric valve consumes considerably less power and has a faster response time over the conventional solenoid-driven valves.

Burkhardt, Wendel M.↗

Models of Saturn's rings which satisfy the optical observations

A theoretical model of Saturn's rings is investigated which includes the shadowing effect and realistic anisotropic phase functions for the ring particles. The effects of multiple scattering and the finite size of the sun, including the penumbra, are rigorously included. The permissible range of the relevant parameters, including optical thickness, single scattering albedo, volume density, and phase function are investigated by comparing the theoretical results to observations of the ring brightness vs phase angle, wavelength, and elevation of the sun and earth. Anisotropic scattering by the ring particles is necessary in order to match the observations. The color dependence of the opposition effect is interpreted in terms of the albedo spectrum of the ring particles.

Kawata, Y.↗

A scattering model for forested area

A forested area is modeled as a volume of randomly oriented and distributed disc-shaped, or needle-shaped leaves shading a distribution of branches modeled as randomly oriented finite-length, dielectric cylinders above an irregular soil surface. Since the radii of branches have a wide range of sizes, the model only requires the length of a branch to be large compared with its radius which may be any size relative to the incident wavelength. In addition, the model also assumes the thickness of a disc-shaped leaf or the radius of a needle-shaped leaf is much smaller than the electromagnetic wavelength. The scattering phase matrices for disc, needle, and cylinder are developed in terms of the scattering amplitudes of the corresponding fields which are computed by the forward scattering theorem. These quantities along with the Kirchoff scattering model for a randomly rough surface are used in the standard radiative transfer formulation to compute the backscattering coefficient. Numerical illustrations for the backscattering coefficient are given as a function of the shading factor, incidence angle, leaf orientation distribution, branch orientation distribution, and the number density of leaves. Also illustrated are the properties of the extinction coefficient as a function of leaf and branch orientation distributions. Comparisons are made with measured backscattering coefficients from forested areas reported in the literature.

Karam, M. A.↗

User's Manual for FEMOM3DS

FEMOM3DS is a computer code written in FORTRAN 77 to compute electromagnetic(EM) scattering characteristics of a three dimensional object with complex materials using combined Finite Element Method (FEM)/Method of Moments (MoM) technique. This code uses the tetrahedral elements, with vector edge basis functions for FEM in the volume of the cavity and the triangular elements with the basis functions similar to that described for MoM at the outer boundary. By virtue of FEM, this code can handle any arbitrarily shaped three-dimensional cavities filled with inhomogeneous lossy materials. The User's Manual is written to make the user acquainted with the operation of the code. The user is assumed to be familiar with the FORTRAN 77 language and the operating environment of the computers on which the code is intended to run.

Reddy, C.J.↗

Viscous hydrodynamic instability theory of the peak and minimum pool boiling heat fluxes

Liquid viscosity was included in the Bellman-Pennington theory of the Taylor wave in a liquid vapor interface. Predictions of the most susceptible wavelength, and of the wave frequency, were made as a function of a liquid viscosity parameter and the Bond number. The stability of a gas jet in a viscous liquid was studied and the result is used to predict the peak heat flux on large horizontal heaters. Experimental measurements of the dominant Taylor wave and its growth rate were made during the film boiling of cyclohexanol on cylindrical heaters. The results bear out the predictions quite well. The thickness of the vapor blanket surrounding a cylindrical heater was measured and a correlation suggested. The effect of large fluxes of vapor volume on the dominant wavelength was also noted. Theoretical results of the peak heat flux are compared with the experimental data, and the effect of finite geometry of flat plate heaters on the peak heat flux is also discussed.

Dhir, V. K.↗

User's manual for CBS3DS, version 1.0

CBS3DS is a computer code written in FORTRAN 77 to compute the backscattering radar cross section of cavity backed apertures in infinite ground plane and slots in thick infinite ground plane. CBS3DS implements the hybrid Finite Element Method (FEM) and Method of Moments (MoM) techniques. This code uses the tetrahedral elements, with vector edge basis functions for FEM in the volume of the cavity/slot and the triangular elements with the basis functions for MoM at the apertures. By virtue of FEM, this code can handle any arbitrarily shaped three-dimensional cavities filled with inhomogeneous lossy materials; due to MoM, the apertures can be of any arbitrary shape. The User's Manual is written to make the user acquainted with the operation of the code. The user is assumed to be familiar with the FORTRAN 77 language and the operating environment of the computer the code is intended to run.

Reddy, C. J.↗

Managing the On-Board Data Storage, Acknowledgment and Retransmission System for Spitzer

The Spitzer Space Telescope has a two-phase downlink system. Data are transmitted during one telecom session. Then commands are sent during the next session to delete those data that were received and to retransmit those data that were missed. We must build sequences that are as efficient as possible to make the best use of our finite supply of liquid helium, One way to improve efficiency is to use only the minimum time needed during telecom sessions to transmit the predicted volume of data. But, we must also not fill the onboard storage and must allow enough time margin to retransmit missed data. We describe tools and procedures that allow us to build observatory sequences that are single-fault tolerant in this regard and that allow us to recover quickly and safely from anomalies that affect the receipt or acknowledgment of data.

Spitzer↗

NASTRAN hydroelastic modal studies. Volume 1: Introduction, theory, and results

The implementation of a general three-dimensional hydroelastic capability in the NASTRAN computer program is described. A method for analyzing the combined mode shapes of arbitrary fluid and structure finite element models was developed. The fluid was modeled with three-dimensional solid elements with options for tetrahedron, wedge, and hexahedron shapes. The elements were then connected to fluid grid points which defined the pressure in the fluid at the specified location. Matrices were formulated to provide efficient solutions for large-order problems. The structure matrices were processed separately and could be reduced using matrix condensation procedures or through a modal formulation using the normal modes of the empty structure as solution coordinates. The fluid/structure interface and the free surface were defined by the user with special NASTRAN boundary elements. A special purpose mesh generator program was provided to generate the actual NASTRAN data cards for the fluid, the structure, and the boundary elements for typical tank-type models.

Source record↗

Research and development program for the development of advanced time-temperature dependent constitutive relationships. Volume 2: Programming manual

The results of a 10-month research and development program for nonlinear structural modeling with advanced time-temperature constitutive relationships are presented. The implementation of the theory in the MARC nonlinear finite element code is discussed, and instructions for the computational application of the theory are provided.

Cassenti, B. N.↗