Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mesh Optimization”

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 451 records · Page 25

On the design of contractions and settling chambers for optimal turbulence manipulations in wind tunnels

The sensitivity of a moderate axisymmetric contraction of turbulence to upstream conditions is parametrically investigated in terms of turbulence intensities and integral length scales. Semi-empirical correlations are derived to characterize the development of turbulence intensities and integral length scales downstream of typical wind-tunnel turbulence manipulators. A new approach is described, permitting wind-tunnel designers to select the characteristic mesh of the turbulence manipulators, a well as the distance from them to the contraction, in order to obtain certain required turbulence characteristics in the test section at the exit of the contraction. The design charts strongly suggest that it may not be possible to achieve small scales and low intensities simultaneously, when using a contraction ratio of nine or larger. The results also confirm that the rapid distortion theory is only valid, with a proper viscous dissipation correction, up to the 'critical point' of a contraction. At that point, the scales of the various components of the velocity are equal and intercomponent transfer from the lateral to the streamwise velocities commences.

Nagib, H. M.↗

Development of aircraft lavatory compartments with improved fire resistance characteristics. Phase 4: Sandwich panel decorative ink development

Five chemically different resin systems with improved fire resistance properties were studied for a possible screenprinting ink application. Fire resistance is hereby defined as the cured ink possessing improvements in flammability, smoke emission, and thermal stability. The developed ink is for application to polyvinyl fluoride film. Only clear inks without pigments were considered. Five formulations were evaluated compared with KC4900 clear acrylic ink, which was used as a baseline. The tests used in the screening evaluation included viscosity, smoke and toxic gas emission, limiting oxygen index (LOI), and polyvinyl fluoride film (PVF) printability. A chlorofluorocarbon resin (FPC461) was selected for optimization studies. The parameters for optimization included screenprinting process performance, quality of coating, and flammability of screenprinted 0.051-mm (0.002-in.) white Tedlar. The quality of the screenprinted coating on Tedlar is dependent on viscosity, curing time, adhesion to polyvinyl fluoride film, drying time (both inscreen and as an applied film), and silk screen mesh material and porosity.

Jayarajan, A.↗

Matrix-free preconditioning for high-order H (curl) discretizations

The greater arithmetic intensity of high-order finite element discretizations makes them attractive for implementation on next-generation hardware, but assembly of high-order finite element operators as matrices is prohibitively expensive. As a result, the development of general algebraic solvers for such operators has been an open research challenge. Fast matrix-free application of high-order operators has received significant attention in the literature in the context of Poisson-type problems, but preconditioners and solvers for inverting more general operators are not very well-developed. In this paper, we consider the problem of preconditioning a definite Maxwell operator at high polynomial order without assembling a matrix. We show that given efficient preconditioners for high-order H 1 finite element problems on the same mesh, efficient H(curl) preconditioners can be constructed in an auxiliary space framework. We demonstrate the resulting preconditioners in a practical setting with tensor-product basis functions on an unstructured mesh of quadrilaterals. Overall, our approach uses a sparsified H 1 solver constructed on a low-order mesh of the nodal points of the underlying high-order space, and we show that the resulting H(curl) preconditioner is effective at very high polynomial orders for two-dimensional model problems with complicated geometry, varying piecewise constant coefficients, and curved elements. The resulting preconditioner scales with nearly optimal O(p d+1 ) floating point operation count and optimal O(p d ) memory transfer requirements, outperforming existing Maxwell preconditioners in the high-order regime.

97 MATHEMATICS AND COMPUTING↗

A neural network‐enhanced reproducing kernel particle method for modeling strain localization

Abstract Modeling the localized intensive deformation in a damaged solid requires highly refined discretization for accurate prediction, which significantly increases the computational cost. Although adaptive model refinement can be employed for enhanced effectiveness, it is cumbersome for the traditional mesh‐based methods to perform while modeling the evolving localizations. In this work, neural network‐enhanced reproducing kernel particle method (NN‐RKPM) is proposed, where the location, orientation, and shape of the solution transition near a localization is automatically captured by the NN approximation via a block‐level neural network (NN) optimization. The weights and biases in the blocked parameterization network control the location and orientation of the localization. The designed basic four‐kernel NN block is capable of capturing a triple junction or a quadruple junction topological pattern, while more complicated localization topological patters are captured by the superposition of multiple four‐kernel NN blocks. The standard RK approximation is then utilized to approximate the smooth part of the solution, which permits a much coarser discretization than the high‐resolution discretization needed to capture sharp solution transitions with the conventional methods. A regularization of the NN approximation is additionally introduced for discretization‐independent material responses. The effectiveness of the proposed NN‐RKPM is verified by a series of numerical verifications.

