Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mesh data”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,117 records · Page 62

Assessment of Coupling One-Dimension and Three-Dimension CFD in Conjugate Heat Transfer Modeling and Simulations for Cooling of Turbine Engine

The CHT (conjugate heat transfer) modeling and simulation work fits under a bigger area in the TTT (Transformational Tools and Technologies) project called “Cooling Tools and Technologies for Combustor and Turbine.” The existing validation data sets such as TTT GRC (Glenn Research Center) THX (Turbulent Heat FLUX) IV experiments will be used to assess the solid/fluid heat transfer development and validation. The experimental configuration was a subsonic convergent nozzle with a square exit blowing a heated exhaust over an effusion-cooled deck. The simulation solutions were obtained and compared to experimental measurements of the profiles of flow velocities and temperatures, as well as deck surface temperatures. Cooling air blowing ratios of 0, 1, and 2 were considered. The three-dimensional cooling tubes were modelled by the medial axis of the cylindrical geometries due to large length over width ratios. It represented a typical coupling of one-dimension and three-dimension simulations of conjugate heat transfer approaches. The surface temperatures on the solid surfaces and the profiles of velocities and temperatures near the surfaces compared to the experimental data were presented.

nonlinear sub-filter stresses↗

SPRUCE Soil Volumetric Water Content in Experimental Plots, Marcell Experimental Forest, Minnesota, 2017-2025

This dataset contains soil volumetric water content (VWC) measurements from the Spruce and Peatland Responses Under Changing Environments (SPRUCE) Experiment. Measurements were made inside SPRUCE experimental plots in the S1 Bog at the US Forest Service Marcell Experimental Forest in northern Minnesota, USA from 2018-2025 (2018-06-14 to 2025-12-31). Observations were made with METER 10HS soil moisture sensors. To reduce bulk density related variability, 10HS sensors were placed inside mesh tubes filled with peat at a standard bulk density. The observations under standard bulk density represent relative differences in water content between hummock and hollow positions, and in response to experimental treatments. Productivity of peatlands, and their keystone species sphagnum, are highly dependent on water availability, which is affected by lateral inputs, precipitation, ground water depth and evapotranspiration. Knowledge of near surface and sphagnum water content is useful to understand peatland function. This dataset contains 8 data files in comma-separate values (.csv) format. Additional metadata are provided: one data dictionary and a file-level metadata file in comma-separate values (*.csv) format and a user guide in PDF (*.pdf) format.

Warren, Jeffrey [ORNL] (ORCID:0000000206804697)↗

Rate of Para-Orthohydrogen Conversion in Cryogenic Vortex Tube

Raising efficiency of the cooling process for cryogenic hydrogen and minimizing hydrogen boil-off during storage and transportation of liquid hydrogen are critically important factors for widening implementation of hydrogen-based clean energy systems. One novel approach involves utilization of relatively simple catalyzed vortex tubes, where hydrogen in the primarily para- nuclear spin isomer form is converted into the ortho- form while consuming a significant amount of heat. To design such systems effectively, better understanding is required about catalyst-assisted para-orthohydrogen conversion rates in high-speed vortical flows. In this study, a computational fluid dynamics simulation has been set up to model an experimental system with cryogenic hydrogen vortex tubes. Mesh-verification and a validation study has been conducted first for a non-catalyzed tube. Then, the rate coefficient of the para-orthohydrogen conversion of cryogenic hydrogen has been determined by matching numerical results with experimental data available for a catalyzed vortex tube. Furthermore, the presented information can help design and optimize cryogenic hydrogen cooling within vortex tubes.

08 HYDROGEN↗

Graph Learning in Physical-informed Mesh-reduced Space for Real-world Dynamic Systems

This Git repository contains codes for the 'Graph Learning in Physical-informed Mesh-reduced Space for Real-world Dynamic Systems' paper that will be published in 2023 SIGKDD. This work uses physical-informed prior (PiP) information to learn and predict fluid dynamics in a reduced mesh space. We propose a two-stage graph-based model for fluid velocity field reconstruction and prediction. In the first stage, we learn a subgraph autoencoder to summarize the information in a mesh-reduced space using physical-informed priors. In the second stage, we learn a dynamics predictor to predict subgraph evolution. We demonstrate the effectiveness of our model on two fluid flow datasets: lid-driven cavity flow data and cylinder flow data. This code is also applicable for any other dynamic systems with corresponding data.

