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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Kinetics of diffusional droplet growth in a liquid/liquid two-phase system

In the case of the diaphragm cell transport equation where the interdiffusion coefficient is a function of concentration, we have derived an integral of the form, t = B(sub 0) + B(sub L)ln(delta(c)) + B(sub 1)(delta(c)) + B(sub 2)(delta(c))(exp 2) +... where t is the time and (delta(c)) is the concentration difference across the frit. The coefficient, B(sub 0), is a constant of integration, while the coefficient, B(sub L), B(sub 1), B(sub 2), ..., depend in general upon the cell constant, the compartment volumes, the interdiffusion coefficient, and various of its concentration derivatives evaluated at the mean concentration for the cell. Explicit formulae for B(sub L), B(sub 1), B(sub 2), ... are given.

Baird, James K.↗
NASA NTRS2002

Spectral (Finite) Volume Method for Conservation Laws on Unstructured Grids II: Extension to Two Dimensional Scalar Equation

The framework for constructing a high-order, conservative Spectral (Finite) Volume (SV) method is presented for two-dimensional scalar hyperbolic conservation laws on unstructured triangular grids. Each triangular grid cell forms a spectral volume (SV), and the SV is further subdivided into polygonal control volumes (CVs) to supported high-order data reconstructions. Cell-averaged solutions from these CVs are used to reconstruct a high order polynomial approximation in the SV. Each CV is then updated independently with a Godunov-type finite volume method and a high-order Runge-Kutta time integration scheme. A universal reconstruction is obtained by partitioning all SVs in a geometrically similar manner. The convergence of the SV method is shown to depend on how a SV is partitioned. A criterion based on the Lebesgue constant has been developed and used successfully to determine the quality of various partitions. Symmetric, stable, and convergent linear, quadratic, and cubic SVs have been obtained, and many different types of partitions have been evaluated. The SV method is tested for both linear and non-linear model problems with and without discontinuities.

Wang, Z. J.↗
DOE OSTI2019

A computational model for nanosecond pulse laser-plasma interactions

A multi-physics numerical model for laser-induced optical breakdown and laser-plasma interaction in a non-equilibrium gas is presented, accounting for: production of priming electrons via multi-photon ionization, energy absorption, cascade ionization, induced hydrodynamic response, and shock formation and propagation. The gas is governed by the Navier-Stokes equations, with non-equilibrium effects taken into account by means of a two-temperature model. The space-time dependence of the laser beam is modeled with a flux-tube formulation for the Radiative Transfer Equation. The flow governing equations are discretized in space using a second-order finite volume method. The semi-discrete equations are marched in time using an implicit-explicit (IMEX) dual time-stepping strategy, where diffusion and chemistry are solved implicitly, whereas convection is explicit. Application to a 20 ns long 50 mJ pulse laser-induced breakdown in quiescent O 2 shows the advantages of this temporal discretization during and just after the laser pulse, while a less-expensive symmetric Strang splitting (with implicit chemistry) is sufficient for the post-breakdown gas dynamics after ≃ 0.1 μs. The integrated model is shown to reproduce key features of corresponding experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

An immersed interface method for the 2D vorticity-velocity Navier-Stokes equations with multiple bodies

We present an immersed interface method for the vorticity-velocity form of the 2D Navier Stokes equations that directly addresses challenges posed by nonconvex immersed bodies, multiply connected domains, and the calculation of force distributions on immersed surfaces. The immersed interface method is re-interpreted as a polynomial extrapolation of flow quantities and boundary conditions into the immersed solid bodies, reducing computational cost and enabling simulations with nonconvex bodies that could not be discretized with previous immersed interface methods. In the flow, the vorticity transport equation is discretized using a conservative finite difference scheme and explicit Runge-Kutta time integration. The velocity reconstruction problem is transformed to a scalar Poisson equation that is discretized with conservative finite differences, and solved using an FFT-accelerated iterative algorithm. The use of conservative differencing throughout leads to exact enforcement of a discrete Kelvin's theorem, allowing for simulations with multiply connected domains and outflow boundaries that have challenged other immersed interface vortex methods. We also explore novel methods for recovering time-dependent pressure distributions on immersed bodies within a vorticity-based method and present a novel control volume formulation for recovering aerodynamic moments from only the vorticity and velocity fields. The method achieves second order spatial accuracy and third order temporal accuracy, and is validated on a variety of 2D flows in internal and free-space domains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Analytic continuation of the relativistic three-particle scattering amplitudes

