Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational Geometry”

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,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

A Combined Geometric Approach for Computational Fluid Dynamics on Dynamic Grids

A combined geometric approach for computational fluid dynamics is presented for the analysis of unsteady flow about mechanisms in which its components are in moderate relative motion. For a CFD analysis, the total dynamics problem involves the dynamics of the aspects of geometry modeling, grid generation, and flow modeling. The interrelationships between these three aspects allow for a more natural formulation of the problem and the sharing of information which can be advantageous to the computation of the dynamics. The approach is applied to planar geometries with the use of an efficient multi-block, structured grid generation method to compute unsteady, two-dimensional and axisymmetric flow. The applications presented include the computation of the unsteady, inviscid flow about a hinged-flap with flap deflections and a high-speed inlet with centerbody motion as part of the unstart / restart operation.

Slater, John W.↗
NASA NTRS1973

Maximally slicing a black hole.

Analytic and computer-derived solutions are presented of the problem of slicing the Schwarzschild geometry into asymptotically flat, asymptotically static, maximal spacelike hypersurfaces. The sequence of hypersurfaces advances forward in time in both halves (u greater than or equal to 0, u less than or equal to 0) of the Kruskal diagram, tending asymptotically to the hypersurface r = 3/2 M and avoiding the singularity at r = 0. Maximality is therefore a potentially useful condition to impose in obtaining computer solutions of Einstein's equations.

Estabrook, F.↗
NASA NTRS2008

Heat Transfer Measurements for a Film Cooled Turbine Vane Cascade

Experimental heat transfer and pressure measurements were obtained on a large scale film cooled turbine vane cascade. The objective was to investigate heat transfer on a commercial high pressure first stage turbine vane at near engine Mach and Reynolds number conditions. Additionally blowing ratios and coolant density were also matched. Numerical computations were made with the Glenn-HT code of the same geometry and compared with the experimental results. A transient thermochromic liquid crystal technique was used to obtain steady state heat transfer data on the mid-span geometry of an instrumented vane with 12 rows of circular and shaped film cooling holes. A mixture of SF6 and Argon gases was used for film coolant to match the coolant-to-gas density ratio of a real engine. The exit Mach number and Reynolds number were 0.725 and 2.7 million respectively. Trends from the experimental heat transfer data matched well with the computational prediction, particularly for the film cooled case.

Poinsatte, Philip E.↗
NASA NTRSDate not supplied

Streamline-Traced, External-Compression Supersonic Inlets for Mach 2

A computational study was performed to explore the aerodynamic design and performance of streamline-traced, external-compression (STEX) inlets for Mach 2.0. The performance metrics included inlet flow rates, total pressure recovery, and total pressure distortion. The study explored the use of round and flattened capture cross-sections and a bleed slot. The design of the inlet and generation of the inlet geometry was performed using the Supersonic Inlet Design and Analysis(SUPIN) Tool. Computational grids were generated, and methods of computational fluid dynamics (CFD) were applied to solve the three-dimensional, turbulent flow through the inlets using the Wind-US CFD flow solver. It was found that the use of a bleed slot with about 5% bleed resulted in an inlet total pressure recovery of 95% at Mach 2 with acceptable radial and circumferential total pressure distortion. It was also found that there was only a slight decrease in performance between a round and flattened capture cross-section. This suggests that the use of streamline tracing offers flexibility in shaping the inlet capture cross-section for a more favorable integration with an aircraft wing or fuselage. This flexibility and the good performance of the inlet provides a novel inlet design for future commercial supersonic aircraft.

Computational Fluid Dynamics↗
NASA NTRS1994

A Cartesian grid approach with hierarchical refinement for compressible flows

