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,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

AI-based design of a nuclear reactor core

The authors developed an artificial intelligence (AI)-based algorithm for the design and optimization of a nuclear reactor core based on a flexible geometry and demonstrated a 3× improvement in the selected performance metric: temperature peaking factor. The rapid development of advanced, and specifically, additive manufacturing (3-D printing) and its introduction into advanced nuclear core design through the Transformational Challenge Reactor program have presented the opportunity to explore the arbitrary geometry design of nuclear-heated structures. The primary challenge is that the arbitrary geometry design space is vast and requires the computational evaluation of many candidate designs, and the multiphysics simulation of nuclear systems is very time-intensive. Therefore, the authors developed a machine learning-based multiphysics emulator and evaluated thousands of candidate geometries on Summit, Oak Ridge National Laboratory’s leadership class supercomputer. The results presented in this work demonstrate temperature distribution smoothing in a nuclear reactor core through the manipulation of the geometry, which is traditionally achieved in light water reactors through variable assembly loading in the axial direction and fuel shuffling during refueling in the radial direction. The conclusions discuss the future implications for nuclear systems design with arbitrary geometry and the potential for AI-based autonomous design algorithms.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Numerical Simulations of Combustion Instabilities in a Combustor with an Augmentor-Like Geometry

With the goal of assessing the capability of Computational Fluid Dynamics (CFD) to simulate combustion instabilities, the present work considers a premixed, bluff-body-stabilized combustor with well-defined inlet and outlet boundary conditions. The present simulations produce flow behaviors in good qualitative agreement with experimental observations. Notably, the flame flapping and standing acoustic waves seen in the experiments are reproduced by the simulations. Moreover, present predictions for the dominant instability frequency have an error of 7% and those of the rmspressure fluctuations show an error of 16%. In addition, an analysis of simulation results for the limit cycle complements previous experimental analyses by supporting the presence of an active frequency-locking mechanism.

Gonzalez-Juez, Esteban↗
NASA NTRS1988

Digital image processing for quantification through full flow field tracing (FFFT) in narrow geometries at low Reynolds numbers

This paper introduces a computer based image processing technique to the field of nonintrusive velocity measurements in fluid mechanics. The method is presented with two alternatives: the first is involving intense interaction between the operator and the computer system; the second is a first generation artificial intelligence based system, where a set of initially imputed rules replaces the operator. The methods are applied to flow in narrow gaps.

Braun, M. J.↗
DOE OSTI2017

Generation of an Activation Map for Decommissioning Planning of the Berlin Experimental Reactor-II

The BER-II is an experimental facility with 10 MW that was operated since 1974. Its operation will end in 2019. To support the decommissioning planning, a map with the overall distribution of relevant radionuclides has to be created according to the state of the art. The activation analysis comprises all load periods until the planned shutdown. A procedure for an activation calculation of the pool structure of the BER-II research reactor including the concrete shielding is presented. With the combination of the options offered by MCNP6 and PARTISN, an effective method to compute the neutron flux for a reactor geometry model with 2,790,000 cells could be achieved. The fact that for calculations with PARTISN some extra effort to generate the XS is needed is more than compensated by the fact that for all volumes of the model, neutron fluxes are produced that can be subsequently used for the activation calculation. With the power history provided together with the material compositions, an activation map was created for the whole structure, even for the most remote parts. The activity can be calculated for all nuclides of the FISPACT library. Depending on material and irradiation history, however, only a limited number of nuclides must be regarded to reduce the size of activation data base. The results for the activation at calibration points give an example that the use of detailed calculation methods allows reliable estimation of the activity of relevant nuclides. The mean and the maximum activity can be determined. For a more detailed analysis, it is desirable to provide some more calibration points and realistic material specifications including impurities where the results of the activation can be validated, especially for more remote regions of the concrete structures.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRSDate not supplied

Comparison of OVERFLOW Computational and Experimental Results for a Blunt Mars Entry Vehicle Concept During Supersonic Retropropulsion

