Search NASA⌕ Search

SEARCH · Search NASA

Results for “Volume integral equation”

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 307 records · Page 17

Computational simulation of intermingled-fiber hybrid composite behavior

Three-dimensional finite-element analysis and a micromechanics based computer code ICAN (Integrated Composite Analyzer) are used to predict the composite properties and microstresses of a unidirectional graphite/epoxy primary composite with varying percentages of S-glass fibers used as hydridizing fibers at a total fiber volume of 0.54. The three-dimensional finite-element model used in the analyses consists of a group of nine fibers, all unidirectional, in a three-by-three unit cell array. There is generally good agreement between the composite properties and microstresses obtained from both methods. The results indicate that the finite-element methods and the micromechanics equations embedded in the ICAN computer code can be used to obtain the properties of intermingled fiber hybrid composites needed for the analysis/design of hybrid composite structures. However, the finite-element model should be big enough to be able to simulate the conditions assumed in the micromechanics equations.

Mital, Subodh K.↗

High accuracy computation of fluid-structure interaction in transonic cascades

A coupling strategy for simulating fluid-structure interaction phenomena is formulated and applied to the prediction of flutter in transonic cascades. The flow is governed by the Euler equations and discretized using a finite volume flux-splitting scheme. The structure is modeled using an isoparametric finite element formulation. The coupling strategy successfully reconciles these two formulations at the fluid-structure interface by enforcing both kinematic and kinetic boundary conditions. In particular, the conservation laws applicable to the combined fluid-structure system are preserved across the interface. Since the primary mechanism driving aeroelastic phenomena involves energy exchange occurring at the interface, this highly accurate coupling mechanism is believed to improve the predictive capability of the scheme. The coupled equations are advanced simultaneously in time using an implicit time integration method. Results obtained using the coupling method are presented for cascade geometries operating in transonic flow.

Boschitsch, Alexander↗

Monte-Carlo Analysis of Minimum Thermocouple Depths using Icarus

Icarus is a three-dimensional, unstructured, finite-volume material response solver developed at NASA Ames Research Center and has been verified against other NASA material response tools like FIAT, which have a long history of successfully designing thermal protection system (TPS). Icarus solves a set of conservation equations for mass and energy and uses Darcy’s Law in place of momentum conservation. An ecosystem of material response tools has been built around a general-purposed Icarus library that in addition to the typical material response analysis also supports TPS sizing (1-D and multi-dimensional), uncertainty quantification, and has been successfully integrated into a multi-physics architecture built around US3D. In this paper, a brief overview of Icarus and its capabilities will be presented using an illustrative Monte Carlo analysis of the one-dimensional, in-depth material response of a representative Dragonfly trajectory.

Material Response↗

Blade loss transient dynamics analysis, volume 2. Task 2: Theoretical and analytical development. Task 3: Experimental verification

The component element method was used to develop a transient dynamic analysis computer program which is essentially based on modal synthesis combined with a central, finite difference, numerical integration scheme. The methodology leads to a modular or building-block technique that is amenable to computer programming. To verify the analytical method, turbine engine transient response analysis (TETRA), was applied to two blade-out test vehicles that had been previously instrumented and tested. Comparison of the time dependent test data with those predicted by TETRA led to recommendations for refinement or extension of the analytical method to improve its accuracy and overcome its shortcomings. The development of working equations, their discretization, numerical solution scheme, the modular concept of engine modelling, the program logical structure and some illustrated results are discussed. The blade-loss test vehicles (rig full engine), the type of measured data, and the engine structural model are described.

Gallardo, V. C.↗

A Comparison of Experimental and Theoretical Results for Labyrinth Gas Seals

The basic equations are derived for a two control volume model for compressible flow in a labyrinth seal. The flow is assumed to be completely turbulent and isoenergetic. The wall friction factors are determined using the Blasius formula. Jet flow theory is used for the calculation of the recirculation velocity in the cavity. Linearized zeroth and first order perturbation equations are developed for small motion about a centered position by an expansion in the eccentricity ratio. The zeroth order pressure distribution is found by satisfying the leakage equation. The circumferential velocity distribution is determined by satisfying the momentum equations. The first order equations are solved by a separation of variable solution. Integration of the resultant pressure distribution along and around the seal defines the reaction force developed by the seal and the corresponding dynamic coefficients. The results of this analysis are compared to experimental test results.

Scharrer, Joseph Kirk↗

3D Material Response Analysis of PICA Pyrolysis Experiments