Many numerical studies of flows that involve complex geometries are limited by the difficulties in generating suitable grids. We present a Cartesian boundary scheme for two-dimensional, compressible flows that is unfettered by the need to generate a computational grid and so it may be used, routinely, even for the most awkward of geometries. In essence, an arbitrary-shaped body is allowed to blank out some region of a background Cartesian mesh and the resultant cut-cells are singled out for special treatment. This is done within a finite-volume framework and so, in principle, any explicit flux-based integration scheme can take advantage of this method for enforcing solid boundary conditions. For best effect, the present Cartesian boundary scheme has been combined with a sophisticated, local mesh refinement scheme, and a number of examples are shown in order to demonstrate the efficacy of the combined algorithm for simulations of shock interaction phenomena.

Quirk, James J.↗
DOE OSTI2022

Demonstrate multi-turbine simulation with hybrid-structured / unstructured-moving-grid software stack running primarily on GPUs and propose improvements for successful KPP-2

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-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, capturing the thin boundary layers, 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.

17 WIND ENERGY↗
DOE OSTI2021

The FMO-DFTB Method

Although the fragment molecular orbital (FMO) method enables electronic structure calculations with near-linear scaling behavior with respect to system size, the computational cost of ab initio methods typically employed in conjunction with FMO is still prohibitive for routine calculations of very large systems or long timescale molecular dynamics simulations. We, therefore, combined the FMO and density-functional tight-binding (DFTB) method, which is one of the emerging semi-empirical quantum chemical methods, and have demonstrated that FMO-DFTB is capable of performing geometry optimizations for systems containing up to one million atoms using limited computational resources. In this chapter, we will review the basics of the DFTB method first before introducing FMO-DFTB, focusing on the relationship with density functional theory and other FMO methodologies. We also demonstrate the latest scalings of FMO-DFTB2 and DFTB3 using three-dimensional water clusters, showing that the most favorable scaling is O(N1.16). Applications of FMO-DFTB to various systems are briefly summarized, and an outlook to future applications is provided.

Nishimoto, Yoshio↗
DOE OSTI2022

Review of Computer-Aided Manufacturing (CAM) strategies for hybrid directed energy deposition

Hybrid additive manufacturing intertwines both additive and subtractive manufacturing layer by layer to digitally fabricate parts with complex geometries, improved surface finish and tight dimensional accuracies, the sum of which is difficult to obtain with any single process. Computer-Aided Manufacturing (CAM) software is required to orchestrate the machine toolpathing for both the deposition as well as the machining processes and is crucial for the successful fabrication of high quality structures. Additionally, CAM requires substantial operator input to account for challenging aspects of each fabricated structure. For example, deciding at which layer of deposition will the machining process continue to maintain access to complex cavities for finishing - internal features that would otherwise be unfinishable due to reach limitations or obstructions. Moreover, of the many commercially-available hybrid systems, each has a unique kinematic environment which can benefit from specific optimization of toolpath planning and substantial research has directly correlated toolpathing with the microstructure evolution, mechanical properties, porosity and residual stress state of the final fabricated part. This review explores the available strategies for CAM in the context of hybrid direct energy deposition, discusses the advantages and disadvantages of each and considers future CAM trends for this transformational digital manufacturing technology.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Demonstrate moving-grid multi-turbine simulations primarily run on GPUs and propose improvements for successful KPP-2

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-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, capturing the thin boundary layers, 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.

17 WIND ENERGY↗
NASA NTRS1991

Analysis and design of planar and non-planar wings for induced drag minimization

The goal of the work was to develop and validate computational tools to be used for the design of planar and non-planar wing geometries for minimum induced drag. Because of the iterative nature of the design problem, it is important that, in addition to being sufficiently accurate for the problem at hand, they are reasonably fast and computationally efficient. Toward this end, a method of predicting induced drag in the presence of a non-rigid wake is coupled with a panel method. The induced drag prediction technique is based on the Kutta-Joukowski law applied at the trailing edge. Until recently, the use of this method has not been fully explored and pressure integration and Trefftz-plane calculations favored. As is shown in this report, however, the Kutta-Joukowski method is able to give better results for a given amount of effort than the more common techniques, particularly when relaxed wakes and non-planar wing geometries are considered. Using these tools, a workable design method is in place which takes into account relaxed wakes and non-planar wing geometries. It is recommended that this method be used to design a wind-tunnel experiment to verify the predicted aerodynamic benefits of non-planar wing geometries.

