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 703 records · Page 39

The Helicopter Antenna Radiation Prediction Code (HARP)

The first nine months effort in the development of a user oriented computer code, referred to as the HARP code, for analyzing the radiation from helicopter antennas is described. The HARP code uses modern computer graphics to aid in the description and display of the helicopter geometry. At low frequencies the helicopter is modeled by polygonal plates, and the method of moments is used to compute the desired patterns. At high frequencies the helicopter is modeled by a composite ellipsoid and flat plates, and computations are made using the geometrical theory of diffraction. The HARP code will provide a user friendly interface, employing modern computer graphics, to aid the user to describe the helicopter geometry, select the method of computation, construct the desired high or low frequency model, and display the results.

Klevenow, F. T.↗

Mechanical properties of triaxially braided composites - Experimental and analytical results

This paper investigates the unnotched tensile properties of two-dimensional (2D) triaxial braid-reinforced composites both experimentally and analytically. The materials are graphite fibers in an epoxy matrix. Three different reinforcing fiber architectures were considered. There were considerable differences in the observed elastic constants from different size strain gage and extensometer readings. Larger strain gages gave more consistent results and correlated better with the extensometer readings. Experimental strains correlated reasonably well with analytical predictions in the longitudinal, 0 deg, fiber direction but not in the transverse direction. Tensile strength results were not always predictable even in reinforcing directions. Minor changes in braid geometry led to disproportionate strength variations. The unit cell structure of the triaxial braid was discussed with the assistance of computer analysis of the microgeometry. Photomicrographs of braid geometry were used to improve upon the computer graphics representations of unit cells. These unit cells were used to predict the elastic moduli with various degrees of sophistication. The simple and the complex analyses were generally in agreement, but none adequately matched the experimental results for all the braids.

Masters, John E.↗

Mechanical properties of triaxially braided composites: Experimental and analytical results

The unnotched tensile properties of 2-D triaxial braid reinforced composites from both an experimental and an analytical viewpoint are studied. The materials are graphite fibers in an epoxy matrix. Three different reinforcing fiber architectures were considered. Specimens were cut from resin transfer molded (RTM) composite panels made from each braid. There were considerable differences in the observed elastic constants from different size strain gage and extensometer reading. Larger strain gages gave more consistent results and correlated better with the extensometer reading. Experimental strains correlated reasonably well with analytical predictions in the longitudinal, 0 degrees, fiber direction but not in the transverse direction. Tensile strength results were not always predictable even in reinforcing directions. Minor changes in braid geometry led to disproportionate strength variations. The unit cell structure of the triaxial braid was discussed with the assistance of computer analysis of the microgeometry. Photomicrographs of braid geometry were used to improve upon the computer graphics representations of unit cells. These unit cells were used to predict the elastic moduli with various degrees of sophistication. The simple and the complex analyses were generally in agreement but none adequately matched the experimental results for all the braids.

Masters, John E.↗

Reduced Navier Stokes Relaxation Procedures for Internal Flows

In spite of significant advancement in the field of high speed computing, flow calculations involving complex geometries and/or flow behavior still require large amounts of CPU time and memory. In order to predict such flows without sacrificing grid convergence and accuracy, adaptive gridding techniques that provide optimal resolution are highly desirable. The present work combines multigrid techniques and domain decomposition concepts to provide local, solution adaptive, grid refinement. Several viscous compressible and incompressible, two and three-dimensional, flows with strong inviscid interaction and/or axial flow reversal, are considered with a segmented multigrid domain decomposition (SMGDD) procedure for which uniform meshes result in each domain. A pressure-based form of flux-vector splitting is applied to the Navier-Stokes equations, which are represented by an implicit lowest-order reduced Navier-Stokes (RNS) system and a purely diffusive, higher-order, deferred-corrector. A trapezoidal or box-like form of discretization insures that all mass conservation properties are satisfied at interfacial and outflow boundaries, even for this primitive-variable non-staggered grid computation. The SMGDD technique presented herein has previously been applied for incompressible two dimensional flows. The present work offers improvement in the gridding strategy, by allowing for disjoint subdomains that provide optimal resolution of disparate flow features. It also extends the SMGDD technique to three dimensional compressible flows. Laminar and turbulent flow in a backward facing step channel is considered; although the procedure is applicable to more severe geometries. The standard K-epsilon model is applied for turbulence closure. For Re greater than 400, differences between two-dimensional theory and experiment are resolved through a three dimensional simulation, which confirms the experimentally observed three dimensionality of the recirculation patterns on the upper and lower surfaces.

