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.

SourceResearch connectionFrom yearThrough year

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

At least 703 records · Page 39

All sources NASA NTRS DOE OSTI

NASA NTRS1974

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.↗

NASA NTRSDate not supplied

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↗

DOE OSTI2016

Comparison between Flownex and Corona thermal-hydraulic models for a single HTR prismatic fuel assembly - 18544

Computational tools and methodologies that can perform thermal-hydraulic analyses to reasonably predict the temperature distribution for the whole core of HTRs are useful for operational and safety analyses. In this paper, Flownex (a system code based on a one-dimensional thermal-fluid network methodology) and CORONA (Core Reliable Optimization and thermo-fluid Network Analysis) thermal-hydraulic models and results for a vertically stacked single fuel assembly will be compared. In Flownex, a basic unit cell containing part of the graphite and the associated fuel elements and coolant channels was developed. CORONA uses a basic unit cell concept for generating computational grids using several types of basic unit cells to model the heat transfer through the hexagonal blocks and a 1-D formulation to model the flow through the coolant channels. In CORONA, computational grids are generated using a combination of pre-generated mesh cells. The coolant channels and fuel rods are modelled explicitly. Flownex in contrast, models fuel and coolant channels implicitly by defining representative coolant channels, fuel rods and coated particle layers associated with each unit cell. Flownex can therefore model temperatures at coated particle level and also through each graphite block. The computational mesh employed by Flownex is much coarser than the mesh employed by CORONA. In both models, the helium gap between the fuel compact and the graphite block is modelled as stagnant. Flownex results were found to be of the same order of magnitude as those obtained using CORONA. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

NASA NTRS2014

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↗

DOE OSTI2024

h5bench: A unified benchmark suite for evaluating HDF5 I/O performance on pre‐exascale platforms

Summary Parallel I/O is a critical technique for moving data between compute and storage subsystems of supercomputers. With massive amounts of data produced or consumed by compute nodes, high‐performant parallel I/O is essential. I/O benchmarks play an important role in this process; however, there is a scarcity of I/O benchmarks representative of current workloads on HPC systems. Toward creating representative I/O kernels from real‐world applications, we have created h5bench , a set of I/O kernels that exercise hierarchical data format version 5 (HDF5) I/O on parallel file systems in numerous dimensions. Our focus on HDF5 is due to the parallel I/O library's heavy usage in various scientific applications running on supercomputing systems. The various tests benchmarked in the h5bench suite include I/O operations (read and write), data locality (arrays of basic data types and arrays of structures), array dimensionality (one‐dimensional arrays, two‐dimensional meshes, three‐dimensional cubes), I/O modes (synchronous and asynchronous). In this paper, we present the observed performance of h5bench executed along several of these dimensions on existing supercomputers (Cori and Summit) and pre‐exascale platforms (Perlmutter, Theta, and Polaris). h5bench measurements can be used to identify performance bottlenecks and their root causes and evaluate I/O optimizations. As the I/O patterns of h5bench are diverse and capture the I/O behaviors of various HPC applications, this study will be helpful to the broader supercomputing and I/O community.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Iteration-based Linearized Distribution-level Locational Marginal Price for Three-phase Unbalanced Distribution Systems

Distributed energy resources (DERs) are rocking the utilities’ business landscape. It calls for competitive market environments that incentivize DERs to form maximum operating efficiency. Among proposed pricing schemes, distribution-level locational marginal price (DLMP) is effective in signaling the marginal generation cost differences driven by energy losses and network constraints. It can be derived from a distribution-level optimal power flow (OPF) framework, as it essentially presents the sensitivity of optimized generation cost towards incremental loads. However, due to the high resistance-to-inductance ratio and unbalanced characteristics of distribution networks, computational affordable DLMPs are highly challenged. This article provides a linear-approximated DLMP that can be solved efficiently and generalized to account for reactive power flow, three-phase unbalanced loads and meshed network structure. The successive linear programming technique is introduced to enhance the model accuracy. Case studies on an IEEE 123-Bus system validate its accuracy against a nonlinear benchmark and capability in offering proper incentives.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

Choosing Transport Events for Initiating Splitting and Rouletting

A study was performed to determine which transport events should be used to initiate a weight window lookup to achieve the best variance reduction performance. A weight window lookup potentially triggers particle splitting (in important regions of phase space) or rouletting (in unimportant regions), thereby optimizing computational effort. Potential initiating transport events include collisions (both pre- and post-collision), geometry surface crossings, traversing a mean-free path, and streaming across a weight window boundary. Permutations of these initiating events were tested on an urban model with background radiation sources and a spent fuel cask with a neutron dose mesh tally. Generally, all methods perform better with finer weight window meshes. Tracking on weight windows performs well for coarse weight window meshes, while a combination of splitting each mean-free path, geometric surface crossing, and before collisions performs well for fine weight window meshes.