Mortara, K.↗
DOE OSTI2020

Edge ICRF simulations in 3D geometry: From MHD equilibrium to coupling determination

We present in this work a consistent numerical scheme that allows the computation of 3D magnetic fields a nd 3D density profiles and their usage in ion cyclotron range of frequencies (ICRF) coupling simulations. We first utilize the PARVMEC code to compute the 3D free-boundary plasma equilibrium in the ideal magnetohydrodynamic (MHD) approximation. Since the PARVMEC solution is only defined within the last closed flux surface (LCFS), the magnetic field domain is extended to the scrape-off layer (SOL) via the BMW code, which computes a divergence-free magnetic field solution arising from the external conductors’ vacuum field a nd t he PARVMEC flux surface currents. This magnetic reconstruction is then used in the EMC3-EIRENE transport code in order to compute 3D density profiles. In the last step, the RAPLICASOL code is utilized to compute the ICRF antenna S-matrices resulting from the 3D density profiles. We exemplify this scheme for the ASDEX Upgrade tokamak. A new implementation of a curved model for the ASDEX Upgrade ICRF 2-strap antenna in RAPLICASOL allows simulations in realistic geometry, without any coordinate transformations.

Lopez, Graham↗
NASA NTRS1997

Efficient Multidisciplinary Analysis Procedure Using Multi-Level Parallelization Approach

Multidisciplinary applications are suitable for parallel computing environment by adopting the domain decomposition method. Immediately, a multidisciplinary application can be parallelized by solving each discipline separately. In order to perform coupled multidisciplinary analysis, coupling of each discipline can be accomplished by exchanging boundary data at the interfaces. This is regarded as discipline-level parallelization. Next level could be a "coarse-grain" parallelization of each discipline, which mainly depends on the physical geometry and nature of each discipline. For example, it is almost impossible for structured-grid based computational fluid dynamics codes to do flow analysis of an aircraft by using a single grid because of the complexity of its configuration. Thus, multi-block grid is commonly used to describe the details of complex geometry. Similarly, in structural analysis, the structure is frequently subdivided into substructures. Thus, the computation of each subdomain can be easily parallelized since each subdomain is solved separately independent of other domains. The parallelization is accomplished by solving each subdomain separately on a separate processor and exchanging the boundary conditions at domain interfaces periodically. However, the physical decomposition of the domain introduces explicit boundary conditions at the domain interfaces. This is not desirable for critical areas such as those containing shock waves or flow separations. Thus, a "fine-grain" parallelization is introduced to overcome this problem. The "fine-grain" parallelization is one that solves exactly the same system of equations of a subdomain by using more than one processors without introducing any explicit boundary conditions. An efficient multidisciplinary analysis procedure can be accomplished by successfully combining the above multi-level parallelism. A multidisciplinary analysis code, ENSAERO developed at NASA Ames Research Center is used in this study to implement the proposed approach. The communication data structure required for the proposed approach will be studied in detail. This work will demonstrate the feasibility of using multi-level parallelization approach in multidisciplinary analysis applications.

Byun, Chansup↗
NASA NTRS1988

Effects of surface pressures and streamline metrics on the calculation of laminar heating rates

The effect of streamline geometry and pressure distributions on surface heating rates is examined for slender, spherically blunted cones. The modifications to the approximate aeroheating code include a curve fit of pressures computed by an Euler solution over a range of Mach numbers and cone angles. The streamline geometry is then found using the surface pressures and inviscid surface properties. Previously, streamlines were determined using the inviscid properties at the edge of the boundary layer when accounting for the effects of entropy-layer swallowing. Streamline calculations are now based on inviscid surface conditions rather than boundary-layer edge properties. However, the heating rates are calculated using inviscid properties at the edge of the boundary layer. Resulting heating rates compare favorably with solutions from the viscous-shock-layer equations.