Rubin, Stanley G.↗

Detection of ferromagnetic domain wall pinning and depinning with a semiconductor device

We demonstrate the detection of a ferromagnetic domain wall using a nanoscale Hall cross. A narrow permalloy wire is defined lithographically on top of a Hall cross fabricated from an InAs quantum well. The width of the Hall cross (500 nm–1 μm) is similar to the width of the ferromagnetic wire (200–500 nm), and a geometric pinning site is fabricated in the ferromagnetic wire to trap a domain wall within the area of the Hall cross. The devices provide a signal that is often the same order of magnitude as the offset Hall voltage when a domain wall is located above the Hall cross, and may be useful for memory applications. Different geometries for the Hall cross and ferromagnetic wire are tested, and radiofrequency pulses are sent into the wire to demonstrate current driven domain wall motion. Further changes to the Hall bar geometry with respect to the wire geometry are investigated by numerical computation. A large gain in signal is seen for Hall bars only slightly wider than the ferromagnetic wires as compared to those twice as wide, as well as a larger sensitivity to the exact position of the domain wall with respect to the center of the Hall cross.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

TU-F-18C-05: Evaluation of a Method to Calculate Patient-Oriented MGD Coefficients Using Estimates of Glandular Tissue Distribution

Purpose: Evaluate mammographic mean glandular dose (MGD) coefficients for particular known tissue distributions using a novel formalism that incorporates the effect of the heterogeneous glandular tissue distribution, by comparing them with MGD coefficients derived from the corresponding anthropomorphic computer breast phantom. Methods: MGD coefficients were obtained using MCNP5 simulations with the currently used homogeneous assumption and the heterogeneously-layered breast (HLB) geometry and compared against those from the computer phantom (ground truth). The tissue distribution for the HLB geometry was estimated using glandularity map image pairs corrected for the presence of non-glandular fibrous tissue. Heterogeneity of tissue distribution was quantified using the glandular tissue distribution index, Idist. The phantom had 5 cm compressed breast thickness (MLO and CC views) and 29% whole breast glandular percentage. Results: Differences as high as 116% were found between the MGD coefficients with the homogeneous breast core assumption and those from the corresponding ground truth. Higher differences were found for cases with more heterogeneous distribution of glandular tissue. The Idist for all cases was in the [−0.8{sup −}+0.3] range. The use of the methods presented in this work results in better agreement with ground truth with an improvement as high as 105 pp. The decrease in difference across all phantom cases was in the [9{sup −}105] pp range, dependent on the distribution of glandular tissue and was larger for the cases with the highest Idist values. Conclusion: Our results suggest that the use of corrected glandularity image pairs, as well as the HLB geometry, improves the estimates of MGD conversion coefficients by accounting for the distribution of glandular tissue within the breast. The accuracy of this approach with respect to ground truth is highly dependent on the particular glandular tissue distribution studied. Predrag Bakic discloses current funding from NIH, NSF, and DoD, former funding from Real Time Tomography, LLC and a current research collaboration with Barco and Hologic.

60 APPLIED LIFE SCIENCES↗

Tree Growth--Hybrid Genetic Algorithm for Predicting the Structure of Small (TiO₂) n , n = 2-13, Nanoclusters

