Engineering Papers⌕ Search

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.

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

At least 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Non-Equilibrium Turbulence Modeling for High Lift Aerodynamics

This phase is discussed in ('Non linear kappa - epsilon - upsilon(sup 2) modeling with application to high lift', Application of the kappa - epsilon -upsilon(sup 2) model to multi-component airfoils'). Further results are presented in 'Non-linear upsilon(sup 2) - f modeling with application to high-lift' The ADI solution method in the initial implementation was very slow to converge on multi-zone chimera meshes. I modified the INS implementation to use GMRES. This provided improved convergence and less need for user intervention in the solution process. There were some difficulties with implementation into the NASA compressible codes, due to their use of approximate factorization. The Helmholtz equation for f is not an evolution equation, so it is not of the form assumed by the approximate factorization method. Although The Kalitzin implementation involved a new solution algorithm ('An implementation of the upsilon(sup 2) - f model with application to transonic flows'). The algorithm involves introducing a relaxation term in the f-equation so that it can be factored. The factorization can be into a plane and a line, with GMRES used in the plane. The NASA code already evaluated coefficients in planes, so no additional memory is required except that associated the the GMRES algorithm. So the scope of this project has expanded via these interactions. . The high-lift work has dovetailed into turbine applications.

Durbin, P. A.↗
NASA NTRS2019

Dynamic Smagorinsky Modeled Large-Eddy Simulations of Turbulence Using Tetrahedral Meshes

Eddy-resolving numerical computations of turbulent flows are emerging as viable alternatives to Reynolds Averaged Navier-Stokes (RANS) calculations for flows with an intrinsically steady mean state due to the advances in large-scale parallel computing. In these computations, medium to large turbulent eddies are resolved by the numerics while the smaller or subgrid scales are either modeled or taken care of by the inherent numerical dissipation. To advance the state of the art of unstructured-mesh turbulence simulation capabilities, large eddy simulations (LES) using the dynamic Smagorinsky model (DSM) on tetrahedral meshes are carried out with the space-time conservation element, solution element (CESE) method. In contrast to what has been reported in the literature, the present implementation of dynamic models allows for active backscattering without any ad-hoc limiting of the eddy viscosity calculated from the subgrid-scale model. For the benchmark problems involving compressible isotropic turbulence decay as well as the shock/turbulent boundary layer interaction benchmark problems, no numerical instability associated with kinetic energy growth is observed and the volume percentage of the backscattering portion accounts for about 38-40% of the simulation domain. A slip-wall model in conjunction with the implemented DSM is used to simulate a relatively high Reynolds number Mach 2.85 turbulent boundary layer over a 30° ramp with several tetrahedral meshes and a wall-normal spacing of either Δγ& = 10 or Δγ& = 20. The computed mean wall pressure distribution, separation region size, mean velocity profiles, and Reynolds stress agree reasonably well with experimental data.

CESE method↗
NASA NTRS1981

Full potential solution of transonic quasi-3-D flow through a cascade using artificial compressability

The three-dimensional flow in a turbomachinery blade row was approximated by correcting for streamtube convergence and radius change in the throughflow direction. The method is a fully conservative solution of the full potential equation incorporating the finite volume technique on body fitted periodic mesh, with an artificial density imposed in the transonic region to insure stability and the capture of shock waves. Comparison of results for several supercritical blades shows good agreement with their hodograph solutions. Other calculations for these profiles as well as standard NACA blade sections indicate that this is a useful scheme analyzing both the design and off-design performance of turbomachinery blading.

Farrell, C.↗
DOE OSTI2024

A B-spline based gradient-enhanced micropolar implicit material point method for large localized inelastic deformations

The quasi-brittle response of cohesive-frictional materials in numerical simulations is commonly represented by softening plasticity or continuum damage models, either individually or in combination. However, classical models, particularly when coupled with non-associated plasticity, often suffer from ill-posedness and a lack of objectivity in numerical simulations. Moreover, the performance of the finite element method significantly degrades in simulations involving finite strains when mesh distortion reaches excessive levels. This represents a challenge for modeling cohesive-frictional materials, given their tendency to experience strongly localized deformations, such as those occurring during shear band dominated failure. Hence, accurate modeling of the response of cohesive-frictional solids is a demanding task. To address these challenges, we present an extension of the material point method (MPM) for the unified gradient-enhanced micropolar continuum, aiming at the analysis of finite localized inelastic deformations in cohesive-frictional materials. The generalized gradient-enhanced micropolar continuum formulation is employed to tackle challenges related to localization and softening material behavior, while the MPM addresses issues arising from excessive deformations. The method utilizes a B-spline formulation for the rigid background mesh to mitigate the well-known cell crossing errors of the MPM. To demonstrate the performance of the method, 2D and 3D numerical studies on localized failure in sandstone in plane strain compression and triaxial extension tests are presented. A comparison with finite element results confirms the suitability of the formulation. Moreover, an efficient numerical implementation of the formulation is presented, and it is demonstrated that the additional MPM specific overhead is negligible.

B-spline↗
DOE OSTI2023

A conservative, interface-resolved, compressible framework for the modeling and simulation of liquid/gas phase change

This paper presents a method for simulating evaporation in a compressible, interface-resolved framework appropriate for modeling problems of engineering interest. In order to achieve robustness and broad applicability, the method has been designed to discretely enforce consistent mass and thermal energy transport at the phase interface, to globally conserve mass, momentum, and energy, and to be capable of modeling compressible and incompressible systems. Verification is performed via the Sod-shock test, one-dimensional heat conduction, evaporation from a planar interface, and evaporation of three-dimensional droplets. Convergence with increasing mesh resolution is demonstrated in all tested configurations, and conservation is maintained near machine precision for a translating droplet. Conservation and accurate phase change rates are preserved at the low numerical resolutions commonly encountered in engineering calculations. Following verification, the method is validated by comparison to an empirical correlation for evaporating droplets in high temperature crossflow, and the presentation concludes with the simulation of an iso-octane spray at conditions representative of gasoline direct injection. In conclusion, successful verification, validation, and demonstrated practical utility suggest the method to be an accurate, efficient, and robust approach for the study of phase change in engineering systems.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1977

Numerical solution of axisymmetric boattail flow fields with plume simulators

Turbulent separating flows over axisymmetric afterbody-boattail configurations with solid sting plume simulators are computed with a time-dependent finite-difference method to solve the compressible Navier-Stokes equations. The Reynolds stress terms are replaced with a two-layer eddy viscosity model including a relaxation formula to model the nonequilibrium effects of the separated flow. The mesh is alined with the boattail body through an analytic transformation which accommodates a wide variety of boattail geometries. Numerical results for a series of boattail geometries over a wide range of Reynolds number (140,000 to 140 million) are presented and, when possible, compared with experimental data or independent numerical results.

Holst, T. L.↗
NASA NTRS1980

Numerical calculation of steady inviscid full potential compressible flow about wind turbine blades

An exact nonlinear mathematical model that accounts for three-dimensional cascade effects about the inner portions of the rotor blades and compressibility effects about the tip regions of the blades was derived. An artificially time dependent version was iteratively solved by a finite volume technique involving an artificial viscosity and a three-level consecutive mesh refinement. The exact boundary conditions were applied by generating a boundary conforming periodic computation mesh.

Dulikravich, D. S.↗
DOE OSTI2026

Scaling kinetic Monte-Carlo simulations of grain growth with combined convolutional and graph neural networks

Graph neural networks (GNN) have emerged as a promising machine learning method for microstructure simulations such as grain growth. However, accurate modeling of realistic grain boundary networks requires large simulation cells, which GNN has difficulty scaling up to. To alleviate the computational costs and memory footprint of GNN, we suggest a hybrid architecture combining a convolutional neural network (CNN) based bijective autoencoder to compress the spatial dimensions, and a GNN that evolves the microstructure in the latent space of reduced spatial sizes. Our results demonstrate that the new design significantly reduces computational costs with using fewer message passing layer (from 12 down to 3) compared with GNN alone. The reduction in computational cost becomes more pronounced as the spatial size increases, indicating strong computational scalability. For the largest mesh evaluated (160 3 ), our method reduces memory usage and runtime in inference by 117× and 115×, respectively, compared with GNN-only baseline. More importantly, it shows higher accuracy and stronger spatiotemporal capability than the GNN-only baseline, especially in long-term testing. Such combination of scalability and accuracy is essential for simulating realistic material microstructures over extended time scales. The improvements can be attributed to the bijective autoencoder’s ability to compress information losslessly from spatial domain into a high dimensional feature space, thereby producing more expressive latent features for the GNN to learn from, while also contributing its own spatiotemporal modeling capability. Training data are generated from stochastic grain growth simulations, providing realistic variability for learning robust microstructure evolution. Comprehensive system validation confirms that the model is accurate, robust, and scalable.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A novel post-processing method for progressive failure analysis of brittle composite compression

Finite element analysis of brittle materials in axial compression typically uses element deletion to allow continued global deformation post-element-failure. However, element deletion produces cyclic load-displacement curves that underestimate energy absorption and are not representative of a continuum system. Two key observations support the conclusion that results from an appropriately discretized model can be an adequate representation of a continuum system. Specifically, the frequency of the oscillations in the load-displacement curve is directly dependent upon element length in the loading direction, and the peak amplitudes of oscillations are mesh size independent. A method of post-processing the analysis results, by connecting the peak amplitudes of oscillations, is proposed and applied to a series of continuous carbon fiber composite crush tubes. The load-displacement curve, stable crushing load, and specific energy absorption of the post-processed results compare well to an experimental study of crush tubes with similar layups.

Materials Science↗
DOE OSTI2016

Trigger Detection for Adaptive Scientific Workflows Using Percentile Sampling

The increasing complexity of both scientific simulations and high-performance computing system architectures are driving the need for adaptive workflows, in which the composition and execution of computational and data manipulation steps dynamically depend on the evolutionary state of the simulation itself. Consider, for example, the frequency of data storage. Critical phases of the simulation should be captured with high frequency and with high fidelity for postanalysis; however, we cannot afford to retain the same frequency for the full simulation due to the high cost of data movement. We can instead look for triggers, indicators that the simulation will be entering a critical phase, and adapt the workflow accordingly. In this paper, we present a methodology for detecting triggers and demonstrate its use in the context of direct numerical simulations of turbulent combustion using S3D. We show that chemical explosive mode analysis (CEMA) can be used to devise a noise-tolerant indicator for rapid increase in heat release. However, exhaustive computation of CEMA values dominates the total simulation, and thus is prohibitively expensive. To overcome this computational bottleneck, we propose a quantile sampling approach. Our sampling-based algorithm comes with provable error/confidence bounds, as a function of the number of samples. Most importantly, the number of samples is independent of the problem size, and thus our proposed sampling algorithm offers perfect scalability. Our experiments on homogeneous charge compression ignition and reactivity controlled compression ignition simulations show that the proposed method can detect rapid increases in heat release, and its computational overhead is negligible. Our results will be used to make dynamic workflow decisions regarding data storage and mesh resolution in future combustion simulations.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1988

Finite element solution of the Euler equations in two and three dimensions

The paper presents a finite element procedure for solving the equations of compressible flow over bodies of arbitrary geometry. The numerical solution algorithm employed is an explicit two-step version of a second order Taylor-Galerkin scheme. The discretization of the computational domain into unstructured meshes of triangles in two dimensions and tetrahedra in three dimensions is performed by an automatic mesh generator. In the approach, the mesh generator is coupled to the finite element solver to produce an adaptive remeshing procedure.

Peiro, J.↗
DOE OSTI2017

A Simulation Framework for Battery Cell Impact Safety Modeling Using LS-DYNA

The development process of electrified vehicles can benefit significantly from computer-aided engineering tools that predict themultiphysics response of batteries during abusive events. A coupled structural, electrical, electrochemical, and thermal model framework has been developed within the commercially available LS-DYNA software. The finite element model leverages a three-dimensional mesh structure that fully resolves the unit cell components. The mechanical solver predicts the distributed stress and strain response with failure thresholds leading to the onset of an internal short circuit. In this implementation, an arbitrary compressive strain criterion is applied locally to each unit cell. A spatially distributed equivalent circuit model provides an empirical representation of the electrochemical responsewith minimal computational complexity.The thermalmodel provides state information to index the electrical model parameters, while simultaneously accepting irreversible and reversible sources of heat generation. The spatially distributed models of the electrical and thermal dynamics allow for the localization of current density and corresponding temperature response. The ability to predict the distributed thermal response of the cell as its stored energy is completely discharged through the short circuit enables an engineering safety assessment. A parametric analysis of an exemplary model is used to demonstrate the simulation capabilities.

25 ENERGY STORAGE↗
NASA NTRS1993

A comprehensive comparison of turbulence models in the far wake

In the present study, the far wake was examined numerically using an implicit, upwind, finite-volume, compressible Navier-Stokes code. The numerical grid started at 500 equivalent circular cylinder diameters in the wave, and extended to 4000 equivalent diameters. By concentrating only on the far wake, the numerical difficulties and fine mesh requirements near the wake-generating body were eliminated. At the time of this writing, results for the K-epsilon and K-omega turbulence models at low Mach number have been completed and show excellent agreement with previous incompressible results and far-wake similarity solutions. The code is presently being used to compare the performance of various other turbulence models, including Reynolds stress models and the new anisotropic two-equation turbulence models being developed at NASA Langley. By increasing our physical understanding of the deficiencies and limits of these models, it is hoped that improvements to the universality of the models can be made. Future plans include examination of two-dimensional momentumless wakes as well.

Cimbala, John M.↗
DOE OSTI2022

Pele: An Exascale-Ready Suite of Combustion Codes

High fidelity simulations of realistic combustion devices are extremely demanding computationally because of the requirements to capture complex fuel chemical decomposition, its intricate interactions with turbulent, often multiphase, flows, and the wide separation of space and time scales between the thin flame and the device boundaries. Software required to carry out such computations tends to be extremely complex, particularly when designed to exploit hardware accelerators, and can be difficult to port and maintain. We present Pele, a performance portable suite of tools for the simulation of combustion systems, including codes to evolve reactive multiphase configurations in the low Mach number and compressible flow regimes, along with a set of inter-compatible post processing and in situ analysis tools. The Pele suite of tools is built on top of the AMReX framework for block-structured adaptive mesh refinement, which provides efficient data structures and algorithms that enable the development of a wide variety of efficient mesh and particle based PDE integration schemes. A hierarchical MPI+X parallelism scheme supports CPU-only and accelerated architectures, where X can be OpenMP, CUDA, and HIP based approaches for intra-node computational work distribution. The algorithms and data structures underlying the Pele simulation and analysis tools are highly scalable and performant across a wide variety of high-performance computing platforms, including DOEs newest exascale-class machines, Frontier and Aurora. The simulation and analysis tools are fully documented and freely distributed as open source via GitHub. We present key algorithmic and software challenges, solution strategies, performance and resulting set of capabilities.

AMReX↗
DOE OSTI2017

In situ imaging during compression of plastic bonded explosives for damage modeling

Here, the microstructure of plastic bonded explosives (PBXs) is known to influence behavior during mechanical deformation, but characterizing the microstructure can be challenging. For example, the explosive crystals and binder in formulations such as PBX 9501 do not have sufficient X-ray contrast to obtain three-dimensional data by in situ, absorption contrast imaging. To address this difficulty, we have formulated a series of PBXs using octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) crystals and low-density binder systems. The binders were hydroxyl-terminated polybutadiene (HTPB) or glycidyl azide polymer (GAP) cured with a commercial blend of acrylic monomers/oligomers. The binder density is approximately half of the HMX, allowing for excellent contrast using in situ X-ray computed tomography (CT) imaging. The samples were imaged during unaxial compression using micro-scale CT in an interrupted in situ modality. The rigidity of the binder was observed to significantly influence fracture, crystal-binder delamination, and flow. Additionally, 2D slices from the segmented 3D images were meshed for finite element simulation of the mesoscale response. At low stiffness, the binder and crystal do not delaminate and the crystals move with the material flow; at high stiffness, marked delamination is noted between the crystals and the binder, leading to very different mechanical properties. Initial model results exhibit qualitatively similar delamination.