Riley, Christopher J.↗
NASA NTRS1996

Research in Parallel Algorithms and Software for Computational Aerosciences

Phase I is complete for the development of a Computational Fluid Dynamics parallel code with automatic grid generation and adaptation for the Euler analysis of flow over complex geometries. SPLITFLOW, an unstructured Cartesian grid code developed at Lockheed Martin Tactical Aircraft Systems, has been modified for a distributed memory/massively parallel computing environment. The parallel code is operational on an SGI network, Cray J90 and C90 vector machines, SGI Power Challenge, and Cray T3D and IBM SP2 massively parallel machines. Parallel Virtual Machine (PVM) is the message passing protocol for portability to various architectures. A domain decomposition technique was developed which enforces dynamic load balancing to improve solution speed and memory requirements. A host/node algorithm distributes the tasks. The solver parallelizes very well, and scales with the number of processors. Partially parallelized and non-parallelized tasks consume most of the wall clock time in a very fine grain environment. Timing comparisons on a Cray C90 demonstrate that Parallel SPLITFLOW runs 2.4 times faster on 8 processors than its non-parallel counterpart autotasked over 8 processors.

Domel, Neal D.↗
NASA NTRSDate not supplied

Aerodynamic and Acoustic Interactions Associated with Inboard Propeller-Wing Configurations

A series of aerodynamic performance and acoustic measurements have been made on a range of inboard propeller-wing interaction configurations in the NASA Langley Low Speed Aeroacoustic Wind Tunnel (LSAWT). The results presented in this paper are part of a more expansive testing campaign encompassing both single propeller-wing and multipropeller-wing interactions, the former of which is discussed in the present work. The primary testing parameters of interest to this study are the axial and vertical positioning of the wing relative to the propeller slipstream under a constant propeller advance ratio. A multi-faceted computational effort was also employed in an effort to identify reflection and scattering effects imposed by both the wing geometry as well as the primary components of the facility test setup. This effort consisted of aerodynamic predictions using high-fidelity computational fluid dynamics (CFD), acoustic predictions using an impermeable Ffowcs Williams and Hawkings (FW-H) solver, and acoustic scattering predictions. Acoustic measurements reveal variations in the acoustic directivity behavior of the propeller blade passage frequency for even modest variations in wing position. CFD-based acoustic predictions reveal discrepancies relative to the experimental data, which is believed to be due to complex acoustic scattering behavior within the test section. Initial attempts at modeling the scattered acoustic field showed functional dependency of the acoustic amplitude variations on the wing position relative to the propeller disk, however discrepancies with experimental data remain.

Nikolas S. Zawodny↗
DOE OSTI2025

A Block-Structured Adaptive Mesh Framework to Solve Radiation Transfer Equation in Irregular Embedded Geometries