The initial structures for the search for the global minimum of TiO₂ nanoclusters were generated by combining a tree growth (TG) algorithm with a hybrid genetic algorithm (HGA). In the TG algorithm, the clusters grow from a small seed to the size of interest stepwise. New atoms are added to the smaller cluster from the previous step, by analogy to new leaves grown by a tree. The addition of the new atoms is controlled by predefined geometry parameters to reduce the computational cost and to provide physically meaningful structures. In each step, the energies for the various generated structures are evaluated, and those with the lowest energies are carried into the next step. The structures that match the formulas of interest are collected as HGA candidates during the various steps. Low energy candidates are fed to the HGA component to search for the global minimum for each formula of interest. The lowest energy structures from the HGA are then optimized by using density functional theory to study the dissociation energies of the clusters and the evolution in the structure as the size of the cluster increases. The optimized geometries of the (TiO₂)n nanoclusters for n = 2-13, do not show the character of a TiO₂ bulk crystal with a hexacoordinate Ti. The average clustering energy (ΔEn) converges slowly to the bulk value for rutile. The TiO₂ dissociation energies for (TiO₂)n clusters approach the bulk value for rutile more quickly but show larger variations. The (TiO₂)12 cluster appears to be quite stable, and the (TiO₂)13 cluster is quite unstable on a relative scale.

Environmental Molecular Sciences Laboratory↗

Computer graphics in aerodynamic analysis

The use of computer graphics and its application to aerodynamic analyses on a routine basis is outlined. The mathematical modelling of the aircraft geometries and the shading technique implemented are discussed. Examples of computer graphics used to display aerodynamic flow field data and aircraft geometries are shown. A future need in computer graphics for aerodynamic analyses is addressed.

Cozzolongo, J. V.↗

Dynamic analysis of astronaut motions in microgravity: Applications for Extravehicular Activity (EVA)

Simulations of astronaut motions during extravehicular activity (EVA) tasks were performed using computational multibody dynamics methods. The application of computational dynamic simulation to EVA was prompted by the realization that physical microgravity simulators have inherent limitations: viscosity in neutral buoyancy tanks; friction in air bearing floors; short duration for parabolic aircraft; and inertia and friction in suspension mechanisms. These limitations can mask critical dynamic effects that later cause problems during actual EVA's performed in space. Methods of formulating dynamic equations of motion for multibody systems are discussed with emphasis on Kane's method, which forms the basis of the simulations presented herein. Formulation of the equations of motion for a two degree of freedom arm is presented as an explicit example. The four basic steps in creating the computational simulations were: system description, in which the geometry, mass properties, and interconnection of system bodies are input to the computer; equation formulation based on the system description; inverse kinematics, in which the angles, velocities, and accelerations of joints are calculated for prescribed motion of the endpoint (hand) of the arm; and inverse dynamics, in which joint torques are calculated for a prescribed motion. A graphical animation and data plotting program, EVADS (EVA Dynamics Simulation), was developed and used to analyze the results of the simulations that were performed on a Silicon Graphics Indigo2 computer. EVA tasks involving manipulation of the Spartan 204 free flying astronomy payload, as performed during Space Shuttle mission STS-63 (February 1995), served as the subject for two dynamic simulations. An EVA crewmember was modeled as a seven segment system with an eighth segment representing the massive payload attached to the hand. For both simulations, the initial configuration of the lower body (trunk, upper leg, and lower leg) was a neutral microgravity posture. In the first simulation, the payload was manipulated around a circular trajectory of 0.15 m radius in 10 seconds. It was found that the wrist joint theoretically exceeded its ulnal deviation limit by as much as 49. 8 deg and was required to exert torques as high as 26 N-m to accomplish the task, well in excess of the wrist physiological limit of 12 N-m. The largest torque in the first simulation, 52 N-m, occurred in the ankle joint. To avoid these problems, the second simulation placed the arm in a more comfortable initial position and the radius and speed of the circular trajectory were reduced by half. As a result, the joint angles and torques were reduced to values well within their physiological limits. In particular, the maximum wrist torque for the second simulation was only 3 N-m and the maximum ankle torque was only 6 N-m.

