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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Nonlinear mesomechanics of composites with periodic microstructure

This work is concerned with modeling the mechanical deformation or constitutive behavior of composites comprised of a periodic microstructure under small displacement conditions at elevated temperature. A mesomechanics approach is adopted which relates the microimechanical behavior of the heterogeneous composite with its in-service macroscopic behavior. Two different methods, one based on a Fourier series approach and the other on a Green's function approach, are used in modeling the micromechanical behavior of the composite material. Although the constitutive formulations are based on a micromechanical approach, it should be stressed that the resulting equations are volume averaged to produce overall effective constitutive relations which relate the bulk, volume averaged, stress increment to the bulk, volume averaged, strain increment. As such, they are macromodels which can be used directly in nonlinear finite element programs such as MARC, ANSYS and ABAQUS or in boundary element programs such as BEST3D. In developing the volume averaged or efective macromodels from the micromechanical models, both approaches will require the evaluation of volume integrals containing the spatially varying strain distributions throughout the composite material. By assuming that the strain distributions are spatially constant within each constituent phase-or within a given subvolume within each constituent phase-of the composite material, the volume integrals can be obtained in closed form. This simplified micromodel can then be volume averaged to obtain an effective macromodel suitable for use in the MARC, ANSYS and ABAQUS nonlinear finite element programs via user constitutive subroutines such as HYPELA and CMUSER. This effective macromodel can be used in a nonlinear finite element structural analysis to obtain the strain-temperature history at those points in the structure where thermomechanical cracking and damage are expected to occur, the so called damage critical points of the structure.

Walker, Kevin P.↗
NASA NTRS1992

Finite element Euler calculations of unsteady transonic cascade flows

A Galerkin finite element procedure incorporating an explicit Runge-Kutta time-stepping scheme has been developed in this work to solve unsteady transonic flow in cascades. The computational domain is discretized by a globally unstructured but locally structured blade-fitted deformable mesh. The Galerkin approximation is applied to the unsteady Euler equations based on a mixed Eulerian-Lagrangian description. The semi-discretized equations are integrated forward in time using a multistage Runge-Kutta scheme. An artificial dissipation operator of the type proposed by Jameson is adapted in the current scheme to capture shocks and suppress nonphysical oscillations. Phase-shifted boundary conditions are used to reduce the computational domain to a single reference passage. Results for both steady and unsteady transonic flows through cascades are presented and compared to existing finite volume solutions.

Hsiao, Chingteng↗
NASA NTRS1995

An approximate Kalman filter for ocean data assimilation: An example with an idealized Gulf Stream model

A practical method of data assimilation for use with large, nonlinear, ocean general circulation models is explored. A Kalman filter based on approximation of the state error covariance matrix is presented, employing a reduction of the effective model dimension, the error's asymptotic steady state limit, and a time-invariant linearization of the dynamic model for the error integration. The approximations lead to dramatic computational savings in applying estimation theory to large complex systems. We examine the utility of the approximate filter in assimilating different measurement types using a twin experiment of an idealized Gulf Stream. A nonlinear primitive equation model of an unstable east-west jet is studied with a state dimension exceeding 170,000 elements. Assimilation of various pseudomeasurements are examined, including velocity, density, and volume transport at localized arrays and realistic distributions of satellite altimetry and acoustic tomography observations. Results are compared in terms of their effects on the accuracies of the estimation. The approximate filter is shown to outperform an empirical nudging scheme used in a previous study. The examples demonstrate that useful approximate estimation errors can be computed in a practical manner for general circulation models.

Fukumori, Ichiro↗
NASA NTRS1983

Computational technique for three-dimensional compressible flow past wings at high angles of attack

A computational technique, based on the integral solution of the full potential equation, has been developed for the solution of three-dimensional high-subsonic flows past wings at high angles of attack. The problem includes two sources of nonlinearities: a boundary oriented nonlinearity (separated flow roll up) and a region oriented nonlinearity (flow compressibility). The former is represented by a nonlinear vortex lattice while the latter is represented by a source distribution inside a finite volume. The solution is obtained by using double iteration cycles; a separated flow (wake) iteration cycle and a compressibility iteration cycle. The computational technique is applied to a delta wing. The results show that the technique is very promising and efficient.