36 MATERIALS SCIENCE↗
DOE OSTI2012

A compatible Lagrangian hydrodynamic scheme for multicomponent flows with mixing

We present a Lagrangian time integration scheme and compatible discretization for total energy conservation in multicomponent mixing simulations. Mixing behavior results from relative motion between species. Species velocities are determined by solving species momentum equations in a Lagrangian manner. Included in the species momentum equations are species artificial viscosity (since each species can undergo compression) and inter-species momentum exchange. Thermal energy for each species is also solved, including compression work and thermal dissipation caused by momentum exchange. The present procedure is applicable to mixing of an arbitrary number of species that may not be in pressure or temperature equilibrium. A traditional staggered stencil has been adopted to describe motion of each species. The computational mesh for the mixture is constructed in a Lagrangian manner using the mass-averaged mixture velocity. Species momentum equations are solved at the vertices of the mesh, and temporary species meshes are constructed and advanced in time using the resulting species velocities. Following the Lagrangian step, species quantities are advected (mapped) from the species meshes to the mixture mesh. Momentum exchange between species introduces work that must be included in an energy-conserving discretization scheme. This work has to be transformed to dissipation in order to effect a net change in species thermal energy. The dissipation between interacting species pairs is obtained by combining the momentum exchange work. The dissipation is then distributed to the species involved using a distribution factor based on species specific heats. The resulting compatible discretization scheme provides total energy conservation of the whole mixture. In addition, the numerical scheme includes conservative local energy exchange between species in mixture. Due to the relatively large species interaction coefficients, both the species momenta and energies are calculated implicitly. Sample calculations have yielded excellent results, including conservation of total energy in Lagrangian steps, symmetry preservation, and correct steady-state behavior.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