Newman, Dava J.↗

ParaStell: parametric modeling and neutronics support for stellarator fusion power plants

The three-dimensional variation inherent to stellarator geometries and fusion sources motivates three-dimensional modeling to obtain accurate results from computational modeling in support of design and analysis of first wall, blanket, and shield (FWBS) systems. Manually constructing stellarator fusion power plant geometries in computer-aided design (CAD) and defining the corresponding fusion source can be cumbersome and challenging. The open-source parametric modeling toolset ParaStell has been developed to automate construction of such geometries in low-fidelity. Low-fidelity modeling is useful during the conceptual phase of engineering design as a means of rapidly exploring the design space of a given device. The modeling capability of ParaStell includes in-vessel components and magnets, for any given stellarator configuration, using a parametric definition and plasma equilibrium data. Furthermore, the toolset automates the generation of detailed, tetrahedral neutron source definitions and DAGMC geometries for use in neutronics modeling. ParaStell assists rapid design iteration, parametric study, and design optimization of stellarator fusion cores. As a demonstration of the design iteration capability, the effect of the three-dimensional parameter space on tritium breeding and magnet shielding is investigated, using the WISTELL-D configuration as a design basis. Blanket and shield thicknesses are varied in three dimensions, using the space available between the plasma edge and magnet coils as a constraint. The corresponding effects on tritium breeding ratio and magnet heating are tallied using the open-source Monte Carlo particle transport code OpenMC. The inclusion of additional and higher-fidelity modeling capabilities is planned for ParaStell’s future, as well as its implementation in machine-driven optimization.

Moreno, Connor↗

A Method for Flow Simulation About Complex Geometries Using Both Structured and Unstructured Grids

A computational fluid dynamics code which utilizes both structured and unstructured grids was developed. The objective of this study was to develop and demonstrate the ability of such a code to achieve solutions about complex geometries in two dimensions. An unstructured grid generator and flow solver were incorporated into the PARC2D structured flow solver. This new unstructured grid generator capability allows for easier generation and manipulation of complex grids. Several examples of the grid generation capabilities are provided. The coupling of different grid topologies and the manipulation of individual grids is shown. Also, grids for realistic geometries, a NACA 0012 airfoil and a wing/nacelle installation, were created. The flow over a NACA 0012 airfoil was used as a test case for the flow solver. Eight separate cases were run. They were both the inviscid and viscous solutions for two freestream Mach numbers and airfoil angle of attacks of 0 to 3.86 degrees. The Mach numbers chosen were for a subsonic case, Mach 0.6, and a case where supersonic regions and a shock wave exists, Mach 0.8. These test case conditions were selected to match experimentally obtained data for code comparison. The results show that the code accurately predicts the flow field for all cases.

Debonis, James R.↗

Maximizing the Performance of a 3D Printed Heat Sink by Accounting for Anisotropic Thermal Conductivity During Filament Deposition

Newly developed high thermal conductivity polymer composite 3D printing filaments are used to characterize the thermal properties as a function of print orientation. The thermal conductivity of a printed part is anisotropic and varies by 2-6 times depending on the print direction - demonstrating higher conductivity in the deposition direction (in-plane) than in the two directions perpendicular to the deposition direction (cross-plane and through-plane). Therefore, deposition path planning greatly affects the overall heat dissipation rate and the performance of the heat sink. Traditionally, 3D printing slicers generate deposition paths based solely on geometric constraints. This work investigates a new approach of deposition path planning assisted by computational predictions of the heat sink thermal performance. The proposed approach uses a thermal simulation of a 3D-printed part, accounting for the anisotropic thermal properties, and the orientation of the local material properties are assigned based on the deposition path in multiple print orientations. The performances predicted via the simulations are compared, and the optimal deposition path is determined. For the highest thermal conductivity 3D printing filament (~ 12 W/m- K in-plane), a heat sink printed with the print direction parallel to the fins z-axis had ~20% improved performance in comparison to a heat sink with print direction perpendicular to the fins z-axis. Moreover, a plastic 3D printed heat sink was able to perform within 7% of an extruded Aluminum heat sink with similar geometry under natural convection. The computational predictions show the same trend as experimental measurements using 3D printed heat sinks.