Kandil, O. A.↗
NASA NTRS1987

The area-time-integral technique to estimate convective rain volumes over areas applied to satellite data - A preliminary investigation

The use of the area-time-integral (ATI) technique, based only on satellite data, to estimate convective rain volume over a moving target is examined. The technique is based on the correlation between the radar echo area coverage integrated over the lifetime of the storm and the radar estimated rain volume. The processing of the GOES and radar data collected in 1981 is described. The radar and satellite parameters for six convective clusters from storm events occurring on June 12 and July 2, 1981 are analyzed and compared in terms of time steps and cluster lifetimes. Rain volume is calculated by first using the regression analysis to generate the regression equation used to obtain the ATI; the ATI versus rain volume relation is then employed to compute rain volume. The data reveal that the ATI technique using satellite data is applicable to the calculation of rain volume.

Doneaud, Andre A.↗
NASA NTRS1994

An Approximate Kalman Filter for Ocean Data Assimilation; An Example with an Idealized Gulf Stream Model

A practical method of data assimilation for use with large, nonlinear, ocean general circulation models is explored. A Kalman filter based on approximations of the state error covariance matrix is presented, employing a reduction of the effective model dimension, the error's asymptotic steady-state limit, and a time-invariant linearization of the dynamic model for the error integration. The approximations lead to dramatic computational savings in applying estimation theory to large complex systems. We examine the utility of the approximate filter in assimilating different measurement types using a twin experiment of an idealized Gulf Stream. A nonlinear primitive equation model of an unstable east-west jet is studied with a state dimension exceeding 170,000 elements. Assimilation of various pseudo measurements is examined, including velocity, density, and volume transport at localized arrays, and realistic distributions of satellite altimetry and acoustic tomography observations. Results are compared in terms of their effects on the accuracies of the estimation. The approximate filter is shown to outperform a previous study that used an empirical nudging scheme. The examples demonstrate that useful approximate estimation errors can be computed in a practical manner for general circulation models.

