SEARCH · Engineering Papers

Results for “Geometry meshing”

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 253 records · Page 14

All sources NASA NTRS DOE OSTI

DOE OSTI2014

Analysis of ballistic transport in nanoscale devices by using an accelerated finite element contact block reduction approach

An accelerated Finite Element Contact Block Reduction (FECBR) approach is presented for computational analysis of ballistic transport in nanoscale electronic devices with arbitrary geometry and unstructured mesh. Finite element formulation is developed for the theoretical CBR/Poisson model. The FECBR approach is accelerated through eigen-pair reduction, lead mode space projection, and component mode synthesis techniques. The accelerated FECBR is applied to perform quantum mechanical ballistic transport analysis of a DG-MOSFET with taper-shaped extensions and a DG-MOSFET with Si/SiO{sub 2} interface roughness. The computed electrical transport properties of the devices obtained from the accelerated FECBR approach and associated computational cost as a function of system degrees of freedom are compared with those obtained from the original CBR and direct inversion methods. The performance of the accelerated FECBR in both its accuracy and efficiency is demonstrated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

The Meshing Complexity of a Solid: An Introduction

This paper explores potential methods for characterizing the meshing complexity of solid geometry. While numerous metrics exist to measure the quality of the finite element, there are currently no metrics that measure the quality of a solid with respect to its meshing complexity. The meshing complexity of a solid is defined by how difficult it is to generate a valid finite element mesh for a given solid. There are many variables that affect meshing complexity. This paper seeks to discuss methods that are decoupled from more subjective variables such as user expertise and software maturity, and it will focus on methods that describe the topological and geometric aspects of a solid. It will present techniques based on: medial axis transformation, wavelets, curvature, proximity, intersection, heuristic topology search, and the measurement of space (volume/area/length) and will analyze their suitability as meshing complexity metrics.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2018

Adaptive wavelet compression of large additive manufacturing experimental and simulation datasets

New manufacturing technologies such as additive manufacturing require research and development to minimize the uncertainties in the produced parts. The research involves experimental measurements and large simulations, which result in huge quantities of data to store and analyze. We address this challenge by alleviating the data storage requirements using lossy data compression. We select wavelet bases as the mathematical tool for compression. Unlike images, additive manufacturing data is often represented on irregular geometries and unstructured meshes. Thus, we use Alpert tree-wavelets as bases for our data compression method. We first analyze different basis functions for the wavelets and find the one that results in maximal compression and miminal error in the reconstructed data. We then devise a new adaptive thresholding method that is data-agnostic and allows a priori estimation of the reconstruction error. Finally, we propose metrics to quantify the global and local errors in the reconstructed data. One of the error metrics addresses the preservation of physical constraints in reconstructed data fields, such as divergence-free stress field in structural simulations. While our compression and decompression method is general, we apply it to both experimental and computational data obtained from measurements and thermal/structural modeling of the sintering of a hollow cylinder from metal powders using a Laser Engineered Net Shape process. The results show that monomials achieve optimal compression performance when used as wavelet bases. The new thresholding method results in compression ratios that are two to seven times larger than the ones obtained with commonly used thresholds. Overall, adaptive Alpert tree-wavelets can achieve compression ratios between one and three orders of magnitude depending on the features in the data that are required to preserve. Furthermore, these results show that Alpert tree-wavelet compression is a viable and promising technique to reduce the size of large data structures found in both experiments and simulations.

42 ENGINEERING↗

DOE OSTI2022

Nanostructural and rheological transitions of pH-responsive supramolecular systems involving a zwitterionic amphiphile and a triamine

