Search NASA⌕ Search

SEARCH · Search NASA

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.

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

At least 235 records · Page 13

Numerical techniques for the solution of the compressible Navier-Stokes equations and implementation of turbulence models

The time-splitting explicit numerical method of MacCormack is applied to separated turbulent boundary layer flow problems. Modifications of this basic method are developed to counter difficulties associated with complicated geometry and severe numerical resolution requirements of turbulence model equations. The accuracy of solutions is investigated by comparison with exact solutions for several simple cases. Procedures are developed for modifying the basic method to improve the accuracy. Numerical solutions of high-Reynolds-number separated flows over an airfoil and shock-separated flows over a flat plate are obtained. A simple mixing length model of turbulence is used for the transonic flow past an airfoil. A nonorthogonal mesh of arbitrary configuration facilitates the description of the flow field. For the simpler geometry associated with the flat plate, a rectangular mesh is used, and solutions are obtained based on a two-equation differential model of turbulence.

Baldwin, B. S.↗

Parallel Adaptive High-Order CFD Simulations Characterizing Cavity Acoustics for the Complete SOFIA Aircraft

This paper presents one-of-a-kind MPI-parallel computational fluid dynamics simulations for the Stratospheric Observatory for Infrared Astronomy (SOFIA). SOFIA is an airborne, 2.5-meter infrared telescope mounted in an open cavity in the aft of a Boeing 747SP. These simulations focus on how the unsteady flow field inside and over the cavity interferes with the optical path and mounting of the telescope. A temporally fourth-order Runge-Kutta, and spatially fifth-order WENO-5Z scheme was used to perform implicit large eddy simulations. An immersed boundary method provides automated gridding for complex geometries and natural coupling to a block-structured Cartesian adaptive mesh refinement framework. Strong scaling studies using NASA's Pleiades supercomputer with up to 32,000 cores and 4 billion cells shows excellent scaling. Dynamic load balancing based on execution time on individual AMR blocks addresses irregularities caused by the highly complex geometry. Limits to scaling beyond 32K cores are identified, and targeted code optimizations are discussed.

Acoustics↗

Rapid Hypersonic Simulations using US3D and Pointwise

For hypersonic simulations, unstructured flow solvers typically have problems predicting surface heat fluxes when strong shocks are present. To address these issues, this paper outlines a workflow that applies best practices developed for structured grids to unstructured meshes. In addition, unstructured grid generation can significantly reduce the time required to create quality grids for complex geometries. Several examples are computed using DPLR, a structured grid flow solver, and US3D flow, an unstructured mesh solver. Results from the two codes are compared, and they show excellent agreement. Overall, the unstructured grid workflow offers a viable and attractive alternative for hypersonic simulations.

C. Tang↗

Rapid Hypersonic Simulations using US3D and Pointwise

For hypersonic simulations, unstructured flow solvers typically have problems predicting surface heat fluxes when strong shocks are present. To address these issues, this paper outlines a workflow that applies best practices developed for structured grids to unstructured meshes. In addition, unstructured grid generation can significantly reduce the time required to create quality grids for complex geometries. Several examples are computed using DPLR, a structured grid flow solver, and US3D flow, an unstructured mesh solver. Results from the two codes are compared, and they show excellent agreement. Overall, the unstructured grid workflow offers a viable and attractive alternative for hypersonic simulations.

Chun Tang↗

CAPRI (Computational Analysis PRogramming Interface): A Solid Modeling Based Infra-Structure for Engineering Analysis and Design Simulations

CAPRI is a CAD-vendor neutral application programming interface designed for the construction of analysis and design systems. By allowing access to the geometry from within all modules (grid generators, solvers and post-processors) such tasks as meshing on the actual surfaces, node enrichment by solvers and defining which mesh faces are boundaries (for the solver and visualization system) become simpler. The overall reliance on file 'standards' is minimized. This 'Geometry Centric' approach makes multi-physics (multi-disciplinary) analysis codes much easier to build. By using the shared (coupled) surface as the foundation, CAPRI provides a single call to interpolate grid-node based data from the surface discretization in one volume to another. Finally, design systems are possible where the results can be brought back into the CAD system (and therefore manufactured) because all geometry construction and modification are performed using the CAD system's geometry kernel.

Haimes, Robert↗

Data, model inputs, and analysis scripts associated with a manuscript on stream intermittency controls across spatial scales in Pacific Northwest watersheds