Simulations of supersonic retropropulsion (SRP) flow over a Hypersonic Inflatable Aerodynamic Decelerator (HIAD) blunt-body vehicle were performed using the Overflow Computational Fluid Dynamics (CFD) solver. Simulation conditions and geometry were designed to match specific test runs in the Descent System Study (DSS) testing campaign. The relative accuracy of simulation predictions are assessed by direct comparison to experimental data. Computational predictions of the SRP flowfield and bow shock shape are compared to experimental schlieren images. Comparisons of the model surface pressure environment are presented for unsteady and static discrete tap data as well as time-averaged Pressure-Sensitive Paint (PSP) data.

Supersonic Retropropulsion↗
NASA NTRS2001

A Combined Experimental/Computational Investigation of a Rocket Based Combined Cycle Inlet

A rocket based combined cycle inlet geometry has undergone wind tunnel testing and computational analysis with Mach 4 flow at the inlet face. Performance parameters obtained from the wind tunnel tests were the mass capture, the maximum back-pressure, and the self-starting characteristics of the inlet. The CFD analysis supplied a confirmation of the mass capture, the inlet efficiency and the details of the flowfield structure. Physical parameters varied during the test program were cowl geometry, cowl position, body-side bleed magnitude and ingested boundary layer thickness. An optimum configuration was determined for the inlet as a result of this work.

Smart, Michael K.↗
NASA NTRS1994

Thermohydrodynamic analysis of cryogenic liquid turbulent flow fluid film bearings, phase 2

The Phase 2 (1994) Annual Progress Report presents two major report sections describing the thermal analysis of tilting- and flexure-pad hybrid bearings, and the unsteady flow and transient response of a point mass rotor supported on fluid film bearings. A literature review on the subject of two-phase flow in fluid film bearings and part of the proposed work for 1995 are also included. The programs delivered at the end of 1994 are named hydroflext and hydrotran. Both codes are fully compatible with the hydrosealt (1993) program. The new programs retain the same calculating options of hydrosealt plus the added bearing geometries, and unsteady flow and transient forced response. Refer to the hydroflext & hydrotran User's Manual and Tutorial for basic information on the analysis and instructions to run the programs. The Examples Handbook contains the test bearing cases along with comparisons with experimental data or published analytical values. The following major tasks were completed in 1994 (Phase 2): (1) extension of the thermohydrodynamic analysis and development of computer program hydroflext to model various bearing geometries, namely, tilting-pad hydrodynamic journal bearings, flexure-pad cylindrical bearings (hydrostatic and hydrodynamic), and cylindrical pad bearings with a simple elastic matrix (ideal foil bearings); (2) improved thermal model including radial heat transfer through the bearing stator; (3) calculation of the unsteady bulk-flow field in fluid film bearings and the transient response of a point mass rotor supported on bearings; and (4) a literature review on the subject of two-phase flows and homogeneous-mixture flows in thin-film geometries.

Sanandres, Luis↗
DOE OSTI2025

Supporting Data: Inverse Trans Influence and Uranium-Arene σ-Bonding Drive Molecular Geometry: Ligand Modification from Hard to Soft Flips the Oxide from Axial to Equatorial

Supporting Computational Data from the manuscript titled "ITI and Uranium-Arene s-Bonding Drive Molecular Geometry: Ligand Modification from Hard to Soft Flips the Oxide from Axial to Equatorial" by Benedikt Kestel, Daniel Pividori, Sabyasachi Roy Chowdhury, Frank W. Heinemann, Bess Vlaisavljevich, and Karsten Meyer published in J. Am. Chem. Soc. (DOI: 10.1021/jacs.5c05029).

