SEARCH · Engineering Papers

Results for “mesh compression”

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

All sources NASA NTRS DOE OSTI

DOE OSTI2015

Numerical insight into the seismic behavior of eight masonry towers in Northern Italy: FE pushover vs non-linear dynamic analyses

This study presents some FE results regarding the behavior under horizontal loads of eight existing masonry towers located in the North-East of Italy. The towers, albeit unique for geometric and architectural features, show some affinities which justify a comparative analysis, as for instance the location and the similar masonry material. Their structural behavior under horizontal loads is therefore influenced by geometrical issues, such as slenderness, walls thickness, perforations, irregularities, presence of internal vaults, etc., all features which may be responsible for a peculiar output. The geometry of the towers is deduced from both existing available documentation and in-situ surveys. On the basis of such geometrical data, a detailed 3D realistic mesh is conceived, with a point by point characterization of each single geometric element. The FE models are analysed under seismic loads acting along geometric axes of the plan section, both under non-linear static (pushover) and non-linear dynamic excitation assumptions. A damage-plasticity material model exhibiting softening in both tension and compression, already available in the commercial code Abaqus, is used for masonry. Pushover analyses are performed with both G1 and G2 horizontal loads distribution, according to Italian code requirements, along X+/− and Y+/− directions. Non-linear dynamic analyses are performed along both X and Y directions with a real accelerogram scaled to different peak ground accelerations. Some few results are presented in this paper. It is found that the results obtained with pushover analyses reasonably well fit expensive non-linear dynamic simulations, with a slightly less conservative trend.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1995

Bending strength model for internal spur gear teeth

Internal spur gear teeth are normally stronger than pinion teeth of the same pitch and face width since external teeth are smaller at the base. However, ring gears which are narrower have an unequal addendum or are made of a material with a lower strength than that of the meshing pinion may be loaded more critically in bending. In this study, a model for the bending strength of an internal gear tooth as a function of the applied load pressure angle is presented which is based on the inscribed Lewis constant strength parabolic beam. The bending model includes a stress concentration factor and an axial compression term which are extensions of the model for an external gear tooth. The geometry of the Lewis factor determination is presented, the iteration to determine the factor is described, and the bending strength J factor is compared to that of an external gear tooth. This strength model will assist optimal design efforts for unequal addendum gears and gears of mixed materials.

Savage, Michael↗

DOE OSTI2022

Fast Solution of Fully Implicit Runge--Kutta and Discontinuous Galerkin in Time for Numerical PDEs, Part II: Nonlinearities and DAEs

Fully implicit Runge--Kutta (IRK) methods have many desirable accuracy and stability properties as time integration schemes, but high-order IRK methods are not commonly used in practice with large-scale numerical PDEs because of the difficulty of solving the stage equations. This paper introduces a theoretical and algorithmic framework for solving the nonlinear equations that arise from IRK methods (and discontinuous Galerkin discretizations in time) applied to nonlinear numerical PDEs, including PDEs with algebraic constraints. Several new linearizations of the nonlinear IRK equations are developed, offering faster and more robust convergence than the often-considered simplified Newton, as well as an effective preconditioner for the true Jacobian if exact Newton iterations are desired. Inverting these linearizations requires solving a set of block 2 x 2 systems. Under quite general assumptions, it is proven that the preconditioned 2 x 2 operator's condition number is bounded by a small constant close to one, independent of the spatial discretization, spatial mesh, and time step, and with only weak dependence on the number of stages or integration accuracy. Moreover, the new method is built using the same preconditioners needed for backward Euler-type time stepping schemes, so can be readily added to existing codes. The new methods are applied to several challenging fluid flow problems, including the compressible Euler and Navier--Stokes equations, and the vorticity-streamfunction formulation of the incompressible Euler and Navier--Stokes equations. Up to 10th-order accuracy is demonstrated using Gauss IRK, while in all cases fourth-order Gauss IRK requires roughly half the number of preconditioner applications as required by standard Singly diagonally implicit Runge--Kutta methods.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2019