The PICA decomposition experiments of Bessire and Minton are investigated using 3D material response analysis. The steady thermoelectric equations have been added to the CHAR code to enable analysis of the Joule-heated experiments and the DAKOTA optimization code is used to define the voltage boundary condition that yields the experimentally observed temperature response. This analysis has identified a potential spatial non-uniformity in the PICA sample temperature driven by the cooled copper electrodes and thermal radiation from the surface of the test article (Figure 1). The non-uniformity leads to a variable heating rate throughout the sample volume that has an effect on the quantitative results of the experiment. Averaging the results of integrating a kinetic reaction mechanism with the heating rates seen across the sample volume yield a shift of peak species production to lower temperatures that is more significant for higher heating rates (Figure 2) when compared to integrating the same mechanism at the reported heating rate. The analysis supporting these conclusions will be presented along with a proposed analysis procedure that permits quantitative use of the existing data. Time permitting, a status on the in-development kinetic decomposition mechanism based on this data will be presented as well.

Oliver, A. Brandon↗

Monte-Carlo Analysis of Minimal Thermocouple Depths using Icarus

Icarus is a three-dimensional, unstructured, finite-volume material response solver developed at NASA Ames Research Center \cite{Schulz_2017} and has been verified against other NASA material response tools like FIAT, which have a long history of successfully designing thermal protection system (TPS). Icarus solves a set of conservation equations for mass and energy and uses Darcy’s Law in place of momentum conservation. An ecosystem of material response tools has been built around a general-purposed Icarus library that in addition to the typical material response analysis also supports TPS sizing (1-D and multi-dimensional), uncertainty quantification, and has been successfully integrated into a multi-physics architecture built around US3D \cite{Schroeder_2021}. In this paper, a brief overview of Icarus and its capabilities will be presented using an illustrative Monte Carlo analysis of the one-dimensional, in-depth material response of a representative Dragonfly trajectory.

Material Response↗

Finite-volume application of high order ENO schemes to multi-dimensional boundary-value problems

The finite volume approach in developing multi-dimensional, high-order accurate essentially non-oscillatory (ENO) schemes is considered. In particular, a two dimensional extension is proposed for the Euler equation of gas dynamics. This requires a spatial reconstruction operator that attains formal high order of accuracy in two dimensions by taking account of cross gradients. Given a set of cell averages in two spatial variables, polynomial interpolation of a two dimensional primitive function is employed in order to extract high-order pointwise values on cell interfaces. These points are appropriately chosen so that correspondingly high-order flux integrals are obtained through each interface by quadrature, at each point having calculated a flux contribution in an upwind fashion. The solution-in-the-small of Riemann's initial value problem (IVP) that is required for this pointwise flux computation is achieved using Roe's approximate Riemann solver. Issues to be considered in this two dimensional extension include the implementation of boundary conditions and application to general curvilinear coordinates. Results of numerical experiments are presented for qualitative and quantitative examination. These results contain the first successful application of ENO schemes to boundary value problems with solid walls.

Casper, Jay↗

Liquid state theory of the structure of model polymerized ionic liquids

We employ polymer integral equation theory to study a simplified model of semiflexible polymerized ionic liquids (PolyILs) that interact via hard core repulsions and short range screened Coulomb interactions. The multi-scale structure in real and Fourier space of PolyILs (ions chosen to mimic Li, Na, K, Br, PF 6 , and TFSI) are determined as a function of melt density, Coulomb interaction strength, and ion size. Comparisons with a homopolymer melt, a neutral polymer–solvent-like athermal mixture, and an atomic ionic liquid are carried out to elucidate the distinct manner that ions mediate changes of polymer packing, the role of excluded volume effects, and the influence of chain connectivity, respectively. The effect of Coulomb strength depends in a rich manner on ion size and density, reflecting the interplay of steric packing, ion adsorption, and charge layering. Ion-mediated bridging of monomers is found, which intensifies for larger ions. Intermediate range charge layering correlations are characterized by a many-body screening length that grows with PolyIL density, cooling, and Coulomb strength, in disagreement with Debye–Hückel theory, but in accord with experiments. Qualitative differences in the collective structure, including an ion-size-dependent bifurcation of the polymer structure factor peak and pair correlation function, are predicted. The monomer cage order parameter increases significantly, but its collective ion counterpart decreases, as ions become smaller. Such behaviors allow one to categorize PolyILs into two broad classes of small and large ions. Furthermore, dynamical implications of the predicted structural results are qualitatively discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

GC/MS Method Development for Separating Lunar Volatile Ice Simulant Headspace Gases

