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 793 records · Page 44

All sources NASA NTRS DOE OSTI

NASA NTRS1998

User Manual for Beta Version of TURBO-GRD: A Software System for Interactive Two-Dimensional Boundary/ Field Grid Generation, Modification, and Refinement

TURBO-GRD is a software system for interactive two-dimensional boundary/field grid generation. modification, and refinement. Its features allow users to explicitly control grid quality locally and globally. The grid control can be achieved interactively by using control points that the user picks and moves on the workstation monitor or by direct stretching and refining. The techniques used in the code are the control point form of algebraic grid generation, a damped cubic spline for edge meshing and parametric mapping between physical and computational domains. It also performs elliptic grid smoothing and free-form boundary control for boundary geometry manipulation. Internal block boundaries are constructed and shaped by using Bezier curve. Because TURBO-GRD is a highly interactive code, users can read in an initial solution, display its solution contour in the background of the grid and control net, and exercise grid modification using the solution contour as a guide. This process can be called an interactive solution-adaptive grid generation.

Choo, Yung K.↗

NASA NTRS2001

Unic Unstructured CFD Methodology Development

Base heating characteristics is crucial to the success and the overall performance of the X-33 engine. Base heating is important throughout the entire flight trajectory due to the aerospike engine design of X-33. The base region is surrounded by the hot-gas plume, which expands, circulates and impinges on the base. An advanced computation fluid dynamics method is employed in an effort to develop a robust, accurate and efficient tool for the X-33 base heating performance predictions. This computational tool is developed based on a Navier-Stokes flow solver, which is suitable for general complex geometry and includes turbulence, finite-rate chemistry, and radiation models. To fulfill the fast turnaround requirement as a design analysis tool, adaptive mesh refinement method and parallel-computing algorithm are also incorporated in the present model. Case study for the X-33 base-region fluid dynamics and heat transfer characteristics are presented.

Chen, Y. S.↗

NASA NTRSDate not supplied

Investigation of Physical Mechanisms for Jet Noise Reduction by Plug Nozzle Porosity

This paper presents the results of a computational study on porous-plug nozzles for their eventual application in reducing takeoff noise levels for supersonic civil transport. The study explores a wide range of porous-plug geometries, encompassing variations in hole size, plug length and porosity. Computational fluid dynamics (CFD) analyses are employed in concert with an advanced automated mesh refinement (AMR) scheme, Sketch-to-Solution (S2S), to investigate porous plug nozzle design strategies that may diminish shock strength, broad-band shock noise (BBSN), flow separation, and other undesirable characteristics associated with high takeoff noise levels. A selected subset of the porous-plug nozzles that were analyzed computationally were manufactured at subscale using stereolithography 3D printing. These nozzles were tested experimentally, obtaining far-field noise spectral data and schlieren flow visualization pictures. Detailed comparisons between CFD predictions and experimental results were conducted. These comparisons focus on gaining a better understanding of the underlying physical mechanisms responsible for noise reduction with the porous plugs relative to corresponding solid plugs.

computational fluid dynamics↗

NASA NTRSDate not supplied

Investigation of Physical Mechanisms for Jet Noise Reduction by Plug Nozzle Porosity

Nozzles↗

DOE OSTI2023

Deep learning closure models for large-eddy simulation of flows around bluff bodies