['Kalman↗
NASA NTRSDate not supplied

Wall-Modeled Large Eddy Simulation Method for Unstructured-Grid Navier-Stokes Solvers

This paper reports on the implementation and assessment of a Wall-Modeled Large-Eddy Simulation (WMLES) methodology in an unstructured-grid, node-centered flow solver, FUN3D that is developed and supported at the NASA Langley Research Center. Finite-volume (FV) and finite-element (FE) discretization schemes considered in the study provide formal second-order spatial accuracy. Large-Eddy Simulations (LES) resolve large-scale turbulent-flow features and filter out small-scale effects using the Vreman subgrid-scale model. At solid-wall boundaries, a shear-stress model is employed to provide a proper boundary-flux closure. The nonlinear equations are integrated in time using either an optimized backward difference formula or an implicit multistage Runge-Kutta temporal scheme. The implicit equations at each time step are solved by strong nonlinear iteration schemes. WMLES demonstrations are shown for two high-lift configurations, namely, the McDonnell Douglas 30P30N multielement airfoil and a NASA High-Lift Common Research Model. Results show that the WMLES approaches implemented in the FV and FE discretization methods produce consistent solutions and are capable of capturing key aerodynamic characteristics and flow structures for high-lift configurations at a wide range of angles of attack including maximum-lift conditions. In the 30P30N example, correct trends in the variations of integrated aerodynamic forces and moments, surface pressure distributions, and boundary-layer profiles are captured as the Reynolds number is increased.

CFD; turbulence modeling; High-Lift flow simulatio↗
NASA NTRS1986

Upwind relaxation algorithms for Euler/Navier-Stokes equations

A description of and results from a solution algorithm for the compressible Navier-Stokes equations are presented. The main features of the algorithm are second or third order accurate upwind discretization of the convection and pressure derivatives and a relaxation scheme for the unfactored implicit backward Euler time method, implemented in a finite-volume formulation. Upwind methods were successfully used to obtain solutions to the Euler equations for flows with strong shock waves. The particular upwind method being used is based on the flux vector splitting technique developed by Van Leer and both second and third order accurate discretizations were developed. Currently, the most widely used implicit solution technique for the Navier-Stokes equations use approximate factorization (AF) methods to treat multidimensional problems. The time integration scheme being used in the present algorithm corresponds to a line Gauss-Seidel relaxation method. This method produces good convergence rates for steady-state flows, and most of the algorithm was vectorized on the NASA Langley VPS 32 computer. The Navier-Stokes algorithm was tested for several two-dimensional flow problems. Solutions for the problems gave excellent results. The presented effort is directed toward the extension of the scheme to the full three-dimensional Navier-Stokes equations.

Thomas, J. L.↗
NASA NTRS1991

Effects of leading-edge flap oscillation on unsteady delta wing flow and rock control

The isolated and interdisciplinary problems of unsteady fluid dynamics and rigid-body dynamics and control of delta wings with and without leading-edge flap oscillation are considered. For the fluid dynamics problem, the unsteady, compressible, thin-layer Navier-Stokes (NS) equations, which are written relative to a moving frame of reference, are solved along with the unsteady, linearized, Navier-displacement (ND) equations. The NS equations are solved for the flowfield using an implicit finite-volume scheme. The ND equations are solved for the grid deformation, if the leading-edge flaps oscillate, using an ADI scheme. For the dynamics and control problem, the Euler equation of rigid-body rolling motion for a wing and its flaps are solved interactively with the fluid dynamics equations for the wing-rock motion and subsequently for its control. A four-stage Runge-Kutta scheme is used to explicitly integrate the dynamics equation.

Kandil, Osama A.↗
NASA NTRS1992

Basic studies in microwave remote sensing

Scattering models were developed in support of microwave remote sensing of earth terrains with particular emphasis on model applications to airborne Synthetic Aperture Radar measurements of forest. Practically useful surface scattering models based on a solution of a pair of integral equations including multiple scattering effects were developed. Comparisons of these models with controlled scattering measurements from statistically known random surfaces indicate that they are valid over a wide range of frequencies. Scattering models treating a forest environment as a two and three layered media were also developed. Extensive testing and comparisons were carried out with the two layered model. Further studies with the three layered model are being carried out. A volume scattering model valid for dense media such as a snow layer was also developed that shows the appropriate trend dependence with the volume fraction of scatterers.

Fung, Adrian K.↗
DOE OSTI2025

A Block-Structured Adaptive Mesh Framework to Solve Radiation Transfer Equation in Irregular Embedded Geometries

Radiation transport arises in various scientific, industrial, and medical fields, and understanding its effect in applications is needed to make accurate predictions, safety assessments and performance optimizations. Solving the Radiation Transport Equation (RTE) is challenging due to its integro-differential nature, which involves both differential and integral terms. The differential term describes the change in radiation intensity due to absorption and emission, while the integral term accounts for scattering. The accurate modeling of radiation is further complicated in many applications due to the complex, irregular geometries. Various methods exist for solving the RTE, including the zonal, Monte Carlo, spherical harmonics, discrete ordinates, and finite volume methods. Traditional mesh-based approaches, which rely on structured or unstructured meshes, struggle with irregular geometries due to: a) the difficulty of conforming structured grids to irregular domains, b) challenges in enforcing boundary conditions correctly, and c) the additional computational cost of unstructured mesh methods. This work presents a second-order accurate method for solving the RTE in irregular geometries. The radiation intensity is discretized using the finite-volume method in both spatial and angular directions on regular Cartesian grid blocks. Leveraging the block-structured adaptive mesh refinement (AMR) framework provided by AMReX, our method refines the grid locally to reduce spatial discretization error, ensuring a converged numerical solution while minimizing computational costs elsewhere. A two-stage deferred correction approach is employed: First, a first-order discretization on grid blocks is solved using an algebraic multigrid method in HYPRE. Second, a correction term is applied explicitly to achieve second-order accuracy. The correction term is calculated by approximating the radiation flux on cell faces using a Total Variation Diminishing (TVD) scheme. This approach ensures quick convergence of the multigrid method while preserving higher-order accuracy of the numerical solution. Irregular geometries are resolved as embedded boundaries (EB), resulting in both cut cells and regular cells. In cut cells, we modify the fluxes using face fractions and incorporate additional contributions from EB boundary conditions. To ensure higher-order convergence near the EB interface, the correction term is modified by interpolating the radiation intensity to fictitious ghost points. The implementation takes advantage of modern supercomputers by leveraging AMReX’sMPI/X parallelization strategy where X can be MPI or a GPU accelerator including CUDA, HIP and DPC++. We validate our solver using classical test cases, both with and without EB, demonstrating accuracy and efficiency. Additionally, we analyze the impact of adaptive mesh refinement on solution accuracy and computational cost, highlighting the advantages of our approach for high-resolution radiation transport simulations.