Supramolecular aqueous complexes of a synthesized zwitterionic surfactant, 2-(dimethyl(octadecyl)ammonio)acetate (stearyl betaine) and diethylenetriamine, have been found to display pH-tunable rheological properties without the need for a co-surfactant. For this supramolecular system, the dynamically changing interplay among the constituents is responsible for the rich phase behavior including hollow nanotubules, bilayers, and vesicles. A further hypothesis is that hollow nanotubular structure is highly desirable for achieving chain entanglement geometries with reduced mesh size. The rheological properties of the suspension were measured using a rotational rheometer. An analysis using zero-shear viscosities at various temperatures as well as differential scanning calorimetry(DSC) helped shed light on thermodynamic and phase transitions taking place. The morphology of the suspension was analyzed by using a combination of small angle x-ray scattering(SAXS) and atomic force microscopy(AFM). Large steady-shear viscosities (up to ~160 Pa.s) and viscoelastic behavior was observed in acidic conditions which diminished upon increasing the pH. A novel analysis consisting of measured thermodynamic quantities such as activation energy(E a ) and enthalpy of transition(ΔH T ) further allude to the rich pH-sensitive nanoarchitecture of the system. Rather than commonly observed wormlike micelles, hollow nanotubules, as long as 1.6 µm were observed at pH 4 which disintegrated to bilayer sheets and vesicles upon increasing pH. Finally, gaining a comprehensive understanding from all the results, a transition mechanism between different suspension architectures was proposed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

TokaMaker: An open-source time-dependent Grad-Shafranov tool for the design and modeling of axisymmetric fusion devices

Here, in this paper, we present a new static and time-dependent MagnetoHydroDynamic (MHD) equilibrium code, TokaMaker, for axisymmetric configurations of magnetized plasmas, based on the well-known Grad-Shafranov equation. This code utilizes finite element methods on an unstructured triangular grid to enable capturing accurate machine geometry and simple mesh generation from engineering-like descriptions of present and future devices. The new code is designed for ease of use without sacrificing capability and speed through a combination of Python, Fortran, and C/C++ components. A detailed description of the numerical methods of the code, including a novel formulation of the boundary conditions for free-boundary equilibria, and validation of the implementation of those methods using both analytic test cases and cross-code validation is shown. Results show expected convergence across tested polynomial degree for analytic and cross-code test cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2018

A Verified Conformal Decomposition Finite Element Method for Implicit, Many-Material Geometries

As computing power rapidly increases, quickly creating a representative and accurate discretization of complex geometries arises as a major hurdle towards achieving a next generation simulation capability. Component definitions may be in the form of solid (CAD) models or derived from 3D computed tomography (CT) data, and creating a surface-conformal discretization may be required to resolve complex interfacial physics. Furthermore, the Conformal Decomposition Finite Element Methods (CDFEM) has been shown to be an efficient algorithm for creating conformal tetrahedral discretizations of these implicit geometries without manual mesh generation.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Mathematical Models and Numerical Methods for High-Fidelity Simulation of Ignition of Reactive Mixtures by Nanosecond Plasma Discharges in Realistic Configurations

We present a newly developed framework for the numerical simulation of ignition of reactive mixtures using single or repeated nanosecond discharge pulses. The framework builds upon the AMReX library, using the existing compressible solver PeleC and low-Mach solver PeleLMeX and allowing for adaptive mesh refinement, complex geometries, and execution on next-generation high-performance computing (HPC) systems. High-fidelity elementary models are adopted for weakly-ionised plasma discharges with significant energy deposition, consistent with nanosecond discharge pulses, and then implemented in the solver. The treatment of non-thermal electrons and charged species, thermodynamics of non-equilbrium species, plasma kinetics, limiting time scales, and boundary conditions for charged species are discussed and addressed for computational efficiency. The framework is demonstrated for three relevant applications: single and multi-pulse discharges in air, single pulse ignition of an ethylene/air mixture, and a three-dimensional plasma discharge in air with temperature stratification. The successful application of the framework demonstrates the feasibility of high-fidelity simulation of ignition of air/hydrocarbon mixtures in three-dimensions with multiple discharge pulses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

The Method of Finite Averages