ECP Milestone Memo WBS 2.3.4.13 ECP/VTK-m FY19Q1 [MS-19/01-03] ZFP / Release / Clip STDA05-17

The STDA05-17 milestone comprises the following 3 deliverables. VTK-m Release 2 We will provide a release of VTK-m software and associated documentation. The source code repository will be tagged at a stable state, and, at a minimum, tarball captures of the source code will be made available from the web site. A version of the VTK-m User's Guide documenting this release will also be made available. Productionize zfp compression The "ZFP: Compressed Floating-Point Arrays" project (WBS 1.3.4.13) is creating an implementation of ZFP compression in VTK-m. Their implementation will be focused on operating in CUDA. The VTK-m project will assist by generalizing the implementation to other devices (such as multi-core CPUs). We will also assist in productionizing the code such that it can be used by external projects and products. Clip Clip operations intersect meshes with implicit functions. It is the foundation of spatial subsetting algorithms, such as "box," and the foundation of data-based subsetting, such as "isovolume." The algorithm requires considering thousands of possible cases, and is thus quite difficult to implement. This milestone will implement clipping to be sufficient for Visit's and ParaView's needs.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1990

Investigation of advanced counterrotation blade configuration concepts for high speed turboprop systems, task 1: Ducted propfan analysis

The time-dependent three-dimensional Euler equations of gas dynamics were solved numerically to study the steady compressible transonic flow about ducted propfan propulsion systems. Aerodynamic calculations were based on a four-stage Runge-Kutta time-marching finite volume solution technique with added numerical dissipation. An implicit residual smoothing operator was used to aid convergence. Two calculation grids were employed in this study. The first grid utilized an H-type mesh network with a branch cut opening to represent the axisymmetric cowl. The second grid utilized a multiple-block mesh system with a C-type grid about the cowl. The individual blocks were numerically coupled in the Euler solver. Grid systems were generated by a combined algebraic/elliptic algortihm developed specifically for ducted propfans. Numerical calculations were initially performed for unducted propfans to verify the accuracy of the three-dimensional Euler formulation. The Euler analyses were then applied for the calculation of ducted propfan flows, and predicted results were compared with experimental data for two cases. The three-dimensional Euler analyses displayed exceptional accuracy, although certain parameters were observed to be very sensitive to geometric deflections. Both solution schemes were found to be very robust and demonstrated nearly equal efficiency and accuracy, although it was observed that the multi-block C-grid formulation provided somewhat better resolution of the cowl leading edge region.

Hall, Edward J.↗

DOE OSTI2023

Data-Driven Smoothers for Extreme-Scale Computing

Patch-based relaxation refers to a family of methods for solving linear systems which partitions the matrix into smaller pieces often corresponding to groups of adjacent degrees of freedom residing within patches of the computational domain. The two most common families of patch-based methods are block-Jacobi and Schwarz methods, where the former typically corresponds to non-overlapping domains and the later implies some overlap. We focus on cases where each patch consists of the degrees of freedom on a finite element method mesh cell. Patch methods often capture complex local physics much more effectively than simpler point-smoothers such as Jacobi; however, forming, inverting, and applying each patch can be prohibitively expensive in terms of both storage and computation time. To this end, we propose several approaches for performing analysis on these patches and constructing a reduced representation. The compression techniques rely on either matrix norm comparisons or unsupervised learning via a clustering approach. We illustrate how it is frequently possible to retain/factor less than 5% of all patches and still develop a method that converges only a little slower than when all patches are stored/factored.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1996

A Genuinely Two-Dimensional Scheme for the Compressible Euler Equations

We present a new genuinely multidimensional discretization for the compressible Euler equations. It is the only high-resolution scheme known to us where Gauss-Seidel relaxation is stable when applied as a smoother directly to the resulting high-resolution scheme. This allows us to construct a very simple and highly efficient multigrid steady-state solver. The scheme is formulated on triangular (possibly unstructured) meshes.

Sidilkover, David↗

NASA NTRS1980

