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 271 records · Page 15

Metal additive manufacturing simulation across length, time, and computing scales

Metal additive manufacturing (AM) offers a unique opportunity for production of advanced materials and complex geometries. However, variability in microstructure and properties challenges conventional approaches to design, process optimization, qualification, and materials selection. Modeling and simulation can improve understanding of AM processing and materials, but also poses major challenges for existing computational methods. Simultaneously, modern scientific computing hardware has become increasingly complex, most notably with the adoption of hybrid architectures such as Graphical Processing Units (GPUs). If appropriately utilized, emerging computational capabilities provide an opportunity to reveal new insight into AM processing and the resulting material structure and properties. In this review we describe the computational AM landscape, identify critical gaps, and highlight opportunities to impact the development and application of AM. First, the requirements and challenges of representative AM problem statements will be defined. Here, these problems range from scientific studies to industrial applications and are designed to capture the breadth of challenges facing the AM community. Next, the current state of AM modeling and simulation is evaluated, broken down by enabling hardware and software, process simulation, microstructure simulation, and property simulation. Each section describes the diversity of simulation approaches and associated trade-offs in physical fidelity and computational expense. Each area is then assessed based on their suitability and readiness for current and developing computational architectures. Lastly, the greatest opportunities for future research and application are highlighted, including gaps in modeling capabilities, opportunities for near-term application, and key scientific challenges.

additive manufacturing↗

Electrocatalytic Reduction of CO 2 to CO over Ag(110) and Cu(211) Modeled by Grand-Canonical Density Functional Theory

We report the results of modeling CO 2 reduction (CO 2 R) to CO over Ag(110) and Cu(211) surfaces at different applied potentials using grand-canonical density functional theory (GC-DFT), a method specifically designed to accurately model electrochemical systems. In addition to demonstrating GC-DFT’s ability to accurately model electrochemical processes, we also compare it with the computational hydrogen electrode (CHE) approach. GC-DFT predicts that the geometries of these reacting systems strongly depend on the applied potential and the Helmholtz free energies vary nonlinearly with the applied potential, which contradicts a central assumption of the CHE approach. The CHE approach neglects the change in the number of electrons on the electrode surface at different applied potentials, which reduces its accuracy as the potential changes from the potential of zero charge. Our results further demonstrate that the grand free energies of the reaction intermediates not only depend on the value of the applied potential but also on the metal surface type, adsorption site, and adsorbate. GC-DFT’s ability to predict the effect of the applied potential on adsorbate geometry enables it to evaluate different possible reaction mechanisms at different applied potentials. For instance, GC-DFT predicts that the first step of CO2R likely switches from proton-coupled electron transfer to sequential electron transfer and then proton transfer at more reducing potentials, a result that cannot be determined by the CHE because it assumes that all electron transfers are coupled to proton transfers and neglects the effect of the applied potential on the adsorbate geometry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Giant Wilson loops and AdS 2 /dCFT 1

The 1/2-BPS Wilson loop in N = 4 supersymmetric Yang-Mills theory is an important and well-studied example of conformal defect. In particular, much work has been done for the correlation functions of operator insertions on the Wilson loop in the fundamental representation. In this paper, we extend such analyses to Wilson loops in the large-rank symmetric and antisymmetric representations, which correspond to probe D3 and D5 branes with AdS 2 × S 2 and AdS 2 × S 4 worldvolume geometries, ending at the AdS 5 boundary along a one-dimensional contour. We first compute the correlation functions of protected scalar insertions from supersymmetric localization, and obtain a representation in terms of multiple integrals that are similar to the eigenvalue integrals of the random matrix, but with important differences. Using ideas from the Fermi Gas formalism and the Clustering method, we evaluate their large N limit exactly as a function of the ’t Hooft coupling. The results are given by simple integrals of polynomials that resemble the Q-functions of the Quantum Spectral Curve, with integration measures depending on the number of insertions. Next, we study the correlation functions of fluctuations on the probe D3 and D5 branes in AdS. We compute a selection of three- and four-point functions from perturbation theory on the D-branes, and show that they agree with the results of localization when restricted to supersymmetric kinematics. We also explain how the difference of the internal geometries of the D3 and D5 branes manifests itself in the localization computation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Serpent transient analyses of GIACINT geometrical change experiments

