Search NASA⌕ Search

SEARCH · Search NASA

Results for “Mesh adaption”

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 379 records · Page 21

Finite element methods for integrated aerodynamic heating analysis

This report gives a description of the work which has been undertaken during the second year of a three year research program. The objectives of the program are to produce finite element based procedures for the solution of the large scale practical problems which are of interest to the Aerothermal Loads Branch (ALB) at NASA Langley Research Establishment. The problems of interest range from Euler simulations of full three dimensional vehicle configurations to local analyses of three dimensional viscous laminar flow. Adaptive meshes produced for both steady state and transient problems are to be considered. An important feature of the work is the provision of specialized techniques which can be used at ALB for the development of an integrated fluid/thermal/structural modeling capability.

Morgan, K.↗

Unstructured grid methods for compressible flows

The implementation of the finite element method on unstructured triangular grids is described and the development of centered finite element schemes for the solution of the compressible Euler equation on general triangular and tetrahedral grids is discussed. Explicit and implicit Lax-Wendroff type methods and a method based upon the use of explicit multistep timestepping are considered. In the latter case, the convergence behavior of the method is accelerated by the incorporation of a fully unstructured multigrid procedure. The advancing front method for generating unstructured grids of triangles and tetrahedra is described and the application of adaptive mesh techniques to both steady and transient flow analysis is illustrated.

Morgan, K.↗

CFD goes to school

The part that universities should play in the future development of CFD, which must be evaluated in light of CFD's pacing elements and challenges, is discussed. Attention is given to CFD pacing items that must be in place before routine aerodynamic simulation can be performed including grid generation and geometry surface definition, solution adaptive meshing, more efficient time-accurate simulation, modeling of real-gas effects, multiple relative body motion, and prediction of transition and turbulence modeling. As universities have contributed to research in CFD from its inception, this research should continue to enhance and motivate teaching, improve CFD as a discipline, and stimulate faculty and students.

Steger, Joseph L.↗

On computational methods for crashworthiness

The evolution of computational methods for crashworthiness and related fields is described and linked with the decreasing cost of computational resources and with improvements in computation methodologies. The latter includes more effective time integration procedures and more efficient elements. Some recent developments in methodologies and future trends are also summarized. These include multi-time step integration (or subcycling), further improvements in elements, adaptive meshes, and the exploitation of parallel computers.

Belytschko, T.↗

The computation of three-dimensional flows using unstructured grids

A general method is described for automatically discretizing, into unstructured assemblies of tetrahedra, the three-dimensional solution domains of complex shape which are of interest in practical computational aerodynamics. An algorithm for the solution of the compressible Euler equations which can be implemented on such general unstructured tetrahedral grids is described. This is an explicit cell-vertex scheme which follows a general Taylor-Galerkin philosophy. The approach is employed to compute a transonic inviscid flow over a standard wing and the results are shown to compare favorably with experimental observations. As a more practical demonstration, the method is then applied to the analysis of inviscid flow over a complete modern fighter configuration. The effect of using mesh adaptivity is illustrated when the method is applied to the solution of high speed flow in an engine inlet.

Morgan, K.↗

AMR on the CM-2

We describe the development of a structured adaptive mesh algorithm (AMR) for the Connection Machine-2 (CM-2). We develop a data layout scheme that preserves locality even for communication between fine and coarse grids. On 8K of a 32K machine we achieve performance slightly less than 1 CPU of the Cray Y-MP. We apply our algorithm to an inviscid compressible flow problem.

Berger, Marsha J.↗

NASA Langley developments in response calculations needed for failure and life prediction

NASA Langley developments in response calculations needed for failure and life predictions are discussed. Topics covered include: structural failure analysis in concurrent engineering; accuracy of independent regional modeling demonstrated on classical example; functional interface method accurately joins incompatible finite element models; interface method for insertion of local detail modeling extended to curve pressurized fuselage window panel; interface concept for joining structural regions; motivation for coupled 2D-3D analysis; compression panel with discontinuous stiffener coupled 2D-3D model and axial surface strains at the middle of the hat stiffener; use of adaptive refinement with multiple methods; adaptive mesh refinement; and studies on quantity effect of bow-type initial imperfections on reliability of stiffened panels.