Baek, Jonghyuk↗

Evaluating pulse-shaping capabilities of next-generation pulsed power architectures

This project evaluated the pulse shaping capabilities of next-generation pulsed power (NGPP) architectures. NGPP architectures share several common attributes including multiple independent pulse-generation lines, a radial water-insulated impedance transformer, and a central vacuum insulated load region. A multi-module circuit model was developed, incorporating independent pulse-generation lines and a 2-D transmission line mesh of the radial impedance transformer to assess the effects of azimuthal asymmetry in pulse-shaped experiments. Circuit model simulations demonstrated that NGPP architectures are able to produce the the desired current pulse shapes for exemplar NGPP experiments. Additionally, the project explored automated methods for experiment design, including derivative -ree optimization and machine learning. Pulse-shaped experiments require designers to determine machine parameters that reliably produce the desired current pulse at the load, a process that typically relies on expert knowledge and iterative adjustments using the Z circuit model. Given the increased complexity of NGPP systems, this manual approach may be impractical. While the evaluated methods do not eliminate the need for manual iteration, they can reduce the time required for experiment design. Derivative-free optimization automates much of the trial-and-error process, providing a close starting point for manual adjustments or making small modifications to near-final designs. Meanwhile, deep neural network methods can generate a good qualitative match to the desired current pulse in under one second without requiring circuit model simulations.

42 ENGINEERING↗

Meeting the challenges with the Douglas Aircraft Company Aeroelastic Design Optimization Program (ADOP)

An overview of the Aeroelastic Design Optimization Program (ADOP) at the Douglas Aircraft Company is given. A pilot test program involving the animation of mode shapes with solid rendering as well as wire frame displays, a complete aircraft model of a high-altitude hypersonic aircraft to test ADOP procedures, a flap model, and an aero-mesh modeler for doublet lattice aerodynamics are discussed.

Rommel, Bruce A.↗

The Design and Evaluation of "CAPTools"--A Computer Aided Parallelization Toolkit

Writing applications for high performance computers is a challenging task. Although writing code by hand still offers the best performance, it is extremely costly and often not very portable. The Computer Aided Parallelization Tools (CAPTools) are a toolkit designed to help automate the mapping of sequential FORTRAN scientific applications onto multiprocessors. CAPTools consists of the following major components: an inter-procedural dependence analysis module that incorporates user knowledge; a 'self-propagating' data partitioning module driven via user guidance; an execution control mask generation and optimization module for the user to fine tune parallel processing of individual partitions; a program transformation/restructuring facility for source code clean up and optimization; a set of browsers through which the user interacts with CAPTools at each stage of the parallelization process; and a code generator supporting multiple programming paradigms on various multiprocessors. Besides describing the rationale behind the architecture of CAPTools, the parallelization process is illustrated via case studies involving structured and unstructured meshes. The programming process and the performance of the generated parallel programs are compared against other programming alternatives based on the NAS Parallel Benchmarks, ARC3D and other scientific applications. Based on these results, a discussion on the feasibility of constructing architectural independent parallel applications is presented.

Yan, Jerry↗

Effects of regenerator geometry on pulse tube refrigerator performance

This paper gives results of the cooling performance of a double-inlet pulse tube refrigerator using various regenerators. The same pulse tube was used for all the experiments and measured 4.76 mm in diameter and 46.2 mm in length. A commercial linear compressor with a swept volume of 4 cm3 was used in these experiments. The operating conditions were held constant at a mean pressure of 2.0 MPa and a frequency of 54 Hz. Using finite difference software called REGEN3.1, developed at NIST, and recent experiment results, we optimized a series of regenerators based on dimensions, materials and screen packing. The values used for calculating the thermal conduction through stacked screens by REGEN3.1 were based on recent experimental results from NIST. The regenerator tubes were designed using 316 stainless steel and titanium materials. The regenerator matrices investigated were 400-mesh and 500-mesh stainless steel screen. The valve settings for both orifices were adjusted to minimize the no-load temperature for all regenerators. A cooling capacity curve from 0 to 3 W was also determined. The performance of the pulse tube refrigerator using the different regenerators is discussed. The experimental results from the various regenerators are evaluated and compared with their corresponding numerically calculated coefficient of performance (COP) and regenerator design as determined by REGEN3.1.