This paper simulates seven transient experiments with geometrical movements performed at the GIACINT facility of Belarus. These experiments include three slow control rods insertion in more than 25 s, one fast control rod insertion in less than 1 s, and three water moderator draining transients. In the fast control rod insertion experiment, the control rods are first pushed by a mechanical spring and then let fall by gravity in the assembly. In this experiment, the control rods speed varies because of the initial spring force, the gravitational acceleration, and the water buoyancy force. In all other transient experiments, the geometry change occurs with a constant speed. The Serpent computer program was utilized to simulate these transient experiments with geometrical changes by performing two separate computations. One time-independent computation is performed in criticality mode and the other time-dependent computation is performed in dynamic-source mode. The first computation writes the neutron population and delayed neutron precursors for starting the transient simulation. The second computation divides the transient time into several time bins, reads the two files from the first computation, and updates them at the end of each time bin. The Serpent simulation results obtained for this set of experiments are in good agreement with the time dependent experimental measurements. In addition, Serpent and MCNP simulations results show an excellent agreement for the water moderator transient

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Computing material volume fractions on a superimposed mesh as applied to Monte Carlo particle transport simulations

Here, we present a newly implemented ray tracing algorithm in OpenMC for efficiently computing material volume fractions on superimposed meshes in complex geometries. By firing rays along each coordinate direction through the geometry, the approach accumulates track-length data in each mesh element, thereby determining the fractional composition of each material. Scaling studies on three different models—a random tetrahedra configuration, the Frascati Neutron Generator ITER dose rate benchmark, and a stellarator design—show excellent parallel performance, with nearly linear speedup on modern multi-threaded and distributed-memory systems. An analysis of the residual error relative to high-resolution reference solutions demonstrated that under optimal conditions it decreases as 1/R, where R is the number of rays fired, making it straightforward to achieve user-prescribed accuracy. This new functionality enables practical, mesh-based approaches for detailed nuclear analyses in production Monte Carlo workflows without resorting to expensive, fully conformal or unstructured meshing.

Monte Carlo↗

Multi-dimensional modeling of mixture preparation in a direct injection engine fueled with gaseous hydrogen

With the recent advances of direct injection (DI) technology, introducing hydrogen into the combustion chamber through DI is being considered as a viable approach to circumvent backfire and pre-ignition encountered in early generations of hydrogen engines. As part of a broader vision to develop a robust numerical model to study hydrogen spark ignition (SI) combustion in internal combustion (IC) engines, the present numerical investigation focuses on mixture preparation in a hydrogen DI SI engine. This study is carried out with a single hole injector with gaseous hydrogen injected at 100 bar injection pressure. Simulations are carried out for high and low tumble configurations and validated against optical data acquired from planar laser induced fluorescence (PLIF) measurements. Varying mesh configurations are investigated for the impact on in-cylinder mixture distribution. A particular emphasis is placed on the effect of nozzle geometry and mesh orientation near the wall. Overall, the computational model is found to predict the mixture distribution in the combustion cylinder reasonably well. The results showed that the alignment of mesh with the flow direction is important to achieve good agreement between numerical analysis and optical measurement data.

08 HYDROGEN↗

Rapid and precise tracking of water influx and efflux across cell membranes induced by a pulsed electric field