Lei, Bo↗

Simulation of Internal Combustion Engines with High-Performance Computing Tools

Traditional Lagrangian spray modeling approaches for internal combustion engines are highly griddependent due to insufficient resolution in the near nozzle region. This is primarily because of inherent restrictions of volume fraction with the Lagrangian assumption together with high computational costs associated with small grid sizes. A state-of-the-art grid-convergent spray modeling approach was developed and implemented by Senecal et al. (ASME-ICEF2012-92043) in the CONVERGE software. The key features of the methodology include Adaptive Mesh Refinement (AMR), advanced liquid-gas momentum coupling, and improved distribution of the liquid phase, which enables use of cell sizes smaller than the nozzle diameter. This modeling approach was rigorously validated against nonevaporating, evaporating, and reacting data from the literature. The current numerical study focusses on further demonstration of the grid-convergent modeling approach for simulating a single-cylinder Cat® compression ignition engine. The simulated injector is characterized with a nominal nozzle exit diameter of 259 μm. Simulations using various minimum grid sizes (ranging from 125 μm to 1000 μm) are compared for engine performance and emissions parameters of interest such as pressure, heat release rate, ignition delay, NOx, HC, and soot emissions. The peak cell-count for the highest resolution simulation was on the order of 34 million. These computationally expensive simulations were facilitated at a high-performance computing facility at Argonne National Laboratory. METIS load-balancing algorithm was developed and implemented in Converge code for the simulations. Scaling studies were also performed. The validity of previously recommended grid settings (ASME-ICEF2012-92043) for accuracy/runtime trade-off is further assessed. Efficacy of a simplified combustion model is also compared against a detailed chemical kinetics based combustion model.

33 ADVANCED PROPULSION SYSTEMS↗

Trinity Phase 2 Open Science: CTH

CTH is an Eulerian hydrocode developed by Sandia National Laboratories (SNL) to solve a wide range of shock wave propagation and material deformation problems. Adaptive mesh refinement is also used to improve efficiency for problems with a wide range of spatial scales. The code has a history of running on a variety of computing platforms ranging from desktops to massively parallel distributed-data systems. For the Trinity Phase 2 Open Science campaign, CTH was used to study mesoscale simulations of the hypervelocity penetration of granular SiC powders. The simulations were compared to experimental data. A scaling study of CTH up to 8192 KNL nodes was also performed, and several improvements were made to the code to improve the scalability.

97 MATHEMATICS AND COMPUTING↗

Developing Light Collection Enhancements and Wire Tensioning Methods for LArTPC Neutrino Detectors

Liquid argon Time Projection Chambers (LArTPCs) are becoming widely used as neutrino detectors because of their image-like event reconstruction which enables precision neutrino measurements. They primarily use ionisation charge to reconstruct neutrino events. It has been shown, however, that the scintillation light emitted by liquid argon could be exploited to improve their performance. As the neutrino measurements planned in the near future require large-scale experiments, their construction presents challenges in terms of both charge and light collection. In this dissertation we present solutions developed to improve the performance in both aspects of these detectors. We present a new wire tensioning measurement method that allows a remote measurement of the tension of the large number wires that constitute the TPC anode. We also discuss the development and installation of WLS-compound covered foils for the SBND neutrino detector at Fermilab, which is a technique proposed t o augment light collection in LArTPCs. This included preparing a SBND-like mesh cathode and testing it in the Run III of LArIAT, a test beam detector also located at Fermilab. Finally, we present a study aimed at understanding late scintillation light emitted by recombining positive argon ions using LArIAT data, which could affect large scale surface detectors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Validation Exercise of a Coarse Finite Element Model of Laser Welds

The objective of this project is to validate low-fidelity models of 304L to 304L stainless steel partial-penetration laser welds for thin sheets. Low-fidelity means that the weld is represented by coarsely meshed element blocks. Here, the hexahedral element size is approx imately half the weld penetration depth. The material behavior of the block is represented by a J2 plasticity model with a Voce hardening function. The source of the data used in this work is an extensive experimental study conducted by Sharlotte Kramer (1528) and published in 2015. Figure 1 shows a cross-section of the weld of interest. The nominal thickness of the sheets is 0.063 in. while the target penetration depth of the weld is in the range of 0.028 to 0.032 in., extending about half the sheet thickness. Uniaxial tension tests provided data for calibration of base material and weld models. Results of two validation geometries were also provided. The principal validation geometry is shown in Fig. 2. It consists of a plate specimen with in-plane dimensions 6 in × 2.875 in loaded in tension. A circular plug with a 1.5 in. diameter was cut from the center of the plate and then welded in place. The details of the welding schedule are given. An important assumption is that the welds in the calibration and validation specimens have similar geometric and material properties as those in the validation tests. The task was to first calibrate models for the base material and the welds and then simulate the validation tests until the point of weld first failure.

