Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mesh Adaptivity”

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 1,225 records · Page 68

Solution verification, goal-oriented adaptive methods for stochastic advection–diffusion problems

A goal-oriented analysis of linear, stochastic advection–diffusion models is presented which provides both a method for solution verification as well as a basis for improving results through adaptation of both the mesh and the way random variables are approximated. A class of model problems with random coefficients and source terms is cast in a variational setting. Specific quantities of interest are specified which are also random variables. A stochastic adjoint problem associated with the quantities of interest is formulated and a posteriori error estimates are derived. These are used to guide an adaptive algorithm which adjusts the sparse probabilistic grid so as to control the approximation error. Numerical examples are given to demonstrate the methodology for a specific model problem.

97 MATHEMATICS AND COMPUTING↗

An assessment of the adaptive unstructured tetrahedral grid, Euler Flow Solver Code FELISA

A three-dimensional solution-adaptive Euler flow solver for unstructured tetrahedral meshes is assessed, and the accuracy and efficiency of the method for predicting sonic boom pressure signatures about simple generic models are demonstrated. Comparison of computational and wind tunnel data and enhancement of numerical solutions by means of grid adaptivity are discussed. The mesh generation is based on the advancing front technique. The FELISA code consists of two solvers, the Taylor-Galerkin and the Runge-Kutta-Galerkin schemes, both of which are spacially discretized by the usual Galerkin weighted residual finite-element methods but with different explicit time-marching schemes to steady state. The solution-adaptive grid procedure is based on either remeshing or mesh refinement techniques. An alternative geometry adaptive procedure is also incorporated.

Djomehri, M. Jahed↗

Comparison of multimesh hp-FEM to interpolation and projection methods for spatial coupling of thermal and neutron diffusion calculations

Multiphysics solution challenges are legion within the field of nuclear reactor design and analysis. One major issue concerns the coupling between heat and neutron flow (neutronics) within the reactor assembly. These phenomena are usually very tightly interdependent, as large amounts of heat are quickly produced with an increase in fission events within the fuel, which raises the temperature that affects the neutron cross section of the fuel. Furthermore, there typically is a large diversity of time and spatial scales between mathematical models of heat and neutronics. Indeed, the different spatial resolution requirements often lead to the use of very different meshes for the two phenomena. As the equations are coupled, one must take care in exchanging solution data between them, or significant error can be introduced into the coupled problem. We propose a novel approach to the discretization of the coupled problem on different meshes based on an adaptive multimesh higher-order finite element method (hp-FEM), and compare it to popular interpolation and projection methods. We show that the multimesh hp-FEM method is significantly more accurate than the interpolation and projection approaches considered in this study.

97 MATHEMATICS AND COMPUTING↗

Emergency core cooling system performance criteria for Multi-Layered Silicon Carbide nuclear fuel cladding

Multi-layered Silicon Carbide (SiC) cladding has been proposed to replace the current Zircaloy (Zr) cladding in existing nuclear power plants. Due to its ceramic property, SiC may have a different failure mechanism from that of metallic Zr clad. This characteristic should be considered in designing the plant’s safety system performance, such that the fuel’s integrity is reasonably assured during accident scenarios. This paper aims at optimizing the Emergency Core Cooling System (ECCS) performance criteria in mitigating LBLOCA accidents at a Pressurized Water Reactor (PWR) using the uranium dioxide fuel and duplex or triplex SiC cladding structures. Clad failures by brittle fractures and thermal degradation were modeled in a fuel analysis code. This fuel model was coupled with the RELAP5 plant model by adapting the fuel’s mesh configuration and formulating the effective heat transfer coefficients. A model to quantify the likelihood of clad functional failure as a combination of individual layer fractures at various axial positions was proposed. Updated SiC material properties were gathered from recent reference studies. Clad reliability estimates using the proposed model were compared with the estimates obtained using other models. The results were used to relax the criterion of required actuation timing of ECCS, such that the risk of the new plant design is less or equal to that of the initial plant configuration. Triplex clad was more reliable than duplex clad and Zr, reducing the conditional core damage probability by three orders of magnitude in the reference LBLOCA scenario. Compressive stresses due to irradiation swelling played the determining role in the vulnerabilities of each clad design. The duplex clad was vulnerable to pressure-loading stresses while triplex clad was prone to fail due to thermal shocks experienced during the LOCA reflood phase. Additionally, it was concluded that the most vulnerable axial fuel rod location in the reference LBLOCA scenario was located at the highest axial power peaking location. A sensitivity study on the ECCS timing density was conducted. The results showed that ECCS actuation timing could be extended from the initial technical specification, without increasing the fuel damage likelihood. The time extension was distributed in the order of several hundreds of seconds, which was comparable with time extensions offered by FeCrAl and Cr-coated-Zr clad designs in LBLOCA scenarios. This time extension may potentially reduce the failure probabilities of EDGs and pumps, hence improving the plant safety further.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Effects of isothermal stratification strength on vorticity dynamics for single-mode compressible Rayleigh-Taylor instability