Radiation transport arises in various scientific, industrial, and medical fields, and understanding its effect in applications is needed to make accurate predictions, safety assessments and performance optimizations. Solving the Radiation Transport Equation (RTE) is challenging due to its integro-differential nature, which involves both differential and integral terms. The differential term describes the change in radiation intensity due to absorption and emission, while the integral term accounts for scattering. The accurate modeling of radiation is further complicated in many applications due to the complex, irregular geometries. Various methods exist for solving the RTE, including the zonal, Monte Carlo, spherical harmonics, discrete ordinates, and finite volume methods. Traditional mesh-based approaches, which rely on structured or unstructured meshes, struggle with irregular geometries due to: a) the difficulty of conforming structured grids to irregular domains, b) challenges in enforcing boundary conditions correctly, and c) the additional computational cost of unstructured mesh methods. This work presents a second-order accurate method for solving the RTE in irregular geometries. The radiation intensity is discretized using the finite-volume method in both spatial and angular directions on regular Cartesian grid blocks. Leveraging the block-structured adaptive mesh refinement (AMR) framework provided by AMReX, our method refines the grid locally to reduce spatial discretization error, ensuring a converged numerical solution while minimizing computational costs elsewhere. A two-stage deferred correction approach is employed: First, a first-order discretization on grid blocks is solved using an algebraic multigrid method in HYPRE. Second, a correction term is applied explicitly to achieve second-order accuracy. The correction term is calculated by approximating the radiation flux on cell faces using a Total Variation Diminishing (TVD) scheme. This approach ensures quick convergence of the multigrid method while preserving higher-order accuracy of the numerical solution. Irregular geometries are resolved as embedded boundaries (EB), resulting in both cut cells and regular cells. In cut cells, we modify the fluxes using face fractions and incorporate additional contributions from EB boundary conditions. To ensure higher-order convergence near the EB interface, the correction term is modified by interpolating the radiation intensity to fictitious ghost points. The implementation takes advantage of modern supercomputers by leveraging AMReX’sMPI/X parallelization strategy where X can be MPI or a GPU accelerator including CUDA, HIP and DPC++. We validate our solver using classical test cases, both with and without EB, demonstrating accuracy and efficiency. Additionally, we analyze the impact of adaptive mesh refinement on solution accuracy and computational cost, highlighting the advantages of our approach for high-resolution radiation transport simulations.

computational fluid dynamics (CFD)↗
DOE OSTI2022

Development and Validation of Resistance-Capacitance Model (RCM) for Phase Change Material (PCM) Embedded in 3D Periodic Structures

The low thermal conductivity of Phase Change Materials (PCM) can be improved with extended surfaces such as additively manufactured 3D periodic lattice structures. Three different aluminum alloy-based lattices (base sizes 10, 20, and 40 mm) with average porosity of 0.95 filled with paraffin wax, with a nominal phase change temperature of 55°C, were experimentally investigated. In this work, a computationally efficient 2D Resistance Capacitance-based model (RCM) was developed for predicting the thermal characteristics of these geometries. Non-uniform porosity in the PCM-metal domain was estimated using image processing and served as model input. The solver does not solve for higher-order physics as in CFD but still can provide a good prediction of thermal resistance and energy storage at a very low computational cost. The simulation-to-real-time factor for this geometry is of the order of 10-4, while CFD simulations typically have a real-time factor greater than 1. The model was validated against the experimental data for melting under three different heat fluxes (6250 W/m2, 12500 W/m2, and 18750 W/m2). The mean deviation of the predicted average PCM temperature was between 1.34 K-2.81 K for different cases. The maximum average temperature deviation of 5.45 K was observed for the 20 mm geometry at the highest heat flux test condition. The effects of natural convection were neglected in the model, but the predicted PCM temperature and energy storage still showed good agreement with the experimental data.

25 ENERGY STORAGE↗
NASA NTRS2000

Thermal Analysis of ISS Service Module Active TCS

ISS Service Module mission must begin in July 2000. The verification of design thermal requirements is mostly due to thermal analysis. The thermal analysis is enough difficult problem because of large number of ISS configurations that had to be investigated and various orbital environments. Besides the ISS structure has articulating parts such as solar arrays and radiators. The presence of articulating parts greatly increases computation times and requires accurate approach to organization of calculations. The varying geometry needs us to calculate the view factors several times during the orbit, while in static geometry case we need do it only once. In this paper we consider the thermal mathematical model of SM that includes the TCS and construction thermal models and discuss the results of calculations for ISS configurations 1R and 9Al. The analysis is based on solving the nodal heat balance equations for ISS structure by Kutta-Merson method and analytical solutions of heat transfer equations for TCS units. The computations were performed using thermal software TERM [1,2] that will be briefly described.

Altov, Vladimir V.↗