The Method of Finite Averages (MoFA) is a rigorous multiscale modeling methodology for efficiently modeling multi-physical phenomena in heterogeneous porous media. The code developed in this project aims to perform the numerical calculations required to formulate, implement, and verify MoFA models for Earth and Energy systems (i.e., model verification refers to performing fully-resolved simulations of the systems and comparing their results to those of the models). In general, MoFA transforms partial differential equations (PDEs) describing the fine-scale physics of a system into coupled ordinary differential equations (ODEs)---in time---that describe the coarse-scale---or "average"---physical behaviors of the system. This transformation significantly expedites system simulation, as the coarse-scale ODEs involve vastly fewer degrees of freedom than the fine-scale PDEs. The code developed under this project will allow users to 1.) generate system geometries and numerical meshes, 2.) solve the PDE and ODE systems required for MoFA model formulation and implementation, 3.) solve the PDE systems required to obtain fully-resolved simulation results for model verification, and 4.) compare and plot results (e.g., the model and fully-resolved simulation solutions, the error between the solutions, etc.).

Pietrzyk, KyleM [Lawrence Livermore National Labor↗

DOE OSTI2024

CurvilinearGrids.jl: A Julia package for curvilinear coordinate transformations

Finite-difference discretizations of partial differential equations are widespread throughout the scientific community. Oftentimes finite-differences are used to compute spatial gradients of fields on a discrete grid, which is typically a uniform or rectilinear Cartesian mesh. Arbitrary multidimensional geometry is difficult to discretize directly with finite differences, however, due to non-uniform grid spacing and non-orthogonality. Curvilinear coordinate transformations can be used as a strategy to enable arbitrary geometry. While these curvilinear transformations are straightforward, the governing PDEs require additional terms (metrics) and must adhere to strict conservation laws; these criteria complicate the application of the transformation and require careful implementation.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2016

Proof-of-Concept Study for Uncertainty Quantification and Sensitivity Analysis using the BRL Shaped-Charge Example

These are the slides for a graduate presentation at Mississippi State University. It covers the following: the BRL Shaped-Charge Geometry in PAGOSA, mesh refinement study, surrogate modeling using a radial basis function network (RBFN), ruling out parameters using sensitivity analysis (equation of state study), uncertainty quantification (UQ) methodology, and sensitivity analysis (SA) methodology. In summary, a mesh convergence study was used to ensure that solutions were numerically stable by comparing PDV data between simulations. A Design of Experiments (DOE) method was used to reduce the simulation space to study the effects of the Jones-Wilkins-Lee (JWL) Parameters for the Composition B main charge. Uncertainty was quantified by computing the 95% data range about the median of simulation output using a brute force Monte Carlo (MC) random sampling method. Parameter sensitivities were quantified using the Fourier Amplitude Sensitivity Test (FAST) spectral analysis method where it was determined that detonation velocity, initial density, C1, and B1 controlled jet tip velocity.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

DOE OSTI2017

Geomechanical Simulation of Bayou Choctaw Strategic Petroleum Reserve - Model Calibration.

A finite element numerical analysis model has been constructed that consists of a realistic mesh capturing the geometries of Bayou Choctaw (BC) Strategic Petroleum Reserve (SPR) site and multi - mechanism deformation ( M - D ) salt constitutive model using the daily data of actual wellhead pressure and oil - brine interface. The salt creep rate is not uniform in the salt dome, and the creep test data for BC salt is limited. Therefore, the model calibration is necessary to simulate the geomechanical behavior of the salt dome. The cavern volumetric closures of SPR caverns calculated from CAVEMAN are used for the field baseline measurement. The structure factor, A 2 , and transient strain limit factor, K 0 , in the M - D constitutive model are used for the calibration. The A 2 value obtained experimentally from the BC salt and K 0 value of Waste Isolation Pilot Plant (WIPP) salt are used for the baseline values. T o adjust the magnitude of A 2 and K 0 , multiplication factors A2F and K0F are defined, respectively. The A2F and K0F values of the salt dome and salt drawdown skins surrounding each SPR cavern have been determined through a number of back fitting analyses. The cavern volumetric closures calculated from this model correspond to the predictions from CAVEMAN for six SPR caverns. Therefore, this model is able to predict past and future geomechanical behaviors of the salt dome, caverns, caprock , and interbed layers. The geological concerns issued in the BC site will be explained from this model in a follow - up report .

58 GEOSCIENCES↗

DOE OSTI2018

Decrease time-to-solution through improved linear-system setup and solve.

17 WIND ENERGY↗

DOE OSTI2018

ECP ALCC Quarterly Report (Oct-Dec 2017)

The scientific goal of ExaWind Exascale Computing Project (ECP) is to advance our fundamental understanding of the flow physics governing whole wind plant performance, including wake formation, complex terrain impacts, and turbine-turbine-interaction effects. Current methods for modeling wind plant performance fall short due to insufficient model fidelity and inadequate treatment of key phenomena, combined with a lack of computational power necessary to address the wide range of relevant length scales associated with wind plants. Thus, our ten-year exascale challenge is the predictive simulation of a wind plant composed of O(100) multi-MW wind turbines sited within a 100 km2 area with complex terrain, involving simulations with O(100) billion grid points. The project plan builds progressively from predictive petascale simulations of a single turbine, where the detailed blade geometry is resolved, meshes rotate and deform with blade motions, and atmospheric turbulence is realistically modeled, to a multi turbine array in complex terrain. The ALCC allocation will be used continually throughout the allocation period. In the first half of the allocation period, small (e.g., for testing Kokkos algorithms) and medium (e.g., 10K cores for highly resolved ABL simulations) sized jobs will be typical. In the second half of the allocation period, we will also have a number of large submittals for our resolved-turbine simulations. A challenge in the latter period is that small time step sizes will require long wall-clock times for statistically meaningful solutions. As such, we expect our allocation-hour burn rate to increase as we move through the allocation period.

17 WIND ENERGY↗

DOE OSTI2018

Decrease time-to-solution through improved linear-system setup and solve (Milestone Report)

The goal of the ExaWind project is to enable predictive simulations of wind farms composed of many MW-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines, and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. We describe in this report our efforts to decrease the setup and solution time for the mass-continuity Poisson system with respect to the benchmark timing results reported in FY18 Q1. In particular, we investigate improving and evaluating two types of algebraic multigrid (AMG) preconditioners: Classical Ruge-Stfiben AMG (C-AMG) and smoothed-aggregation AMG (SA-AMG), which are implemented in the Hypre and Trilinos/MueLu software stacks, respectively.

17 WIND ENERGY↗

DOE OSTI2019

Geomechanical Simulation of Big Hill Strategic Petroleum Reserve - Calibration of Model Containing Shear Zone

A finite element numerical analysis model, that consists of a realistic mesh capturing the geometries of Big Hill (BH) Strategic Petroleum Reserve (SPR) site using the multi-mechanism deformation (MD) salt constitutive model and including data taken daily of the wellhead pressure and level of the oil-brine interface, has been upgraded. The upgraded model contains the shear zone to examine the interbed behavior in a realistic manner. The salt creep rate is not uniform in the salt dome, and creep test data for BH salt is limited. Therefore, a model calibration is necessary to simulate the geomechanical behavior of the salt dome. Cavern volumetric closures of SPR caverns calculated from sonar survey reports are used for the field baseline measurement. The structure factor, A 2 , and transient strain limit factor, K 0 , in the M-D constitutive model are used for model calibration. An A 2 value obtained experimentally from the BH salt and K 0 value of WIPP salt are used as the baseline values. To adjust the magnitude of A 2 and K 0 , multiplication factors A 2 F and K 0 F are defined, respectively. The A 2 F and K 0 F values of the salt dome and salt drawdown layer of elements surrounding each SPR cavern have been determined through a number of back fitting analyses. The trendlines of the predictions and sonar data match up well for BH 101, 103, 104, 106, 110, 111, 112, and 113. The prediction curves are close to the sonar data for BH 102 and 114. However, the prediction curves for BH 105, 107, 108, and 109 are not close to the sonar data. An inconsistency was found in the sonar data, i.e. the volume measured later is larger than that before in some time intervals, even if the leached volume is taken into account, for BH 101, 104, 106, 107, and 112. Project discussions are needed to determine possibilities on how to resolve the issues and determine the best path forward for future computer modeling attempts.

02 PETROLEUM↗

DOE OSTI2019

Integrated Simulation Capabilities for Analysis of Experiments in the Versatile Test Reactor

The DOE-initiated Versatile Test Reactor (VTR) program aims to establish a domestic fast-neutron irradiation testing capability that meets a variety of domestic and international nuclear energy needs. Currently, integrated tools capable of mechanistic modeling of VTR and a test vehicle do not exist. To address this modeling and analysis deficiency, Argonne's SAS4A/SYSSYS-1 safety analysis code has been coupled with SAM (System Analysis Module) to provide a novel modeling capability that supports VTR and other fast test facilities. This report documents the design, implementation, and testing of this integrated capability. A coupling boundary has been identified at the test vehicle and primary coolant interface, where SAS4A/SASSYS-1 treats primary coolant thermal hydraulics outside the test vehicle, while SAM treats all thermal hydraulic behavior within the test vehicle, including the vehicle walls. Essential to this integrated tool is its newly developed capability to properly model the conjugate heat transfer process which ensures equality of temperatures and heat fluxes at the vehicle wall interface while ensuring energy conservation. A range of testing has been completed to support demonstration of this interface. The testing scope includes steady-state verification using a series of increasingly complex analytical solutions and transient demonstration for a range of design basis accident scenarios using prototypic VTR and test vehicle (i.e. cartridge loop) configurations. Results of this testing confirm verification of the steady-state solutions and provide reasonable and expected transient behavior. The coupling interface has been developed to be robust yet flexible. The implementation within SAS4A/SASSYS-1 supports coupling to any external software or module. The use of decomposed domains and the existing in-core thermal-hydraulic calculational framework within SAS4A/SASSYS-1 enable treatment of arbitrary test vehicle geometries with nonuniform meshing between SAS4A/SASSYS-1 and the coupled code. Modeling of multiple test locations is also supported. Furthermore, computational efficiency is enhanced by utilizing named pipes (or FIFOs), an interface for interprocess communication (IPC).

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2020

Watertight Tensor-Product Spline Reconstruction (KCNSC PDRD Technical Report)

The standard for representing geometric models in Computer Aided Design (CAD) systems is boundary representation with trimmed surfaces. However, the approximate nature of surface-surface intersections in CAD systems leaves many small gaps between mating surfaces. These gaps lead to problems for downstream analysis and manufacturing applications, requiring time consuming repair operations and further approximation of the geometry as polygonal meshes, which do not support design iteration. In this study, nVariate, Inc. evaluated its watertight reconstruction technology on a test part. The resulting watertight model has the same absolute model accuracy of the original model zero gaps, providing a representation that requires no repair and supports further design iterations in the CAD system.

42 ENGINEERING↗

DOE OSTI2018

Geomechanical Simulation of Big Hill Strategic Petroleum Reserve (Model Calibration)

A finite element numerical analysis model has been constructed that consists of a realistic mesh capturing the geometries of Big Hill (BH) Strategic Petroleum Reserve (SPR) site using the multi-mechanism deformation (M-D) salt constitutive model and including data taken daily of the wellhead pressure and level of the oil-brine interface. The salt creep rate is not uniform in the salt dome, and creep test data for BH salt is limited. Therefore, a model calibration is necessary to simulate the geomechanical behavior of the salt dome. Cavern volumetric closures of SPR caverns calculated from sonar survey reports are used for the field baseline measurement. The structure factor, A 2 , and transient strain limit factor, K o , in the M-D constitutive model are used for model calibration. An A 2 value obtained experimentally from the BH salt and K o value of WIPP salt are used as the baseline values. To adjust the magnitude of A 2 and K o , multiplication factors A2F and KOF are defined, respectively. The A2F and KOF values of the salt dome and salt drawdown layer of elements surrounding each SPR cavern have been determined through a number of back fitting analyses. The trendlines of the predictions and sonar data match up well for BH 101, 103, 104, 106, 110, 111, and 113. The prediction curves are close to the sonar data for BH 102 and 114. However, the prediction curves for BH 105, 107, 108, 109, and 112 are not close to the sonar data. An inconsistency was found in the sonar data, i.e. the sonar measurements of cavern volumes increase with time, during some periods for BH 101, 104, 106, 107, and 112. A follow-up report in 2019 will provide a resolution for these issues.

02 PETROLEUM↗