The effects of isothermal initial stratification on the dynamics of the vorticity for single-mode Rayleigh-Taylor instability (RTI) are examined using two-dimensional fully compressible wavelet-based direct numerical simulations. Here, the simulations model low Atwood number ( A = 0.04) RTI development for four different initial stratification strengths, corresponding to Mach numbers from 0.3 (weakly stratified) to 1.2 (strongly stratified), and for three different Reynolds numbers, from 25 500 to 102 000. Here, the Mach number is based on the Atwood-independent gravity wave speed and characterizes the strength of the initial stratification. All simulations use adaptive wavelet-based mesh refinement to achieve very fine spatial resolutions at relatively low computational cost. For all stratifications, the RTI bubble and spike go through the exponential growth regime, followed by a slowing of the RTI evolution.

42 ENGINEERING↗

PVmesh [SWR-24-123]

PVmesh is an open-source Python tool for generating high-fidelity, adaptive finite-element meshes for photovoltaic modules. It supports multilayer panel geometries, frame details, and mounting zones that are often difficult to mesh robustly with generic workflows. Built on Gmsh, PVmesh automates geometry construction, partitioning, domain/boundary tagging, and export for downstream finite-element solvers.

He, Xin [National Laboratory of the Rockies (NLR),↗

Stability and convergence of nuclear detonations in white dwarf collisions

We investigate the numerical stability of thermonuclear detonations in 1D accelerated reactive shocks and 2D binary collisions of equal-mass magnetized and unmagnetized white dwarf stars. To achieve high resolution at initiation sites, we devised geometric gridding and mesh velocity strategies specially adapted to the unique requirements of head-on collisional geometries, scenarios in which one expects maximum production of iron-group products. We study the effects of grid resolution and the limiting of temperature, energy, and reactants for different stellar masses, separations, magnetic fields, inigenerationtial compositions, detonation mechanisms, and limiter parameters across a range of cell sizes from 1 to 100 km. Our results set bounds on the parameter space of limiter amplitudes for which both temperature- and energy-limiting procedures yield consistent and monotonically convergent solutions. Within these bounds, we find that grid resolutions of 5 km or better are necessary for uncertainties in total released energy and iron-group products to drop below 10%. Intermediate-mass products (e.g., calcium) exhibit similar convergence trends but with somewhat greater uncertainty. These conclusions apply equally to pure C/O white dwarfs, multispecies compositions (including helium shells), magnetized and unmagnetized cores, and either single or multiple detonation scenarios.

79 ASTRONOMY AND ASTROPHYSICS↗

A Navier-Stokes solver for high speed equilibrium flows and application to blunt bodies

This paper presents a finite element method for the solution of Navier-Stokes equations with the assumption of thermodynamic and chemical equilibrium. The method employs an upwind finite element technique with an implicit time-marching scheme for the solution, and uses an adaptively generated unstructured triangular mesh with several layers of quadrilateral elements near solid walls. The complexity associated with the assumption that the flow is in equilibrium is treated consistently, and the inviscid flux Jacobian matrices are derived. Several problems involving inviscid and viscous hypersonic flow past blunt are solved. Results are compared with other numerical results and experimental data, and generally good agreement is observed.

Prabhu, Ramadas K.↗

Efficient Fluid Dynamic Design Optimization Using Cartesian Grids

This report is subdivided in three parts. The first one reviews a new approach to the computation of inviscid flows using Cartesian grid methods. The crux of the method is the curvature-corrected symmetry technique (CCST) developed by the present authors for body-fitted grids. The method introduces ghost cells near the boundaries whose values are developed from an assumed flow-field model in vicinity of the wall consisting of a vortex flow, which satisfies the normal momentum equation and the non-penetration condition. The CCST boundary condition was shown to be substantially more accurate than traditional boundary condition approaches. This improved boundary condition is adapted to a Cartesian mesh formulation, which we call the Ghost Body-Cell Method (GBCM). In this approach, all cell centers exterior to the body are computed with fluxes at the four surrounding cell edges. There is no need for special treatment corresponding to cut cells which complicate other Cartesian mesh methods.

Dadone, A.↗

Aerothermodynamic CFD Analysis

This presentation provides an introduction to two Computational Fluid Dynamics (CFD) codes, LAURA and FUN3D, which have been developed and widely used at NASA Langley Research Center in Entry, Descent, and Landing applications. A workflow using the LAURA and FUN3D CFD codes to predict aerodynamic and aerothermodynamic engineering quantities is presented. The presentation will cover the basics of setting up and running simulations using both CFD codes, as well as how to post-process results obtained from each code. More advanced topics will also be presented, including a state-of-the-art uncertainty quantification approach to be included in the next LAURA release, and a new geometry-based workflow used with FUN3D focused on replacing user-defined meshing with an adaptation-based approach.

Kyle B Thompson↗

Comparison of multimesh hp-FEM to interpolation and projection methods for spatial coupling of reactor thermal and neutron diffusion calculations

Multiphysics solution challenges are legion within the ?eld of nuclear reactor design and analysis. One major issue concerns the coupling between heat and neutron ?ow (neutronics) within the reactor assembly. These phenomena are usually very tightly interdependent, as large amounts of heat are quickly produced with an increase in ?ssion events within the fuel, which raises the temperature that a?ects the neutron cross section of the fuel. Furthermore, there typically is a large diversity of time and spatial scales between mathematical models of heat and neutronics. Indeed, the di?erent spatial resolution requirements often lead to the use of very di?erent meshes for the two phenomena. As the equations are coupled, one must take care in exchanging solution data between them, or signi?cant error can be introduced into the coupled problem. We propose a novel approach to the discretization of the coupled problem on di?erent meshes based on an adaptive multimesh higher-order ?nite element method (hp-FEM), and compare it to popular interpolation and projection methods. We show that the multimesh hp-FEM method is signi?cantly more accurate than the interpolation and projection approaches considered in this study.

97 MATHEMATICS AND COMPUTING↗

Arbitrary-level hanging nodes for adaptive hphp-FEM approximations in 3D

In this paper we discuss constrained approximation with arbitrary-level hanging nodes in adaptive higher-order finite element methods (hphp-FEM) for three-dimensional problems. This technique enables using highly irregular meshes, and it greatly simplifies the design of adaptive algorithms as it prevents refinements from propagating recursively through the finite element mesh. The technique makes it possible to design efficient adaptive algorithms for purely hexahedral meshes. We present a detailed mathematical description of the method and illustrate it with numerical examples.

99 GENERAL AND MISCELLANEOUS↗

An efficient adaptive patched grid gas dynamic solver for complex flows

The patched mesh CSCM flow solver is used with a fast adaptive grid technique to compute complex two dimensional and axisymmetric, inviscid and viscous flows accurately and efficiently. The patched grid technique allows the user to generate the initial grid around complex regions with minimal effort. The adaptive grid strategy permits the grid points to adjust to the solution in accordance with the flow structure and user specified grid constraints and thus improve the accuracy of the solution. Solutions to challenging flow problems are obtained with adapted grids and the numerical examples and the results presented here show clearly the effectiveness of the adaptive grid solution strategy to obtain highly accurate numerical solutions with minimal human effort and computational resources.

Venkatapathy, E.↗

Adaptive Control Parameters for Dispersal of Multi-Agent Mobile Ad Hoc Network (MANET) Swarms

A mobile ad hoc network is a collection of independent nodes that communicate wirelessly with one another. This paper investigates nodes that are swarm robots with communications and sensing capabilities. Each robot in the swarm may operate in a distributed and decentralized manner to achieve some goal. This paper presents a novel approach to dynamically adapting control parameters to achieve mesh configuration stability. The presented approach to robot interaction is based on spring force laws (attraction and repulsion laws) to create near-optimal mesh like configurations. In prior work, we presented the extended virtual spring mesh (EVSM) algorithm for the dispersion of robot swarms. This paper extends the EVSM framework by providing the first known study on the effects of adaptive versus static control parameters on robot swarm stability. The EVSM algorithm provides the following novelties: 1) improved performance with adaptive control parameters and 2) accelerated convergence with high formation effectiveness. Simulation results show that 120 robots reach convergence using adaptive control parameters more than twice as fast as with static control parameters in a multiple obstacle environment.