42 ENGINEERING↗

NASA NTRS2021

Applying the Cognitive Space Gateway to Swarm Topologies

NASA's future vision for interplanetary networking includes a lunar network, Cube Satellite (CubeSat) constellations, and deep space robotic missions, comprising what could be viewed as a network of networks. Delay-tolerant networking (DTN) architecture and protocols provide a standard network layer among these varying scenarios and mitigate many challenges of the space environment, such as long delays, unplanned service interruptions, and asymmetric links. The Cognitive Space Gateway (CSG) is a routing method in a DTN architecture that uses spiking neural networks as the learning element to optimize routing decisions in a complex environment. This work aims to further develop cognitive networking technologies in several critical areas, including DTN, the CSG algorithm, CubeSat swarm topologies, and cloud services. To test the algorithm in a realistic scenario, the emulated network topology is based on a CubeSat swarm. The swarm may function as a mesh of nodes or as a hub-and-spoke network. An emulation environment will be built upon a commercial cloud service, such as Amazon Web Services (AWS) Elastic Compute Cloud. The cloud environment may enable a flexible, lower maintenance approach versus a multi-hop network based in a physical laboratory. The cloud platform will provide a secure environment allowing for collaboration among government and academic entities.

Ricardo Lent↗

DOE OSTI2021

Method and apparatus for shape-based energy analysis of solids

A computer-readable medium stores instructions including storing a model of a physical structure and defining a mesh for the model. The mesh includes finite elements, each defined by a respective set of edges. The instructions include, for each finite element, identifying a governing differential equation and a set of complementary functions that exactly satisfy the governing differential equation. The instructions include determining an applied physical stimulus for the physical structure. The instructions include generating an energy optimization model that minimizes a difference between internal energy and external energy of the finite elements in response to the applied physical stimulus. The model includes a matrix of respective scalar multipliers for the complementary functions for each finite element. The instructions include transforming the matrix and calculating a physical parameter of interest. The instructions include, in response to the physical parameter not satisfying a design parameter, updating the model of the physical structure.

Spears, Robert E.↗

NASA NTRS1995

Multigrid solution strategies for adaptive meshing problems

This paper discusses the issues which arise when combining multigrid strategies with adaptive meshing techniques for solving steady-state problems on unstructured meshes. A basic strategy is described, and demonstrated by solving several inviscid and viscous flow cases. Potential inefficiencies in this basic strategy are exposed, and various alternate approaches are discussed, some of which are demonstrated with an example. Although each particular approach exhibits certain advantages, all methods have particular drawbacks, and the formulation of a completely optimal strategy is considered to be an open problem.

Mavriplis, Dimitri J.↗

NASA NTRS1995

Multigrid solution strategies for adaptive meshing problems

Mavriplis, Dimitri J.↗

DOE OSTI2010

EMGeo 3.0

EMGeo 3.0 an algorithm that improves both the computational capabilities of joint 3D electromagnetic, EM and magnetotelluric and acoustic seismic wavefield simulation and inverse modeling. This seismic wavefield module requires that data be transformed first to the Laplace-Fourier Domain before processing. The inverse modeling is based upon non-liner conjugated gradients and 3D finite differences approximations of the seismic, EM and MT wavefields. Improved modeling efficiencies are obtain by separating the modeling grid from the simulation grid. This grid separation method allows for treatment of very large data sets and imaging volumes. Further computational efficiency is obtained by combining different levels of parallezation using the Message Passing Interface (MPI). Bound constraints are employed in the imaging process to insure stability. Additional acceleration of the inverse modeling is achieved by preconditioning the conjugate gradient optimizer using an approximate Hessian. The algorithm includes improved capabilities to accurately treat models that exhibit transverse anisotropy in the electrical conductivity. However, the seismic acoustic velocity is assumed to be isotropic. Models exhibiting topography and bathymetry can be treated. The software also includes a set of pre and post processing tools for designing input model meshes and data plotting.

Gregory Newman, Michael Commer↗

DOE OSTI2025

Analysis of the SiMPL Method for Density-Based Topology Optimization

We present a rigorous convergence analysis of a new method for density-based topology optimization that provides pointwise bound-preserving design updates and faster convergence than other popular first-order topology optimization methods. Due to its strong bound preservation, the method is exceptionally robust, as demonstrated in numerous examples here and in the companion article [D. Kim et al., Struct. Multidiscip. Optim., 68 (2025), 74]. Furthermore, it is easy to implement with clear structure and analytical expressions for the updates. Our analysis covers two versions of the method, characterized by the employed line search strategies. We consider a modified Armijo backtracking line search and a Bregman backtracking line search. For both line search algorithms, our algorithm delivers a strict monotone decrease in the objective function and further intuitive convergence properties, e.g., strong and pointwise convergence of the density variables on the active sets, norm convergence to zero of the increments, convergence of the Lagrange multipliers, and more. In addition, the numerical experiments demonstrate apparent mesh-independent convergence of the algorithm. Here, we refer to the new algorithm as the SiMPL method (pronounced “simple”), which stands for Sigmoidal Mirror descent with a Projected Latent variable.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Multi-plane moment-of-fluid interface reconstruction in 3D