36 MATERIALS SCIENCE↗

PASS Spacecraft Antenna Technology Assessment

The purpose was to generate estimates of mechanical performance for the classes of spacecraft antenna under construction for application to the Personal Access Satellite System (PASS). These performance data are needed for the support of trade studies involving antenna system development. The classes of antenna considered included: (1) rigid non-deployable antenna structures; (2) mechanical deployable antenna concepts; (3) inflatable deployable antenna concepts; and (4) mesh deployable antenna concepts. The estimates of mechanical performance are presented in terms of structural weight and cost as a function of the reflector size. Estimates of aperture surface precision are presented for a few discrete antenna sizes. The range of reflector size is 1 to 4 meters for non-deployable structures and 2 to 8 meters for deployable structures. The range of reflector surface precision is lambda/30 to lambda/50 for 20 and 30 GHz, respectively.

Freeland, R. E.↗

Experience with parametric binary dissection

Parametric Binary Dissection (PBD) is a new algorithm that can be used for partitioning graphs embedded in 2- or 3-dimensional space. It partitions explicitly on the basis of nodes + (lambda)x(edges cut), where lambda is the ratio of time to communicate over an edge to the time to compute at a node. The new algorithm is faster than the original binary dissection algorithm and attempts to obtain better partitions than the older algorithm, which only takes nodes into account. The performance of parametric dissection with plain binary dissection on 3 large unstructured 3-d meshes obtained from computational fluid dynamics and on 2 random graphs were compared. It was showm that the new algorithm can usually yield partitions that are substantially superior, but that its performance is heavily dependent on the input data.

Bokhari, Shahid H.↗

A New Method for Reduction of Photomultiplier Signal-Induced Noise

For lidar measurements of ozone, photomultiplier tube (PMT) detector signal-induced noise represents a fundamental problem that complicates the extraction of information from lidar data. A new method is developed to significantly reduce signal-induced noise in lidar receiver PMT detectors. The electron optics of the lidar photomultiplier detector is modified to filter the source of signal-induced noise. A mesh electrode external to the PMT is utilized to control photoemission and disorient electron trajectories from the photocathode to the first dynode. Experiments were taken both with simulated and actual lidar return signals at Langley Research Center. Results show at least 40 percent more accurate ozone number density values with a mesh voltage of 60 V applied than with no voltage applied.

Koble, Andrea↗

Aeroheating Predictions for a Hypersonic, Turbulent Near-Wake

The accuracy of heating predictions using various turbulence models is examined for an axisymmetric near-wake at Mach 6. The CFD predictions are compared with experimental data collected under AGARD Working Group 18 on the wake of a 70-degree sphere-cone. The impact of grid resolution and discretization error is estimated, which allows a comparison of stacked-block and conventional structured meshes. The accuracy of steady Reynolds-averaged Navier-Stokes (RANS) models is contrasted with that of a hybrid RANS/Large-Eddy Simulation model. The predictions are made with three different CFD codes (LAURA, FUN3D, and HyperSolve), to demonstrate the code-to-code variation in the results. Steady SST models capture the qualitative nature of the heating in the wake, including the increasing peak heating with increasing Reynolds number. Spalart-Allmaras models, including SA-Catris, under-predicted the peak heating in the wake. Hybrid RANS/LES models improve upon the SA results but have their own modeling difficulties near the shear layer impingement. These results are generally consistent across solvers and grid topologies.

RANS↗

Aeroheating Predictions for a Hypersonic, Turbulent Near-Wake