Various investigators propose the lunar surface contains widely distributed volatiles, especially water- like species, i.e. OH and H2O. Surface volatiles are theorized to exist as a hydrated regolith layer, concentrated in extremely cold polar permanently shadowed regions (PSR), and/or solar wind implantation reservoirs in lunar glasses. The proposed sources of lunar surface volatiles range from cometary impacts, solar wind, or a supply present during moon formation. Future Artemis missions aim to collect and return the samples containing volatiles collected near lunar polar craters or PSRs. We, as advanced curation scientists, are responsible for developing techniques and methodologies for preserving returned sample integrity as much as possible. Pristine volatile-bearing samples are invaluable to the scientific community seeking to unravel the history of the solar system. Realistically, a sample will experience alteration during collection, transportation back to earth, and storage. The Planetary Exploration and Astromaterials Research Lab (PEARL) seeks to understand temperature and pressure effects on high-fidelity volatile-containing regolith simulants, the foundation for the future of cold curation. This abstract outlines the separation, identification, and quantification of headspace gases over volatile ice feed stock material using gas chromatography/mass spectrometry (GC/MS). Preliminary objectives concentrated on sample handling, reproducibility, and understanding the elution characteristics for each analyte. Initial GC/MS method development experiments utilized diluted static headspace sample preparation. Diluted samples were used because sampling headspace gases directly from a vial containing liquid analyte resulted in overloading of the column and detector. Overloading is evident based on chromatogram peak shapes and instrument contamination, or carry over, between experiments. A mixture of three alcohols were used for a majority of the sample handling and reproducibility studies. Reproducibility was tested via multiple users, calibration curves, and check standards. Stock solutions of condensed lunar volatile analytes included methanol, ammonia in methanol, hydrogen sulfide in water, and an equal volume mixture of methanol, ethanol, and isopropanol. Current samples use room air as the headspace sample matrix, however future experiments will incorporate an inert purge gas, such as argon or nitrogen. Three mL of each analyte solution were capped in separate 20 mL crimp top GC vials. Dilutions were carried out by removing an aliquot of headspace gases with a calibrated 1 mL gastight syringe and immediately transferring to a 20 mL capped crimp top vial. The GC/MS is a Thermo Fisher Trace 1310/ISQ 7000 with a TriPlus RSH autosampler and split/splitless injector module. The experiments outlined in this abstract use the following hardware: a 2.5 mL gastight headspace syringe tool, 1 mm ID x 78.5 mm length ultra-inert straight injection liner, and a TG-BondQ 30 m × 0.32 mm × 10 μm column. Various parameters, such as hardware selection and the temperature, pressure, and split ratio set points, continue to evolve as the overall experiment is refined. Diluted headspace chromatograms were collected for the individual stock solutions. Retention times, peak shapes, and mass spectra were evaluated and added to the data processing method for each molecule of interest. Figure 1 shows the total ion chromatograms for the three major lunar volatile simulant stock solutions: methanol, 7 N ammonia in methanol, and 0.4% hydrogen sulfide in water. Tailing peak shapes for ammonia (2.98 min rt) and water (4.06 min rt) indicate the molecules are not properly eluting from the selected column with the current separation method. Additionally, hydrogen sulfide and ammonia have overlapping peak windows, which could impact quantification. Ongoing experiments aim to address the peak shape and overlapping via the separation method and hardware selection. Sample preparation reproducibility experiments used stock solution containing equal volumes of a non- interactive mixture of methanol, ethanol, and isopropanol. Mass spectrum ion traces were used to identify and quantify all three alcohols. Peaks were automatically detected, identified, and integrated through the mass spectra detection and processing parameters. Calibration response curves and check standards were used to evaluate the validity of the sample preparation procedure. Figure 2 shows the methanol chromatogram peak area versus total headspace dilution volume transferred from the alcohol mixture vial. The calibration response curves and check standards validate sample preparation procedure. Continuing data analysis efforts are working towards correlating the peak area and instrument response factor to the headspace analyte concentration and condensed phase composition. Static headspace gas chromatography theory relies on Dalton’s law, Raoult’s law, Henry’s Law, and the Kolb and Ettre equation to associate peak area to the analyte composition in a non-ideal solution. Equation 1 is a simplified expression derived from the aforementioned theories. Future experiments involve liquid injections of the individual stock solutions, liquid and headspace analysis of various stock solution combinations, and the addition of regolith simulants to the mixtures. Temperature is another variable expected to affect reaction rates and will be explored.