Quantitative measurements of water content within a single cell are notoriously difficult. In this work, we introduce a single-shot optical method for tracking the intracellular water content, by mass and volume, of a single cell at video rate. We utilize quantitative phase imaging and a priori knowledge of a spherical cellular geometry, leveraging a two-component mixture model to compute the intracellular water content. We apply this technique to study CHO-K1 cells responding to a pulsed electric field, which induces membrane permeabilization and rapid water influx or efflux depending upon the osmotic environment. The effects of mercury and gadolinium on water uptake in Jurkat cells following electropermeabilization are also examined.

59 BASIC BIOLOGICAL SCIENCES↗

Performance Evaluation of LPBF Manufactured 316H Components

This work represents the continuation of a benchmark study that includes modeling, fabrication and characterization as demonstration to support industry’s adoption of advanced manufacturing processes in a variety of structures. This comprehensive study investigated the feasibility of using additive manufacturing (AM) technologies, specifically Laser Powder Direct Energy Deposition (LP-DED) and Laser Powder Bed Fusion (LPBF), to produce complex nuclear microreactor components using 316H stainless steel. The research focused on manufacturing an expanded elbow pipe component with transitioning sections, which are traditionally difficult and costly to produce through conventional manufacturing methods. The overall study’s primary objectives are therefore demonstrating AM viability for nuclear applications, optimizing process parameters, developing comprehensive material characterization protocols, validating computational modeling approaches, and establishing manufacturing guidelines for complex geometries. Although the initial work included the phased approach of cubical, upscaled cylindrical components, it is to enable to obtain more knowledge for the printing of the expanded elbow structure. The project achieved significant progress in process development by successfully optimizing LP-DED parameters to achieve 99.16-99.97% relative density in 316H stainless steel components. Through systematic evaluation of sixteen cube samples with varied laser powers (400-700W) and scan speeds (600-900 mm/min), optimal processing windows were identified at 500-550W with 600-700 mm/min or 650-700W with 650-900 mm/min scan speeds. The DED manufactured 316H demonstrated mechanical properties comparable or superior to wrought materials, with Young's modulus ranging from 153-208 GPa and controlled microstructural characteristics including greater than 95% face-centered cubic (FCC) phases and engineered cellular structures with sizes between 3.23-6.17 µm.

36 MATERIALS SCIENCE↗

Reconstruction of Charged Particle Tracks in Realistic Detector Geometry Using a Vectorized and Parallelized Kalman Filter Algorithm

One of the most computationally challenging problems expected for the High-Luminosity Large Hadron Collider (HL-LHC) is finding and fitting particle tracks during event reconstruction. Algorithms used at the LHC today rely on Kalman filtering, which builds physical trajectories incrementally while incorporating material e ects and error estimation. Recognizing the need for faster computational throughput, we have adapted Kalman-filterbased methods for highly parallel, many-core SIMD and SIMT architectures that are now prevalent in high-performance hardware. Previously we observed significant parallel speedups, with physics performance comparable to CMS standard tracking, on Intel Xeon, Intel Xeon Phi, and (to a limited extent) NVIDIA GPUs. While early tests were based on artificial events occurring inside an idealized barrel detector, we showed subsequently that our mkFit software builds tracks successfully from complex simulated events (including detector pileup) occurring inside a geometrically accurate representation of the CMS-2017 tracker. Here, we report on advances in both the computational and physics performance of mkFit, as well as progress toward integration with CMS production software. Recently we have improved the overall eciency of the algorithm by preserving short track candidates at a relatively early stage rather than attempting to extend them over many layers. Moreover, mkFit formerly produced an excess of duplicate tracks; these are now explicitly removed in an additional processing step. We demonstrate that with these enhancements, mkFit becomes a suitable choice for the first iteration of CMS tracking, and eventually for later iterations as well. We plan to test this capability in the CMS High Level Trigger during Run 3 of the LHC, with an ultimate goal of using it in both the CMS HLT and oine reconstruction for the HL-LHC CMS tracker.

Cerati, Giuseppe↗

Experimental measurement of energy release from an initiating layer in an insensitive explosive