Chowdhury, Sabyasachi Roy [Department of Chemistry↗
NASA NTRS1981

Simulation of charge exchange plasma propagation near an ion thruster propelled spacecraft

A model describing the charge exchange plasma and its propagation is discussed, along with a computer code based on the model. The geometry of an idealized spacecraft having an ion thruster is outlined, with attention given to the assumptions used in modeling the ion beam. Also presented is the distribution function describing charge exchange production. The barometric equation is used in relating the variation in plasma potential to the variation in plasma density. The numerical methods and approximations employed in the calculations are discussed, and comparisons are made between the computer simulation and experimental data. An analytical solution of a simple configuration is also used in verifying the model.

Robinson, R. S.↗
NASA NTRS2010

Analysis of a Channeled Centerbody Supersonic Inlet for F-15B Flight Research

The Propulsion Flight Test Fixture at the NASA Dryden Flight Research Center is a unique test platform available for use on the NASA F-15B airplane, tail number 836, as a modular host for a variety of aerodynamics and propulsion research. The first experiment that is to be flown on the test fixture is the Channeled Centerbody Inlet Experiment. The objectives of this project at Dryden are twofold: 1) flight evaluation of an innovative new approach to variable geometry for high-speed inlets, and 2) flight validation of channeled inlet performance prediction by complex computational fluid dynamics codes. The inlet itself is a fixed-geometry version of a mixed-compression, variable-geometry, supersonic in- let developed by TechLand Research, Inc. (North Olmsted, Ohio) to improve the efficiency of supersonic flight at off-nominal conditions. The concept utilizes variable channels in the centerbody section to vary the mass flow of the inlet, enabling efficient operation at a range of flight conditions. This study is particularly concerned with the starting characteristics of the inlet. Computational fluid dynamics studies were shown to align well with analytical predictions, showing the inlet to remain unstarted as designed at the primary test point of Mach 1.5 at an equivalent pressure altitude of 29,500 ft local conditions. Mass-flow-related concerns such as the inlet start problem, as well as inlet efficiency in terms of total pressure loss, are assessed using the flight test geometry.

Ratnayake, Nalin A.↗
NASA NTRS1985

A computer-assisted process for supersonic aircraft conceptual design

Design methodology was developed and existing major computer codes were selected to carry out the conceptual design of supersonic aircraft. A computer-assisted design process resulted from linking the codes together in a logical manner to implement the design methodology. The process does not perform the conceptual design of a supersonic aircraft but it does provide the designer with increased flexibility, especially in geometry generation and manipulation. Use of the computer-assisted process for the conceptual design of an advanced technology Mach 3.5 interceptor showed the principal benefit of the process to be the ability to use a computerized geometry generator and then directly convert the geometry between formats used in the geometry code and the aerodynamics codes. Results from the interceptor study showed that a Mach 3.5 standoff interceptor with a 1000 nautical-mile mission radius and a payload of eight Phoenix missiles appears to be feasible with the advanced technologies considered. A sensitivity study showed that technologies affecting the empty weight and propulsion system would be critical in the final configuration characteristics with aerodynamics having a lesser effect for small perturbations around the baseline.

Johnson, V. S.↗
NASA NTRS1976

Prediction of sonic boom at a focus

The behavior of sonic boom at a focus has been reviewed for the purpose of extending present sonic boom computational methods to include focal zones. The geometry of a focal zone - whether a smooth caustic, a cusped caustic, or a perfect focus to a point - determines the character of focused signatures. The seeming contradiction of various experimental data can be resolved by noting these differences. A ray acoustic analysis has been developed for quantitative determination of caustic geometry. The only reliable theory presently available for signatures at a focus is for a smooth caustic. There has been some controversy between theoretical and experimental values of a constant in the scaling law for this case. It has been found that this discrepancy can be resolved by accounting for the finite thickness of real sonic boom shock waves. These findings have been incorporated into an existing sonic boom computer program.

Plotkin, K. J.↗
DOE OSTI2014

High Energy Electromagnetic Particle Transportation on the GPU

We present massively parallel high energy electromagnetic particle transportation through a finely segmented detector on a Graphics Processing Unit (GPU). Simulating events of energetic particle decay in a general-purpose high energy physics (HEP) detector requires intensive computing resources, due to the complexity of the geometry as well as physics processes applied to particles copiously produced by primary collisions and secondary interactions. The recent advent of hardware architectures of many-core or accelerated processors provides the variety of concurrent programming models applicable not only for the high performance parallel computing, but also for the conventional computing intensive application such as the HEP detector simulation. The components of our prototype are a transportation process under a non-uniform magnetic field, geometry navigation with a set of solid shapes and materials, electromagnetic physics processes for electrons and photons, and an interface to a framework that dispatches bundles of tracks in a highly vectorized manner optimizing for spatial locality and throughput. Core algorithms and methods are excerpted from the Geant4 toolkit, and are modified and optimized for the GPU application. Program kernels written in C/C++ are designed to be compatible with CUDA and OpenCL and with the aim to be generic enough for easy porting to future programming models and hardware architectures. To improve throughput by overlapping data transfers with kernel execution, multiple CUDA streams are used. Issues with floating point accuracy, random numbers generation, data structure, kernel divergences and register spills are also considered. Performance evaluation for the relative speedup compared to the corresponding sequential execution on CPU is presented as well.

