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,153 records · Page 64

Three-dimensional potential energy surfaces of ArNO ($\tilde{X}$ 2 Π)

Until now, the potential energy surfaces (PESs) of the ArNO complex found in the literature were two-dimensional, with the NO interatomic distance being fixed. In this work, we present the first accurate three-dimensional ground state X 2 Π PESs (both A ' and A ") of ArNO computed at the CCSD(T)/CBS level of theory. The equilibrium geometries and the well depths ( D e ) are compared to several other electronic structure methods. We found that using the multireference method, MRCI-F12 makes the surfaces much shallower (by 25%) and the depth of the surfaces does not agree with experimental data. The explicitly correlated coupled-cluster method underestimates the well depth as well. Analytic representations for both A ' and A " surfaces were fit to 4380 ab initio points to within 2.71 cm -1 .A three-dimensional Numerov propagator method in Delves coordinates is used to compute the bound state spectrum up to J tot = 6.5. The recommended dissociation energies are D 0 = 97.2 cm -1 for the adiabatic ground state and D e = 133.7 (128.1) cm -1 for A ' ( A ")

36 MATERIALS SCIENCE↗

TranAir: Recent advances and applications

TranAir is a computer code which solves the full-potential equation for transonic flow about very general and complex configurations. Piecewise flat surface panels are used to describe the surface geometry. This paneled definition is then embedded in an unstructured cartesian flow field grid. Finite elements are used in the discretization of the flow field grid in a manner which is fully conservative and second-order accurate. Since geometries may be defined with relative ease, and since the user is not involved in the generation of the flow field grid, computational results may be generated rather quickly for a wide range of geometries. For transonic cases in the cruise angle-of-attack range, TranAir has generated results which are in generally good agreement with both Euler results and wind tunnel data. A typical transonic case runs in 1 to 2 CPU hours on a Cray X-MP. For subcritical cases, the code runs in 15 to 30 CPU minutes, even for geometries in which several thousand surface panels are used in the definition. This ability to rapidly and accurately provide both subsonic and transonic predictions about very complex aircraft configurations gives TranAir the potential of being a very powerful and widely used design tool.

Madson, Michael D.↗

Computational Fluid Dynamic Modeling of Rocket Based Combined Cycle Engine Flowfields

Computational Fluid Dynamic techniques are used to study the flowfield of a fixed geometry Rocket Based Combined Cycle engine operating in rocket ejector mode. Heat addition resulting from the combustion of injected fuel causes the subsonic engine flow to choke and go supersonic in the slightly divergent combustor-mixer section. Reacting flow computations are undertaken to predict the characteristics of solutions where the heat addition is determined by the flowfield. Here, adaptive gridding is used to improve resolution in the shear layers. Results show that the sonic speed is reached in the unheated portions of the flow first, while the heated portions become supersonic later. Comparison with results from another code show reasonable agreement. The coupled solutions show that the character of the combustion-based thermal choking phenomenon can be controlled reasonably well such that there is opportunity to optimize the length and expansion ratio of the combustor-mixer.

Daines, Russell L.↗

Status of Mars Retropropulsion Testing in the Langley Unitary Plan Wind Tunnel

Future Mars human landings will be enabled by a powered descent phase starting at supersonic conditions, something which has never been done before on a Mars mission. Significant aerosciences challenges exist due to jet interactions between the retrorocket engine plumes, freestream flow, and vehicle that will affect the aerodynamic behavior during powered descent. Historically, wind tunnel tests have been used to study the interactions with inert gas exhaust simulants in place of rocket engines. On the computational side, flowfield simulations have been completed at full-scale conditions, but the available ground and flight data are not appropriate for calibrating computational uncertainties for aerodynamic interference on proposed Mars descent vehicles, due to insufficient data, dissimilar vehicle geometries, and disparate operating conditions. A wind tunnel test has been designed to begin addressing powered descent aerodynamics risks for large-scale human Mars entry concepts and to identify gaps in computational predictive capabilities. The test will be conducted in the NASA Langley Unitary Plan Wind Tunnel and is designed with improvements in model design and data products over past tests. The test campaign will be run using sub-scale model geometries derived from NASA powered descent reference vehicles: a blunt low lift-to-drag vehicle and a more slender geometry that generates higher unpowered lift. Both models have been fabricated and are ready for testing. The blunt model is equipped with the flexibility to examine the effects of nozzle pointing direction, number, location, size, and area ratio. The main measurements are heatshield aerodynamic interference forces and moments with a custom flow-through balance, discrete and distributed heatshield pressure, and high-speed flowfield visualization. This paper covers the test objectives, facility, models and instrumentation, and planned test matrix.

Supersonic Retropropulsion↗