NOTE: The manuscript associated with this data package is currently in review. The data may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final dataset and additional metadata. This data package is associated with the manuscript "Hydroclimatic Memory and Watershed Template Shape Stream Intermittency: Multi-scale Attribution Using Process-based Simulation and Explainable ML" by Niroula et al. (2026), submitted to Water Resources Research (WRR). The study investigates the dominant controls on stream intermittency across local, reach, and watershed scales using a coupled process-based simulation and explainable machine-learning framework. Long-term daily simulations from the Advanced Terrestrial Simulator (ATS) were used to generate wetness states and ponded-depth responses over river-corridor cells. These ATS outputs were then aggregated across scales and used to train XGBoost (eXtreme Gradient Boosting) models. SHAP (SHapley Additive exPlanations) was applied to quantify the relative importance of hydroclimatic forcings, watershed template attributes, and antecedent-memory effects in shaping intermittency behavior. The analysis is carried out for three contrasting Pacific Northwest watersheds: Oak Creek (OCW), American River Watershed (ARW), and H.J. Andrews (HJA). Across these testbeds, the package contains ATS-ready watershed inputs, ATS run configuration and selected output files, model-evaluation data products, intermittency-analysis datasets, machine-learning target-feature tables, SHAP outputs, and notebooks used to organize, analyze, and visualize results. At a high level, the package documents a workflow in which ATS provides the physically based simulation backbone and explainable machine learning is used as a post-processing attribution tool. The contents are intended to support interpretation of the manuscript figures and results, provide context for how intermittency metrics were generated at multiple scales, and preserve the key artifacts needed to understand and reuse the analysis workflow. The package contains a high-level directory summary file (`summary.txt`) and four main content folders (1) `evaluation_plots` contains evaluation figures and supporting evaluation datasets; (2) `intermittency_plots` contains intermittency-focused analysis notebook and prepared datasets; (3) `ml-training-and-shap_values_plots` contains ML training inputs, SHAP outputs, and figure-generation notebooks; and (4) `watershed_mesh_and_ats_input` contains ATS model setup materials, forcing inputs, geometry, and selected run files. More specifically, the `evaluation_plots` folder contains the notebook used for ATS evaluation plotting and site-specific evaluation datasets. These include evapotranspiration and water-balance products for three watersheds, as well as an Oak Creek field-measurement discharge file. The `intermittency_plots` folder contains the notebook used for intermittency analysis and the prepared datasets used to analyze intermittent and non-intermittent wetness behavior across the study watersheds. The `ml-training-and-shap_values_plots` folder contains notebooks and outputs for the machine-learning and explainability workflow. This includes the main XGBoost and SHAP notebook(s), a beeswarm plotting notebook, target-feature tables for machine-learning training, SHAP summary tables, and per-sample SHAP value archives. The `watershed_mesh_and_ats_input` folder contains ATS-related watershed inputs and supporting materials. This includes mesh and shape products, ATS-readable LAI and meteorological forcing inputs, selected ATS spinup and transient-run files, and a watershed workflow example notebook. Subdirectories are organized by watershed where applicable.All files are .cpg (codepage files), .csv (comma-separated values), .dbf (database files), .exo (Exodus mesh format), .h5 (HDF5 format), .ipynb (Jupyter notebooks), .pkl (Python pickle), .prj (projection files), .sh (shell scripts), .shp (shapefile geometry), .shx (shapefile index), .txt (text files), or .xml (markup data).

Advanced Terrestrial Simulator↗

Computation of turbulent near wake for asymmetric airfoils

A numerical procedure for studying the turbulent near wake of two dimensional airfoil sections is presented. The Reynolds Navier-Stokes equations were written for flow about bodies of arbitrary geometry and solved on an arbitrary nonuniform curvilinear computational mesh. Eddy viscosity and Reynolds stress turbulence transport models are considered. Specific examples are shown for airfoil section by using an algebraic viscosity model with streamwise relaxation and the interactive Reynolds stress model.

Deiwert, G. S.↗

Computational Issues Associated with Temporally Deforming Geometries Such as Thrust Vectoring Nozzles

During the past decade, computational simulation of fluid flow around complex configurations has progressed significantly and many notable successes have been reported, however, unsteady time-dependent solutions are not easily obtainable. The present effort involves unsteady time dependent simulation of temporally deforming geometries. Grid generation for a complex configuration can be a time consuming process and temporally varying geometries necessitate the regeneration of such grids for every time step. Traditional grid generation techniques have been tried and demonstrated to be inadequate to such simulations. Non-Uniform Rational B-splines (NURBS) based techniques provide a compact and accurate representation of the geometry. This definition can be coupled with a distribution mesh for a user defined spacing. The present method greatly reduces cpu requirements for time dependent remeshing, facilitating the simulation of more complex unsteady problems. A thrust vectoring nozzle has been chosen to demonstrate the capability as it is of current interest in the aerospace industry for better maneuverability of fighter aircraft in close combat and in post stall regimes. This current effort is the first step towards multidisciplinary design optimization which involves coupling the aerodynamic heat transfer and structural analysis techniques. Applications include simulation of temporally deforming bodies and aeroelastic problems.

