SEARCH · Search NASA

Results for “Integral Equations Finite Element”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

NASA NTRS1972

Assessment of advances in structural analysis for space shuttle

The latest technology for structural analysis in relation to the design tasks that lie ahead for the space shuttle is reviewed. For shell-of-revolution structures, the analysis can be formulated as a one-dimensional problem which is readily solved by using finite-difference or numerical-integration techniques. For more general asymmetric shells, a two-dimensional formulation is required. However, the governing equations are readily formulated and are amenable to solution by finite-difference techniques. For a completely general structural arrangement, such as structural frameworks, recourse is usually made to discretized formulations using finite elements. Of course, the finite-element programs could be used for shell structures, but at a loss in accuracy and increase in computer time compared with the special purpose programs.

Dixon, S. C.↗

NASA NTRS1988

MHOST: An efficient finite element program for inelastic analysis of solids and structures

An efficient finite element program for 3-D inelastic analysis of gas turbine hot section components was constructed and validated. A novel mixed iterative solution strategy is derived from the augmented Hu-Washizu variational principle in order to nodally interpolate coordinates, displacements, deformation, strains, stresses and material properties. A series of increasingly sophisticated material models incorporated in MHOST include elasticity, secant plasticity, infinitesimal and finite deformation plasticity, creep and unified viscoplastic constitutive model proposed by Walker. A library of high performance elements is built into this computer program utilizing the concepts of selective reduced integrations and independent strain interpolations. A family of efficient solution algorithms is implemented in MHOST for linear and nonlinear equation solution including the classical Newton-Raphson, modified, quasi and secant Newton methods with optional line search and the conjugate gradient method.

Nakazawa, S.↗

DOE OSTI2022

Integration of Electromagnetic Geophysics Forward Simulation in Coupled Flow and Geomechanics for Monitoring a Gas Hydrate Deposit Located in the Ulleung Basin, East Sea, Korea

We investigate the feasibility of electromagnetic (EM) geophysics methods to detect the dissociation of gas hydrate specifically from a gas hydrate deposit located in the Ulleung Basin, East Sea, Korea via an integrated flow-geomechanics-EM geophysics simulation. To this end, coupled flow and geomechanics simulation is first performed with the multiple porosity model employed, where a mixed formulation with the finite volume (FV) and finite element (FE) methods are taken for the flow and geomechanics, respectively. From the saturation and porosity fields obtained from the coupled flow and geomechanics, the electrical conductivity model is established for the EM simulation. Solving the partial differential equation of electrical diffusion which is linearized using the 3D finite element method (FEM), the EM fields are then computed. For numerical experiments, particularly two approaches in the configuration for the EM methods are compared in this contribution: the surface-to-surface and the surface-to-borehole methods. When the surface-to-surface EM method is employed, the EM is found to be less sensitive, implying low detectability. Especially for the short term of production, the low detectability is attributed to the similarity of electrical resistivity between the dissociated gas (CH4) and hydrate as well as the specific dissociation pattern within the intercalated composites of the field. On the other hand, when the surface-to-borehole EM method is employed, its sensitivity to capture the produced gas flow is improved, confirming its detectability in monitoring gas flow. Hence, the EM geophysics simulation integrated with coupled flow and geomechanics can be a potential tool for monitoring gas hydrate deposits.

depressurization↗

NASA NTRS1989

Analysis of finite phased arrays of circular microstrip patches

A method is presented for analyzing a finite planar array of circular microstrip patches fed by coaxial probes. The self- and mutual impedances between array elements are calculated using the method of moments with the dyadic Green's function for a dielectric layer on a ground plane. The patch circuits are determined by using the reaction integral equation. The active input impedance as well as the active element pattern of the array are computed from a knowledge of the resultant patch currents. The calculated results for two-element and eight-element linear arrays are in good agreement with experimental data. The active reflection coefficient and element pattern for the center and edge elements of a two-dimensional array as a function of scan angle are also presented.

Deshpande, Manohar D.↗

NASA NTRS1999

Development of Analysis Methods for Designing with Composites