NASA Discipline Life Sciences Technologies↗

A Comprehensive Structural Dynamic Analysis Approach for Multi Mission Earth Entry Vehicle (MMEEV) Development

The anticipated NASA Mars Sample Return Mission (MSR) requires a simple and reliable method in which to return collected Martian samples back to earth for scientific analysis. The Multi-Mission Earth Entry Vehicle (MMEEV) is NASA's proposed solution to this MSR requirement. Key aspects of the MMEEV are its reliable and passive operation, energy absorbing foam-composite structure, and modular impact sphere (IS) design. To aid in the development of an EEV design that can be modified for various missions requirements, two fully parametric finite element models were developed. The first model was developed in an explicit finite element code and was designed to evaluate the impact response of the vehicle and payload during the final stage of the vehicle's return to earth. The second model was developed in an explicit code and was designed to evaluate the static and dynamic structural response of the vehicle during launch and reentry. In contrast to most other FE models, built through a Graphical User Interface (GUI) pre-processor, the current model was developed using a coding technique that allows the analyst to quickly change nearly all aspects of the model including: geometric dimensions, material properties, load and boundary conditions, mesh properties, and analysis controls. Using the developed design tool, a full range of proposed designs can quickly be analyzed numerically and thus the design trade space for the EEV can be fully understood. An engineer can then quickly reach the best design for a specific mission and also adapt and optimize the general design for different missions.

Perino, Scott↗

Machine Learning Assisted Safety Modeling and Analysis of Advanced Reactors

With the advances in computational power and numerical methods, analysts can now rely on first-principle simulations to predict ultra-fine details in a variety of applications. Advances in machine learning (ML) have produced algorithms that can now learn high-level abstractions via hierarchical models. This project aims to leverage advances in ML techniques and the available high-resolution simulation data to develop a novel modeling and simulation (M\&S) methodology for reactor safety analysis. While application-agnostic ML techniques are available, complex physics constraints need to be incorporated into ML techniques to build ML-based closures for computationally efficient predictive simulations. This project intends to develop a physics-guided data-driven multi-scale methodology for M\&S of advanced reactors. The project focuses on thermal fluid (T/F) phenomena, which play major roles in advanced reactor safety. Specifically, we propose a data-driven coarse-mesh turbulence model based on local flow features for the transient analysis of thermal mixing and stratification in a sodium-cooled fast reactor (SFR). The model has a coarse-mesh setup to ensure computational efficiency, while it is trained by fine-mesh computational fluid dynamics (CFD) data with Reynolds-averaged Navier-Stokes (RANS) turbulence model to ensure accuracy. Three different neural networks are developed and tested for loss-of-flow transients in the hot pool of SFR, i.e. the densely connected convolutional neural network (DCNN), long-short-term-memory network based on proper orthogonal decomposition (POD-LSTM), and the DCNN informed by LSTM (DCNN-LSTM). The performances of these three neural networks are evaluated based on baseline models. The DCNN-LSTM model has been chosen for further hyperparameter optimization. Furthermore, based on a simplified two-dimensional case, uncertainty quantification (UQ) of the developed ML-based closure are investigated with three methods, i.e. Monte Carlo dropout, deep ensemble, and Bayesian neural network. The developed ML-based turbulent viscosity closure relation based on deep ensemble is then integrated into the system analysis module SAM and serves as a term in the conservation equations. Such a SAM-ML based procedure guarantees that the obtained results are consistent with the physical constraints of the thermal-fluid system. The SAM-ML simulation on the same loss-of-flow transient showed comparable accuracy with the CFD simulation but with a much coarser mesh setup. Last but not least, the ML-based closure improvement with the support of higher-fidelity data from large eddy simulation (LES) is discussed. As a first step towards this direction, a baseline LES simulation is performed to obtain comparable data with RANS results. Based on the early results, future investigation on further improving the ML-based closure is discussed. We believe the developed approach that combines scientific machine learning with nuclear system analysis code can benefit the advanced reactor community as more accurate safety analyses will better characterize reactor safety margins and reduce licensing efforts.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Foundations for high-order, conservative cut-cell methods: Stable discretizations on degenerate meshes

