Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geometry meshing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 865 records · Page 48

Viscous Aerodynamic Shape Optimization with Installed Propulsion Effects

Aerodynamic shape optimization is demonstrated to tailor the under-track pressure signature of a conceptual low-boom supersonic aircraft. Primarily, the optimization reduces nearfield pressure waveforms induced by propulsion integration effects. For computational efficiency, gradient-based optimization is used and coupled to the discrete adjoint formulation of the Reynolds-averaged Navier Stokes equations. The engine outer nacelle, nozzle, and vertical tail fairing are axi-symmetrically parameterized, while the horizontal tail is shaped using a wing-based parameterization. Overall, 48 design variables are coupled to the geometry and used to deform the outer mold line. During the design process, an inequality drag constraint is enforced to avoid major compromise in aerodynamic performance. Linear elastic mesh morphing is used to deform volume grids between design iterations. The optimization is performed at Mach 1.6 cruise, assuming standard day altitude conditions at 51,707-ft. To reduce uncertainty, a coupled thermodynamic engine cycle model is employed that captures installed inlet performance effects on engine operation.

computational fluid dynamics↗

Finite element solution of the Euler equations in two and three dimensions

The paper presents a finite element procedure for solving the equations of compressible flow over bodies of arbitrary geometry. The numerical solution algorithm employed is an explicit two-step version of a second order Taylor-Galerkin scheme. The discretization of the computational domain into unstructured meshes of triangles in two dimensions and tetrahedra in three dimensions is performed by an automatic mesh generator. In the approach, the mesh generator is coupled to the finite element solver to produce an adaptive remeshing procedure.

Peiro, J.↗

First benchmark of the Unstructured Grid Adaptation Working Group

Unstructured grid adaptation is a technology that holds the potential to improve the automation and accuracy of computational fluid dynamics and other computational disciplines. Difficulty producing the highly anisotropic elements necessary for simulation on complex curved geometries that satisfies a resolution request has limited this technology's widespread adoption. The Unstructured Grid Adaptation Working Group is an open gathering of researchers working on adapting simplicial meshes to conform to a metric field. Current members span a wide range of institutions including academia, industry, and national laboratories. The purpose of this group is to create a common basis for understanding and improving mesh adaptation. We present our first major contribution: a common set of benchmark cases, including input meshes and analytic metric specifications, that are publicly available to be used for evaluating any mesh adaptation code. We also present the results of several existing codes on these benchmark cases, to illustrate their utility in identifying key challenges common to all codes and important differences between available codes. Future directions are defined to expand this benchmark to mature the technology necessary to impact practical simulation workflows.

Ibanez, Daniel↗

Formulation and Implementation of Inflow/Outflow Boundary Conditions to Simulate Propulsive Effects

Boundary conditions appropriate for simulating flow entering or exiting the computational domain to mimic propulsion effects have been implemented in an adaptive Cartesian simulation package. A robust iterative algorithm to control mass flow rate through an outflow boundary surface is presented, along with a formulation to explicitly specify mass flow rate through an inflow boundary surface. The boundary conditions have been applied within a mesh adaptation framework based on the method of adjoint-weighted residuals. This allows for proper adaptive mesh refinement when modeling propulsion systems. The new boundary conditions are demonstrated on several notional propulsion systems operating in flow regimes ranging from low subsonic to hypersonic. The examples show that the prescribed boundary state is more properly imposed as the mesh is refined. The mass-flowrate steering algorithm is shown to be an efficient approach in each example. To demonstrate the boundary conditions on a realistic complex aircraft geometry, two of the new boundary conditions are also applied to a modern low-boom supersonic demonstrator design with multiple flow inlets and outlets.

Boundary↗

Sub-Grid Model Improving Shrinking Core Simulations of Redox-Sensitive Technetium Release From Slag Cement Waste Forms - 17295

Numerical simulation of Tc transport in reducing cementitious materials undergoing oxidation is challenging because the effective solid-solute partition coefficient (K{sub d}) varies greatly between the reduced and oxidized states, assuming the Tc mass loading far exceeds that required to reach solubility in the aqueous phase, as is typical in waste disposal. In a conventional finite-volume formulation, the local aqueous concentration increases orders of magnitude when a computational cell switches from a reducing to oxidizing condition, which is an artifact of numerical discretization. In principal such artifacts can be reduced to an insignificant level by refining the computational grid. However, this approach is not feasible in typical facility scale applications due to computer runtime and storage limitations. In this paper an alternative sub-grid Tc release model is derived to approximately capture redox processes occurring at scales smaller than the computational mesh. This alternative approach introduces minimal additional computational expense, while significantly reducing numerical artifacts in the form of Tc concentration / flux spikes. The method is illustrated for two simple geometries and the Savannah River Site Saltstone waste Disposal Facility. The sub-grid model is shown to significantly reduce artificial spikes in Tc-99 release, resulting in more physically realistic numerical simulations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Experimental Investigation on Additively Manufactured Transpiration and Film Cooling Structures