Near-wall flow simulation remains a central challenge in aerodynamics modelling: Reynolds-averaged Navier–Stokes predictions of separated flows are often inaccurate, and large-eddy simulation (LES) can require prohibitively small near-wall mesh sizes. A deep learning (DL) closure model for LES is developed by introducing untrained neural networks into the governing equations and training in situ for incompressible flows around rectangular prisms at moderate Reynolds numbers. The DL-LES models are trained using adjoint partial differential equation (PDE) optimization methods to match, as closely as possible, direct numerical simulation (DNS) data. They are then evaluated out-of-sample – for aspect ratios, Reynolds numbers and bluff-body geometries not included in the training data – and compared with standard LES models. The DL-LES models outperform these models and are able to achieve accurate LES predictions on a relatively coarse mesh (downsampled from the DNS mesh by factors of four or eight in each Cartesian direction). We study the accuracy of the DL-LES model for predicting the drag coefficient, near-wall and far-field mean flow, and resolved Reynolds stress. A crucial challenge is that the LES quantities of interest are the steady-state flow statistics; for example, a time-averaged velocity component $\langle {u}_i\rangle (x) = \lim _{t \rightarrow \infty } ({1}/{t}) \int _0^t u_i(s,x)\, {\rm d}s$ . Calculating the steady-state flow statistics therefore requires simulating the DL-LES equations over a large number of flow times through the domain. It is a non-trivial question whether an unsteady PDE model with a functional form defined by a deep neural network can remain stable and accurate on $t \in [0, \infty )$ , especially when trained over comparatively short time intervals. Our results demonstrate that the DL-LES models are accurate and stable over long time horizons, which enables the estimation of the steady-state mean velocity, fluctuations and drag coefficient of turbulent flows around bluff bodies relevant to aerodynamics applications.

Mechanics↗

DOE OSTI2020

High Order Implicit Residual-Based Spatial Discretization Error Estimation for S N Neutron Transport

This work demonstrates our novel residual source spatial discretization error estimator (LeR/TEAD) for a DGFEM-1 discretization and assesses it along with two contemporary estimators, Ragusa and Wang's h -refinement estimator (RW) and Duo, Azmy, and Zikatanov's explicit residual-based estimator (DAZ), on a suite of Method of Manufactured Solutions (MMS) 2D problems and three realistic problem geometries. LeR/TE-AD is attractive because it directly estimates the local error in the angular flux, as opposed to a mere indicator of the error's behavior, on the same mesh and method order as the original numerical solution, thus typically being less computationally intensive than a refinement-based method. On the MMS suite, LeR/TE-AD consistently displayed a reduced performance versus its DGFEM-0 results in terms of accuracy and precision metrics, though it was not typically grossly inaccurate. This is attributed to the irregularities in the true solution across singular characteristics limiting the local accuracy of the numerical flux solution, leading to poor derivative approximations used in the residual approximations. The error transport problem then spreads the error in the residual to nearby cells, causing a greater degree of imprecision that did not afflict DAZ or RW. In testing the estimators on realistic problem geometries, however, LeR/TE-AD fared better. In practice, the true error is much larger in non-idealized geometries like in MMS, and a superlinear true solution means that RW and DAZ are not beneficially biased for DGFEM-1 error estimation. LeR/TE-AD was typically first or second in accuracy, primarily competing with RW, but the latter usually consumed 2-4 times the computational time as LeR/TE-AD, and requires a solution with four times as many unknowns. Furthermore, RW and LeR/TE-AD can be used to compute direct estimates of the error in any quantity of interest that is based on the angular ux solution, such as the fission rate density in a fuel pin, whereas DAZ requires a heuristic extension due to its norm-based nature.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2019

Separability of Mesh Bias and Parametric Uncertainty for a Full System Thermal Analysis

When making computational simulation predictions of multiphysics engineering systems, sources of uncertainty in the prediction need to be acknowledged and included in the analysis within the current paradigm of striving for simulation credibility. A thermal analysis of an aerospace geometry was performed at Sandia National Laboratories. Here, for this analysis, a verification, validation, and uncertainty quantification (VVUQ) workflow provided structure for the analysis, resulting in the quantification of significant uncertainty sources including spatial numerical error and material property parametric uncertainty. It was hypothesized that the parametric uncertainty and numerical errors were independent and separable for this application. This hypothesis was supported by performing uncertainty quantification (UQ) simulations at multiple mesh resolutions, while being limited by resources to minimize the number of medium and high resolution simulations. In conclusion, based on this supported hypothesis, a prediction including parametric uncertainty and a systematic mesh bias is used to make a margin assessment that avoids unnecessary uncertainty obscuring the results and optimizes use of computing resources.