Cut-cell methods for unsteady flow problems can greatly simplify the grid generation process and allow for high-fidelity simulations on complex geometries. However, cut-cell methods have been limited to low orders of accuracy. This is driven, largely, by the variety of procedures typically introduced to evaluate derivatives in a stable manner near the highly irregular embedded geometry. Here, a completely new approach, termed TEMO (truncation error matching and optimization), is taken to solve this problem. The approach is based on two simple and intuitive design principles. These principles directly allow for the construction of stable 8th To the best of the authors' knowledge, these are the highest orders ever achieved for a cut-cell discretization by a significant margin. This is done for both explicit and compact finite differences and is accomplished without any geometric transformations or artificial stabilization procedures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Physical Model for Understanding the Activation of MoS 2 Basal-plane Sulfur Atoms for the Hydrogen Evolution Reaction

All the DFT calculations are done with the Vienna Ab Initio Simulation Package (VASP) using the projector augmented wave method. The Bayesian error estimation exchange-correlation functionals (BEEF) with van der Waals interactions are employed. A plane-wave cutoff energy of 400 eV is used together with PAW-PBE potentials where semi core p states are treated as valence. All the calculations allow for spin-polarization. The structures are relaxed until the force is converged to < 0.01 eV/Å. The lattice parameter of MoS 2 unit cell, optimized with this functional, is 3.19 Å. A (4×4) supercell is used to model all the transition metal doped MoS 2 systems studied here, including those with S vacancies. For calculations in the initial dopant structure exploration, the Brillouin zone is sampled with a 3x3x1 Monkhorst-Pack k-point mesh. A 6×6×1 Monkhorst-Pack k-point mesh is used for the H binding energy and density of states calculations. In all calculations, the vacuum layer is set as 15 Å to eliminate periodic interaction perpendicular to the basal plane.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Modeling Gravitational Radiation Waveforms from Black Hole Mergers

Gravitational radiation from merging binary black hole systems is anticipated as a key source for gravitational wave observations. Ground-based instruments, such as the Laser Interferometer Gravitational-wave Observatory (LIGO) may observe mergers of stellar-scale black holes, while the space-based Laser Interferometer Space Antenna (LISA) observatory will be sensitive to mergers of massive galactic-center black holes over a broad range of mass scales. These cataclysmic events may emit an enormous amount of energy in a brief time. Gravitational waves from comparable mass mergers carry away a few percent of the system's mass-energy in just a few wave cycles, with peak gravitational wave luminosities on the order of 10^23 L_Sun. Optimal analysis and interpretation of merger observation data will depend on developing a detailed understanding, based on general relativistic modeling, of the radiation waveforms. We discuss recent progress in modeling radiation from equal mass mergers using numerical simulations of Einstein's gravitational field equations, known as numerical relativity. Our simulations utilize Adaptive Mesh Refinement (AMR) to allow high-resolution near the black holes while simultaneously keeping the outer boundary of the computational domain far from the black holes, and making it possible to read out gravitational radiation waveforms in the weak-field wave zone. We discuss the results from simulations beginning with the black holes orbiting near the system's innermost stable orbit, comparing the recent simulations with earlier "Lazarus" waveform estimates based on an approximate hybrid numerical/perturbative technique.

Baker, J. G.↗

A mixed-integer PDE-constrained optimization formulation for constructing electromagnetic cloaks with multiple materials

We study the design of an electromagnetic cloak from multiple materials with an additional constraint on the mass of the cloak. Our problem is an example of a topology optimization problem, and we formulate this problem as a mixed-integer partial-differential equation constrained optimization (MIPDECO) problem, where Maxwell’s equation models the propagation of the wave through the cloak and surrounding medium. We use binary variables to model the assignment of the different materials, and their relevant properties (permittivity and density). The mass constraint adds a nontrivial constraint to this problem. We propose a two-phase strategy to solve this problem. In the first phase, we solve a continuous relaxation, and then propose a new variant of the feasibility pump that exploits the structure of the PDE to obtain an initial integral solution candidate. In the second phase, we use a trust-region approach to improve this incumbent. We also consider a continuation or mesh-sequencing approach to find better solutions faster on consecutively finer meshes. We present detailed numerical results to illustrate the effectiveness of our approaches for constructing multi-material cloaks with a mass constraint.