Housner, Jerrold M.↗

On the dynamics of a shock-bubble interaction

We present a detailed numerical study of the interaction of a weak shock wave with an isolated cylindrical gas inhomogenity. Such interactions have been studied experimentally in an attempt to elucidate the mechanisms whereby shock waves propagating through random media enhance mixing. Our study concentrates on the early phases of the interaction process which are dominated by repeated refractions of acoustic fronts at the bubble interface. Specifically, we have reproduced two of the experiments performed by Haas and Sturtevant : M(sub s) = 1.22 planar shock wave, moving through air, impinges on a cylindrical bubble which contains either helium or Refrigerant 22. These flows are modelled using the two-dimensional, compressible Euler equations for a two component fluid (air-helium or air-Refrigerant 22). Although simulations of shock wave phenomena are now fairly commonplace, they are mostly restricted to single component flows. Unfortunately, multi-component extensions of successful single component schemes often suffer from spurious oscillations which are generated at material interfaces. Here we avoid such problems by employing a novel, nonconservative shock-capturing scheme. In addition, we have utilized a sophisticated adaptive mesh refinement algorithm which enables extremely high resolution simulations to be performed relatively cheaply. Thus we have been able to reproduce numerically all the intricate mechanisms that were observed experimentally (e.g., transitions from regular to irregular refraction, cusp formation and shock wave focusing, multi-shock and Mach shock structures, jet formation, etc.), and we can now present an updated description for the dynamics of a shock-bubble interaction.

Quirk, James J.↗

On the hydrodynamic interaction of shock waves with interstellar clouds. 1: Nonradiative shocks in small clouds

The interstellar medium (ISM) is inhomogeneous, with clouds of various temperatures and densities embedded in a tenuous intercloud medium. Shocks propagating through the ISM can ablate or destroy the clouds, at the same time significantly altering the properties of the intercloud medium. This paper presents a comprehensive numerical study of the simplest case of the interaction between a shock wave and a spherical cloud, in which the shock far from the cloud is steady and planar, and in which radiative losses, thermal conduction, magnetic fields, and gravitational forces are all neglected. As a result, the problem is completely specified by two numbers: the Mach number of the shock, M, and the ratio of the density of the cloud to that of the intercloud medium, Chi. For strong shocks we show that the dependence on M scales out, so the primary independent parameter is Chi. Variations from this simple case are also considered: the potential effect of radiative losses is assessed by calculations in which the ratio of specific heats in the cloud is 1.1 instead of 5/3; the effect of the initial shape of the cloud is studied by using a cylindrical cloud instead of a spherical one; and the role of the initial shock is determined by considering the case of a cloud embedded in a wind. Local adaptive mesh refinement techniques with a second-order, two-fluid, two-dimensional Godunov hydrodynamic scheme are used to address these problems, allowing heretofore unobtainable numerical resolution. Convergence studies to be described in a subsequent paper demonstrate that about 100 zones per cloud radius are needed for accurate results; previous calculations have generally used about a third of this number. The results of the calculations are analyzed in terms of global quantities which provide an overall description of te shocked cloud: the size and shape of the cloud, the mean density, the mean pressure, the mean velocity, the velocity dispersion, and the total circulation.

Klein, Richard I.↗

An Adaptive Semi-Implicit Scheme for Simulations of Unsteady Viscous Compressible Flows

A numerical scheme for simulation of unsteady, viscous, compressible flows is considered. The scheme employs an explicit discretization of the inviscid terms of the Navier-Stokes equations and an implicit discretization of the viscous terms. The discretization is second order accurate in both space and time. Under appropriate assumptions, the implicit system of equations can be decoupled into two linear systems of reduced rank. These are solved efficiently using a Gauss-Seidel method with multigrid convergence acceleration. When coupled with a solution-adaptive mesh refinement technique, the hybrid explicit-implicit scheme provides an effective methodology for accurate simulations of unsteady viscous flows. The methodology is demonstrated for both body-fitted structured grids and for rectangular (Cartesian) grids.

Steinthorsson, Erlendur↗

Adaptively-refined overlapping grids for the numerical solution of systems of hyperbolic conservation laws