42 ENGINEERING↗

DOE OSTI2018

Separability of Mesh Bias and Parametric Uncertainty for a Full System Thermal Analysis

When making computational simulation predictions of multi-physics engineering systems, sources of uncertainty in the prediction need to be acknowledged and included in the analysis within the current paradigm of striving for simulation credibility. A thermal analysis of an aerospace geometry was performed at Sandia National Laboratories. For this analysis a verification, validation and uncertainty quantification workflow provided structure for the analysis, resulting in the quantification of significant uncertainty sources including spatial numerical error and material property parametric uncertainty. It was hypothesized that the parametric uncertainty and numerical errors were independent and separable for this application. This hypothesis was supported by performing uncertainty quantification simulations at multiple mesh resolutions, while being limited by resources to minimize the number of medium and high resolution simulations. Based on this supported hypothesis, a prediction including parametric uncertainty and a systematic mesh bias are used to make a margin assessment that avoids unnecessary uncertainty obscuring the results and optimizes computing resources.

42 ENGINEERING↗

DOE OSTI2022

Time-explicit Darwin PIC algorithm

A new approach to Darwin particle-in-cell plasma simulation is described. Using a finite-element approach, the vector potential is assured to be exactly solenoidal. This allows writing the system of multi-species particles as an action-at-a-distance Hamiltonian system. Applying recently-developed explicit symplectic methods to this gives a time-explicit algorithm with all desired properties. Elliptic systems for the electrostatic and magneto-static portions of the problem are inverted efficiently by an algebraic multi-grid algorithm. The algorithm is implemented in a two-dimensional Cartesian-geometry code and tested by application to Weibel instability and to Alfvén wave dynamics. Results show the effectiveness of this approach. Extensions to arbitrary meshes and to a partially time-implicit scheme are briefly discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2014

SU-E-T-565: RAdiation Resistance of Cancer CElls Using GEANT4 DNA: RACE

Purpose: The objective of the RACE project is to develop a comparison between Monte Carlo simulation using the Geant4-DNA toolkit and measurements of radiation damage on 3D melanoma and chondrosarcoma culture cells coupled with gadolinium nanoparticles. We currently expose the status of the developments regarding simulations. Methods: Monte Carlo studies are driven using the Geant4 toolkit and the Geant4-DNA extension. In order to model the geometry of a cell population, the opensource CPOP++ program is being developed for the geometrical representation of 3D cell populations including a specific cell mesh coupled with a multi-agent system. Each cell includes cytoplasm and nucleus. The correct modeling of the cell population has been validated with confocal microscopy images of spheroids. The Geant4 Livermore physics models are used to simulate the interactions of a 250 keV X-ray beam and the production of secondaries from gadolinium nanoparticles supposed to be fixed on the cell membranes. Geant4-DNA processes are used to simulate the interactions of charged particles with the cells. An atomistic description of the DNA molecule, from PDB (Protein Data Bank) files, is provided by the so-called PDB4DNA Geant4 user application we developed to score energy depositions in DNA base pairs and sugar-phosphate groups. Results: At the microscopic level, our simulations enable assessing microscopic energy distribution in each cell compartment of a realistic 3D cell population. Dose enhancement factors due to the presence of gadolinium nanoparticles can be estimated. At the nanometer scale, direct damages on nuclear DNA are also estimated. Conclusion: We successfully simulated the impact of direct radiations on a realistic 3D cell population model compatible with microdosimetry calculations using the Geant4-DNA toolkit. Upcoming validation and the future integration of the radiochemistry module of Geant4-DNA will propose to correlate clusters of ionizations with in vitro experiments. All those developments will be released publicly. This work was supported by grants from Plan Cancer 2009-2013 French national initiative managed by INSERM (Institut National de la Sante et de la Recherche Medicale)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2014