Forebody and base region real-gas flow in severe planetary entry by a factored implicit numerical method. I - Computational fluid dynamics

A new code for the simulation of full (forebody and base region) flowfields about bluff bodies in the hypersonic regime of severe planetary entry is described. The present 'maximally conservative, maximally differenced' formulation of the unsteady compressible Navier-Stokes equations for 2-D axisymmetric 3-D flow is contrasted for stability with previous formulations of Viviand, Kutler, et al, and Thomas and Lombard. Discrete metric relations peculiar to the axisymmetric finite volume formulation are presented along with a general discussion of their relations to and consequences of failure to close computational cells. A computational mesh of curvilinear coordinate topology singular in the flow regime is presented that permits aligned capturing of the major physical features of the complex flowfield.

Lombard, C. K.↗

NASA NTRS1985

A fourth-order scheme for the unsteady compressible Navier-Stokes equations

A computational scheme is described which is second-order accurate in time and fourth-order accurate in space (2-4). This method is applied to study the stability of compressible boundary layers. The laminar compressible Navier-Stokes equations are solved with a time harmonic inflow superimposed on the steady state solution. This results in spatially unstable modes. It is shown that the second-order methods are inefficient for calculating the growth rates and phases of the unstable modes. In contrast the fourth-order method yields accurate results on relatively coarse meshes.

Bayliss, A.↗

NASA NTRS1985

A fourth-order scheme for the unsteady compressible Navier-Stokes equations

Bayliss, A.↗

NASA NTRS2008

Compressible Boundary Layer Predictions at High Reynolds Number using Hybrid LES/RANS Methods