Adaptive mesh refinement (AMR) in conjunction with higher-order upwind finite-difference methods have been used effectively on a variety of problems in two and three dimensions. In this paper we introduce an approach for resolving problems that involve complex geometries in which resolution of boundary geometry is important. The complex geometry is represented by using the method of overlapping grids, while local resolution is obtained by refining each component grid with the AMR algorithm, appropriately generalized for this situation. The CMPGRD algorithm introduced by Chesshire and Henshaw is used to automatically generate the overlapping grid structure for the underlying mesh.

Brislawn, Kristi D.↗

COMET-AR User's Manual: COmputational MEchanics Testbed with Adaptive Refinement

The COMET-AR User's Manual provides a reference manual for the Computational Structural Mechanics Testbed with Adaptive Refinement (COMET-AR), a software system developed jointly by Lockheed Palo Alto Research Laboratory and NASA Langley Research Center under contract NAS1-18444. The COMET-AR system is an extended version of an earlier finite element based structural analysis system called COMET, also developed by Lockheed and NASA. The primary extensions are the adaptive mesh refinement capabilities and a new "object-like" database interface that makes COMET-AR easier to extend further. This User's Manual provides a detailed description of the user interface to COMET-AR from the viewpoint of a structural analyst.

Moas, E.↗

LOOPREF: A Fluid Code for the Simulation of Coronal Loops

This report documents the code LOOPREF. LOOPREF is a semi-one dimensional finite element code that is especially well suited to simulate coronal-loop phenomena. It has a full implementation of adaptive mesh refinement (AMR), which is crucial for this type of simulation. The AMR routines are an improved version of AMR1D. LOOPREF's versatility makes is suitable to simulate a wide variety of problems. In addition to efficiently providing very high resolution in rapidly changing regions of the domain, it is equipped to treat loops of variable cross section, any non-linear form of heat conduction, shocks, gravitational effects, and radiative loss.

deFainchtein, Rosalinda↗

Studies of Tenuous Planetary Atmospheres

The final report includes an overall project overview as well as scientific background summaries of dust and sodium in comets, and tenuous atmospheres of Jupiter's natural satellites. Progress and continuing work related to dust coma and tenuous atmospheric studies are presented. Also included are published articles written during the course of the report period. These are entitled: (1) On Europa's Magnetospheric Interaction: An MHD Simulation; (2) Dust-Gas Interrelations in Comets: Observations and Theory; and (3) Io's Plasma Environment During the Galileo Flyby: Global Three Dimensional MHD Modeling with Adaptive Mesh Refinement.

Combi, Michael R.↗

Development of Numerical Extended Hydrodynamics for Transition-Regime Non-Equilibrium Flows Encountered in Semiconductor Manufacturing Processes

Six months of funding was received for the proposed three year research program (funding for the period from March 1, 1997 to August 31, 1997). Although the official starting date for the project was March 1, 1997, no funding for the project was received until July 1997. In the funded research period, considerable progress was made on Phase I of the proposed research program. The initial research efforts concentrated on applying the 10-, 20-, and 35-moment Gaussian-based closures to a series of standard two-dimensional non-reacting single species test flow problems, such as the flat plate, couette, channel, and rearward facing step flows, and to some other two-dimensional flows having geometries similar to those encountered in chemical-vapor deposition (CVD) reactors. Eigensystem analyses for these systems for the case of two spatial dimensions was carried out and efficient formulations of approximate Riemann solvers have been formulated using these eigenstructures. Formulations to include rotational non-equilibrium effects into the moment closure models for the treatment of polyatomic gases were explored, as the original formulations of the closure models were developed strictly for gases composed of monatomic molecules. The development of a software library and computer code for solving relaxing hyperbolic systems in two spatial dimensions of the type arising from the closure models was also initiated. The software makes use of high-resolution upwind finite-volumes schemes, multi-stage point implicit time stepping, and automatic adaptive mesh refinement (AMR) to solve the governing conservation equations for the moment closures. The initial phase of the code development was completed and a numerical investigation of the solutions of the 10-moment closure model for the simple two-dimensional test cases mentioned above was initiated. Predictions of the 10-moment model were compared to available theoretical solutions and the results of direct-simulation Monte Carlo (DSMC) calculations. The first results of this study were presented at a meeting last year.