Cecilia L. Amick↗

Analysis and conceptual design of a lunar radiator parabolic shade

On the moon, the available heat sink temperature for a vertical unshaded radiator at the equator is 322 K. A method of reducing this heat sink temperature using a parabolic trough shading device was investigated. A steady state heat balance was performed to predict the available heat sink temperature. The effect of optical surface properties on system performance was investigated. Various geometric configurations were also evaluated. A flexible shade conceptual design is presented which greatly reduces the weight and stowed volume of the system. The concept makes use of the natural catenary shape assumed by a flexible material when supported at two points. The catenary shape is very near parabolic. The lunar radiator parabolic shade design presented integrates the energy collection and rejection of a solar dynamic power cycle with the moderate temperature waste heat rejection of a lunar habitat.

Ewert, Michael K.↗

Unsteady hybrid vortex technique for transonic vortex flows and flutter applications

The well developed surface panel method is extended by adding a volume integral term allowing calculation of the full effect of compressibility. The full effect of compressibility is calculated by using mixed type finite difference scheme consistent with the mixed nature of transonic flow. The solution is obtained through successive iteration cycles for subcritical flows and for critical flows the solution is obtained through satisfying the Rankine-Hugoniot relations across the captured shock in addition to the captured shock. The method is applied to airfoils in subcritical and critical flows and the results are in good agreement with the experimental data and finite-difference solutions of the full-potential and Euler equations.

Kandil, Osama A.↗

CFD Analysis of Modular Thrusters Performance

The effective performance of modular thrusters in an aerospike configuration is difficult to determine. Standard analytical tools are applicable to conventional nozzle shapes, but are limited when applied to an aerospike nozzle (An aerospike nozzle is an altitude compensating external nozzle). Three baseline nozzle shapes are derived using standard analytical procedures. The baseline nozzle sizes are restricted to fill a volume envelope. The three shapes are an axi-symmetric round nozzle, a two dimensional planar square exit nozzle, and a super elliptic round to nearly square nozzle. The integrated (thruster/aerospike) performance of the three nozzles is determined through the use of three dimensional viscous computational fluid dynamic (CFD) calculation where complex features of the flow field can be accurately captured. The resulting installed performance is then used to evaluate the efficiency of these nozzle shapes for aerospike applications. The determination of effective performance of a thruster nozzle integrated into an aerospike nozzle require the solution of the three dimensional turbulent Navier-Stokes equations. The model used in this study consisted of two zones; one of the upstream thruster cowl surface so freestream conditions can be accurately predicted, and two, the aerospike surface beginning with the thruster outflow and extending to the end of the aerospike surface. The numerical grid consisted of over 120,000 nodes and used symmetry on the thruster centerline and edge. A two species non-reacting chemistry model was used to capture the variation of fluid properties between the hot plume base and freestream air. From the results of the three baseline nozzle aerospike calculations, the effictive performance of the nozzle was determined. The flow fields of these calculations do show some variation between the cases. Recirculation zones on the cowl surface is predicted for the two dimensional planar nozzle and a smaller one for the super elliptic nozzle. The recirculation is caused by the strong pressure gradient between the plume and freestream flows. The axi-symmetric nozzle results indicate recirculation zones on the thruster face. These recirculation zones smooth the pressure gradient between the plume and freestream flow limiting the formation of recirculation on the cowl surface. Thruster to thruster interaction is evident for the axi-symmetric and supper elliptic calculation while the two dimensional planar nozzle did not have any lateral expansion in the nozzle, so thruster to thruster interaction is limited. The integrated performance results, at the altitude choosen, show very little variation between the three thruster shapes. This result allows for nozzle shape determination based on additional considerations (thermal, structural, weight) besides performance.

Ungewitter, Ronald J.↗

Numerical computation of shock wave-turbulent boundary layer interaction in transonic flow over an axisymmetric curved hill

A control-volume based finite difference computation of a turbulent transonic flow over an axisymmetric curved hill is presented. The numerical method is based on the SIMPLE algorithm, and hence the conservation of mass equation is replaced by a pressure correction equation for compressible flows. The turbulence is described by a k-epsilon turbulence model supplemented by a near-wall turbulence model. In the method, the dissipation rate in the region very close to the wall is obtained from an algebraic equation and that for the rest of the flow domain is obtained by solving a partial differential equation for the dissipation rate. The other flow equations are integrated up to the wall. It is shown that the present turbulence model yields the correct location of the compression shock. The other computational results are also in good agreement with experimental data.