TU-C-17A-03: An Integrated Contour Evaluation Software Tool Using Supervised Pattern Recognition for Radiotherapy

Purpose: Radiotherapy (RT) contours delineated either manually or semiautomatically require verification before clinical usage. Manual evaluation is very time consuming. A new integrated software tool using supervised pattern contour recognition was thus developed to facilitate this process. Methods: The contouring tool was developed using an object-oriented programming language C# and application programming interfaces, e.g. visualization toolkit (VTK). The C# language served as the tool design basis. The Accord.Net scientific computing libraries were utilized for the required statistical data processing and pattern recognition, while the VTK was used to build and render 3-D mesh models from critical RT structures in real-time and 360° visualization. Principal component analysis (PCA) was used for system self-updating geometry variations of normal structures based on physician-approved RT contours as a training dataset. The inhouse design of supervised PCA-based contour recognition method was used for automatically evaluating contour normality/abnormality. The function for reporting the contour evaluation results was implemented by using C# and Windows Form Designer. Results: The software input was RT simulation images and RT structures from commercial clinical treatment planning systems. Several abilities were demonstrated: automatic assessment of RT contours, file loading/saving of various modality medical images and RT contours, and generation/visualization of 3-D images and anatomical models. Moreover, it supported the 360° rendering of the RT structures in a multi-slice view, which allows physicians to visually check and edit abnormally contoured structures. Conclusion: This new software integrates the supervised learning framework with image processing and graphical visualization modules for RT contour verification. This tool has great potential for facilitating treatment planning with the assistance of an automatic contour evaluation module in avoiding unnecessary manual verification for physicians/dosimetrists. In addition, its nature as a compact and stand-alone tool allows for future extensibility to include additional functions for physicians’ clinical needs.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2025

Q3 Report for FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment

This report describes the work and activities carried out towards the completion of each of the following milestones in FY25 Q3: 1. Demonstrate workflow for generating self-consistent CESOL plasma profiles + first wall and divertor loading prediction and generate the CAT plasma and neutron loading needed for further engineering analysis. • Benchmark between two first wall heat flux mapping methods, identify importance of various heat flux sources and physics impact of using fully coupled CESOL vs post-analysis evaluation. 2. Generate medium fidelity parametrized CAD. • Generate parametrized CAD components for the CAT example case via either user-defined modules called within the geometry generation or by defeatured/parametrized CAD, including DCLL blanket matched to divertor boundary and magnets. Define materials, labels, and boundary conditions for passing the mesh to CFD tools. 3. Demonstrate multiphysics magnet analysis. • Demonstrate magnet analysis workflow called from the FREDA workflow, and 4. Demonstrate nuclear analysis. • Add model to OpenFOAM and/or other codes possibly including Diablo to account for tritium diffusion in solids.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1997

RIACS

Topics considered include: high-performance computing; cognitive and perceptual prostheses (computational aids designed to leverage human abilities); autonomous systems. Also included: development of a 3D unstructured grid code based on a finite volume formulation and applied to the Navier-stokes equations; Cartesian grid methods for complex geometry; multigrid methods for solving elliptic problems on unstructured grids; algebraic non-overlapping domain decomposition methods for compressible fluid flow problems on unstructured meshes; numerical methods for the compressible navier-stokes equations with application to aerodynamic flows; research in aerodynamic shape optimization; S-HARP: a parallel dynamic spectral partitioner; numerical schemes for the Hamilton-Jacobi and level set equations on triangulated domains; application of high-order shock capturing schemes to direct simulation of turbulence; multicast technology; network testbeds; supercomputer consolidation project.

Oliger, Joseph↗

NASA NTRS1985

Embedded mesh solutions of the Euler equation using a multiple-grid method