The last 50 years has witnessed significant improvement in film cooling technologies while transpiration cooling is still not implemented in turbine airfoil cooling. Although transpiration cooling could provide higher cooling efficiency with less coolant consumption compared to film cooling, the fine pore structure and high porosity in transpiration cooling metal media always raised difficulties in conventional manufacturing. Recently, the rapid development of additive manufacturing (AM) has provided a new perspective to address such challenge. With the capability of the innovative powder bed selective laser metal sintering (SLMS) AM technology, the complex geometries of transpiration cooling part could be precisely fabricated and endued with improved mechanical strength. This study utilized the SLMS AM technology to fabricate the transpiration cooling and film cooling structures with Inconel 718 superalloy. Five different types of porous media including two perforated plates with different hole pitches, metal sphere packing, metal wire mesh, and blood vessel shaped passages for transpiration cooling were fabricated by EOS M290 system. One laidback fan-shaped film cooling coupon was also fabricated with the same printing process as the control group. Heat transfer tests under three different coolant mass flow rates and four different mainstream temperatures were conducted to evaluate the cooling performance of the printed coupons. The effects of geometry parameters including porosity, surface outlet area ratio, and internal solid–fluid interface area ratio were investigated as well. The results showed that the transpiration cooling structures generally had higher cooling effectiveness than film cooling structure. The overall average cooling effectiveness of blood vessel-shaped transpiration cooling reached 0.35, 0.5, and 0.57, respectively, with low (1.2%), medium (2.4%), and high (3.6%) coolant injection ratios. The morphological parameters analysis showed the major factor that affected the cooling effectiveness most was the internal solid–fluid interface area ratio for transpiration cooling. This study showed that additive manufactured transpiration cooling could be a promising alternative method for turbine blade cooling and worthwhile for further investigations.

Engineering↗

Estimation of metal strength at very high rates using free-surface Richtmyer–Meshkov Instabilities

Recently, Richtmyer–Meshkov Instabilities (RMI) have been proposed for studying the average strength at strain rates up to at least 10 7 /s. RMI experiments involve shocking a metal interface that has initial sinusoidal perturbations. The perturbations invert and grow subsequent to shock and may arrest because of strength effects. In this work we present new RMI experiments and data on a copper target that had five regions with different perturbation amplitudes on the free surface opposite the shock. We estimate the high-rate, low-pressure copper strength by comparing experimental data with Lagrangian numerical simulations. From a detailed computational study we find that mesh convergence must be carefully addressed to accurately compare with experiments, and numerical viscosity has a strong influence on convergence. We also find that modeling the as-built perturbation geometry rather than the nominal makes a significant difference. Because of the confounding effect of tensile damage on total spike growth, which has previously been used as the metric for estimating strength, we instead use a new strength metric: the peak velocity during spike growth. Furthermore, this new metric also allows us to analyze a broader set of experimental results that are sensitive to strength because some larger initial perturbations grow unstably to failure and so do not have a finite total spike growth.

36 MATERIALS SCIENCE↗

An automated data generator for NASTRAN

A modeling system based on construction-in-context of geometry and physical properties data is adapted for NASTRAN finite element data generation. This system was originally developed for generating three-dimensional finite data for composite structures in PATCHES-III using parametric cubic models. These finite line, surface and volume models constructed using simple data directives will be subdivided into NASTRAN finite elements and their supporting data in the new data generator. Any geometry can be constructed as illustrated by a compressor fan blade with camber, twist and changing airfoil section and the shape then subdivided using a uniform or non-uniform mesh. Special emphasis is given to the adaption of construction-in-context to composite material property modeling. A PMAT series of directives is described for any one, two or three-dimensionally reinforced composite that allows the phase properties, as well as composite properties, to be generated.

Stanton, E. L.↗

Modeling of Fillets in Thin-Walled Structures for Dynamic Analysis

Relatively new developments in manufacturing methods have made it possible to produce machined parts with a wall thickness of less than 0.010". While the parts are being machined. Fillets of relatively large radius to thickness ratios are created by the milling tool. This paper discusses an accurate new technique for finite-element modeling of fillets for dynamic and stiffness analysis using efficient plate elements rather than more computationally intensive, high density solid meshes, A simple filleted cantilever beam of 0.040" thickness, 1.6" length, 1" depth, and 0.250" radius in the corner was modeled using both solid elements and several different plate element geometries. The finite element results were then compared with static and modal testing of a machined sample. The highest-scoring plate element technique uses a bridge of elements that is constructed through the tangent point of the Fillet radius such that the thickness of the element matches tile volume of the fillet. This model produces errors of less than 4.6 percent for static loading and less than 7.4 percent for modal analysis. This simple plate element technique will prove critical for efficient, timely, and accurate dynamic analysis of complex thin-walled structures.