Kim, S.-W.↗

Prediction and verification of creep behavior in metallic materials and components for the space shuttle thermal protection system. Volume 2: Phase 2 subsize panel cyclic creep predictions

A method for predicting permanent cyclic creep deflections in stiffened panel structures was developed. The resulting computer program may be applied to either the time-hardening or strain-hardening theories of creep accumulation. Iterative techniques were used to determine structural rotations, creep strains, and stresses as a function of time. Deflections were determined by numerical integration of structural rotations along the panel length. The analytical approach was developed for analyzing thin-gage entry vehicle metallic-thermal-protection system panels subjected to cyclic bending loads at high temperatures, but may be applied to any panel subjected to bending loads. Predicted panel creep deflections were compared with results from cyclic tests of subsize corrugation and rib-stiffened panels. Empirical equations were developed for each material based on correlation with tensile cyclic creep data and both the subsize panels and tensile specimens were fabricated from the same sheet material. For Vol. 1, see N75-21431.

Cramer, B. A.↗

Derivation of the fundamental equation of sound generated by moving aerodynamic surfaces

Ffowcs Williams and Hawkings (1969) based their derivation of the fundamental equation of the sound generated by arbitrarily moving aerodynamic surfaces on the study of mass and momentum balance of a control volume imbedding a mathematical surface(s) exactly corresponding to real surface(s). These investigators also sketched an alternative method, employing generalized functions, for its derivation. This latter method, which was later developed by Farassat (1975), is purely mathematical and formal. Goldstein (1976) used the free-space Green function to produce an implicit derivation of the Ffowcs Williams and Hawkings equation. In the study presented here, Lowson's (1965) concept of moving point singularities is generalized to moving surface singularities, and a new derivation is given of the fundamental equation. The derivation is based on topological considerations of the underlying space, the fluid medium, and the integral properties of the Dirac delta function.

Aggarwal, H. R.↗

Toward Higher-order Accuracy in Self-gravitating Hydrodynamics

High-order algorithms have emerged in numerical astrophysics as a promising avenue to reduce truncation error (proportional to a power of the linear resolution Δ x ) with only a moderate increase to computational expense. Significant effort has been placed in the development of finite-volume algorithms for (magneto)hydrodynamics; however, state-of-the-art astrophysical simulations tightly couple a plenitude of physics, additionally including gravity, photon transport, cosmic-ray transport, chemistry, and/or diffusion, to name a few. Algorithms frequently operator-split this additional physics (often a first-order error in time) and/or adopt a model wherein their evaluation is limited to second-order accuracy in space. In this work, we present a fourth-order-accurate finite-volume scheme for self-gravitating hydrodynamics on a uniform Cartesian grid. The method supplies source terms for the gravitational acceleration ( ρ g ) and gravitational energy release ( ρ v · g ) associated with fourth-order-accurate solutions to the Poisson equation. Our scheme (1) guarantees the conservation of total linear momentum while (2) decreasing (in proportion to Δ x 4 ) the effects of spurious heating and/or cooling associated with truncation error in the gravity. We demonstrate expected convergence rates for the algorithm by measuring errors in test problems evolving self-gravity modified linear waves and 3D polytropic equilibria. We test robustness of the algorithm by integrating an induced “inside-out” adiabatic collapse. We also discuss a method to smoothly downgrade the solution to second-order spatial accuracy to avoid spurious overshoots near steep density and/or pressure gradients.

79 ASTRONOMY AND ASTROPHYSICS↗

Numerical computation of viscous flow about unconventional airfoil shapes

A new two-dimensional computer code was developed to analyze the viscous flow around unconventional airfoils at various Mach numbers and angles of attack. The Navier-Stokes equations are solved using an implicit, upwind, finite-volume scheme. Both laminar and turbulent flows can be computed. A new nonequilibrium turbulence closure model was developed for computing turbulent flows. This two-layer eddy viscosity model was motivated by the success of the Johnson-King model in separated flow regions. The influence of history effects are described by an ordinary differential equation developed from the turbulent kinetic energy equation. The performance of the present code was evaluated by solving the flow around three airfoils using the Reynolds time-averaged Navier-Stokes equations. Excellent results were obtained for both attached and separated flows about the NACA 0012 airfoil, the RAE 2822 airfoil, and the Integrated Technology A 153W airfoil. Based on the comparison of the numerical solutions with the available experimental data, it is concluded that the present code in conjunction with the new nonequilibrium turbulence model gives excellent results.

Ahmed, S.↗