When subjected to a shock of insufficient strength to trigger prompt reaction, heterogenous condensed phase explosives can form regions where significant amounts of the explosive remain unreacted for times much greater than the reaction time of the detonating explosive. This phenomena is observed for the explosive PBX 9502 (95 wt% TATB) both for planar and oblique input shocks. Here, we build on previous results by performing cylinder expansion (CYLEX) tests where the explosive charge is comprised of a faster core of PBX 9501 (95 wt% HMX) inside a slower annulus of PBX 9502. The detonation in the faster PBX 9501 drives an oblique shock into the adjacent PBX 9502, and an annular transverse initiating layer (IL) results. In the test geometry, the IL travels steadily down the length of the test after a short run distance. At radial positions beyond the IL, an annular region of detonating PBX 9502 is observed. Using standard CYLEX test diagnostics, we infer the total energy release of this experiment. By making the assumptions that (1) the combined energy release is comprised of contributions from detonating PBX 9501, detonating PBX 9502, and the IL in the PBX 9502 and (2) mass-specific energy release for the detonating explosives is approximately the same as typically observed for each explosive, the IL energy release and reaction efficiency can be computed. Results are compared to prior results for a similar geometry, and indicate that while shock deadened PBX 9502 does not detonate promptly, it does eventually release a significant portion of its chemical potential energy over longer timescales on the order of 10 µs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Analytical Study and Experimental Verification of Shear-Thinning Ink Flow in Direct Ink Writing Process

Abstract Direct ink writing (DIW) process is a facile additive manufacturing technology to fabricate three-dimensional (3D) objects with various materials. Its versatility has attracted considerable interest in academia and industry in recent years. As such, upsurging endeavors are invested in advancing the ink flow behaviors in order to optimize the process resolution and the printing quality. However, so far, the physical phenomena during the DIW process are not revealed in detail, leaving a research gap between the physical experiments and its underlying theories. Here, we present a comprehensive analytical study of non-Newtonian ink flow behavior during the DIW process. Different syringe-nozzle geometries are modeled for the comparative case studies. By using the computational fluid dynamics (CFD) simulation method, we reveal the shear-thinning property during the ink extrusion process. Besides, we study the viscosity, shear stress, and velocity fields, and analyze the advantages and drawbacks of each syringe-nozzle model. On the basis of these investigations and analyses, we propose an improved syringe-nozzle geometry for stable extrusion and high printing quality. A set of DIW printing experiments and rheological characterizations are carried out to verify the simulation studies. The results developed in this work offer an in-depth understanding of the ink flow behavior in the DIW process, providing valuable guidelines for optimizing the physical DIW configuration toward high-resolution printing and, consequently, improving the performance of DIW-printed objects.

Engineering↗

Integrated Computational Materials and Mechanical Modeling for Additive Manufacturing of Alloys with Graded Structure Used in Fossil Fuel Power Plants

