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

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.↗

LS-DYNA User-Defined Internal Ballistic Modeling

- LS-DYNA Explicit Finite Element Analysis software and CADPROG Internal Ballistic Analysis code are integrated as a single modeling tool to analyze pyro-mechanical devices. - Combining these two codes allows the ballistic calculations to be performed without modeling the propellant in an FEA mesh. At each LS-DYNA time step, CADPROG is analytically run using the kinematic data (stroke, velocity, chamber volume) fed back from the LS-DYNA simulation. Calculated pressure is then applied to a piston or any applicable surface mesh in LS-DYNA model interactively. -All external loads such as friction, damping loads, assisting/resisting loads, locking/unlocking loads are modeled in LS-DYNA and coupled with pressure loads from CADPROG to create a fully defined equation of motion -Both codes are written in Fortran 77. CADPROG is inserted in LS-DYNA source code as a user-defined subroutine and re-compiled to generate a custom executable.

pyrotechnics↗

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

Finite-volume applications of high-order accurate ENO schemes to two-dimensional boundary-value problems are studied. These schemes achieve high-order spatial accuracy, in smooth regions, by a piecewise polynomial approximation of the solution from cell averages. In addition, this spatial operation involves an adaptive stencil algorithm in order to avoid the oscillatory behavior that is associated with interpolation across steep gradients. High-order TVD Runge-Kutta methods are employed for time integration, thus making these schemes best suited for unsteady problems. Fifth- and sixth-order accurate applications are validated through a grid refinement study involving the solutions of scalar hyperbolic equations. A previously proposed extension for the Euler equations of gas dynamics is tested, including its application to solutions of boundary-value problems involving solid walls and curvilinear coordinates.

Casper, Jay↗

SAM-ML: Integrating data-driven closure with nuclear system code SAM for improved modeling capability

Advanced reactors often involve complicated thermal-fluid (T-F) phenomena. Modeling such phenomena with the traditional one-dimensional (1-D) system code is a challenging task. The System Analysis Module (SAM), a modern nuclear system code, has developed a coarse mesh multi-dimensional (multi-D) flow model to capture the spatial effect of T-F phenomena in advanced reactors. As a coarse mesh solver, constitutive relations are required for SAM's multi-D model for unresolved fine-scale physics, such as turbulence. Here this work presents a novel approach that integrates neural networks as data-driven closure for SAM's multi-D flow model. The data-driven closure is trained with fine-resolution data to ensure its accuracy while maintaining a coarse mesh setup to ensure its efficiency and consistency with SAM. We demonstrate the applicability of this SAM-ML capability in an open volume thermal stratification problem, where a neural network model serves as the eddy viscosity closure. A customized interface between the neural network and SAM is developed to ensure flexible and efficient data exchange. The SAM-ML results demonstrate superior performance compared to SAM's built-in zero-equation eddy viscosity closure. The case study shows that although the generalization capability of the data-driven closure still needs to be improved for different transient case or different geometric setup, SAM -ML demonstrates good potential for challenging simulation problems with improved accuracy and computational efficiency.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Determination of constant-volume balloon capabilities for aeronautical research

The proper application of constant-volume balloons (CVB) for measurement of atmospheric phenomena was determined. And with the proper interpretation of the resulting data. A literature survey covering 176 references is included. the governing equations describing the three-dimensional motion of a CVB immersed in a flow field are developed. The flowfield model is periodic, three-dimensional, and nonhomogeneous, with mean translational motion. The balloon motion and flow field equations are cast into dimensionless form for greater generality, and certain significant dimensionless groups are identified. An alternate treatment of the balloon motion, based on first-order perturbation analysis, is also presented. A description of the digital computer program, BALLOON, used for numerically integrating the governing equations is provided.

Tatom, F. B.↗

Slow manifolds of classical Pauli particle enable structure-preserving geometric algorithms for guiding center dynamics

Since variational symplectic integrators for the guiding center was proposed, structure-preserving geometric algorithms have become an active research field in plasma physics. We found that the slow manifolds of the classical Pauli particle enable a family of structure-preserving geometric algorithms for guiding center dynamics with long-term stability and accuracy. This discovery overcomes the difficulty associated with the unstable parasitic modes for variational symplectic integrators when applied to the degenerate guiding center Lagrangian. It is a pleasant surprise that Pauli's Hamiltonian for electrons, which predated the Dirac equation and marks the beginning of particle physics, reappears in classical physics as an effective algorithm for solving an important plasma physics problem. This technique is applicable to other degenerate Lagrangians reduced from regular Lagrangians.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Elasto visco-plastic flow with special attention to boundary conditions

A simple but nontrivial steady-state creeping elasto visco-plastic (Maxwell fluid) radial flow problem is analyzed, with special attention given to the effects of the boundary conditions. Solutions are obtained through integration of a governing equation on stress using the Runge-Kutta method for initial value problems and finite differences for boundary value problems. A more general approach through the finite element method, an approach that solves for the velocity field rather than the stress field and that is applicable to a wide range of problems, is presented and tested using the radial flow example. It is found that steady-state flows of elasto visco-plastic materials are strongly influenced by the state of stress of material as it enters the region of interest. The importance of this boundary or initial condition in analyses involving materials coming into control volumes from unusual stress environments is emphasized.

Shimazaki, Y.↗