We investigate the relativistic scattering of three identical scalar bosons interacting via pair-wise interactions. Extending techniques from the nonrelativistic three-body scattering theory, we provide a detailed and general prescription for solving and analytically continuing integral equations describing the three-body reactions. We use these techniques to study a system with zero angular momenta described by a single scattering length leading to a bound state in a two-body subchannel. We obtain bound-state-particle and three-particle amplitudes in the previously unexplored kinematical regime; in particular, for real energies below elastic thresholds and complex energies in the physical and unphysical Riemann sheets. We extract positions of three-particle bound-states that agree with previous finite-volume studies, providing further evidence for the consistency of the relativistic finite-volume three-body quantization conditions. We also determine previously unobserved virtual bound states in this theory. Lastly, we find numerical evidence of the breakdown of the two-body finite-volume formalism in the vicinity of the left-hand cuts and argue for the generalization of the existing formalism.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1991

A general multiblock Euler code for propulsion integration. Volume 1: Theory document

A general multiblock Euler solver was developed for the analysis of flow fields over geometrically complex configurations either in free air or in a wind tunnel. In this approach, the external space around a complex configuration was divided into a number of topologically simple blocks, so that surface-fitted grids and an efficient flow solution algorithm could be easily applied in each block. The computational grid in each block is generated using a combination of algebraic and elliptic methods. A grid generation/flow solver interface program was developed to facilitate the establishment of block-to-block relations and the boundary conditions for each block. The flow solver utilizes a finite volume formulation and an explicit time stepping scheme to solve the Euler equations. A multiblock version of the multigrid method was developed to accelerate the convergence of the calculations. The generality of the method was demonstrated through the analysis of two complex configurations at various flow conditions. Results were compared to available test data. Two accompanying volumes, user manuals for the preparation of multi-block grids (vol. 2) and for the Euler flow solver (vol. 3), provide information on input data format and program execution.

Chen, H. C.↗
DOE OSTI2021

A weighted Shifted Boundary Method for free surface flow problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods and was recently introduced for the Poisson, linear advection/diffusion, Stokes, Navier-Stokes, acoustics, and shallow-water equations. By reformulating the original boundary value problem over a surrogate (approximate) computational domain, the SBM avoids integration over cut cells and the associated problematic issues regarding numerical stability and matrix conditioning. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we extend the SBM to the simulation of incompressible Navier-Stokes flows with moving free-surfaces, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach prevents spurious pressure oscillations in time, which would otherwise be produced if the total active fluid volume were to change abruptly over a time step. In fact, the proposed weighted SBM method induces small mass (i.e., volume) conservation errors, which converge quadratically in the case of piecewise-linear finite element interpolations, as the grid is refined. Finally, we present an extensive set of two- and three-dimensional tests to demonstrate the robustness and accuracy of the method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Shock-ramp analysis test problem

Quasi-isentropic (ramp) compression is now a well-established experimental method and so are the analysis techniques to give Lagrangian sound speed, pressure, and density along the sample material's isentrope. A shock followed by ramp compression is a natural extension to investigate, for example, shock melt and refreeze on compression, or isentropes of states off the Hugoniot or principal isentrope. In practice, graded-density impactors produce initial shocks, compression by shaped laser pulses may be unable to produce a smooth pressure increase from zero, and incidental perturbations on the drive pulse may also give rise to shocks, so robust shock-ramp analysis methods will be needed. Appropriate analysis methods are needed for shock-ramp experiments, based on those for quasi-isentropic compression, and these require validation. This paper describes three different analyses of a shock-ramp test problem, including an assessment of their estimated errors. The methods tested were based on hydrodynamic characteristics or integration backward in space. All methods gave the known Lagrangian sound speed to within ~1%, and pressure and volume to within less than 2% and 1%, demonstrating that the analysis methods of isentropic compression experiments can be confidently extended to the analysis of shock and ramp compression.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2001

Accurate Time-Dependent Traveling-Wave Tube Model Developed for Computational Bit-Error-Rate Testing