42 ENGINEERING↗

Algebraic grid generation with boundary orthogonality control

This paper presents a new approach in applying the blending function method to the generation of structured finite-difference meshes. It is shown that the surface equation developed by Coons (1967) can be easily adapted to generate two-dimensional meshes, and its extension to hyperspace allows the construction of three-dimensional meshes. The robust method in its basic form blends among the boundary curves/surfaces conforming to their intrinsic slopes. When parameter-controlled boundary-slope corrections are introduced, it is possible to establish varying degrees of orthogonality along the boundaries. Additional controls are made possible by moving points in the parameter space, which affects the movement of points in the physical space. The combination of these controls provides a simple and yet powerful means for exploiting the blending-function method to satisfy the needs of practical CFD applications. Results from a validation program for two-dimensions are shown for a variety of geometries.

Luh, Raymond Ching-Chung↗

A Parallel Ghosting Algorithm for The Flexible Distributed Mesh Database

Critical to the scalability of parallel adaptive simulations are parallel control functions including load balancing, reduced inter-process communication and optimal data decomposition. In distributed meshes, many mesh-based applications frequently access neighborhood information for computational purposes which must be transmitted efficiently to avoid parallel performance degradation when the neighbors are on different processors. This article presents a parallel algorithm of creating and deleting data copies, referred to as ghost copies, which localize neighborhood data for computation purposes while minimizing inter-process communication. The key characteristics of the algorithm are: (1) It can create ghost copies of any permissible topological order in a 1D, 2D or 3D mesh based on selected adjacencies. (2) It exploits neighborhood communication patterns during the ghost creation process thus eliminating all-to-all communication. (3) For applications that need neighbors of neighbors, the algorithm can create n number of ghost layers up to a point where the whole partitioned mesh can be ghosted. Strong and weak scaling results are presented for the IBM BG/P and Cray XE6 architectures up to a core count of 32,768 processors. The algorithm also leads to scalable results when used in a parallel super-convergent patch recovery error estimator, an application that frequently accesses neighborhood data to carry out computation.