Simulations of compressible boundary layer flow at three different Reynolds numbers (Re(sub delta) = 5.59x10(exp 4), 1.78x10(exp 5), and 1.58x10(exp 6) are performed using a hybrid large-eddy/Reynolds-averaged Navier-Stokes method. Variations in the recycling/rescaling method, the higher-order extension, the choice of primitive variables, the RANS/LES transition parameters, and the mesh resolution are considered in order to assess the model. The results indicate that the present model can provide good predictions of the mean flow properties and second-moment statistics of the boundary layers considered. Normalized Reynolds stresses in the outer layer are found to be independent of Reynolds number, similar to incompressible turbulent boundary layers.

Choi, Jung-Il↗

NASA NTRS1984

The solution of non-linear hyperbolic equation systems by the finite element method

A finite-element method for the solution of nonlinear hyperbolic systems of equations, such as those encountered in non-self-adjoint problems of transient phenomena in convection-diffusion or in the mixed representation of wave problems, is developed and demonstrated. The problem is rewritten in moving coordinates and reinterpolated to the original mesh by a Taylor expansion prior to a standard Galerkin spatial discretization, and it is shown that this procedure is equivalent to the time-discretization approach of Donea (1984). Numerical results for sample problems are presented graphically, including such shallow-water problems as the breaking of a dam, the shoaling of a wave, and the outflow of a river; compressible flows such as the isothermal flow in a nozzle and the Riemann shock-tube problem; and the two-dimensional scalar-advection, nonlinear-shallow-water, and Euler equations.

Loehner, R.↗

NASA NTRS2001

A Highly Accurate Technique for the Treatment of Flow Equations at the Polar Axis in Cylindrical Coordinates using Series Expansions

Numerical methods for solving the flow equations in cylindrical or spherical coordinates should be able to capture the behavior of the exact solution near the regions where the particular form of the governing equations is singular. In this work we focus on the treatment of these numerical singularities for finite-differences methods by reinterpreting the regularity conditions developed in the context of pseudo-spectral methods. A generally applicable numerical method for treating the singularities present at the polar axis, when nonaxisymmetric flows are solved in cylindrical, coordinates using highly accurate finite differences schemes (e.g., Pade schemes) on non-staggered grids, is presented. Governing equations for the flow at the polar axis are derived using series expansions near r=0. The only information needed to calculate the coefficients in these equations are the values of the flow variables and their radial derivatives at the previous iteration (or time) level. These derivatives, which are multi-valued at the polar axis, are calculated without dropping the accuracy of the numerical method using a mapping of the flow domain from (0,R)*(0,2pi) to (-R,R)*(0,pi), where R is the radius of the computational domain. This allows the radial derivatives to be evaluated using high-order differencing schemes (e.g., compact schemes) at points located on the polar axis. The proposed technique is illustrated by results from simulations of laminar-forced jets and turbulent compressible jets using large eddy simulation (LES) methods. In term of the general robustness of the numerical method and smoothness of the solution close to the polar axis, the present results compare very favorably to similar calculations in which the equations are solved in Cartesian coordinates at the polar axis, or in which the singularity is removed by employing a staggered mesh in the radial direction without a mesh point at r=0, following the method proposed recently by Mohseni and Colonius (1). Extension of the method described here for incompressible flows or for any other set of equations that are solved on a non-staggered mesh in cylindrical or spherical coordinates with finite-differences schemes of various level of accuracy is immediate.

Constantinescu, George S.↗

DOE OSTI2016

A critical appraisal of the zero-multipole method: Structural, thermodynamic, dielectric, and dynamical properties of a water system

We performed extensive and strict tests for the reliability of the zero-multipole (summation) method (ZMM), which is a method for estimating the electrostatic interactions among charged particles in a classical physical system, by investigating a set of various physical quantities. This set covers a broad range of water properties, including the thermodynamic properties (pressure, excess chemical potential, constant volume/pressure heat capacity, isothermal compressibility, and thermal expansion coefficient), dielectric properties (dielectric constant and Kirkwood-G factor), dynamical properties (diffusion constant and viscosity), and the structural property (radial distribution function). We selected a bulk water system, the most important solvent, and applied the widely used TIP3P model to this test. In result, the ZMM works well for almost all cases, compared with the smooth particle mesh Ewald (SPME) method that was carefully optimized. In particular, at cut-off radius of 1.2 nm, the recommended choices of ZMM parameters for the TIP3P system are α ≤ 1 nm{sup −1} for the splitting parameter and l = 2 or l = 3 for the order of the multipole moment. We discussed the origin of the deviations of the ZMM and found that they are intimately related to the deviations of the equilibrated densities between the ZMM and SPME, while the magnitude of the density deviations is very small.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Deployment of the Finite Volume Method in Pronghorn for Gas and Salt cooled Pebble Bed Reactors

This report summarizes the activities related to ”Complete FHR and HTGR pebble bed simulator, including initial validation” funded by the NEAMS thermal-hydraulics focus area. The activity revolves around the coarse-mesh thermal-hydraulics code Pronghorn and its application to gas and salt cooled Pebble bed reactor (PBR). The main difference between gas and salt cooled PBR from a thermal-hydraulics perspective is the fluid. To address the difference in fluid behavior, two separate approaches are implemented in the MOOSE Navier Stokes module: 1) a Boussinesq approximation and 2) a fully compressible formulation. The developed finite volume method capabilities are used for improving pre-existing gas-cooled and salt-cooled pebble-bed reactor models. A steady-state, multiphysics gas-cooled pebble-bed reactor model is created that couples the equilibrium core depletion capability developed in previous work, and the finite volume method capability developed for this report. The salt-cooled pebble-bed reactor model is upgraded to use the incompressible finite volume method capability and then extended to three spatial dimensions. Finally, several verification-and-validation exercises performed with Pronghorn are documented using the verification-and-validation report of the MooseDoc system. The goal of this effort to document the verification-and-validation level of Pronghorn and improve stakeholder confidence in the results obtained with Pronghorn.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1975

On the calculation of supersonic separating and reattaching flows