Wire-arc additive manufacturing (WAAM) has demonstrated its unique capability of producing large-size alloy components with a significantly reduced fabrication time and enhanced geometry design freedom. In this project, the team has developed an ICME (Integrated Computational Materials Engineering) modeling framework, which supports the WAAM of the AUSC (Advanced Ultra-Supercritical) power plant components. The manufacturing design has been applied to Inconel 740H, steel P91, as well as the dissimilar alloy components between steel P91 and Inconel 740H. The ICME model framework is developed by considering two types of modeling. First, mechanistic modeling has been applied to control the printing quality and understand the sequence of the dissimilar printing of the wall structure. The following models have been included in the developed ICME framework: finite element thermal model, grain structure model, residual stress simulation, crystal plasticity model, CALPHAD-based precipitation kinetic model, phase stability prediction, thermal expansion predictive model, and heuristic creep model. Secondary, a physics-based machine learning model has also been developed based on the ICME model structure. The machine learning model development is based on the ICME model prediction with calibration of the experiments. In addition, the WAAM has been utilized as a high-throughput experimental tool rapidly generating a gradient of alloy composition to facilitate experimental database generation for process-structure-property relationships. Such a database directly supported the ICME-enhanced machine learning, which further assisted in intermediate composition block design between P91 and 740H. A high-throughput screening study of the oxidation resistance has been performed based on such high-throughput experimentation. Based on the computational design, several dissimilar alloy manufacturing with post-heat treatment have been performed with a comprehensive evaluation of mechanical performance, including hardness mapping, yield strength, creep resistance. In this project, the single component of P91 and 740H processed by WAAM after heat treatment designed by ICME has demonstrated higher performance in yield strength and creep resistance than the wrought materials. The P91 sample prepared by WAAM with ICME-designed heat treatment performs better than P92 in creep resistance. The designed graded alloy printing with intermediate block shows a promising performance that exceeds the traditional welding. Moreover, the current research indicates the high need for location-specific design analysis with uncertainty quantification, an important topic that deserves more dedicated research. The achievement of this project demonstrated the promising future of WAAM in structural alloy manufacturing for energy power plant development. Successful printing requires synergetic efforts made by manufacturing, mechanical, and materials sciences.

20 FOSSIL-FUELED POWER PLANTS↗

Influence of intermolecular interactions on the infrared complex indices of refraction for binary liquid mixtures

This paper investigates the accuracy of deriving the composite optical constants of binary mixtures from only the complex indices of refraction of the neat materials. These optical constants enable the reflectance spectra of the binary mixtures to be modeled for multiple scenarios (e.g., different substrates, thicknesses, volume ratios), which is important for contact and standoff chemical detection. Using volume fractions, each mixture's complex index of refraction was approximated via three different mixing rules. To explore the impact of intermolecular interactions, these predictions are tested by experimental measurements for two representative sets of binary mixtures: (1) tributyl phosphate combined with n-dodecane, a non-polar medium, to represent mixtures which primarily interact via dispersion forces and (2) tributyl phosphate and 1-butanol to represent mixtures with polar functional groups that can also interact via dipole–dipole interactions, including hydrogen bonding. The residuals and the root-mean-square error between the experimental and calculated index values are computed and demonstrate that for miscible liquids in which the average geometry of the cross-interactions can be considered isotropic (e.g., dispersion), the refractive indices of the mixtures can be modeled using composite n and k values derived from volume fractions of the neat liquids. Conversely, in spectral regions where the geometry of the cross-interactions is more restricted and anisotropic (e.g., hydrogen bonding), the calculated n and k values vary from the measured values. The impact of these interactions on the reflectance spectra are then compared by modeling a thin film of the binary mixtures on an aluminum substrate using both the measured and the mathematically computed indices of refraction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Coordinate Transformations

This whitepaper provides details for the background, conventions, definitions, transformation equations, examples, and implementation used in the Solid Body Geometry (SBG) transformations in Monte Carlo Application ToolKit (MCATK).

97 MATHEMATICS AND COMPUTING↗

The Page curve of Hawking radiation from semiclassical geometry

We consider a gravity theory coupled to matter, where the matter has a higher-dimensional holographic dual. In such a theory, finding quantum extremal surfaces becomes equivalent to finding the RT/HRT surfaces in the higher-dimensional theory. Using this we compute the entropy of Hawking radiation and argue that it follows the Page curve, as suggested by recent computations of the entropy and entanglement wedges for old black holes. The higher-dimensional geometry connects the radiation to the black hole interior in the spirit of ER=EPR. The black hole interior then becomes part of the entanglement wedge of the radiation. Inspired by this, we propose a new rule for computing the entropy of quantum systems entangled with gravitational systems which involves searching for “islands” in determining the entanglement wedge.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Particle-in-cell simulation of laser wakefield accelerators with oblique lasers in quasicylindrical geometry