Seugling, Richard M.↗

Computation of External Aerodynamics for a Canard Rotor/Wing Aircraft

The aerodynamic loads on a Canard Rotor/Wing vehicle are investigated using inviscid numerical simulations to understand the aerodynamic characteristics of the vehicle during conversion from rotorcraft to fixed wing flight. Steady numerical simulations at five azimuthal rotor indices are presented over a quarter turn of the rotor, producing 19 points during a single rotation due to symmetry. A Cartesian mesh approach is used to compute the steady flow field. All computations are done with a faired over engine inlet and exit to be consistent with the wind tunnel model geometry. Modification to the geometry is suggested and the aerodynamic effect of the modification is discussed.

Pandya, S.↗

Modeling of Biological Fluid Systems

This study presents a detailed methodology for modeling and validating three- dimensional fluid flow in realistic reconstructed vasculature geometry. Optical scanned data from a sectioned stenotic lumen having -65% area constriction within an arterial vessel of 9 mm diameter is processed and converted into a finite-element mesh. Mass and momentum transfer equations within the reconstructed stenotic vessel with a-flow rate of 200 ml/rnin is solved using the Galarkin Finite Element formulation. The calculated velocities at several axial locations along the length of the vessel are compared with the Particle Image Velocimetry (PIV) measurements. The computation results are in agreement with the initial steady state PIV measurements. In addition, initial results from pulsatile flow computations at different time points of a pulse cycle will be presented.

Banerjee, Rupak↗

High-Fidelity Multidisciplinary Design Optimization of Aircraft Configurations

To evaluate new airframe technologies we need design tools based on high-fidelity models that consider multidisciplinary interactions early in the design process. The overarching goal of this NRA is to develop tools that enable high-fidelity multidisciplinary design optimization of aircraft configurations, and to apply these tools to the design of high aspect ratio flexible wings. We develop a geometry engine that is capable of quickly generating conventional and unconventional aircraft configurations including the internal structure. This geometry engine features adjoint derivative computation for efficient gradient-based optimization. We also added overset capability to a computational fluid dynamics solver, complete with an adjoint implementation and semiautomatic mesh generation. We also developed an approach to constraining buffet and started the development of an approach for constraining utter. On the applications side, we developed a new common high-fidelity model for aeroelastic studies of high aspect ratio wings. We performed optimal design trade-o s between fuel burn and aircraft weight for metal, conventional composite, and carbon nanotube composite wings. We also assessed a continuous morphing trailing edge technology applied to high aspect ratio wings. This research resulted in the publication of 26 manuscripts so far, and the developed methodologies were used in two other NRAs. 1

Martins, Joaquim R. R. A.↗

ExaWind: Predictive Wind Energy Simulations

This presentation describes the ExaWind project and the team's progress in creating a suite of performance-portable codes designed for predictive simulations of wind farms on next-generation exascale-class supercomputers. Such simulations will require the resolution of scales spanning many orders of magnitude, from blade boundary layers to wind farm flow structures. In the U.S., the first exascale systems will be GPU accelerated, and different GPU manufacturers have been chosen for the different systems. At the heart of the ExaWind software is a hybrid-solver approach based on the codes Nalu-Wind and AMR-Wind, which are computational fluid dynamics solvers for the incompressible Navier-Stokes equations. Nalu-Wind is an unstructured-grid code used to resolve wind turbine geometry and blade boundary layers, whereas AMR-Wind is a structured-grid background solver for atmospheric turbulent flow and turbine wake propagation. The models are coupled with overset meshes and global linear systems are approximated through a loose-coupling algorithm. Results will include validation-quality high-fidelity simulations and strong/weak scaling results from the Summit supercomputer.