Smith, Matthew↗

Open quantum systems and Schwinger-Keldysh holograms

We initiate the study of open quantum field theories using holographic methods. Specifically, we consider a quantum field theory (the system) coupled to a holographic field theory at finite temperature (the environment). We investigate the effects of integrating out the holographic environment with an aim of obtaining an effective dynamics for the resulting open quantum field theory. The influence functionals which enter this open effective action are determined by the real-time (Schwinger-Keldysh) correlation functions of the holographic thermal environment. To evaluate the latter, we exploit recent developments, wherein the semiclassical gravitational Schwinger-Keldysh saddle geometries were identified as complexified black hole spacetimes. We compute real-time correlation functions using holographic methods in these geometries, and argue that they lead to a sensible open effective quantum dynamics for the system in question, a question that hitherto had been left unanswered. In addition to shedding light on open quantum systems coupled to strongly correlated thermal environments, our results also provide a principled computation of Schwinger-Keldysh observables in gravity and holography. In particular, these influence functionals we compute capture both the dissipative physics of black hole quasi- normal modes, as well as that of the fluctuations encoded in outgoing Hawking quanta, and interactions between them. We obtain results for these observables at leading order in a low frequency and momentum expansion in general dimensions, in addition to determining explicit results for two dimensional holographic CFT environments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

PROTEUS two-dimensional Navier-Stokes computer code, version 1.0. Volume 1: Analysis description

A new computer code was developed to solve the two-dimensional or axisymmetric, Reynolds averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The thin-layer or Euler equations may also be solved. Turbulence is modeled using an algebraic eddy viscosity model. The objective was to develop a code for aerospace applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The equations are written in nonorthogonal body-fitted coordinates, and solved by marching in time using a fully-coupled alternating direction-implicit procedure with generalized first- or second-order time differencing. All terms are linearized using second-order Taylor series. The boundary conditions are treated implicitly, and may be steady, unsteady, or spatially periodic. Simple Cartesian or polar grids may be generated internally by the program. More complex geometries require an externally generated computational coordinate system. The documentation is divided into three volumes. Volume 1 is the Analysis Description, and describes in detail the governing equations, the turbulence model, the linearization of the equations and boundary conditions, the time and space differencing formulas, the ADI solution procedure, and the artificial viscosity models.

Towne, Charles E.↗

PROTEUS two-dimensional Navier-Stokes computer code, version 1.0. Volume 2: User's guide

A new computer code was developed to solve the two-dimensional or axisymmetric, Reynolds averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The thin-layer or Euler equations may also be solved. Turbulence is modeled using an algebraic eddy viscosity model. The objective was to develop a code for aerospace applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The equations are written in nonorthogonal body-fitted coordinates, and solved by marching in time using a fully-coupled alternating direction-implicit procedure with generalized first- or second-order time differencing. All terms are linearized using second-order Taylor series. The boundary conditions are treated implicitly, and may be steady, unsteady, or spatially periodic. Simple Cartesian or polar grids may be generated internally by the program. More complex geometries require an externally generated computational coordinate system. The documentation is divided into three volumes. Volume 2 is the User's Guide, and describes the program's general features, the input and output, the procedure for setting up initial conditions, the computer resource requirements, the diagnostic messages that may be generated, the job control language used to run the program, and several test cases.

