SEARCH · Search NASA

Results for “non-matching grids”

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.

6 records

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Coupling Approaches with Non-matching Grids for Classical Linear Elasticity and Bond-based Peridynamic Models in 1D

Local-nonlocal coupling approaches provide a means to combine the computational efficiency of local models and the accuracy of nonlocal models. To facilitate the coupling of the two models, non-matching grids are often desirable as nonlocal grids usually require a finer resolution than local grids. In that case, it is often convenient to resort to interpolation operators so that models can exchange information in the overlap regions when nodes from the two grids do not coincide. This paper studies three existing coupling approaches, namely 1) a method that enforces matching displacements in an overlap region, 2) a variant that enforces a constraint on the stresses instead, and 3) a method that considers a variable horizon in the vicinity of the interfaces. Further, the effect of the interpolation order and of the grid ratio on the performance of the three coupling methods with non-matching grids is carefully studied on one-dimensional examples using polynomial manufactured solutions. The numerical results show that the degree of the interpolants should be chosen with care to avoid introducing additional modeling errors, or simply minimize these errors, in the coupling approach.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1987

A multidomain global modeling technique for analysis of Space Shuttle main engine

A multidomain global modeling technique is presented. The purpose of this multidomain methodology is to simultaneously analyze the three-dimensional compressible flow field in Space Shuttle main engine components. This technique provides flexibility and modularity for individual domain and also simplifies the grid generation task for each domain. A two-way coupling procedure is adopted for exchanging interface boundary conditions. This algorithm which requires that the adjacent domains share an overlapping grid slab, is fully conservative and implicit. For non-matching grid systems and for non-overlapping interfaces, an interpolation and/or extrapolation procedure is used when transferring the interface boundary conditions at the adjacent domains. To maintain overall conservation, a flux correction method is subsequently used. The coupling procedure was verified by modeling several geometries with features similar to SSME components. Satisfactory results for all test cases were obtained. Finally, this algorithm has been applied to the SSME flow simulation.

Tam, L. T.↗

DOE OSTI2021

Continuum gyrokinetic simulations of edge plasmas in single-null geometries

The first continuum gyrokinetic calculations of electrostatic ion scale turbulence are presented for the case of a diverted tokamak geometry. The simulation model solves the long-wavelength limit of the full-F gyrokinetic equation for ion species coupled to the quasi-neutrality equation for electrostatic potential variations, where a fluid model is used for an electron response. Furthermore, in order to facilitate simulations of highly-anisotropic microturbulence in the presence of strong magnetic shear and a magnetic X-point, a numerical algorithm utilizing a locally field-aligned multiblock coordinate system has been developed and implemented in the 5D finite-volume code COGENT. In this approach, the toroidal direction is divided into blocks, such that within each block, the cells are field-aligned and a non-matching grid interface is allowed at block boundaries. The toroidal angle corresponds to the “coarse” field-aligned coordinate, whereas the poloidal cross section, comprised of the radial and poloidal directions, is finely gridded to resolve short-scale perpendicular turbulence structures and to support accurate re-mapping (interpolation) at block boundaries. The 5D simulations explore cross-separatrix ion scale turbulence in the presence of a self-consistent radial electric field and address the effects of magnetic-shear stabilization in the X-point region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2019

Domain decomposition and multiscale mortar mixed finite element methods for linear elasticity with weak stress symmetry

Two non-overlapping domain decomposition methods are presented for the mixed finite element formulation of linear elasticity with weakly enforced stress symmetry. The methods utilize either displacement or normal stress Lagrange multiplier to impose interface continuity of normal stress or displacement, respectively. By eliminating the interior subdomain variables, the global problem is reduced to an interface problem, which is then solved by an iterative procedure. The condition number of the resulting algebraic interface problem is analyzed for both methods. A multiscale mortar mixed finite element method for the problem of interest on non-matching multiblock grids is also studied. It uses a coarse scale mortar finite element space on the non-matching interfaces to approximate the trace of the displacement and impose weakly the continuity of normal stress. A priori error analysis is performed. It is shown that, with appropriate choice of the mortar space, optimal convergence on the fine scale is obtained for the stress, displacement, and rotation, as well as some superconvergence for the displacement. Computational results are presented in confirmation of the theory of all proposed methods.

Mathematics↗

DOE OSTI2017

Asynchronous space–time domain decomposition method with localized uncertainty quantification

The computational cost associated with uncertainty quantification of engineering problems featuring localized phenomenon can be reduced by confining the random variability of the model parameters within a region of interest. Here in this case, a localized treatment of mesh and time resolutions is required to capture the effect of the confined material uncertainty on the global response. We present a computational approach for localized uncertainty quantification with the capability of asynchronous treatment of mesh and time resolutions. In particular, we allow each subdomain to have its local uncertainty representation and the corresponding mesh and time resolutions. As a result, computing resources can be directed toward a small region of interest where a model with high spatial and temporal resolutions is required. To verify the numerical implementation, we consider elastic wave propagation in an axially loaded beam. Moreover, we perform convergence studies with respect to the spatial and temporal discretizations as well as the size of an uncertain subdomain. A projectile impacting a composite sandwich plate is considered as an engineering application for the proposed method.

42 ENGINEERING↗

DOE OSTI2020

Progress with the 5D full-F continuum gyrokinetic code COGENT

COGENT is an Eulerian gyrokinetic code being developed for edge plasma modelling. The code is distinguished by the use of a high-order finite-volume (conservative) discretization combined with mapped multi-block grid technology. Our recent work is focused on the development of a 5D full-F COGENT version.Anumerical algorithm utilizing locally a field-aligned multi-block coordinate system is implemented to facilitate simulations of highly anisotropic microturbulence in the presence of a strong magnetic shear. In this approach, the toroidal direction is divided into blocks such that,within each block, the cells are field-aligned and a non-matching (non-conformal) grid interface is allowed at the block boundaries. In this paper we report on details of the numerical implementation and present preliminary results of verification studies performed for the case of the ion temperature gradient (ITG) instability in a sheared toroidal annulus geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