The project involved the development of new analysis methods to achieve efficient design of composite structures. We developed a complex variational formulation to analyze the in-plane and bending coupling response of an unsymmetrically laminated plate with an elliptical cutout subjected to arbitrary edge loading as shown in Figure 1. This formulation utilizes four independent complex potentials that satisfy the coupled in-plane and bending equilibrium equations, thus eliminating the area integrals from the strain energy expression. The solution to a finite geometry laminate under arbitrary loading is obtained by minimizing the total potential energy function and solving for the unknown coefficients of the complex potentials. The validity of this approach is demonstrated by comparison with finite element analysis predictions for a laminate with an inclined elliptical cutout under bi-axial loading.The geometry and loading of this laminate with a lay-up of [-45/45] are shown in Figure 2. The deformed configuration shown in Figure 3 reflects the presence of bending-stretching coupling. The validity of the present method is established by comparing the out-of-plane deflections along the boundary of the elliptical cutout from the present approach with those of the finite element method. The comparison shown in Figure 4 indicates remarkable agreement. The details of this method are described in a manuscript by Madenci et al. (1998).

Madenci, E.↗

NASA NTRS2001

Computational Aeroacoustics by the Space-time CE/SE Method

In recent years, a new numerical methodology for conservation laws-the Space-Time Conservation Element and Solution Element Method (CE/SE), was developed by Dr. Chang of NASA Glenn Research Center and collaborators. In nature, the new method may be categorized as a finite volume method, where the conservation element (CE) is equivalent to a finite control volume (or cell) and the solution element (SE) can be understood as the cell interface. However, due to its rigorous treatment of the fluxes and geometry, it is different from the existing schemes. The CE/SE scheme features: (1) space and time treated on the same footing, the integral equations of conservation laws are solve( for with second order accuracy, (2) high resolution, low dispersion and low dissipation, (3) novel, truly multi-dimensional, simple but effective non-reflecting boundary condition, (4) effortless implementation of computation, no numerical fix or parameter choice is needed, an( (5) robust enough to cover a wide spectrum of compressible flow: from weak linear acoustic waves to strong, discontinuous waves (shocks) appropriate for linear and nonlinear aeroacoustics. Currently, the CE/SE scheme has been developed to such a stage that a 3-13 unstructured CE/SE Navier-Stokes solver is already available. However, in the present paper, as a general introduction to the CE/SE method, only the 2-D unstructured Euler CE/SE solver is chosen as a prototype and is sketched in Section 2. Then applications of the CE/SE scheme to linear, nonlinear aeroacoustics and airframe noise are depicted in Sections 3, 4, and 5 respectively to demonstrate its robustness and capability.

Loh, Ching Y.↗

DOE OSTI2024

A large deformation multiphase continuum mechanics model for shock loading of soft porous materials

A large deformation, coupled finite-element (FE) model is developed to simulate the multiphase response of soft porous materials subjected to high strain-rate loading. The approach is based on the theory of porous media (TPM) at large deformations. Simplifications to the one-dimensional regime studied in the numerical simulations follow. An overview of several different time integration schemes is presented for the purpose of solving the nonlinear dynamic coupled balance of momenta (mixture and fluid) and balance of mass of the mixture equations. Numerical examples are presented for (i) verification against closed-form analytical solutions assuming small loads, (ii) demonstrating large deformation effects at high strain-rate, and (iii) showing differences in deformations between a single-phase elastodynamics model with occluded compressible pore fluid and a multiphase poroelastodynamics model at high strain-rate. The multiphase model shows that the relative motion of the pore fluid significantly dampens the deformation response of the solid skeleton as compared to the single-phase model, and makes it possible to extract quantitative values for the stresses of the different constituents, thereby allowing one to form preliminary conclusions about the onset of damage in the solid skeleton. The novelty of the current work is developing a multiphase, large deformation, mixture theory numerical model for high strain-rate loading of soft porous materials. It was discovered that explicit, adaptive time-stepping Runge–Kutta schemes offer high accuracy at relatively low cost when compared to traditional implicit or explicit central difference time-stepping schemes for shock-like loadings. Here, shock viscosity is added to the mixture momentum balance equation to regularize the shock front, and a stabilization term is added to the mixture mass balance equation to stabilize equal order interpolation finite elements for the coupled finite element solution of multiphase materials.

Engineering↗

NASA NTRS1996

Simulation of Spiral Slot Antennas on Composite Platforms

The project goals, plan and accomplishments up to this point are summarized in the viewgraphs. Among the various accomplishments, the most important have been: the development of the prismatic finite element code for doubly curved platforms and its validation with many different antenna configurations; the design and fabrication of a new slot spiral antennas suitable for automobile cellular, GPS and PCs communications; the investigation and development of various mesh truncation schemes, including the perfectly matched absorber and various fast integral equation methods; and the introduction of a frequency domain extrapolation technique (AWE) for predicting broadband responses using only a few samples of the response. This report contains several individual reports most of which have been submitted for publication to referred journals. For a report on the frequency extrapolation technique, the reader is referred to the UM Radiation Laboratory report A total of 14 papers have been published or accepted for publication with the full or partial support of this grant. Several more papers are in preparation.

Volakis, John L.↗

NASA NTRS1998

Closed Form Equations for the Preliminary Design of a Heat-Pipe-Cooled Leading Edge

A set of closed form equations for the preliminary evaluation and design of a heat-pipe-cooled leading edge is presented. The set of equations can provide a leading-edge designer with a quick evaluation of the feasibility of using heat-pipe cooling. The heat pipes can be embedded in a metallic or composite structure. The maximum heat flux, total integrated heat load, and thermal properties of the structure and heat-pipe container are required input. The heat-pipe operating temperature, maximum surface temperature, heat-pipe length, and heat pipe-spacing can be estimated. Results using the design equations compared well with those from a 3-D finite element analysis for both a large and small radius leading edge.

Glass, David E.↗

NASA NTRS1982

On the dynamic collapse of a column impacting a rigid surface

Results are presented of an analytical investigation on the dynamic collapse of an elastic periodically supported column. The column has an attached mass at one end and impacts a rigid surface with prescribed velocity and angle of incidence at the other end. A first-order approximate nonlinear analysis is developed in which it is assumed that only first-order nonlinear terms need be retained. Differential equations are derived and solved numerically by explicit time integration in order to determine specific ranges of the basic nondimensional parameters for which the response characteristics are markedly different. To assess modeling sophistication relative to prediction accuracy, a comparison is made of the results from this approximate analysis with those of a finite element computer code based on a convected coordinate formulation.

Housner, J. M.↗

NASA NTRS1982

Subsonic aerodynamic and flutter characteristics of several wings calculated by the SOUSSA P1.1 panel method

Several applications of the steady, oscillatory, and unsteady subsonic and supersonic aerodynamics (SOUSSA) computer program to wings with steady and oscillatory motion, including flutter, are discussed. The program employs a generalized Green's function to the full, time-dependent potential-flow equation to obtain an integral equation for the velocity potential at any point in a flow, even points on a body or whole bodies in a flow. Aerodynamic calculations are provided for two rectangular wings, a clipped-tip delta wing, and two swept wings with and without a fuselage. The number and distribution of the finite element panels are varied in order to demonstrate the convergence of the results. The results are shown to be close to those of lifting-surface theory, and further applications with bodies having deformities, arbitrary shapes, motions, and deformations are indicated.

Yates, E. C., Jr.↗

NASA NTRS1989

A new uniformly valid asymptotic integration algorithm for elasto-plastic-creep and unified viscoplastic theories including continuum damage

A new scheme to integrate a system of stiff differential equations for both the elasto-plastic creep and the unified viscoplastic theories is presented. The method has high stability, allows large time increments, and is implicit and iterative. It is suitable for use with continuum damage theories. The scheme was incorporated into MARC, a commercial finite element code through a user subroutine called HYPELA. Results from numerical problems under complex loading histories are presented for both small and large scale analysis. To demonstrate the scheme's accuracy and efficiency, comparisons to a self-adaptive forward Euler method are made.

Chulya, A.↗

NASA NTRS1991

A new uniformly valid asymptotic integration algorithm for elasto-plastic creep and unified viscoplastic theories including continuum damage

Chulya, Abhisak↗

NASA NTRS2009

Thermal Transport Model for Heat Sink Design

A document discusses the development of a finite element model for describing thermal transport through microcalorimeter arrays in order to assist in heat-sinking design. A fabricated multi-absorber transition edge sensor (PoST) was designed in order to reduce device wiring density by a factor of four. The finite element model consists of breaking the microcalorimeter array into separate elements, including the transition edge sensor (TES) and the silicon substrate on which the sensor is deposited. Each element is then broken up into subelements, whose surface area subtends 10 10 microns. The heat capacity per unit temperature, thermal conductance, and thermal diffusivity of each subelement are the model inputs, as are the temperatures of each subelement. Numerical integration using the Finite in Time Centered in Space algorithm of the thermal diffusion equation is then performed in order to obtain a temporal evolution of the subelement temperature. Thermal transport across interfaces is modeled using a thermal boundary resistance obtained using the acoustic mismatch model. The document concludes with a discussion of the PoST fabrication. PoSTs are novel because they enable incident x-ray position sensitivity with good energy resolution and low wiring density.

Chervenak, James A.↗

NASA NTRS1975

On an integro-differential equation model for the study of the response of an acoustically coupled panel

The response of a clamped panel to supersonically convected turbulence is considered. A theoretical model in the form of an integro-differential equation is employed that takes into account the coupling between the panel motion and the surrounding acoustic medium. The kernels of the integrals, which represent induced pressures due to the panel motion, are Green's functions for sound radiations under various moving and stationary sources. An approximate analysis is made by following a finite-element Ritz-Galerkin procedure. Preliminary numerical results, in agreement with experimental findings, indicate that the acoustic damping is the controlling mechanism of the response.

Yen, D. H. Y.↗

NASA NTRS1981

Accuracy and convergence of a finite element algorithm for laminar boundary layer flow

The Galerkin-weighted residuals formulation is employed to derive an implicit finite element solution algorithm for a generally non-linear initial-boundary value problem. Solution accuracy and convergence with discretization refinement are quantized in several error norms, for the non-linear parabolic partial differential equation system governing laminar boundary layer flow, using linear, quadratic and cubic functions. Richardson extrapolation is used to isolate integration truncation error in all norms, and Newton iteration is employed for all equation solutions performed in double-precision. The mathematical theory supporting accuracy and convergence concepts for linear elliptic equations appears extensible to the non-linear equations characteristic of laminar boundary layer flow.

Soliman, M. O.↗

DOE OSTI2026

A Full-Induction Magnetohydrodynamics Solver for Liquid Metal Fusion Blankets in Vertex-CFD

Multiphysics modeling of liquid metal fusion blankets, which produce tritium and convert energy of neutrons created via fusion reactions into heat, is crucial for predicting performance, ensuring structural integrity, and optimizing energy production. While traditional blanket modeling of liquid metal flows during normal steady operating conditions commonly employs the inductionless approximation of the magnetohydrodynamics (MHD) equations, transient scenarios, when the plasma-confining magnetic field varies on millisecond time scales, require a full-induction MHD approach that dynamically evolves the magnetic field via the time-dependent induction equation. This paper presents the formulation, implementation, and initial verification of a full-induction MHD solver integrated within the open-source Vertex-CFD framework, which aims to achieve tight multiphysics coupling, a flexible software design enabling easy extension and addition of physics models, and performance portability across computing platforms. The solver utilizes finite element spatial discretization, implicit Runge–Kutta time integration, and an inexact Newton method to solve the resulting discrete nonlinear system, leveraging Trilinos packages for efficient computation. Verification against selected benchmark problems demonstrates accuracy and robustness of the solver. Furthermore, when the solver is applied to an idealized blanket model in 2.5D and full 3D, results obtained with Vertex-CFD are in good agreement with recently published quasi-2D simulations. These findings establish a computational foundation for future simulations of transient MHD phenomena in liquid metal blankets with Vertex-CFD, and open avenues for future extensions and performance optimizations.

Endeve, Eirik [ORNL] (ORCID:0000000312519507)↗

NASA NTRS1991

Microfracture in high temperature metal matrix laminates

Computational simulation procedures are described to evaluate the composite microfracture behavior, establish the hierarchy/sequence of fracture modes, and the influence of compliant layers and partial debonding on composite properties and microfracture initiation. These procedures are based upon three-dimensional finite element analysis and composite micromechanics equations. Typical results for the effects of compliant layers and partial debonding, microfracture initiation, and propagation and the thermomechanical cyclic loading on a SiC/Ti15 composite system are presented and discussed. The results show that interfacial debonding follows fiber or matrix fracture, and the thermomechanical cyclic loading severely degrades the composite integrity.

Mital, Subodh K.↗