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,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Exploration of Failure and Potential Damage Markers in Ti-6Al-4V

To improve prognostics for predicting component life, coupon tests were conducted on Ti-6Al-4V to investigate the development of damage. Various test types were employed (tension, creep, stress relaxation, and fatigue) under a wide range of temperatures and loading rates to engage various amounts of time-dependent and damage behavior. Both stiffness degradation and Poisson’s ratio were monitored to see if they were useful in indicating the initiation and accumulation of damage and subsequently the onset of failure. Stiffness exhibited large decreases with increasing strain when calculated using engineering stress and strain. However, minimal stiffness change was observed when using the true-stress and true- strain values. Detailed microscopy and image analysis were performed to document the actual physical damage that occurred within the various specimens. Two types of damage were observed: shallow surface cracks and internal pores. The severity of damage was greater in the region of localized necking and diminished further away from the neck. The area fraction of pores was higher in the necked region, yet even there it did not exceed more than a few percent and, in most cases, was significantly lower. It was observed that conditions that induced more time-dependent deformation (e.g., slower strain rate testing) resulted in a larger number of pores and larger pore sizes. The pores were always found to occur in the softer β-phase. The severity of observed damage was inconsistent with the large drop in engineering stiffness, but was more aligned with stiffness changes from the true-stress and -strain values. Poisson’s ratio was monitored on every test and showed promise as a damage indicator. It was found to be more sensitive than either the ultimate tensile strength or the inception of tertiary creep. However, calculation of Poisson’s ratio was found to be more involved than expected because of anisotropy and the presence of strain gradients along the specimen at large applied strains. First, the material had a moderate texture and resulted in anisotropy of the strains. Second, the evolution of Poisson’s ratio as a function of axial strain showed strange and unexpected behaviors, which necessitated more extensive analyses of its development. This led to an investigation on the use of small versus large strain approaches and the importance of making consistent measurements; that is, axial and transverse strain measurements over the same volume of material. Additionally, finite element analysis of a tensile sample indicated that non- uniform strain occurred in the sample prior to attaining the maximum strength. This, along with the location at which transverse strain was measured, had a significant influence on the resulting values of Poisson’s ratio. Based on this, lessons learned are offered for future testing.

Ti-6Al-4V↗
NASA NTRS1996

Spray Combustion Modeling with VOF and Finite-Rate Chemistry

A spray atomization and combustion model is developed based on the volume-of-fluid (VOF) transport equation with finite-rate chemistry model. The gas-liquid interface mass, momentum and energy conservation laws are modeled by continuum surface force mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed range flows. The objectives of the present study are: (1) to develop and verify the fractional volume-of-fluid (VOF) cell partitioning approach into a predictor-corrector algorithm to deal with multiphase (gas-liquid) free surface flow problems; (2) to implement the developed unified algorithm in a general purpose computational fluid dynamics (CFD) code, Finite Difference Navier-Stokes (FDNS), with droplet dynamics and finite-rate chemistry models; and (3) to demonstrate the effectiveness of the present approach by simulating benchmark problems of jet breakup/spray atomization and combustion. Modeling multiphase fluid flows poses a significant challenge because a required boundary must be applied to a transient, irregular surface that is discontinuous, and the flow regimes considered can range from incompressible to highspeed compressible flows. The flow-process modeling is further complicated by surface tension, interfacial heat and mass transfer, spray formation and turbulence, and their interactions. The major contribution of the present method is to combine the novel feature of the Volume of Fluid (VOF) method and the Eulerian/Lagrangian method into a unified algorithm for efficient noniterative, time-accurate calculations of multiphase free surface flows valid at all speeds. The proposed method reformulated the VOF equation to strongly couple two distinct phases (liquid and gas), and tracks droplets on a Lagrangian frame when spray model is required, using a unified predictor-corrector technique to account for the non-linear linkages through the convective contributions of VOF. The discontinuities within the sharp interface will be modeled as a volume force to avoid stiffness. Formations of droplets, tracking of droplet dynamics and modeling of the droplet breakup/evaporation, are handled through the same unified predictor-corrector procedure. Thus the new algorithm is non-iterative and is flexible for general geometries with arbitrarily complex topology in free surfaces. The FDNS finite-difference Navier-Stokes code is employed as the baseline of the current development. Benchmark test cases of shear coaxial LOX/H2 liquid jet with atomization/combustion and impinging jet test cases are investigated in the present work. Preliminary data comparisons show good qualitative agreement between data and the present analysis. It is indicative from these results that the present method has great potential to become a general engineering design analysis and diagnostics tool for problems involving spray combustion.