The phenomenal growth of the satellite communications industry has created a large demand for traveling-wave tubes (TWT's) operating with unprecedented specifications requiring the design and production of many novel devices in record time. To achieve this, the TWT industry heavily relies on computational modeling. However, the TWT industry's computational modeling capabilities need to be improved because there are often discrepancies between measured TWT data and that predicted by conventional two-dimensional helical TWT interaction codes. This limits the analysis and design of novel devices or TWT's with parameters differing from what is conventionally manufactured. In addition, the inaccuracy of current computational tools limits achievable TWT performance because optimized designs require highly accurate models. To address these concerns, a fully three-dimensional, time-dependent, helical TWT interaction model was developed using the electromagnetic particle-in-cell code MAFIA (Solution of MAxwell's equations by the Finite-Integration-Algorithm). The model includes a short section of helical slow-wave circuit with excitation fed by radiofrequency input/output couplers, and an electron beam contained by periodic permanent magnet focusing. A cutaway view of several turns of the three-dimensional helical slow-wave circuit with input/output couplers is shown. This has been shown to be more accurate than conventionally used two-dimensional models. The growth of the communications industry has also imposed a demand for increased data rates for the transmission of large volumes of data. To achieve increased data rates, complex modulation and multiple access techniques are employed requiring minimum distortion of the signal as it is passed through the TWT. Thus, intersymbol interference (ISI) becomes a major consideration, as well as suspected causes such as reflections within the TWT. To experimentally investigate effects of the physical TWT on ISI would be prohibitively expensive, as it would require manufacturing numerous amplifiers, in addition to acquiring the required digital hardware. As an alternative, the time-domain TWT interaction model developed here provides the capability to establish a computational test bench where ISI or bit error rate can be simulated as a function of TWT operating parameters and component geometries. Intermodulation products, harmonic generation, and backward waves can also be monitored with the model for similar correlations. The advancements in computational capabilities and corresponding potential improvements in TWT performance may prove to be the enabling technologies for realizing unprecedented data rates for near real time transmission of the increasingly larger volumes of data demanded by planned commercial and Government satellite communications applications. This work is in support of the Cross Enterprise Technology Development Program in Headquarters' Advanced Technology & Mission Studies Division and the Air Force Office of Scientific Research Small Business Technology Transfer programs.

Kory, Carol L.↗
NASA NTRS2013

Zero-gravity Mean Free Surface Curvature of a Confined Liquid in a Radially-Vaned Container

A variety of increasingly intricate container geometries are under consideration for the passive manipulation of liquids aboard spacecraft where the impact of gravity may be neglected. In this study we examine the mean curvature of a liquid volume confined within a radial array of disconnected vanes of infinite extent. This particular geometry possesses a number of desirable characteristics relevant to waste water treatment aboard spacecraft for life support. It is observed that under certain conditions the slender shape of the free surface approaches an asymptote, which can be predicted analytically using new hybrid boundary conditions proposed herein. This contribution represents possibly the final extension of what has been referred to as the method of de Lazzer et al. (1996). The method enables the integration of the Young-Laplace equation over a domain with its boundaries, including the wetted portion of the solid boundaries, symmetry planes, and circular arcs representing free surfaces at the center plane of the liquid body. Asymptotic solutions at several limits are obtained and the analysis is confirmed with numerical computations.

Chen, Yongkang↗
NASA NTRSDate not supplied

Molecular flow venting of a volume with an outgassing or desorbing source

A venting equation commonly used to describe the transient evolution of pressure within a volume containing an outgassing or offgassing source and a restrictive vent conductance under conditions of molecular flow has been solved analytically. Solutions are found for sources of finite thickness characterized by classical diffusion-limited behavior (proportional to inverse square root of time), as well as responses for thick material sources often observed in testing that are characterized by a more general form of power-law decay, up to inverse time behavior associated with surface desorption. Solutions involve evaluating integrals where both numerators and denominators of the integrands diverge with time, making wide-ranging transient solutions difficult to directly compute numerically. Usually, one can avoid evaluating these integrals by assuming quasistatic conditions at long time scales. A novel approach is used in this work to analytically produce solutions by generating bespoke mathematical functions, some of which solve integrals that have apparently had no previous analytical solution.

Michael Woronowicz↗
NASA NTRSDate not supplied

Molecular flow venting of a volume with an outgassing or desorbing source

"A venting equation commonly used to describe the transient evolution of pressure within a volume containing an outgassing or offgassing source and a restrictive vent conductance under conditions of molecular flow has been solved analytically. Solutions are found for sources of finite thickness characterized by classical diffusion-limited behavior (proportional to inverse square root of time), as well as responses for thick material sources often observed in testing that are characterized by a more general form of power-law decay, up to inverse time behavior associated with surface desorption. Solutions involve evaluating integrals where both numerators and denominators of the integrands diverge with time, making wide-ranging transient solutions difficult to directly compute numerically. Usually, one can avoid evaluating these integrals by assuming quasistatic conditions at long time scales. A novel approach is used in this work to analytically produce solutions by generating bespoke mathematical functions, some of which solve integrals that have apparently had no previous analytical solution."

Michael Woronowicz↗
NASA NTRS1992

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

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↗
NASA NTRS2022

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↗
NASA NTRS1981

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

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↗
NASA NTRS2017

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↗