A method is developed for solving the laminar and turbulent compressible boundary layer equations for separating and reattaching flows. Results of this method are compared with experimental data for two laminar and three turbulent layer, shock wave interactions. Several Navier-Stokes solutions are obtained for each of the laminar boundary layer, shock wave interactions considered. Comparison of these solutions indicates a first order sensitivity in C sub f to the computational mesh selected in both the viscous and inviscid portions of the flow. Comparison of the present boundary layer solutions with the Navier-Stokes solutions and with data for a given Mach number indicates that as long as the separation bubble is small, the boundary layer approximation yields solutions whose accuracy is comparable to the Navier-Stokes solutions.

Murphy, J. D.↗

NASA NTRS1991

Numerical methods in laminar and turbulent flow; Proceedings of the 7th International Conference, Stanford Univ., CA, July 15-19, 1991. Vol. 7, pts. 1 & 2

Consideration is given to the impulse response of a laminar boundary layer and receptivity; numerical transition to turbulence in plane Poiseuille flow; large eddy simulation of turbulent wake flow; a viscous model and loss calculation of a multisplitter cascade; vortex initiation during dynamic stall of an airfoil; a numerical analysis of isothermal flow in a combustion chamber; and compressible flow calculations with a two-equation turbulence model and unstructured grids. Attention is also given to a 2D calculation of a buoyant flow around a burning sphere, a fast multigrid method for 3D turbulent incompressible flows, a streaming flow induced by an oscillating cascade of circular cylinders, an algebraic multigrid scheme for solving the Navier-Stokes equations on unstructured meshes; and nonlinear coupled multigrid solutions to thermal problems employing different nodal grid arrangements and convective transport approximations.

Taylor, C.↗

NASA NTRS2006

Hybrid Deployable Foam Antennas and Reflectors

Hybrid deployable radio antennas and reflectors of a proposed type would feature rigid narrower apertures plus wider adjoining apertures comprising reflective surfaces supported by open-cell polymeric foam structures (see figure). The open-cell foam structure of such an antenna would be compressed for compact stowage during transport. To initiate deployment of the antenna, the foam structure would simply be released from its stowage mechanical restraint. The elasticity of the foam would drive the expansion of the foam structure to its full size and shape. There are several alternatives for fabricating a reflective surface supported by a polymeric foam structure. One approach would be to coat the foam with a metal. Another approach would be to attach a metal film or a metal-coated polymeric membrane to the foam. Yet another approach would be to attach a metal mesh to the foam. The hybrid antenna design and deployment concept as proposed offers significant advantages over other concepts for deployable antennas: 1) In the unlikely event of failure to deploy, the rigid narrow portion of the antenna would still function, providing a minimum level of assured performance. In contrast, most other concepts for deploying a large antenna from compact stowage are of an "all or nothing" nature: the antenna is not useful at all until and unless it is fully deployed. 2) Stowage and deployment would not depend on complex mechanisms or actuators, nor would it involve the use of inflatable structures. Therefore, relative to antennas deployed by use of mechanisms, actuators, or inflation systems, this antenna could be lighter, cheaper, amenable to stowage in a smaller volume, and more reliable. An open-cell polymeric (e.g., polyurethane) foam offers several advantages for use as a compressible/expandable structural material to support a large antenna or reflector aperture. A few of these advantages are the following: 3) The open cellular structure is amenable to compression to a very small volume - typically to 1/20 of its full size in one dimension. 4) At a temperature above its glass-transition temperature (T(sub g)), the foam strongly damps vibrations. Even at a temperature below T(sub g), the damping should exceed that of other materials. 5) In its macroscopic mechanical properties, an open-cell foam is isotropic. This isotropy facilitates computational modeling of antenna structures. 6) Through chemical formulation, the T(sub g) of an open-cell polyurethane foam can be set at a desired value between about - 100 and about 0 C. Depending on the application, it may or may not be necessary to rigidify a foam structure after deployment. If rigidification is necessary, then the T(sub g) of the foam can be tailored to exceed the temperature of the deployment environment, in conjunction with providing a heater to elasticize the foam for deployment. Once deployed, the foam would become rigidified by cooling to below T(sub g). 7) Techniques for molding or machining polymeric foams (especially including open-cell polyurethane foams) to desired sizes and shapes are well developed.

Rivellini, Tommaso↗