Boyalakuntla, Kishore↗

Computations of Aerodynamic Performance Databases Using Output-Based Refinement

Objectives: Handle complex geometry problems; Control discretization errors via solution-adaptive mesh refinement; Focus on aerodynamic databases of parametric and optimization studies: 1. Accuracy: satisfy prescribed error bounds 2. Robustness and speed: may require over 105 mesh generations 3. Automation: avoid user supervision Obtain "expert meshes" independent of user skill; and Run every case adaptively in production settings.

Nemec, Marian↗

Handbook on Face Gear Drives with a Spur Involute Pinion

The use of face gears in power transmission and drive systems has a significant number of benefits. Face gears allow a variety of new transmission arrangements as well as high reduction ratio capability. This leads to drive system weight reduction and improvements in performance. In this work, basic information about the design and analysis of face gear drives is presented. The work considers face gears in mesh with spur involute pinions for both intersecting axes and offset drives. Tooth geometry, kinematics, generation of face gears with localized bearing contact by cutting and grinding, avoidance of tooth undercutting, avoidance of tooth pointing, tooth contact analysis, and algorithms for the simulation of meshing and contact arc all topics which are discussed. In addition, applications of face gear drives are presented. Included are design uses in aerospace applications such as helicopter transmissions, split-torque face gear arrangements, comparisons of face gears with bevel gears, and general design considerations.

Litvin, F. L.↗

FORTRAN program for generating a two-dimensional orthogonal mesh between two arbitrary boundaries

A FORTRAN 4 program is presented which computes and plots coordinates for a two-dimensional orthogonal mesh in the region between the walls of a flow channel. The program is designed for a channel containing a body about which flow passes and which spans the channel from one wall to the other. However, the condition that the channel contain an immersed body can be easily removed from the program. Input to the program consists of spline points of the channel walls and the body geometry. Output includes printed and plotted coordinates of the generated orthogonal mesh and angles of the mesh with the horizontal plane.

Mcnally, W. D.↗

Modeling Geometry and Progressive Failure of Material Interfaces in Plain Weave Composites

A procedure combining a geometrically nonlinear, explicit-dynamics contact analysis, computer aided design techniques, and elasticity-based mesh adjustment is proposed to efficiently generate realistic finite element models for meso-mechanical analysis of progressive failure in textile composites. In the procedure, the geometry of fiber tows is obtained by imposing a fictitious expansion on the tows. Meshes resulting from the procedure are conformal with the computed tow-tow and tow-matrix interfaces but are incongruent at the interfaces. The mesh interfaces are treated as cohesive contact surfaces not only to resolve the incongruence but also to simulate progressive failure. The method is employed to simulate debonding at the material interfaces in a ceramic-matrix plain weave composite with matrix porosity and in a polymeric matrix plain weave composite without matrix porosity, both subject to uniaxial cyclic loading. The numerical results indicate progression of the interfacial damage during every loading and reverse loading event in a constant strain amplitude cyclic process. However, the composites show different patterns of damage advancement.

Hsu, Su-Yuen↗

An Overset Grid Motion Driver for Resolving Entry Vehicle CFD Simulations with Large Changes in Wake Orientation

Time-accurate computational fluid dynamics (CFD) simulations using overset meshes can enable reduced computational cost compared to a rigid, monolithic mesh for blunt-body atmospheric entry vehicles. If the vehicle has a significant change in attitude, the high-density region of the overset mesh designed for resolving the subsonic wake can be kept in position as the wake orientation changes relative to the body. Current simulations generally keep the wake mesh fixed in attitude while the body rotates, assuming the wake orientation remains relatively constant. However, changes in the relative velocity over a long-duration simulation or effects of the body geometry can invalidate this assumption. This work presents a method for calculating the desired mesh orientation to resolve a blunt body wake using the atmosphere-relative velocity of the vehicle. A motion driver using these calculations is implemented in a CFD-in-the-loop flight model bridging the FUN3D flow solver and POST2 trajectory propagator. Overset meshing is performed using the Yoga extension for FUN3D. Results are presented for simulations with both overset and monolithic meshes under large changes in relative velocity (and therefore wake orientation). The results demonstrate that the motion driver successfully orients an overset wake mesh to capture the subsonic wake with sufficient resolution despite variable body rotation and body-relative freestream velocity.

CFD↗

An Overset Grid Motion Driver for Resolving Entry Vehicle CFD Simulations with Large Changes in Wake Orientation