Groth, Clinton P. T.↗

Studies for the Loss of Atomic and Molecular Species from IO

The general objective of this project has been to advance our theoretical understanding of lo's atmosphere and how various atomic and molecular species are lost from this atmosphere and are distributed in the circumplanetary environment of Jupiter. This grant has provided support for the activities of Dr. Michael Combi at the University of Michigan to serve as a small part in collaboration with a larger project awarded to Atmospheric & Environmental Research, Inc., with primary principal investigator Dr. William H. Smyth. Dr. Combi is the Principal Investigator and Project Manager for the Michigan grant NAG5-6187. This Michigan grant has provided for a continuation of a collaboration between Drs. Smyth and Combi in related efforts beginning in 1981, and with the object to develop and apply sophisticated theoretical models to interpret and to relate a number of new and exciting observations for the atmospheric gases of the satellite. The ability to interpret and then to relate through the theoretical fabric a number of these otherwise independent observations are a central strength of this program. This comprehensive approach provides a collective power, extracting more from the sum of the parts and seeing beyond various limitations that are inherent in any one observation. Although the approach is designed to unify, the program is divided into well-defined studies for the likely dominant atmospheric gases involving species of the SO2 family (SO2, SO, O2, S and O) and for the trace atmospheric gas atomic sodium and a likely escaping molecular ion NaX(+) (where Na(X) is the atmospheric molecule and X represents one or more atoms).Attachments: IO's sodium corona and spatially cloud: a consistent flux speed distribution. and Io's plasma environment during the Galileo flyby: global three-dimensional MHD modeling with adaptive mesh refinement.

Combi, Michael R.↗

Stress Recovery and Error Estimation for Shell Structures

The Penalized Discrete Least-Squares (PDLS) stress recovery (smoothing) technique developed for two dimensional linear elliptic problems is adapted here to three-dimensional shell structures. The surfaces are restricted to those which have a 2-D parametric representation, or which can be built-up of such surfaces. The proposed strategy involves mapping the finite element results to the 2-D parametric space which describes the geometry, and smoothing is carried out in the parametric space using the PDLS-based Smoothing Element Analysis (SEA). Numerical results for two well-known shell problems are presented to illustrate the performance of SEA/PDLS for these problems. The recovered stresses are used in the Zienkiewicz-Zhu a posteriori error estimator. The estimated errors are used to demonstrate the performance of SEA-recovered stresses in automated adaptive mesh refinement of shell structures. The numerical results are encouraging. Further testing involving more complex, practical structures is necessary.

Yazdani, A. A.↗

Crack Turning and Arrest Mechanisms for Integral Structure

In the course of several years of research efforts to predict crack turning and flapping in aircraft fuselage structures and other problems related to crack turning, the 2nd order maximum tangential stress theory has been identified as the theory most capable of predicting the observed test results. This theory requires knowledge of a material specific characteristic length, and also a computation of the stress intensity factors and the T-stress, or second order term in the asymptotic stress field in the vicinity of the crack tip. A characteristic length, r(sub c), is proposed for ductile materials pertaining to the onset of plastic instability, as opposed to the void spacing theories espoused by previous investigators. For the plane stress case, an approximate estimate of r(sub c), is obtained from the asymptotic field for strain hardening materials given by Hutchinson, Rice and Rosengren (HRR). A previous study using of high order finite element methods to calculate T-stresses by contour integrals resulted in extremely high accuracy values obtained for selected test specimen geometries, and a theoretical error estimation parameter was defined. In the present study, it is shown that a large portion of the error in finite element computations of both K and T are systematic, and can be corrected after the initial solution if the finite element implementation utilizes a similar crack tip discretization scheme for all problems. This scheme is applied for two-dimensional problems to a both a p-version finite element code, showing that sufficiently accurate values of both K(sub I) and T can be obtained with fairly low order elements if correction is used. T-stress correction coefficients are also developed for the singular crack tip rosette utilized in the adaptive mesh finite element code FRANC2D, and shown to reduce the error in the computed T-stress significantly. Stress intensity factor correction was not attempted for FRANC2D because it employs a highly accurate quarter-point scheme to obtain stress intensity factors.

Pettit, Richard↗