Hypercomplex Automatic Differentiation in the Eulerian Hydrocode PAGOSA

Enabling the computation of partial derivatives or sensitivities in production hydrocodes is beneficial for design, optimization, sensitivity analysis, and uncertainty quantification. Traditional finite difference approximations of these sensitivities are inefficient since convergence studies of the step size is required for each parameter of interest. For these reasons, HYPercomplex Automatic Differentiation (HYPAD) was implemented in the Eulerian hydrocode PAGOSA. HYPAD is analogous to forward-mode automatic differentiation except hypercomplex numbers (numbers with multiple imaginary parts) are used instead of dual numbers. Accurate partial derivatives can be computed of all state variables with respect to multiple input variables in a single run. The method was implemented using operator overloading to handle hypercomplex algebra. HYPAD was demonstrated and verified on Sod’s shock tube problem to compute derivatives of the state variables with respect to a material parameter, initial conditions, and geometry.

97 MATHEMATICS AND COMPUTING↗

A computational study of incipient leading-edge separation on a 65-deg delta wing at M = 1.60

A computational study on a 65-deg delta wing at a freestream Mach number of 1.60 has been conducted by obtaining conical Reynolds-averaged Navier-Stokes solutions on a parametric series of geometries which varied in leading-edge radius and/or circular-arc camber. The computational results showed that increasing leading-edge radius or camber can delay the onset of leading-edge separation on the leeside of a delta wing at a specific angle of attack. Reynolds number was varied from 1 x 10 to the 6th to 5 x 10 to the 6th for a turbulent boundary-layer and was shown to have a minor effect on the effectiveness of leading-edge radius and/or camber in delaying the onset of leading-edge separation. Both laminar and turbulent boundary-layer models were investigated at a Reynolds number of 1 x 10 to the 6th, and the predicted flow pattern was found to change from attached flow for the turbulent boundary-layer model to separated flow for the laminar boundary-layer model. Based upon these results, three wind-tunnel models have been designed to be tested in the Langley Unitary Plan Wind Tunnel.

Mcmillin, S. Naomi↗

Approach to nonlinear magnetohydrodynamic simulations in stellarator geometry

The capability to model the nonlinear magnetohydrodynamic (MHD) evolution of stellarator plasmas is developed by extending the M3D-C 1 code to allow non-axisymmetric domain geometry. We introduce a set of logical coordinates, in which the computational domain is axisymmetric, to utilize the existing finite-element framework of M3D-C 1 . A C 1 coordinate mapping connects the logical domain to the non-axisymmetric physical domain, where we use the M3D-C 1 extended MHD models essentially without modifications. We present several numerical verifications on the implementation of this approach, including simulations of the heating, destabilization, and equilibration of a stellarator plasma with strongly anisotropic thermal conductivity, and of the relaxation of stellarator equilibria to integrable and non-integrable magnetic field configurations in realistic geometries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An Investigation of the Interrelationships Between Linear and NonLinear Polarizabilities and Bond Length Alternation on Connugated Organic Molecules

A computational method was devised to explore the relationship between charge separation, geometry, molecular dipole moment, polarizability, and hyperpolarizabilities in comjugated organic molecules. It is shown that bond length alternation (the average difference in length between single and double bonds in the molecule) is a key structurally observable parameter that can be correlated with hyperpolarizabilities and is thus relevant to the optimizaton of molecules and materials.

nonlinear optics computation electric field cyanin↗

Performance Assessment of OVERFLOW on Distributed Computing Environment

The aerodynamic computer code, OVERFLOW, with a multi-zone overset grid feature, has been parallelized to enhance its performance on distributed and shared memory paradigms. Practical application benchmarks have been set to assess the efficiency of code's parallelism on high-performance architectures. The code's performance has also been experimented with in the context of the distributed computing paradigm on distant computer resources using the Information Power Grid (IPG) toolkit, Globus. Two parallel versions of the code, namely OVERFLOW-MPI and -MLP, have developed around the natural coarse grained parallelism inherent in a multi-zonal domain decomposition paradigm. The algorithm invokes a strategy that forms a number of groups, each consisting of a zone, a cluster of zones and/or a partition of a large zone. Each group can be thought of as a process with one or multithreads assigned to it and that all groups run in parallel. The -MPI version of the code uses explicit message-passing based on the standard MPI library for sending and receiving interzonal boundary data across processors. The -MLP version employs no message-passing paradigm; the boundary data is transferred through the shared memory. The -MPI code is suited for both distributed and shared memory architectures, while the -MLP code can only be used on shared memory platforms. The IPG applications are implemented by the -MPI code using the Globus toolkit. While a computational task is distributed across multiple computer resources, the parallelism can be explored on each resource alone. Performance studies are achieved with some practical aerodynamic problems with complex geometries, consisting of 2.5 up to 33 million grid points and a large number of zonal blocks. The computations were executed primarily on SGI Origin 2000 multiprocessors and on the Cray T3E. OVERFLOW's IPG applications are carried out on NASA homogeneous metacomputing machines located at three sites, Ames, Langley and Glenn. Plans for the future will exploit the distributed parallel computing capability on various homogeneous and heterogeneous resources and large scale benchmarks. Alternative IPG toolkits will be used along with sophisticated zonal grouping strategies to minimize the communication time across the computer resources.