49 EE - Wind and Water Power Program - Wind (EE-4W↗

Novel Solver Algorithms for Nearly Singular Linear Systems Arising in Combustion Modelling

Direct Numerical Simulations of realistic combustion devices are extremely challenging due to the wide separation of scales in the simulation, for example an internal combustion (IC) engine chamber, and the flame thickness of a high-pressure flame. The PeleLMeX solver uses adaptive mesh refinement (AMR) to evolve multi-species reacting flows in the low Mach number limit at the Exascale and relies on an embedded boundary (EB) approach to represent complex geometries. In that framework, the EB geometries often give rise to very small cut-cells along the boundary, which translate into extreme ill-conditioning of the pressure-projection, with eigenvalues that span 15-16 orders of magnitude. In this talk, we focus on the case of a typical IC piston bowl geometry for which we present on a novel approach towards solving these nearly singular linear systems with ILU-based, C-AMG smoothers on massively parallel architectures. In particular, we use scaling and equilibration algorithms to handle the non-normality of the upper triangular factors. This enables us to approximate the highly sequential triangular solve algorithm, embedded in the AMG smoothing-solve phase, with Jacobi iterations. This approximation can be written as a convergent Neumann series whose terms are composed of highly parallel sparse matrix vector multiplications. The result is an algorithm that substantially decreases setup and solve time, compared to state-of-the-art, for these challenging linear systems.

combustion modelling↗

A consistent and conservative scheme for MHD flows with complex boundaries on an unstructured Cartesian adaptive system

The numerical simulation of Magnetohydrodynamics (MHD) flows with complex boundaries has been a topic of great interest in the development of a fusion reactor blanket for the difficulty to accurately simulate the Hartmann layers and side layers along arbitrary geometries. An adaptive version of a consistent and conservative scheme has been developed for simulating the MHD flows. Besides, the present study forms the first attempt to apply the cut-cell approach for irregular wall-bounded MHD flows, which is more flexible and conveniently implemented under adaptive mesh refinement (AMR) technique. It employs a Volume-of-Fluid (VOF) approach to represent the fluid–conducting wall interface that makes it possible to solve the fluid–solid coupling magnetic problems, emphasizing at how electric field solver is implemented when conductivity is discontinuous in cut-cell. For the irregular cut-cells, the conservative interpolation technique is applied to calculate the Lorentz force at cell-center. On the other hand, it will be shown how consistent and conservative scheme is implemented on fine/coarse mesh boundaries when using AMR technique. Then, the applied numerical schemes are validated by five test simulations and excellent agreement was obtained for all the cases considered, simultaneously showed good consistency and conservative properties.

42 ENGINEERING↗

String flash-boiling in gasoline direct injection simulations with transient needle motion

A computational study was performed to investigate the influence of transient needle motion on gasoline direct injection (GDI) internal nozzle flow and near-field sprays. Simulations were conducted with a compressible Eulerian flow solver modeling liquid, vapor, and non-condensable gas phases with a diffuse interface. Variable rate generation and condensation of fuel vapor were captured using the homogeneous relaxation model (HRM). The non-flashing (spray G) and flashing (spray G2) conditions specified by the Engine Combustion Network were modeled using the nominal spray G nozzle geometry and transient needle lift and wobble were based upon ensemble averaged x-ray imaging preformed at Argonne National Lab. The minimum needle lift simulated was 5 μm and dynamic mesh motion was achieved with Laplacian smoothing. The results were qualitatively validated against experimental imaging and the experimental rate of injection profile was captured accurately using pressure boundary conditions and needle motion to actu- ate the injection. Needle wobble was found to have no measurable effect on the flow. Low needle lift is shown to result in vapor generation as fuel rushes past the needle. In conclusion, the internal injector flow is shown to contain many transient and interacting vortices which cause perturbations in the spray angle, fluctuations in the mass flux, and frequently result in string flash-boiling.

33 ADVANCED PROPULSION SYSTEMS↗

Steady RANS methodology for calculating pressure drop in an in-line molten salt compact crossflow heat exchanger

We report the development of molten salt cooled reactors (MSR) and fluoride-salt cooled high temperature reactors (FHR) requires the use of advanced design tools for the primary heat exchanger design. Due to geometric and flow characteristics, compact (pitch to diameter ratios equal to or less than 1.25) heat exchangers with a crossflow flow arrangement can become desirable for these reactors. Unfortunately, the available experimental data is limited for compact tube bundles or banks in crossflow. Computational Fluid Dynamics can be used to alleviate the lack of experimental data in these tube banks. Previous computational efforts have been primarily focused on large S/D ratios (larger than 1.4) using unsteady Reynolds averaged Navier-Stokes and Large Eddy Simulation frameworks. These approaches are useful, but have large computational requirements that make comprehensive design studies impractical. A CFD study was conducted with steady RANS in an effort to provide a starting point for future design work. The study was performed for an in-line tube bank geometry with FLiBe (LiF-BeF2), a frequently selected molten salt, as the working fluid. Based on the estimated pressure drops, the pressure and velocity distributions in the domain, an appropriate meshing strategy was determined and presented. Periodic boundaries in the spanwise direction transverse flow were determined to be an appropriate boundary condition for reduced computational domains. The domain size was investigated and a minimum of 2-flow channels for a domain is recommended to ensure the behavior is accounted for. Finally, the standard low Re κ-ε (Lien) turbulence model was determined to be the most appropriate for steady RANS of this case at the time of writing.

42 ENGINEERING↗