In the studies of optical injections in laser or beam-driven wakefield accelerators, there is a frequent demand for using obliquely propagating assistant lasers in particle-in-cell simulations. In conventional methods, this is only possible in either two- or three-dimensional Cartesian geometries, which have the drawbacks of either lack of fidelity for the actual situation or requiring a huge amount of computational resources. In this work, we develop a new method that uses an expression-defined oblique laser field to simulate such situations in quasicylindrical geometry particle-in-cell simulations, having the advantages of both maintaining good fidelity and saving computational resources. As an example, we use this method in the scissor-cross ionization injection scheme for the optimization of the injected beam quality. This method is widely applicable to particle injections with assistant lasers in wakefield accelerators, as long as the assistant lasers only influence the injected particles during the injection process. Published by the American Physical Society 2025

Ma, Minghao (ORCID:0009000602317775)↗

Spectral Element Applications in Complex Nuclear Reactor Geometries: Tet-to-Hex Meshing

The spectral element code Nek5000 is an open-source, higher-order computational fluid dynamics code developed at Argonne National Laboratory. It is designed to solve incompressible Navier-Stokes equations, but it also has a low-Mach-number approximation feature available. Large eddy simulation is approached by explicit filtering of the velocity field (and other fields) to mimic the effect of dissipation due to the unresolved scale. The computational domain is decomposed into second-order hexahedral elements that conform to the boundaries. However, generating a high-quality pure-hexahedral mesh can be challenging for some problems. For simple geometries, traditional blocking methods can be used to decompose the domain into smaller blocks to generate a so-called structural mesh. A structural mesh can maintain good orthogonality but can have a highly skewed mesh to conform to the geometry, as well as unnecessary refinement in the far field. Moreover, for geometries with relative complexity, blocking the geometry becomes impossible. To address these issues, we adopted a let-to-hex strategy to generate a pure hexahedral mesh for Nek5000. First, we generate a pure tetrahedral mesh for the geometry; then we divide one tetrahedral element into four hexahedral elements. A pure tetrahedral mesh could be easily generated for complex geometries by using many current meshing codes. In this paper, we use the commercial codes ANSYS meshing and ANSYSICEM to generate the pure tetrahedral mesh and then convert it to a pure hexahedral mesh. Boundary layers are extruded in ANSYSICEM to maintain near-wall resolution.

ANSYS ICEM↗

A Bezier Curve Informed Melt Pool Geometry to Model Additive Manufacturing Microstructures Using SPPARKS

Additive manufacturing is a transformative technology with the potential to manufacture designs which traditional subtractive machining methods cannot. Additive manufacturing offers fast builds at near final desired geometry; however, material properties and variability from part to part remain a challenge for certification and qualification of metallic components. AM induced metallic microstructures are spatially heterogeneous and highly process dependent. Engineering properties such as strength and toughness are significantly affected by microstructure morphologies resulting from the manufacturing process Linking process parameters to microstructures and ultimately to the dynamic response of AM materials is critical to certifying and qualifying AM built parts and components and improving the performance of AM materials. The AM fabrication process is characterized by building parts layer by layer using a selective laser melt process guided by a computer. A laser selectively scans and melts metal according to a designated geometry. As the laser scans, metal melts, fuses, and solidifies forming the final geometry in a layerwise fashion. As the laser heat source moves away, the metal cools and solidifies forming metallic microstructures. This work describes a microstructure modeling application implemented in the SPPARKS kinetic Monte Carlo computational framework for simulating the resulting microstructures. The application uses Bzier curves and surfaces to model the melt pool surface and spatial temperature profile induced by moving the laser heat source; it simulates the melting and fusing of metal at the laser hot spot and microstructure formation and evolution when the laser moves away. The geometry of the melt pool is quite flexible and we explore effects of variances in model parameters on simulated microstructures.

36 MATERIALS SCIENCE↗