Canal, P.↗
DOE OSTI2016

Localized Radial Basis Functions for Radiation Transport in Slab Geometry

Meshing the spatial domain for the neutron transport equation can be computationally expensive and labor-intensive. Mesh-free methods, such as radial basis function (RBF) approximations, may alleviate these concerns by using a pointwise solution of the transport equation instead of the standard volume-integrated solution using cells or elements. This allows for the use of any geometric description that can provide the material at solution points and the boundary normal directions for boundary points. RBFs give spectral convergence when used to solve differential equations, but only if they are infinitely smooth. This occurs when the RBFs are effectively flat, creating an ill-conditioned, dense system of equations. One simple way to alleviate this issue is to use localized RBFs. Global RBFs have been applied previously to the monoenergetic radiative transport equation in X-Y geometry using source iteration. In this summary, localized RBFs are applied to the multigroup neutron transport equation in slab geometry using Krylov iteration. The results show agreement with benchmark solutions, provided the coefficients for the RBF basis are properly chosen. (authors)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

An Approach to Quad Meshing Based On Cross Valued Maps and the Ginzburg-Landau Theory

A generalization of vector fields, referred to as N-direction fields or cross fields when N=4, has been recently introduced and studied for geometry processing, with applications in quadrilateral (quad) meshing, texture mapping, and parameterization. We make the observation that cross field design for two-dimensional quad meshing is related to the well-known Ginzburg-Landau problem from mathematical physics. This identification yields a variety of theoretical tools for efficiently computing boundary-aligned quad meshes, with provable guarantees on the resulting mesh, for example, the number of mesh defects and bounds on the defect locations. The procedure for generating the quad mesh is to (i) find a complex-valued "representation" field that minimizes the Dirichlet energy subject to a boundary constraint, (ii) convert the representation field into a boundary-aligned, smooth cross field, (iii) use separatrices of the cross field to partition the domain into four sided regions, and (iv) mesh each of these four-sided regions using standard techniques. Under certain assumptions on the geometry of the domain, we prove that this procedure can be used to produce a cross field whose separatrices partition the domain into four sided regions. To solve the energy minimization problem for the representation field, we use an extension of the Merriman-Bence-Osher (MBO) threshold dynamics method, originally conceived as an algorithm to simulate motion by mean curvature, to minimize the Ginzburg-Landau energy for the optimal representation field. Lastly, we demonstrate the method on a variety of test domains.

Computational Geometry↗
NASA NTRS1984

Computation of three-dimensional viscous flows using a space-marching method

A space-marching method, developed to compute three-dimensional flows for internal geometries, has been utilized to predict viscous flows through a curved duct and over a swept wing. The Navier-Stokes equations have been posed as an initial value problem by neglecting the streamwise viscous diffusion terms and by treating the pressure gradient as a known source term. The resulting equations have been solved by a non-iterative (single pass) algorithm at each streamwise step. The results are compared with earlier computations (based on iterative methods) and the experimental data. The agreement between the present predictions, the experimental data, and the earlier predictions is good for the cases computed. The computation time is only a fraction of the iterative methods.

Murthy, K. N. S.↗
NASA NTRS2002

Progress Toward Generation of a Navier-Stokes Database for a Harrier in Ground Effect