computational fluid dynamics (CFD)↗
NASA NTRS1967

Program for calculating optimum dimensions of alpha radioisotope capsules exposed to varying stress and temperature

A method and computer program were developed for calculating the creep and optimizing the dimensions of capsules filled with alpha-emitting radioisotopes. The method solves an integral equation that was developed assuming linear accumulation of partial creep lives and relating life to time-dependent stress and temperature using the Larson-Miller parameter. The computer program, CAPSUL, is written in Fortran language for the IBM 360/75 computer. The program makes a least squares fit of the creep life function using conventional constant stress, constant temperature creep data. Dimensions of capsules having maximum thermal power per unit of weight, volume, or area are calculated for a given creep life and pressure-temperature history using a numerical Lagrange Multiplier formulation. The program also calculates the life to a prescribed strain for capsules of given dimensions and pressure-temperature history. The method has been used to analyze creep data for the alloys 304 stainless steel, Hastelloy N, Cb-1% Zr, FS-85, and T-222.

J. P. Nichols↗
NASA NTRS1959

A Numerical Method for Obtaining Monoenergetic Neutron Flux Distributions and Transmissions in Multiple-Region Slabs

This method is investigated for semi-infinite multiple-slab configurations of arbitrary width, composition, and source distribution. Isotropic scattering in the laboratory system is assumed. Isotropic scattering implies that the fraction of neutrons scattered in the i(sup th) volume element or subregion that will make their next collision in the j(sup th) volume element or subregion is the same for all collisions. These so-called "transfer probabilities" between subregions are calculated and used to obtain successive-collision densities from which the flux and transmission probabilities directly follow. For a thick slab with little or no absorption, a successive-collisions technique proves impractical because an unreasonably large number of collisions must be followed in order to obtain the flux. Here the appropriate integral equation is converted into a set of linear simultaneous algebraic equations that are solved for the average total flux in each subregion. When ordinary diffusion theory applies with satisfactory precision in a portion of the multiple-slab configuration, the problem is solved by ordinary diffusion theory, but the flux is plotted only in the region of validity. The angular distribution of neutrons entering the remaining portion is determined from the known diffusion flux and the remaining region is solved by higher order theory. Several procedures for applying the numerical method are presented and discussed. To illustrate the calculational procedure, a symmetrical slab ia vacuum is worked by the numerical, Monte Carlo, and P(sub 3) spherical harmonics methods. In addition, an unsymmetrical double-slab problem is solved by the numerical and Monte Carlo methods. The numerical approach proved faster and more accurate in these examples. Adaptation of the method to anisotropic scattering in slabs is indicated, although no example is included in this paper.

Schneider, Harold↗
NASA NTRS2015

Numerical Calibration of Mass Flow Plug for Inlet Testing

A simple control volume model has been developed to calculate the discharge coefficient through a mass flow plug (MFP) and validated with a calibration experiment. The maximum error of the model within the operating region of the MFP is 0.54%. The control volume analysis developed work is comprised of a sequence of flow calculations through the MFP. The model uses the MFP geometry and operating pressure and temperature to couple continuity, momentum, energy, an equation of state, and wall shear. The discharge coefficient calculation also includes the effects of boundary layer growth, including the reduction in cross-sectional flow area as characterized by the boundary layer displacement thickness. The last calculation in the sequence uses an integral method to calculate the growth of the boundary layer, from which the displacement thickness is then determined. The result of these successive calculations is an accurate one-dimension model of the velocity, pressure, and temperature through the MFP. For comparison, a computational fluid dynamic (CFD) calibration is shown, which when compared to the presented numerical model, had a lower accuracy with a maximum error of 1.35% in addition to being slower by a factor of 100."

Flow Measurement↗
NASA NTRS2009

Spacecraft Internal Acoustic Environment Modeling