Chen, Yen-Sen↗
DOE OSTI2026

Advancing SiC clad fuel performance model: bridging micro- and macro-scale models and experimental validation

This report presents a workflow for advancing fuel performance modeling of SiC composite cladding for light-water reactors by linking microscale, experimental data-informed finite element analysis with rod-scale fuel performance codes such as BISON. The workflow uses X-ray computed tomography (XCT) to capture the actual geometry and processing-induced defects of as-fabricated SiC composite tube specimens, particularly porosity and wall-thickness variations, and converts the segmented XCT volumes into image-based finite element meshes for high fidelity structural analysis.

Koyanagi, Takaaki [Oak Ridge National Laboratory (↗
NASA NTRS1977

DESAP 1: A structural design program with stress and displacement constraints. Volume 1: Theoretical and user's manual

A finite element program is presented for computer-automated, minimum weight design of elastic structures with constraints on stresses (including local instability criteria) and displacements. Volume 1 of the report contains the theoretical and user's manual of the program. Sample problems and the listing of the program are included in Volumes 2 and 3. The element subroutines are organized so as to facilitate additions and changes by the user. As a result, a relatively minor programming effort would be required to make DESAP 1 into a special purpose program to handle the user's specific design requirements and failure criteria.

Kiusalaas, J.↗
NASA NTRS1995

ANSYS duplicate finite-element checker routine

An ANSYS finite-element code routine to check for duplicated elements within the volume of a three-dimensional (3D) finite-element mesh was developed. The routine developed is used for checking floating elements within a mesh, identically duplicated elements, and intersecting elements with a common face. A space shuttle main engine alternate turbopump development high pressure oxidizer turbopump finite-element model check using the developed subroutine is discussed. Finally, recommendations are provided for duplicate element checking of 3D finite-element models.

Ortega, R.↗
DOE OSTI2025

Shear viscosity from perturbative quantum chromodynamics to the hadron resonance gas at finite baryon, strangeness, and electric charge densities

Through model-to-data comparisons from heavy-ion collisions, it has been shown that the quark gluon plasma has an extremely small shear viscosity at vanishing densities. At large baryon densities, significantly less is known about the nature of the shear viscosity from quantum chromodynamics (QCD). Within heavy-ion collisions, there are three conserved charges: baryon number (B), strangeness (S), and electric charge (Q). Here we calculate the shear viscosity in two limits using perturbative QCD (pQCD) and an excluded-volume hadron resonance gas at finite BSQ densities. We then develop a framework that interpolates between these two limits such that shear viscosity is possible to calculate across a wide range of finite BSQ densities. We find that the pQCD and hadron resonance gas calculations have different BSQ density dependencies such that a rather nontrivial shear viscosity appears at finite densities.

Phenomenology↗
NASA NTRS1993

Analytical evaluation of a surface integral expressing the coupling between interior and exterior volumes in a FE-IE approach

Recently an approach which combines the finite element technique and an integral equation to determine the fields scattered by inhomogeneous bodies of complicated shape has been proposed. Basically, a mathematical surface which encloses the scatterers is introduced, thus dividing the space into an interior and an exterior volume, in which the finite element technique and an integral equation for EM scattering, respectively, are applied. The integral equation is set up for the tangential components of the fields at the surface, while the interior volume the unknowns are the total fields. Continuity of the tangential fields at the boundary, as required by Maxwell's equations, is imposed, thus coupling the two methods to obtain a consistent solution. The coupling term is expressed by a surface integral formed by the dot product of a FE basis function and an IE testing function, or viceversa. By choosing the boundary to be a surface of revolution and by making a convenient selection of IE basis (testing) functions, it is possible to evaluate the integrals analytically on surfaces such as curved triangles, curved quadrilaterals and curved pentagons. We will illustrate the salient steps involved in setting up and carrying out these integrals and discuss what class of basis (testing) functions and analytic surfaces of revolution they are applicable to. Analytic calculations offer the advantage of better accuracy than purely numerical ones, and, when combined with them, often shed light on issues of numerical convergence and limiting values. Furthermore, they may reduce computation time and storage requirements.

Zuffada, C.↗
NASA NTRS1997

Binary Oscillatory Crossflow Electrophoresis

Electrophoresis has long been recognized as an effective analytic technique for the separation of proteins and other charged species, however attempts at scaling up to accommodate commercial volumes have met with limited success. In this report we describe a novel electrophoretic separation technique - Binary Oscillatory Crossflow Electrophoresis (BOCE). Numerical simulations indicate that the technique has the potential for preparative scale throughputs with high resolution, while simultaneously avoiding many problems common to conventional electrophoresis. The technique utilizes the interaction of an oscillatory electric field and a transverse oscillatory shear flow to create an active binary filter for the separation of charged protein species. An oscillatory electric field is applied across the narrow gap of a rectangular channel inducing a periodic motion of charged protein species. The amplitude of this motion depends on the dimensionless electrophoretic mobility, alpha = E(sub o)mu/(omega)d, where E(sub o) is the amplitude of the electric field oscillations, mu is the dimensional mobility, omega is the angular frequency of oscillation and d is the channel gap width. An oscillatory shear flow is induced along the length of the channel resulting in the separation of species with different mobilities. We present a model that predicts the oscillatory behavior of charged species and allows estimation of both the magnitude of the induced convective velocity and the effective diffusivity as a function of a in infinitely long channels. Numerical results indicate that in addition to the mobility dependence, the steady state behavior of solute species may be strongly affected by oscillating fluid into and out of the active electric field region at the ends of the cell. The effect is most pronounced using time dependent shear flows of the same frequency (cos((omega)t)) flow mode) as the electric field oscillations. Under such conditions, experiments indicate that solute is drawn into the cell from reservoirs at both ends of the cell leading to a large mass build up. As a consequence, any initially induced mass flux will vanish after short times. This effect was not captured by the infinite channel model and hence numerical and experimental results deviated significantly. The revised model including finite cell lengths and reservoir volumes allowed quantitative predictions of the time history of the concentration profile throughout the system. This latter model accurately describes the fluxes observed for both oscillatory flow modes in experiments using single protein species. Based on the results obtained from research funded under NASA grant NAG-8-1080.S, we conclude that binary separations are not possible using purely oscillatory flow modes because of end effects associated with the cos((omega)t) mode. Our research shows, however, that a combination of cos(2(omega)t) and steady flow should lead to efficient separation free of end effects. This possibility is currently under investigation.

Molloy, Richard F.↗
NASA NTRS2010

Acoustic-Modal Testing of the Ares I Launch Abort System Attitude Control Motor Valve

The Attitude Control Motor (ACM) is being developed for use in the Launch Abort System (LAS) of NASA's Ares I launch vehicle. The ACM consists of a small solid rocket motor and eight actuated pintle valves that directionally allocate.thrust_- 1t.has-been- predicted-that significant unsteady. pressure.fluctuations.will.exist. inside the-valves during operation. The dominant frequencies of these oscillations correspond to the lowest several acoustic natural frequencies of the individual valves. An acoustic finite element model of the fluid volume inside the valve has been critical to the prediction of these frequencies and their associated mode shapes. This work describes an effort to experimentally validate the acoustic finite model of the valve with an acoustic modal test. The modal test involved instrumenting a flight-like valve with six microphones and then exciting the enclosed air with a loudspeaker. The loudspeaker was configured to deliver broadband noise at relatively high sound pressure levels. The aquired microphone signals were post-processed and compared to results generated from the acoustic finite element model. Initial comparisons between the test data and the model results revealed that additional model refinement was necessary. Specifically, the model was updated to implement a complex impedance boundary condition at the entrance to the valve supply tube. This boundary condition models the frequency-dependent impedance that an acoustic wave will encounter as it reaches the end of the supply tube. Upon invoking this boundary condition, significantly improved agreement between the test data and the model was realized.

Davis, R. Benjamin↗
DOE OSTI2025

2-loop scattering on superstring and supermembrane in flat space

We consider the S-matrix of transverse scalar excitations on an infinite D = 10 GS superstring and D = 11 supermembrane in flat target space. We compute the 4-particle scattering amplitude in the 2-loop approximation and demonstrate that, like in the Nambu string case, the D = 10 GS string S-matrix does not contain non-trivial 2d UV divergences (UV pole not accompanied by terms with logarithms of momenta is an artifact of dimensional regularization). This is consistent with underlying integrability of this model which is maintaining by adding appropriate local counterterms. In the supermembrane case there are no 1-loop divergences but we find a genuine 2-loop UV pole. This demonstrates non-finiteness of the world-volume S-matrix of the M2 brane theory.

AdS-CFT Correspondence↗
DOE OSTI2020

Thermodynamics of spin-orbit-coupled bosons in two dimensions from the complex Langevin method

We investigate the thermal properties of interacting spin-orbit-coupled bosons with contact interactions in two spatial dimensions. To that end, we implement the complex Langevin method, motivated by the appearance of a sign problem, on a square lattice with periodic boundary conditions. We calculate the density equation of state nonperturbatively in a range of spin-orbit couplings and chemical potentials. In this work, our results show that mean-field solutions tend to underestimate the average density, especially for stronger values of the spin-orbit coupling. Additionally, the finite nature of the simulation volume induces the formation of pseudocondensates. These have been observed to be destroyed by the spin-orbit interactions.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

An $\bar{F}$ Meshfree Treatment for Nearly Incompressible Materials

The computational modeling of nearly incompressible materials is a difficult task for many numerical methods, and even after several decades of investigation, it is still an active research area. This report seeks to address the treatment of incompressible materials in meshfree methods using a synergistic combination of two treatments. The first treatment is an $\bar{F}$ method, where the decomposed dilatational and deviatoric parts are calculated over different smoothing domains. The second treatment “activates” additional nodes throughout the domain to increase the flexibility of the model. We implement this synergistic combination in the context of the reproducing kernel particle method (RKPM) and present results for the Cook’s membrane benchmark problem. The results are compared with those using the composite tet10 finite element with a volume-averaged J formulation. We show that the combined treatment is an effective way to deal with nearly incompressible materials in a meshfree framework and compares well with other highly-effective treatments.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1987

Diffusion length measurements in bulk and epitaxially grown 3-5 semiconductors using charge collection microscopy

Diffusion lengths and surface recombination velocities were measured in GaAs diodes and InP finished solar cells. The basic techniques used was charge collection microscopy also known as electron beam induced current (EBIC). The normalized currents and distances from the pn junction were read directly from the calibrated curves obtained while using the line scan mode in an SEM. These values were then equated to integral and infinite series expressions resulting from the solution of the diffusion equation with both extended generation and point generation functions. This expands previous work by examining both thin and thick samples. The surface recombination velocity was either treated as an unknown in a system of two equations, or measured directly using low e(-) beam accelerating voltages. These techniques give accurate results by accounting for the effects of surface recombination and the finite size of the generation volume.

Leon, R. P.↗
NASA NTRS1987

Diffusion length measurement in bulk and epitaxially grown III-V semiconductors using charge collection microscopy

Diffusion lengths and surface recombination velocities were measured in GaAs diodes and InP finished solar cells. The basic technique used was charge collection microscopy, also known as electron beam induced current (EBIC). The normalized currents and distances from the pn junction were read directly from the calibrated curves obtained while using the line-scan mode in an SEM. These values were then equated to integral and infinite series expressions resulting from the solution of the diffusion equation with both extended-generation and point-generation functions. This expands previous work by examining both thin and thick samples. The surface recombination velocity was either treated as an unknown in a system of two equations or measured directly using low e(-) beam accelerating voltages. These techniques give accurate results by accounting for the effects of surface recombination and the finite size of the generation volume.

Leon, R. P.↗
NASA NTRS2002

Investigation of Venus Surface Properties

Strong localized radar echoes have been observed at decimeter wavelengths from the highlands of Venus since the earliest radar maps were obtained over 30 years ago. These echoes are some five to ten times stronger than those from the presumably basaltic rocks seen at lower altitudes elsewhere on Venus. Observations of thermal emission from the visible disk of Venus at wavelengths corresponding to those used in the radar mapping confirm that the regions of high reflectivity also exhibit low emissivity, as expected from considerations of detailed thermodynamic balance. Two possibilities have been put forward to explain this unexpected aspect of the Venus highlands: 1) surface materials of high effective dielectric constant, probably associated with finite electrical conductivity, and 2) volume scattering associated with multiple scattering from a layer of very-low-loss material containing voids and extending down a few hundred wavelengths beneath the surface. Analogs to these two mechanisms are found elsewhere in the solar system, and each is capable of explaining the basic observations. as of the early 1990's. In 1993, however, it became possible to carry out a bistatic observation of the anomalous highland regions using the Magellan spacecraft, then in orbit about Venus. In this experiment the on-board telemetry transmitter was aimed at the planet's surface with its linear S-band polarization vector oriented at 45 deg to the spacecraft-Venus-Earth scattering plane. The pointing of the transmitting antenna was adjusted so that the spacecraft-to-illuminated-surface incidence angle equalled the Earth-to-Venus-surface incidence angle. In this way, the experiment emphasized the specular scattering component. A full Stokes-Vector analysis of the reflected signal as received on Earth was carried out as the illuminated region scanned across the highland regions of Venus. From the observed position angle of the echo, it was possible to calculate the Fresnel reflectivity of a number of regions on Venus, including the high-altitude Maxwell Montes, and to show unequivocally that the complex dielectric constant of that anomalous area was 100 +/- 50i, as compared to 4.0 (real) for the typical lowland area. Thus this experiment was able to rule out definitively an hypothesis of volume scattering as the explanation for the anomalous scattering in the Venus highlands.