The accuracy of heating predictions using various turbulence models is examined for an axisymmetric near-wake at Mach 6. The CFD predictions are compared with experimental data collected under AGARD Working Group 18 on the wake of a 70-degree sphere-cone. The impact of grid resolution and discretization error is estimated, which allows a comparison of stacked-block and conventional structured meshes. The accuracy of steady Reynolds-averaged Navier-Stokes (RANS) models is contrasted with that of a hybrid RANS/Large-Eddy Simulation model. The predictions are made with three different CFD codes (LAURA, FUN3D, and HyperSolve), to demonstrate the code-to-code variation in the results. Steady SST models capture the qualitative nature of the heating in the wake, including the increasing peak heating with increasing Reynolds number. Spalart-Allmaras models, including SA-Catris, under-predicted the peak heating in the wake. Hybrid RANS/LES models improve upon the SA results but have their own modeling difficulties near the shear layer impingement. These results are generally consistent across solvers and grid topologies.

RANS↗

ISSM: Ice Sheet System Model

In order to have the capability to use satellite data from its own missions to inform future sea-level rise projections, JPL needed a full-fledged ice-sheet/iceshelf flow model, capable of modeling the mass balance of Antarctica and Greenland into the near future. ISSM was developed with such a goal in mind, as a massively parallelized, multi-purpose finite-element framework dedicated to ice-sheet modeling. ISSM features unstructured meshes (Tria in 2D, and Penta in 3D) along with corresponding finite elements for both types of meshes. Each finite element can carry out diagnostic, prognostic, transient, thermal 3D, surface, and bed slope simulations. Anisotropic meshing enables adaptation of meshes to a certain metric, and the 2D Shelfy-Stream, 3D Blatter/Pattyn, and 3D Full-Stokes formulations capture the bulk of the ice-flow physics. These elements can be coupled together, based on the Arlequin method, so that on a large scale model such as Antarctica, each type of finite element is used in the most efficient manner. For each finite element referenced above, ISSM implements an adjoint. This adjoint can be used to carry out model inversions of unknown model parameters, typically ice rheology and basal drag at the ice/bedrock interface, using a metric such as the observed InSAR surface velocity. This data assimilation capability is crucial to allow spinning up of ice flow models using available satellite data. ISSM relies on the PETSc library for its vectors, matrices, and solvers. This allows ISSM to run efficiently on any parallel platform, whether shared or distrib- ISSM: Ice Sheet System Model NASA's Jet Propulsion Laboratory, Pasadena, California uted. It can run on the largest clusters, and is fully scalable. This allows ISSM to tackle models the size of continents. ISSM is embedded into MATLAB and Python, both open scientific platforms. This improves its outreach within the science community. It is entirely written in C/C++, which gives it flexibility in its design, and the power/speed that C/C++ allows. ISSM is svn (subversion) hosted, on a JPL repository, to facilitate its development and maintenance. ISSM can also model propagation of rifts using contact mechanics and mesh splitting, and can interface to the Dakota software. To carry out sensitivity analysis, mesh partitioning algorithms are available, based on the Scotch, Chaco, and Metis partitioners that ensure equal area mesh partitions can be done, which are then usable for sampling and local reliability methods.

Larour, Eric↗

Study of the Correlation between the Performances of Lunar Vehicle Wheels Predicted by the Nepean Wheeled Vehicle Performance Model and Test Data

This paper describes the results of a study of the correlation between the performances of wheels for lunar vehicles predicted using the Nepean wheeled vehicle performance model (NWVPM), developed under the auspices of Vehicle Systems Development Corporation, Ottawa, Canada, and the corresponding test data presented in Performance evaluation of wheels for lunar vehicles , Technical Report M-70-2, prepared for George C. Marshall Space Flight Center, National Aeronautics and Space Administration (NASA), USA, by the US Army Engineer Waterways Experiment Station (WES). The NWVPM was originally developed for design and performance evaluation of terrestrial off-road wheeled vehicles. The purpose of this study is to assess the potential of the NWVPM for evaluating wheel candidates for the new generation of extra-terrestrial vehicles. Two versions of a wire-mesh wheel and a hoop-spring wheel, which were considered as candidates for lunar roving vehicles for the NASA Apollo program in the late 1960s, together with a pneumatic wheel were examined in this study. The tractive performances of these wheels and of a 464 test vehicle with the pneumatic wheels on air-dry sand were predicted using the NWVPM and compared with the corresponding test data obtained under Earth s gravity and previously documented in the above-named report. While test data on wheel or vehicle performances obtained under Earth s gravity may not necessarily be representative of those on extra-terrestrial bodies, because of the differences in gravity and in environmental conditions, such as atmospheric pressure, it is still a valid approach to use test data obtained under Earth s gravity to evaluate the predictive capability of the NWVPM and its potential applications to predicting wheel or wheeled rover performances on extra-terrestrial bodies. Results of this study show that, using the ratio (P20/W) of the drawbar pull to normal load at 20 per cent slip as a performance indicator, there is a reasonable correlation between the predictions and experimental data. This indicates that the NWVPM has the potential as an engineering tool for evaluating wheel candidates for a future generation of extra-terrestrial vehicles, provided that appropriate input data are available.