Calculus of Variations and Optimization↗

Rate of Para-Orthohydrogen Conversion in Cryogenic Vortex Tube

Raising efficiency of the cooling process for cryogenic hydrogen and minimizing hydrogen boil-off during storage and transportation of liquid hydrogen are critically important factors for widening implementation of hydrogen-based clean energy systems. One novel approach involves utilization of relatively simple catalyzed vortex tubes, where hydrogen in the primarily para- nuclear spin isomer form is converted into the ortho- form while consuming a significant amount of heat. To design such systems effectively, better understanding is required about catalyst-assisted para-orthohydrogen conversion rates in high-speed vortical flows. In this study, a computational fluid dynamics simulation has been set up to model an experimental system with cryogenic hydrogen vortex tubes. Mesh-verification and a validation study has been conducted first for a non-catalyzed tube. Then, the rate coefficient of the para-orthohydrogen conversion of cryogenic hydrogen has been determined by matching numerical results with experimental data available for a catalyzed vortex tube. Furthermore, the presented information can help design and optimize cryogenic hydrogen cooling within vortex tubes.

08 HYDROGEN↗

Eight cm technology thruster development

A structural integrated ion thruster with 8-cm beam diameter (SIT-8) was developed for attitude control and stationkeeping of synchronous satellites. As optimized, the system demonstrates a thrust T=1.14 mlb (not corrected for beam V sub B = 1200 V (I sub sp = 2200 sec) total propellant utilization efficiency nu sub u = 59.8% (is approximately 72% without auxiliary pulse-igniter electrode), and electrical efficiency n sub E 61.9%. The thruster incorporates a wire-mesh anode and tantalum cover surfaces to control discharge chamber flake formation and employs an auxiliary pulse-igniter electrode for hollow-cathode ignition. When the SIT-8 is integrated with the compatible SIT-5 propellant tankage, the system envelope is 35 cm long by 13 cm flange bolt circle with a mass of 9.8 kg including 6.8 kg of mercury propellant. Two thrust vectoring systems which generate beam deflections in two orthogonal directions were also developed under the program and tested with the 8-cm thruster. One system vectors the beam over + or - 10 degrees by gimbaling of the entire thruster (not including tankage), while the other system vectors the beam over + or - 7 degrees by translating the accel electrode relative to the screen electrode.

Hyman, J., Jr.↗

Parametric Mixer Design Methods for Acoustic Explorations of Internally Mixed Nozzles

As part of a jet noise reduction effort under the Commercial Supersonics Technology project, a process for the design and optimization of a parametric dual-stream lobe mixer is presented. The lobe mixer is fully defined from five independent parameters, including mixer length, lobe count, lobe penetration, lobe bias, and trailing edge slope. The mixers are inserted into an academic, constant-area nozzle with an external plug. The full geometry is meshed in Pointwise v18.4, and the RANS equations are solved in FUN3D v13.6. The results of a DoE study, involving ~80 different mixer designs, are presented, along with an associated parametric model for various performance response quantities. With a design space fully explored, observations and recommendations are given, and future work is discussed to further improve the design.

Dual-stream Mixer↗

Parallel Adaptive High-Order CFD Simulations Characterizing Cavity Acoustics for the Complete SOFIA Aircraft

This paper presents one-of-a-kind MPI-parallel computational fluid dynamics simulations for the Stratospheric Observatory for Infrared Astronomy (SOFIA). SOFIA is an airborne, 2.5-meter infrared telescope mounted in an open cavity in the aft of a Boeing 747SP. These simulations focus on how the unsteady flow field inside and over the cavity interferes with the optical path and mounting of the telescope. A temporally fourth-order Runge-Kutta, and spatially fifth-order WENO-5Z scheme was used to perform implicit large eddy simulations. An immersed boundary method provides automated gridding for complex geometries and natural coupling to a block-structured Cartesian adaptive mesh refinement framework. Strong scaling studies using NASA's Pleiades supercomputer with up to 32,000 cores and 4 billion cells shows excellent scaling. Dynamic load balancing based on execution time on individual AMR blocks addresses irregularities caused by the highly complex geometry. Limits to scaling beyond 32K cores are identified, and targeted code optimizations are discussed.

Acoustics↗