Moment-of-fluid (MOF) methods for interface reconstruction approximate the region occupied by material in each mesh element only through reference to its geometric moments. Here, we present a 3D MOF method that represents the material (POM) in each cell as the convex intersection of the cell and multiple half-spaces, each selected to minimize the least-squares error between computed moments of the approximated material and provided reference moments. This optimization problem is highly non-linear and non-convex, making the numerical result very sensitive to the initial guess. To create an effective initial guess in each cell, we construct an ellipsoid from 0th–2nd order reference moments such that its shape corresponds with that of the POM. Within this ellipsoid we inscribe a polyhedron, and initialize the minimization problem with the half-spaces defined by each of its faces. The inscribed polyhedron has minimally 4 faces, and using up to 3rd order moments permits optimization over up to 20 unknown values. We therefore define MOF methods that utilize 4, 5, or 6 half-spaces, correspondingly initialized with the faces of a single inscribed tetrahedron, triangular prism, or hexahedron. Stability of the non-linear optimization is further improved with a prepossessing step that normalizes the reference moments according to the axes of the reference ellipsoid. Using this approach, the non-linear least-squares solver reliably converges to a near-global minimum from a single initial guess. We demonstrate accuracy and robustness using single-cell and multi-cell examples over a wide spectrum of geometry. In particular, we demonstrate our ability to exactly reproduce several important and complex features defined by up to four half-spaces, such as corners, filaments, filament tips, and embedded material in the cell.

3D interface reconstruction↗

DOE OSTI2011

ePLAS Development for Jet Modeling and Applications

Plasma jets provide an alternate approach to the creation of high energy density laboratory plasmas (HEDLP). For the Plasma Liner Experiment (PLX), typically 30 partially ionized argon jets, produced with mini-rail guns, will be focused into a central volume for subsequent magnetic compression into high density plasma liners that can reach high (0.1 Mbar) peak pressures upon stagnation. The jets are typically 2.5 cm in radius traveling at Mach number 30. Ultimate success will require optimized tuning of the rail configurations, the nozzles injecting the gases, and the careful implementation of pre-ionization. The modeling of plasma jet transport is particularly challenging, due the large space (100 sq cm) and time scales (microseconds) involved. Even traditional implicit methods are insufficient, due to the usual need to track electrons explicitly on the mesh. Wall emission and chemistry must be managed, as must ionization of the jet plasma. Ions in the jets are best followed as particles to account properly for collisions upon jet merger. This Phase I Project developed the code ePLAS to attack and successfully surmount many of these challenges. It invented a new 'super implicit' electromagnetic scheme, using implicit electron moment currents that allowed for modeling of jets over multi-cm and multi-picoseconds on standard, single processor 2 GHz PCs. It enabled merger studies of two jets, in preparation for the multi-jet merger problem. The Project explored particle modeling for the ions, and prepared for the future addition of a grid-base jet ion collision model. Access was added to tabular equations of state for the study of ionization effects in merging jets. The improved code was discussed at the primary plasma meetings (IEEE and APS) during the Project period. Collaborations with National Laboratory and industrial partners were nurtured. Code improvements were made to facilitate code use. See: http://www.researchapplicationscorp.com. The ePLAS code enjoys EAR99 export control treatment, permitting distribution to most foreign countries without a license.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2011

Recovery Act: Federspiel Controls (now Vigilent) and State of California Department of General Services Data Center Energy Efficient Cooling Control Demonstration. Final technical project report