Mubarak, Misbah↗

Mesh refinement in finite element analysis by minimization of the stiffness matrix trace

Most finite element packages provide means to generate meshes automatically. However, the user is usually confronted with the problem of not knowing whether the mesh generated is appropriate for the problem at hand. Since the accuracy of the finite element results is mesh dependent, mesh selection forms a very important step in the analysis. Indeed, in accurate analyses, meshes need to be refined or rezoned until the solution converges to a value so that the error is below a predetermined tolerance. A-posteriori methods use error indicators, developed by using the theory of interpolation and approximation theory, for mesh refinements. Some use other criterions, such as strain energy density variation and stress contours for example, to obtain near optimal meshes. Although these methods are adaptive, they are expensive. Alternatively, a priori methods, until now available, use geometrical parameters, for example, element aspect ratio. Therefore, they are not adaptive by nature. An adaptive a-priori method is developed. The criterion is that the minimization of the trace of the stiffness matrix with respect to the nodal coordinates, leads to a minimization of the potential energy, and as a consequence provide a good starting mesh. In a few examples the method is shown to provide the optimal mesh. The method is also shown to be relatively simple and amenable to development of computer algorithms. When the procedure is used in conjunction with a-posteriori methods of grid refinement, it is shown that fewer refinement iterations and fewer degrees of freedom are required for convergence as opposed to when the procedure is not used. The mesh obtained is shown to have uniform distribution of stiffness among the nodes and elements which, as a consequence, leads to uniform error distribution. Thus the mesh obtained meets the optimality criterion of uniform error distribution.

Kittur, Madan G.↗