Acoustic modeling can be used to identify key noise sources, determine/analyze sub-allocated requirements, keep track of the accumulation of minor noise sources, and to predict vehicle noise levels at various stages in vehicle development, first with estimates of noise sources, later with experimental data. In FY09, the physical mockup developed in FY08, with interior geometric shape similar to Orion CM (Crew Module) IML (Interior Mode Line), was used to validate SEA (Statistical Energy Analysis) acoustic model development with realistic ventilation fan sources. The sound power levels of these sources were unknown a priori, as opposed to previous studies that RSS (Reference Sound Source) with known sound power level was used. The modeling results were evaluated based on comparisons to measurements of sound pressure levels over a wide frequency range, including the frequency range where SEA gives good results. Sound intensity measurement was performed over a rectangular-shaped grid system enclosing the ventilation fan source. Sound intensities were measured at the top, front, back, right, and left surfaces of the and system. Sound intensity at the bottom surface was not measured, but sound blocking material was placed tinder the bottom surface to reflect most of the incident sound energy back to the remaining measured surfaces. Integrating measured sound intensities over measured surfaces renders estimated sound power of the source. The reverberation time T6o of the mockup interior had been modified to match reverberation levels of ISS US Lab interior for speech frequency bands, i.e., 0.5k, 1k, 2k, 4 kHz, by attaching appropriately sized Thinsulate sound absorption material to the interior wall of the mockup. Sound absorption of Thinsulate was modeled in three methods: Sabine equation with measured mockup interior reverberation time T60, layup model based on past impedance tube testing, and layup model plus air absorption correction. The evaluation/validation was carried out by acquiring octave band microphone data simultaneously at ten fixed locations throughout the mockup. SPLs (Sound Pressure Levels) predicted by our SEA model match well with measurements for our CM mockup, with a more complicated shape. Additionally in FY09, background NC noise (Noise Criterion) simulation and MRT (Modified Rhyme Test) were developed and performed in the mockup to determine the maximum noise level in CM habitable volume for fair crew voice communications. Numerous demonstrations of simulated noise environment in the mockup and associated SIL (Speech Interference Level) via MRT were performed for various communities, including members from NASA and Orion prime-/sub-contractors. Also, a new HSIR (Human-Systems Integration Requirement) for limiting pre- and post-landing SIL was proposed.

Chu, SShao-sheng R.↗
DOE OSTI2025

Implicit full-F simulations of neoclassical ion transport

The development of implicit time integration capabilities for axisymmetric full-F continuum simulations of ion neoclassical transport is reported. The approach involves the implicit treatment of the gyrokinetic Vlasov equation coupled to the nonlinear Fokker–Planck collision model in the long-wavelength limit approximation. To facilitate implicit simulations, advanced preconditioning of individual physics operators is developed, and a global multi-physics preconditioner is constructed by adopting an operator splitting methodology. The algorithm is implemented in the finite-volume code COGENT and is applied to study neoclassical transport properties for both the main ion species and the lithium impurity species in the closed-field-line region of the LTX- β tokamak. The implicit COGENT simulations elucidate the role of non-local transport effects, while demonstrating substantial speedup over the corresponding explicit approach.

Dorf, Mikhail [Lawrence Livermore National Laborat↗
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 NTRS1994

Design of 3-D Nacelle near Flat-Plate Wing Using Multiblock Sensitivity Analysis (ADOS)

One of the major design tasks involved in reducing aircraft drag is the integration of the engine nacelles and airframe. With this impetus, nacelle shapes with and without the presence of a flat-plate wing nearby were optimized. This also served as a demonstration of the 3-D version of the recently developed aerodynamic design optimization methodology using sensitivity analysis, ADOS. The required flow analyses were obtained by solving the three-dimensional, compressible, thin-layer Navier-Stokes equations using an implicit, upwind-biased, finite volume scheme. The sensitivity analyses were performed using the preconditioned version of the SADD scheme (sensitivity analysis on domain decomposition). In addition to demonstrating the present method's capability for automatic optimization, the results offered some insight into two important issues related to optimizing the shapes of multicomponent configurations in close proximity. First, inclusion of the mutual interference between the components resulted in a different shape as opposed to shaping an isolated component. Secondly, exclusion of the viscous effects compromised not only the flow physics but also the optimized shapes even for isolated components.

Eleshaky, Mohamed E.↗