Time-accurate computational fluid dynamics (CFD) simulations using overset meshes can enable reduced computational cost compared to a rigid, monolithic mesh for blunt-body atmospheric entry vehicles. If the vehicle has a significant change in attitude, the high-density region of the overset mesh designed for resolving the subsonic wake can be kept in position as the wake orientation changes relative to the body. Current simulations generally keep the wake mesh fixed in attitude while the body rotates, assuming the wake orientation remains relatively constant. However, changes in the relative velocity over a long-duration simulation or effects of the body geometry can invalidate this assumption. This work presents a method for calculating the desired mesh orientation to resolve a blunt body wake using the atmosphere-relative velocity of the vehicle. A motion driver using these calculations is implemented in a CFD-in-the-loop flight model bridging the FUN3D flow solver and POST2 trajectory propagator. Overset meshing is performed using the Yoga extension for FUN3D. Results are presented for simulations with both overset and monolithic meshes under large changes in relative velocity (and therefore wake orientation). The results demonstrate that the motion driver successfully orients an overset wake mesh to capture the subsonic wake with sufficient resolution despite variable body rotation and body-relative freestream velocity.

Computational Fluid Dynamics↗

Investigation of advanced counterrotation blade configuration concepts for high speed turboprop systems. Task 4: Advanced fan section aerodynamic analysis computer program user's manual

The computer program user's manual for the ADPACAPES (Advanced Ducted Propfan Analysis Code-Average Passage Engine Simulation) program is included. The objective of the computer program is development of a three-dimensional Euler/Navier-Stokes flow analysis for fan section/engine geometries containing multiple blade rows and multiple spanwise flow splitters. An existing procedure developed by Dr. J. J. Adamczyk and associates at the NASA Lewis Research Center was modified to accept multiple spanwise splitter geometries and simulate engine core conditions. The numerical solution is based upon a finite volume technique with a four stage Runge-Kutta time marching procedure. Multiple blade row solutions are based upon the average-passage system of equations. The numerical solutions are performed on an H-type grid system, with meshes meeting the requirement of maintaining a common axisymmetric mesh for each blade row grid. The analysis was run on several geometry configurations ranging from one to five blade rows and from one to four radial flow splitters. The efficiency of the solution procedure was shown to be the same as the original analysis.

Crook, Andrew J.↗

Accelerating Large Data Analysis By Exploiting Regularities

We present techniques for discovering and exploiting regularity in large curvilinear data sets. The data can be based on a single mesh or a mesh composed of multiple submeshes (also known as zones). Multi-zone data are typical to Computational Fluid Dynamics (CFD) simulations. Regularities include axis-aligned rectilinear and cylindrical meshes as well as cases where one zone is equivalent to a rigid-body transformation of another. Our algorithms can also discover rigid-body motion of meshes in time-series data. Next, we describe a data model where we can utilize the results from the discovery process in order to accelerate large data visualizations. Where possible, we replace general curvilinear zones with rectilinear or cylindrical zones. In rigid-body motion cases we replace a time-series of meshes with a transformed mesh object where a reference mesh is dynamically transformed based on a given time value in order to satisfy geometry requests, on demand. The data model enables us to make these substitutions and dynamic transformations transparently with respect to the visualization algorithms. We present results with large data sets where we combine our mesh replacement and transformation techniques with out-of-core paging in order to achieve significant speed-ups in analysis.

Moran, Patrick J.↗

Aerodynamics of Engine-Airframe Interaction

The report describes progress in research directed towards the efficient solution of the inviscid Euler and Reynolds-averaged Navier-Stokes equations for transonic flows through engine inlets, and past complete aircraft configurations, with emphasis on the flowfields in the vicinity of engine inlets. The research focusses upon the development of solution-adaptive grid procedures for these problems, and the development of multi-grid algorithms in conjunction with both, implicit and explicit time-stepping schemes for the solution of three-dimensional problems. The work includes further development of mesh systems suitable for inlet and wing-fuselage-inlet geometries using a variational approach. Work during this reporting period concentrated upon two-dimensional problems, and has been in two general areas: (1) the development of solution-adaptive procedures to cluster the grid cells in regions of high (truncation) error;and (2) the development of a multigrid scheme for solution of the two-dimensional Euler equations using a diagonalized alternating direction implicit (ADI) smoothing algorithm.

Caughey, D. A.↗

COMGEN-BEM: Boundary element model generation for composite materials micromechanical analysis

Composite Model Generation-Boundary Element Method (COMGEN-BEM) is a program developed in PATRAN command language (PCL) which generates boundary element models of continuous fiber composites at the micromechanical (constituent) scale. Based on the entry of a few simple parameters such as fiber volume fraction and fiber diameter, the model geometry and boundary element model are generated. In addition, various mesh densities, material properties, fiber orientation angles, loads, and boundary conditions can be specified. The generated model can then be translated to a format consistent with a boundary element analysis code such as BEST-CMS.

Goldberg, Robert K.↗