Eight State of California data centers were equipped with an intelligent energy management system to evaluate the effectiveness, energy savings, dollar savings and benefits that arise when powerful artificial intelligence-based technology measures, monitors and actively controls cooling operations. Control software, wireless sensors and mesh networks were used at all sites. Most sites used variable frequency drives as well. The system dynamically adjusts temperature and airflow on the fly by analyzing real-time demands, thermal behavior and historical data collected on site. Taking into account the chaotic interrelationships of hundreds to thousands of variables in a data center, the system optimizes the temperature distribution across a facility while also intelligently balancing loads, outputs, and airflow. The overall project will provide a reduction in energy consumption of more than 2.3 million kWh each year, which translates to $240,000 saved and a reduction of 1.58 million pounds of carbon emissions. Across all sites, the cooling energy consumption was reduced by 41%. The average reduction in energy savings across all the sites that use VFDs is higher at 58%. Before this case study, all eight data centers ran the cooling fans at 100% capacity all of the time. Because of the new technology, cooling fans run at the optimum fan speed maintaining stable air equilibrium while also expending the least amount of electricity. With lower fan speeds, the life of the capital investment made on cooling equipment improves, and the cooling capacity of the data center increases. This case study depicts a rare technological feat: The same process and technology worked cost effectively in eight very different environments. The results show that savings were achieved in centers with diverse specifications for the sizes, ages and types of cooling equipment. The percentage of cooling energy reduction ranged from 19% to 78% while keeping temperatures substantially within the limits recommended by the American Society of Heating, Refrigeration, and Air-Conditioning Engineers (ASHRAE) for data center facilities.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2012

Zgoubi Users' Guide

The computer code zgoubi calculates trajectories of charged particles in magnetic and electric fields. At the origin specially adapted to the definition and adjustment of beam lines and magnetic spectrometers, it has so evolved that it allows the study of systems including complex sequences of optical elements such as dipoles, quadrupoles, arbitrary multipoles, FFAG magnets and other magnetic or electric devices, and is able as well to handle periodic structures. Compared to other codes, it presents several peculiarities, as follows - a non-exhaustive list : • a numerical method for integrating the Lorentz equation, based on Taylor series, which optimizes computing time and provides high accuracy and strong symplecticity, • spin tracking, using the same numerical method as for the Lorentz equation, • account of stochastic photon emission, and its effects on particle dynamics, • calculation of the synchrotron radiation electric field and spectra in arbitrary magnetic fields, from the ray-tracing outcomes, • the possibility of using a mesh, which allows ray-tracing from simulated or measured (1-D, 2-D or 3-D) electric and magnetic field maps, • a number Monte Carlo procedures : unlimited number of trajectories, in-flight decay, stochastic radiation, etc., • built-in fitting procedures allowing arbitrary variables and a large variety of constraints, easily expandable, • multi-turn tracking in circular accelerators including features proper to machine parameter calculation and survey, • simulation of time-varying power supplies, and • simulation of arbitrary radio-frequency programs. The initial version of zgoubi was dedicated to ray-tracing in magnetic elements, beam lines, spectrometers. It was perfected for the purpose of studying, and operating, the four spectrometers SPES I, II, III, IV at the Laboratoire National SATURNE (CEA-Saclay, France), and, later, SPEG at GANIL (Caen, France). Developments regarding spin tracking and ray-tracing in combined electric and magnetic fields were implemented, in the late 1980s and early 1990s respectively. In the mid-1990s, the computation of synchrotron radiation electromagnetic impulse and spectra was introduced, for the purpose of synchrotron radiation diagnostic R&D at LEP, and further applied to the design of the SR diagnostics installations at LHC in the early 2000s. In the mean time, several new optical elements were added, such as electro-magnetic and other electrostatic lenses. Used since several years for special studies in periodic machines (e.g., SATURNE at Saclay, COSY at Julich, LEP and LHC at CERN), zgoubi has also undergone extensive developments regarding storage ring related features. Many developments have been accomplished since the early 2000s in the frame of a number of project design and beam dynamics studies, as the neutrino factory, lepton and hadron colliders, spin studies at AGS and RHIC, etc. As a consequence the list of optical elements and the compendium of numerical methods, so-called “Glossary of Keywords” list, pp. 7 and 169, has stretched with new simulation and computing procedures, ranging from fitting (FIT2 procedure ; additional constraints) to overlapping magnetic field capabilities (DIPOLES, FFAG), spin manipulation (SPINR), optical elements (BEAMBEAM), radiation damping tools (SRLOSS) and many others (FAISTORE, MAP2D-E, OPTICS, PICKUPS, TWISS, etc.). The graphic interface to zgoubi (zpop, Part D) has undergone extensive developments, making it a convenient companion tool to the use of zgoubi.

43 PARTICLE ACCELERATORS↗

NASA NTRS2012

Flow Simulation of N2B Hybrid Wing Body Configuration

The N2B hybrid wing body aircraft was conceptually designed to meet environmental and performance goals for the N+2 generation transport set by the subsonic fixed wing project. In this study, flow fields around the N2B configuration is simulated using a Reynolds-averaged Navier-Stokes flow solver using unstructured meshes. Boundary conditions at engine fan face and nozzle exhaust planes are provided by response surfaces of the NPSS thermodynamic engine cycle model. The present flow simulations reveal challenging design issues arising from boundary layer ingestion offset inlet and nacelle-airframe interference. The N2B configuration can be a good test bed for application of multidisciplinary design optimization technology.

Kim, Hyoungjin↗