Ford, Peter G.↗
NASA NTRS1987

Component mode synthesis and large deflection vibration of complex structures. Volume 2: Single-mode large deflection vibrations of beams and plates using finite element method

A finite element method is presented for the large amplitude vibrations of complex structures that can be modelled with beam and rectangular plate elements subjected to harmonic excitation. Both inplane deformation and inertia are considered in the formulation. Derivation of the harmonic force and nonlinear stiffness matrices for a beam and a rectangular plate element are presented. Solution procedures and convergence characteristics of the finite element method are described. Nonlinear response to uniform and concentrated harmonic loadings and improved nonlinear free vibration results are presented for beams and rectangular plates of various boundary conditions.

Mei, Chuh↗
NASA NTRS1987

Component mode synthesis and large deflection vibration of complex structures. Volume 3: Multiple-mode nonlinear free and forced vibrations of beams using finite element method

Multiple-mode nonlinear forced vibration of a beam was analyzed by the finite element method. Inplane (longitudinal) displacement and inertia (IDI) are considered in the formulation. By combining the finite element method and nonlinear theory, more realistic models of structural response are obtained more easily and faster.

Mei, Chuh↗
NASA NTRS1989

Space Shuttle Main Engine structural analysis and data reduction/evaluation. Volume 1: Aft Skirt analysis

Using the ANSYS finite element program, a global model of the aft skirt and a detailed nonlinear model of the failure region was made. The analysis confirmed the area of failure in both STA-2B and STA-3 tests as the forging heat affected zone (HAZ) at the aft ring centerline. The highest hoop strain in the HAZ occurs in this area. However, the analysis does not predict failure as defined by ultimate elongation of the material equal to 3.5 percent total strain. The analysis correlates well with the strain gage data from both the Wyle influence test of the original design aft sjirt and the STA-3 test of the redesigned aft skirt. it is suggested that the sensitivity of the failure area material strength and stress/strain state to material properties and therefore to small manufacturing or processing variables is the most likely cause of failure below the expected material ultimate properties.

Berry, David M.↗