SWinzip v. 1.0

SWinzip is a Matlab and C++ library for scientific lossy data compression and reconstruction using compressed sensing and tree-wavelets transforms. These methods are known for their large compression and usefulness in data analytics such as features extraction. Compressed sensing and wavelets methods rely heavily on sparse and dense linear algebra operations implemented through the Boost codes in our library. SWinzip accommodates data represented on both regular grids (e.g. image data) and point-clouds (e.g. unstructured meshes).

Salloum, Maher↗
DOE OSTI2020

Progressive collapse analysis of stiffened panel structure with consideration of actual welding distortion and residual stresses

To date, inclusion of welding distortion and residual stresses in determining the strength of marine structures is based on empirical formulations due to welding complexity. In this study, a more feasible numerical model was developed from the perspective of inherent strain theory. The accuracy of this simplified numerical model was confirmed by transient thermo-elastic–plastic finite element analysis using a novel dynamic mesh refining and iterative substructure method which is 50 times faster than conventional finite element method. Welding distortion by the two numerical approaches was validated with experimental measurement. The collapse analyses indicated that the inherent strain model predicted a compressive strength and load-displacement curve with high accuracy, as validated by the transient thermo-mechanical model. The global stiffness and the ultimate strength of the welded structure would be overestimated by about 22% if no welding effects were considered. The proposed computation procedure for determining the ultimate strength of welded structures is accurate, simple, and practical and enables improvement of ultimate strength through optimization of welding joint design to reduce the welding residual stress and distortion.

Huang, Hui↗