Wong, J. Y.↗

Nonlinear elasticity with the Shifted Boundary Method

Here, we propose a new unfitted/immersed computational framework for nonlinear solid mechanics, which bypasses the complexities associated with the generation of CAD representations and subsequent body-fitted meshing. This approach allows to speed up the cycle of design and analysis in complex geometry and requires relatively simple computer graphics representations of the surface geometries to be simulated, such as the Standard Tessellation Language (STL format). Complex data structures and integration on cut elements are avoided by means of an approximate boundary representation and a modification (shifting) of the boundary conditions to maintain optimal accuracy. An extensive set of computational experiments in two and three dimensions is included.

97 MATHEMATICS AND COMPUTING↗

Multilevel summation with B-spline interpolation for pairwise interactions in molecular dynamics simulations

The multilevel summation method for calculating electrostatic interactions in molecular dynamics simulations constructs an approximation to a pairwise interaction kernel and its gradient, which can be evaluated at a cost that scales linearly with the number of atoms. The method smoothly splits the kernel into a sum of partial kernels of increasing range and decreasing variability with the longer-range parts interpolated from grids of increasing coarseness. Multilevel summation is especially appropriate in the context of dynamics and minimization, because it can produce continuous gradients. This article explores the use of B-splines to increase the accuracy of the multilevel summation method (for nonperiodic boundaries) without incurring additional computation other than a preprocessing step (whose cost also scales linearly). To obtain accurate results efficiently involves technical difficulties, which are overcome by a novel preprocessing algorithm. Numerical experiments demonstrate that the resulting method offers substantial improvements in accuracy and that its performance is competitive with an implementation of the fast multipole method in general and markedly better for Hamiltonian formulations of molecular dynamics. The improvement is great enough to establish multilevel summation as a serious contender for calculating pairwise interactions in molecular dynamics simulations. In particular, the method appears to be uniquely capable for molecular dynamics in two situations, nonperiodic boundary conditions and massively parallel computation, where the fast Fourier transform employed in the particle–mesh Ewald method falls short.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development of a Method for Shape Optimization for a Gas Turbine Fuel Injector Design Using Metal-Additive Manufacturing

Adjoint shape optimization has enabled physics-based optimal designs for aerodynamic surfaces. Additive manufacturing (AM) makes it possible to manufacture complex shapes. However, there has been a gap between optimal and manufacturable surfaces due to the inherent limitations of commercial computational fluid dynamics (CFD) codes to implement geometric constraints during adjoint computation. In such cases, the design sensitivities are exported and used to perform constrained shape modifications using parametric information stored in computer aided design (CAD) files to satisfy manufacturability constraints. However, modifying the design using adjoint methods in CFD solvers and performing constrained shape modification in CAD can lead to inconsistencies due to different shape parameterization schemes. This paper describes a method to enable the simultaneous optimization of the fluid domain and impose AM manufacturability constraints, resolving one of the key issues of geometry definition for isogeometric analysis. Similar to a grid convergence study, the proposed method verifies the consistencies between shape parameterization techniques present within commercial CAD and CFD software during mesh movement as a part of the adjoint shape optimization routine. By identifying the appropriate parameters essential to a shape optimization study, the error metric between the different parameterization techniques converges to demonstrate sufficient consistencies for justifiable exchange of data between CAD and CFD. For the identified shape optimization parameters, the error metric to measure the deviation between the two parameterization schemes lies within the AM laser-powder bed fusion (L-PBF) process tolerance. Additionally, comparison for subsequent objective function calculations between iterations of the optimization loop showed acceptable differences within 1% variation between the modified geometries obtained using the two parameterization schemes. This method provides justification for the use of multiphysics guided adjoint design sensitivities computed in CFD software to perform shape modifications in CAD to incorporate AM manufacturability constraints during the shape optimization loop such that optimal designs are also additively manufacturable.

33 ADVANCED PROPULSION SYSTEMS↗