Djomehri, M. Jahed↗

Colors of reflection nebulae. I - Phase function effects in the Merope nebula

The subject of color differences between reflection nebulae and their illuminating stars is reexamined in the light of developments of observational techniques, permitting accurate surface-brightness photometry over an expanded spectral region from the UV to the IR. Color-color diagrams for reflection nebulae can yield useful information about the wavelength dependence of the scattering properties of nebular dust without excessive sensitivity to the specific nebular geometry or the presence of multiple scattering, resulting in considerable savings in computational efforts. As an illustration, the color-difference method was applied to existing data for the Merope nebula, covering the spectral region 1550-5500 A. Strong evidence for a monotonically changing phase function of scattering at wavelengths less than or equal to 3500 A is found. The result is interpreted in the context of a plausible geometry for the Merope environment as providing support for a bimodal size distribution of nebular dust grains.

Witt, A. N.↗

Impacts of Biomass Feedstock Pre-Processing on Heat and Mass Transfer During Pyrolysis Using X-Ray Computed Tomography and Multiscale Modeling

Knowledge of the transport properties of biomass particles such as porosity, tortuosity, and permeability is paramount for high-fidelity modeling of biomass pyrolysis due to the heat and mass transfer limitations imposed by particle microstructure. X-ray computed tomography (XCT) is a non-destructive imaging method that enables full 3D reconstructions of the biomass particle microstructure with high resolution, permitting direct calculation of porosity, tortuosity, and permeability from real particle geometries. In this study, XCT imaging revealed the 3D microstructures of particles and chars from pyrolytic conversion of cylindrically cut or milled/pelletized loblolly pine samples. The porosity, tortuosity, and permeability were calculated directly from the XCT geometries via open-source microstructural analysis tool MATBOX+TauFactor (https://github.com/NREL/MATBOX_Microstructure_analysis_toolbox) and computational fluid dynamics (CFD) simulations using our solver, Mesoflow (https://github.com/NREL/mesoflow). These properties were used in a reactor scale model developed in COMSOL of the single particle reactor at NREL to investigate the impact of feedstock pre-processing on biomass conversion during pyrolysis with rigorous experimental validation.

biomass↗

Spur gears: Optimal geometry, methods for generation and Tooth Contact Analysis (TCA) program

The contents of this report include the following: (1) development of optimal geometry for crowned spur gears; (2) methods for their generation; and (3) tooth contact analysis (TCA) computer programs for the analysis of meshing and bearing contact on the crowned spur gears. The method developed for synthesis is used for the determination of the optimal geometry for crowned pinion surface and is directed to reduce the sensitivity of the gears to misalignment, localize the bearing contact, and guarantee the favorable shape and low level of the transmission errors. A new method for the generation of the crowned pinion surface has been proposed. This method is based on application of the tool with a surface of revolution that slightly deviates from a regular cone surface. The tool can be used as a grinding wheel or as a shaver. The crowned pinion surface can also be generated by a generating plane whose motion is provided by an automatic grinding machine controlled by a computer. The TCA program simulates the meshing and bearing contact of the misaligned gears. The transmission errors are also determined.

Litvin, Faydor L.↗

Lightning Radio Source Retrieval Using Advanced Lightning Direction Finder (ALDF) Networks

A linear algebraic solution is provided for the problem of retrieving the location and time of occurrence of lightning ground strikes from an Advanced Lightning Direction Finder (ALDF) network. The ALDF network measures field strength, magnetic bearing and arrival time of lightning radio emissions. Solutions for the plane (i.e., no Earth curvature) are provided that implement all of tile measurements mentioned above. Tests of the retrieval method are provided using computer-simulated data sets. We also introduce a quadratic planar solution that is useful when only three arrival time measurements are available. The algebra of the quadratic root results are examined in detail to clarify what portions of the analysis region lead to fundamental ambiguities in source location. Complex root results are shown to be associated with the presence of measurement errors when the lightning source lies near an outer sensor baseline of the ALDF network. In the absence of measurement errors, quadratic root degeneracy (no source location ambiguity) is shown to exist exactly on the outer sensor baselines for arbitrary non-collinear network geometries. The accuracy of the quadratic planar method is tested with computer generated data sets. The results are generally better than those obtained from the three station linear planar method when bearing errors are about 2 deg. We also note some of the advantages and disadvantages of these methods over the nonlinear method of chi(sup 2) minimization employed by the National Lightning Detection Network (NLDN) and discussed in Cummins et al.(1993, 1995, 1998).

Koshak, William J.↗

Hybrid Curve Fitting For Reducing Motion Commands In Object Construction

Existing slicing software for additive manufacturing typically requires a triangulated mesh as its input. Triangulated meshes are approximate representations of exact CAD models. Despite the loss in dimensional accuracy, triangulated meshes are used because they are computationally easier to cross-section and offset than the exact geometry in CAD format. When a triangulated object is prepared, the resulting machine instructions include only linear motion commands. Numerous modern motion controllers can move in arc and spline motions; however, the absence of slicing software that supports curvature prevents these commands from being leveraged. To address this limitation, this paper presents a method for the hybrid reconstruction of arcs and splines as a post-processing step to traditional slicing. This method can greatly reduce the number of motion commands required to construct an object by printing smooth curved surfaces. This concise representation of tool-pathing allows for more even extrusion and is computed without a major impact on slicing time.

Wade, Charles↗

An efficient numerical model for predicting residual stress and strain in parts manufactured by laser powder bed fusion

Abstract Computational modeling of additively manufactured structures plays an increasingly important role in product design and optimization. For laser powder bed fusion processes, the accurate modeling of stress and distortion requires large amount of computational cost due to very localized heat input and evolving complex geometries. The current study takes advantage of a graphics processing unit accelerated explicit finite element analysis code and approximated heat conduction analysis to predict the macroscopic thermo-mechanical behavior in laser selective melting. Adjacent layers and tracks were lumped to reduce the number of time steps and elements in the finite element model. The effects of track and layer grouping on prediction accuracy and solution efficiency are investigated to provide a guidance for a cost-effective simulation. Thin-wall builds from Inconel alloy 625 (IN625) powders were simulated by applying the developed modeling approach to get the detailed residual stress and distortion at a computational speed 50 times higher than conventional approach. Under repeated heating and cooling cycles, a high tensile stress was produced near surfaces of a build due to a larger shrinkage on surface than that in central area. It is also shown that horizontal stresses concentrate near the root and top layers of the IN625 build. The predicted residual elastic strain distribution was validated by the experimental measurement using x-ray synchrotron diffraction.

36 MATERIALS SCIENCE↗

Scale-Resolving Simulations of Laminar-to-Turbulent Transition in Swept-Wing Flow

Aerodynamic shape optimization of swept Natural Laminar Flow (NLF) wing design plays an important role for the development of future energy-efficient commercial transport aircraft. Lack of predictive laminar-to-turbulent transition modeling capabilities is still limiting the extent to which design can rely on computational tools, while making costly wind tunnel and flight tests a necessity. Prediction of transition on practical 3D wing geometries using the e (sup N) method coupled with laminar boundary layer and linear stability codes or local statistical correlation CFD (Computational Fluid Dynamics) models such as the gamma - Re (sub theta) Langtry-Menter model have proven challenging, whereas transitional hybrid RANS-LES solvers have shown encouraging but still mitigated results.

Scale-resolving↗

Two MCNP Models for Computational Performance Benchmarking

The purpose of this report is to describe two non-trivial MCNP models and execution configurations that are suitable to characterizing computational performance. To that end, this report uses a constructive solid geometry (CSG) representation of the Oak Ridge National Laboratory (ORNL) Pool Critical Assembly (PCA) [1–3] to perform a k-eigenvalue calculation and an unstructured mesh (UM) representation of the International Commission on Radiological Protection (ICRP) publication 145 (ICRP145) male human phantom [4, 5] to perform a fixed-source calculation assuming that a 1 MeV photon source is distributed throughout the phantom’s liver.

97 MATHEMATICS AND COMPUTING↗

Corrections for Convective Heat Flux Gauges Subjected to a Surface Temperature Discontinuity

A two-dimensional Navier-Stokes computational fluid dynamics (CFD) analysis has been carried out in an effort to determine the convective heat transfer corrections for circular heat flux gauges subjected to a surface temperature discontinuity. Solutions were obtained at a Reynolds number Of 1 x 10(exp 6) and a Mach number of 4. The CFD results are compared with the existing correlations for the correction factors. In general, the CFD corrections exceed those provided by the correlations. The discrepancy increases with increasing upstream surface temperature, thus indicating the role of property variations, which are not accounted for in the correlations. A quasi-two-dimensional analysis is also performed to treat the cylindrical geometry of the heat flux gauges by area-averaging the computed two-dimensional results from CFD.

Kandula, M.↗