Towne, Charles E.↗

PROTEUS two-dimensional Navier-Stokes computer code, version 1.0. Volume 3: Programmer's reference

A new computer code was developed to solve the 2-D or axisymmetric, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The thin-layer or Euler equations may also be solved. Turbulence is modeled using an algebraic eddy viscosity model. The objective was to develop a code for aerospace applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The equations are written in nonorthogonal body-fitted coordinates, and solved by marching in time using a fully-coupled alternating-direction-implicit procedure with generalized first- or second-order time differencing. All terms are linearized using second-order Taylor series. The boundary conditions are treated implicitly, and may be steady, unsteady, or spatially periodic. Simple Cartesian or polar grids may be generated internally by the program. More complex geometries require an externally generated computational coordinate system. The documentation is divided into three volumes. Volume 3 is the Programmer's Reference, and describes the program structure, the FORTRAN variables stored in common blocks, and the details of each subprogram.

Towne, Charles E.↗

Comparison Study of Conventional and Adaptive Mesh Refinement in Organic Material Decomposition Models

This study compares conventional mesh refinement techniques, specifically Uniform Mesh Refinement (UMR), with a new Adaptive Mesh Refinement (AMR) method, applied to Organic Material Decomposition (OMD) models. The proposed benefit of AMR is that only areas that require refinement, based on minimizing a specific field gradient, are refined thus decreasing model wall time compared to conventional UMR methods. This work specifically focuses on comparing UMR and AMR methods on decomposing (both No-Flow and Porous-Flow material models) Polymeric Methylene Diisocyanate (PMDI) polyurethane foam. Throughout the work, the geometry increased in complexity to assess the refinement methods performance at varying levels geometric intricacy. While AMR has been shown to work well in a variety of applications, the UMR approach proved to be computationally faster, for many of the geometries and foam decomposition models, than AMR. However, it was observed that at higher levels of refinement, greater than 3 UMR, AMR begins to be computationally better. Additionally, the settings used to perform AMR greatly impact its performance, and lessons learned, in terms of OMD models, are shared. Due to physics involved in material decomposition, specifically the evolution of state variables, these problems don’t fully benefit from the advantages of AMR.

36 MATERIALS SCIENCE↗

Computational Study of Droplet Dynamics in Dispersed Flow Film Boiling Regime in PWR Geometries

Dispersed flow film boiling (DFFB) regime is exhibited during low-flooding rate conditions following the postulated Loss-of-Coolant Accident (LOCA) in nuclear reactor flow channels. LOCA is a type of design-basis accident scenario, which is crucial for reactor design and licensing. Accurate models to predict the heat transfer in DFFB regime is critical for the reactor to maintain safety limits. The nominal operating linear heat rate is also set based on the maximum linear heat rate determined for a LOCA transient. It is understood that the behavior of entrained droplets in DFFB regime greatly affects the heat transfer by evaporating and desuperheating the surrounding steam. Furthermore, interaction with flow structures also aids in increasing the evaporation rate of entrained droplets. In order to quantify this enhancement in heat transfer, the size of droplets downstream of flow structures (Spacer Grids and Mixing Vanes) is of crucial importance. Thermal hydraulic analysis codes like COBRA-TF, CATHARE-3 and adopt a three-field equation model that allows them to model the droplet field separately. However, to ensure the accuracy of these field equations, various correlations and physical models are required. This paper proposes a method for characterizing the mean Sauter Diameter of droplets across spacer-grids and mixing vanes in PWR sub-channels using PHASTA - a high fidelity flow solver with multiphase capabilities using the level-set method. The PHASTA code is a finite element based computational fluid dynamics code designed for solving compressible and incompressible turbulent flows. It has been demonstrated to be massively scalable (up to 786,432 cores), which is crucial for simulating large realistic domains. The code incorporates the level-set method for interface tracking of the gas-liquid interface.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