New developments make it now possible to obtain solutions to Euler equations for many problems on the basis of acceptable computing times. The possibility of an extension of these methods to more difficult flows and increasingly complex geometries can depend on obtaining suitable grid structures. Thus, an attractive approach to complex geometries involves the recasting of the problem into the framework of a multiple-grid structure. The discrete equations represented on the grid are solved as a simultaneous system. In the present investigation, this is accomplished with the aid of Ni's multiple-grid algorithm. The multiple-grid structure provided in this study represents the first step towards a completely general modular approach to the solution of the complete compressible flow problem. In the presented cases, a fine mesh resolution is obtained within the embedded mesh regions, while preserving the global coarse mesh convergence rates.

Usab, W. J., Jr.↗

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↗

DOE OSTI2022

Neutron diffusion calculation in heterogeneous geometry based on local/global iteration using proper orthogonal decomposition

This study newly proposes a heterogeneous core calculation method based on local/global iteration using proper orthogonal decomposition (POD). By using the singular value decomposition (SVD) and the low-rank approximation, appropriate POD bases for expanding the neutron flux can be obtained from snapshot data of the neutron flux obtained by fine mesh calculations. By projection using the POD bases, the dimension of the target equation (e.g., discretized neutron diffusion equation) can be dramatically reduced. In the proposed method, POD is effectively applied to each single assembly calculation (local calculation). Furthermore, using the local/global iteration, the effective neutron multiplication factor and the neutron flux distribution in the whole core geometry can be obtained by combining the numerical results of the local calculation for each fuel assembly and the global calculation for the whole core. As a feasibility study, the proposed method is applied to a one-dimensional heterogeneous core analysis, and the accuracy is investigated by changing the total number of POD bases. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

An unstructured mesh based neutronics optimization workflow

We have developed a fully automated workflow to optimize the neutronics performance of the Second Target Station (STS) at the Oak Ridge National Laboratory’s Spallation Neutron Source. The optimization workflow starts with the parametrized solid CAD engineering models and converts them into the unstructured mesh (UM) models for the neutronics calculations with MCNP6.2. Calculations are executed and their results are loaded into the Dakota optimization toolkit. Dakota analyzes the results and proposes new geometry parameters for the next design iteration. The cycle repeats until the optimal parameters are found. The automated CAD to MCNP conversion, the use of high-fidelity UM models, and the use of modern optimizer are the key elements that advance the entire optimization workflow in comparison with the original workflow. The original workflow was based on a simplified constructive solid geometry (CSG) modeling with MCNPX, mcnp_pstudy tool, and an in-house optimizer. Herein to demonstrate the new workflow, we present a case of neutronics optimization of the moderator–reflector assembly (MRA). Apart from the MRA, the workflow can optimize other major STS components, such as the spallation target, neutron beamlines, radiation shielding, and various accelerator components. Importantly, the new workflow opens the door to the advanced multi-physics multi-parameter optimization and has the potential for use in other nuclear physics and accelerator applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

NASA NTRS2001

Thermal Analysis of a High-Speed Aircraft Wing Using p-Version Finite Elements

This paper presents the results of conceptual level thermal analyses of a High Speed Civil Transport (HSCT) wing using p-version finite elements. The work was motivated by a thermal analysis of a HSCT wing structure which showed the importance of radiation heat transfer throughout the structure. The analysis also showed that refining a traditional finite element mesh to accurately capture the temperature distribution on the internal structure led to very large meshes with unacceptably long execution times. Further study indicated using p-version finite elements might improve computation performance for this class of problem. Methods for determining internal radiation heat transfer were then developed and demonstrated on test problems representative of the geometry found in an aircraft wing structure. This paper presents the results of the application of these new methods to the analysis of a high speed aircraft wing. Results for both a wing box model as well as a full wing model are presented. 'Me reduced wing box model allows for a comparison of the traditional finite element method with mesh refinement (h-refinement) to the new p-version finite elements while the full wing model demonstrates the applicability and efficiency of p-version finite elements for large models.

Gould, Dana C.↗