The Harrier YAV-8B aircraft is capable of vertical and short-field take-off and landing (V/STOL) by directing its four exhaust nozzles toward the ground, or conventional flight by rotating its nozzles into a horizontal position. The British Royal Air Force and the United States Marine Corps have used this aircraft for more than 30 years to provide a quick reaction time for troop support, and reduce the need for long runways. The success of this powered-lift (PL) vehicle has also prompted the more recent design of the Joint Strike Fighter (JSF). However there are significant safety issues that must be addressed when operating a PL vehicle in close proximity to the ground. Hot Gas Ingestion (HGI) by the inlets can result in a rapid loss of powered lift; and high-speed jet flows along the ground plane can induce low pressures underneath the vehicle, causing a 'suck-down' effect. Under these conditions, departure from controlled flight may occur. Moreover, unsteady ground vortices and jet fountains can affect the aircraft,s controllability and its proximity to ground troops. The viscous, time-dependent flow fields of PL vehicles are difficult to accurately and efficiently predict using Computational Fluid Dynamics (CFD). A number of researchers have used the time-dependent Reynolds-averaged Navier-Stokes (RANS) equations to compute flows for single and multiple jets in a cross-flow. A few have added some geometric complexity to the problem by computing flows for jet-augmented delta wings near a ground plane. Smith et.al. computed for the first time a single RANS solution about a simplified Harrier. This geometry included a fuselage, wing, leading edge root extension (LERX), inlets, and exhaust nozzles. All of these investigations cite two practical problems with computing these flows: 1) the need for improved solution accuracy; and, 2) the need for faster solution methods. We view the need for faster solution methods as key to improving the solution accuracy and making this class of computation more routine. One can hardly refine grids, explore the use of advanced turbulence models, and generate databases when it takes weeks of dedicated computer time for a single solution. Chaderjian, Ahmad, Pandya, and Murman have focused on reducing the time-to-solution for this very difficult and complex problem through process automation and exploitation of parallel computing. They began with the Harrier geometry reported, and added a deflected wing flap and empennage for greater realism. To date more than 80 solutions have been carried out. This paper will describe this process and progress made in reducing the time required to generate a simple longitudinal force and moment database for a Harrier in ground effect. It shows a typical snap-shot from an unsteady streakline animation, where fluid particles are colored by temperature. The ground vortex and a jet-fountain vortex are highlighted. It also shows a similar streakline image, where HGI occurs due to the vehicle in close proximity to the ground. It is show the mean lift coefficient as a function of angle of attack and height. The angle of attack range was 4 deg less than or = alpha less than or = 10 deg with an increment of 1 degree, and the height range was 10 ft less than or = h less than or = 30ft with an increment of 5 feet. This 35 solution database was extended to over 2500 cases using a monotone cubic-spline interpolation procedure. The suck-down effect (reduction of lift near the ground) is highlighted in the figure. The "cushion effect," the conventional reduction of lift as the vehicle moves out of ground effect, is also indicated. All 35 RANS solutions were obtained using 952 Silicon Graphics Origin 2000 and 3000 processors in dedicated mode for one week. Typically, 112 processors were assigned to each case. Some other cases used fewer processors to utilize all available CPUS. The final paper will report on the automation of the solution process, including: grid generation, job monitoring, solution completion criteria, and post processing. Moreover, improvements in parallel efficiency for a dual time-step algorithm for the RANS equations will also be presented. Results will be discussed in detail using unsteady streakline flow visualization to correlate unsteady flow structures with dominant aerodynamic frequencies. The stability derivatives, CL, and CL, will also be presented.

Chaderjian, Neal M.↗
NASA NTRS1994

Sonic boom predictions using a solution-adaptive full-potential code

The TranAir full-potential code, which utilizes a solution-adaptive, unstructured Cartesian grid, was used to compute near-field pressure signatures of three generic geometries: (1) a cone-cylinder; (2) a low aspect ratio rectangular wing; and (3) a delta-wing/body. Pressure signatures computed by TranAir are extrapolated to near- and midfield distances using a quasilinear extrapolation technique. Results are compared with wind-tunnel data and, where possible, with extrapolated wind-tunnel data and Whitham-method results. In general, extrapolated TranAir data is in good agreement with experimental and Whitham-method pressure signatures, but peak pressures are slightly underpredicted. It is felt that data computed too close to the model caused some of the discrepancies. The flowfield grid structure was also found to be an important issue in obtaining accurate off-body data. Results were encouraging enough to proceed with the application of TranAir to conceptual configurations in support of the NASA high speed research program (HSRP).

Madson, Michael D.↗