Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Modeling of resistive plasma response in toroidal geometry using an asymptotic matching approach

The method of solving linear resistive plasma response, based on the asymptotic matching approach, is developed for full toroidal tokamaks by upgrading the Resistive DCON code [A.H. Glasser, Z.R. Wang and J.-K. Park, Physics of Plasmas, \textbf{23}, 112506 (2016)]. The derived matching matrix, asymptotically matching the outer and inner regions, indicates that the applied three dimension (3-D) magnetic perturbations contribute additional small solutions at each resonant surface due to the toroidal coupling of poloidal modes. In contrast, the resonant harmonic only affects the corresponding resonant surface in the cylindrical plasma. Since the solution of ideal outer region is critical to the asymptotic matching and is challenging to be solved in the toroidal geometry due to the singular power series solution at the resonant surfaces, systematic verification of the outer region $\Delta^\prime$ matrix is made by reproducing the well known analytical $\Delta^{\prime}$ result in [H.P. Furth, P.H. Rutherford and H. Selberg, The Physics of Fluids, \textbf{16}, 1054-1063 (1073)] as well as by making a quantitative benchmark with the PEST3 code [A. Pletzer and R.L. Dewar, J. Plasma Physics, \textbf{45}, 427-451 (1991)]. Finally, the reconstructed numerical solution of resistive plasma response from the toroidal matching matrix is presented. Comparing with the ideal plasma response, the global structure of the response can be affected by the small finite island at the resonant surfaces.

MHD simulation↗
DOE OSTI2025

Geometry Mesher

This program implements tetrahedral meshing for implicitly-defined geometries, focusing on meshing unions of geometric primitives like spheres and cylinders. This enables users to quickly generate quality meshes of truss lattice and filament structures.

Watts, SethE [Lawrence Livermore National Laborato↗
DOE OSTI2024

Learning Only on Boundaries: A Physics-Informed Neural Operator for Solving Parametric Partial Differential Equations in Complex Geometries

Recently, deep learning surrogates and neural operators have shown promise in solving partial differential equations (PDEs). However, they often require a large amount of training data and are limited to bounded domains. In this work, we present a novel physics-informed neural operator method to solve parameterized boundary value problems without labeled data. By reformulating the PDEs into boundary integral equations (BIEs), we can train the operator network solely on the boundary of the domain. This approach reduces the number of required sample points from $O(N^d)$ to $O(N^{d-1}$), where $d$ is the domain’s dimension, leading to a significant acceleration of the training process. Additionally, our method can handle unbounded problems, which are unattainable for existing physics-informed neural networks (PINNs) and neural operators. Finally, our numerical experiments show the effectiveness of parameterized complex geometries and unbounded problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Comparison of nested geometry treatments within GPU-based Monte Carlo neutron transport simulations of fission reactors

Monte Carlo (MC) neutron transport provides detailed estimates of radiological quantities within fission reactors. This involves tracking individual neutrons through a computational geometry. CPU-based MC codes use multiple polymorphic tracker types with different tracking algorithms to exploit the repeated configurations of reactors, but virtual function calls have high overhead on the GPU. The Shift MC code was modified to support GPU-based tracking with three strategies: dynamic polymorphism with virtual functions, static polymorphism, and a single tracker type with tree-based acceleration. On the Frontier supercomputer these methods achieve 77.8%, 91.2%, and 83.4%, respectively, of the tracking rate obtained using a specialized tracker optimized for rectilinear-grid-based reactors. This indicates that all three methods are suitable for typical reactor problems in which tracking does not dominate runtime. The flexibility of the single tracker method is highlighted with a hexagonal-grid microreactor problem, performed without hexagonal-grid-specific tracking routines, providing a 2.19× speedup over CPU execution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Decay Dose Shielding Analysis with Hybrid Unstructured Mesh/Constructive Solid Geometry Monte Carlo Calculation and ADVANTG Acceleration

The Oak Ridge National Laboratory (ORNL) Second Target Station (STS) neutron production facility is an accelerator driven pulsed neutron source that is currently being actively developed at ORNL. The neutrons are produced by proton-induced spallation reactions. A proton beam of 700 kW power is delivered to a spallation target in short, less than 1 µs long pulses, with 15 Hz repetition rate. The spallation target of ORNL STS is a rotating water-cooled tungsten target with tantalum cladding housed in a stainless-steel shroud. It is divided into 21 segments. These segments become highly activated due to spallation reactions or nuclei transmutation by the emitted neutrons. The radioactive nuclides continue to decay after ceasing operation. The decay dose rates generated from the target segments once they are removed from their operational location within the core vessel must be accurately quantified to determine the shielding configurations of remote handling tools and transport casks and to aid in planning maintenance events. To determine the shielding configurations needed for an activated target segment after ceasing operation, both the hybrid unstructured mesh (UM)/constructive solid geometry (CSG) approach that was previously utilized for STS analyses [1] and the ADVANTG code [2] were used. Even though the ADVANTG code does not include UM capability, the utilization of its advanced variance reduction technique was crucial to accelerate the extremely difficult final photon transport calculation in this analysis. This paper also describes the procedures taken to mitigate the convergence issues that occur when ADVANTG uses a source definition that does not match the source of the final Monte Carlo (MC) calculation. These convergence issues often occur because of inconsistencies between source and transport biasing parameters.

Ibrahim, Ahmad↗
DOE OSTI2022

Temporal and spatial laser intensification within nodular defects overcoated with multilayer dielectric mirrors over a wide range of defect geometries

The nodular defect shape and the laser incidence angle have a dramatic impact on the spatial distribution of light intensification within the nodule as well as how the laser light is drained from the defect. Nodular defect geometries unique to ion beam sputtering, ion-assisted deposition, and electron-beam (e-beam) deposition, respectively, are modeled in this parametric study over a wide range of nodular inclusion diameters and layer count for optical interference mirror coatings constructed with quarter-wave thicknesses and capped with a half wave of the low index material. It was found for hafnia ($n$=1.9) and silica ($n$=1.45) multilayer mirrors that the light intensification in nodular defects with a $C$ factor of 8, typical of e-beam deposited coatings deposited with a wide range of deposition angles, was maximized for a 24-layer design. For intermediate size inclusion diameters, increasing the layer count for normal incidence multilayer mirrors reduced the light intensification within the nodular defect. A second parametric study explored the impact of the nodule shape on the light intensification for a fixed number of layers. In this case, there is a strong temporal trend for the different nodule shapes. Narrow nodules tend to drain more laser energy through the bottom of the nodule into the substrate while wide nodules tend to drain more laser energy through the top of the nodule when irradiated at normal incidence. Further, at a 45° incidence angle, waveguiding is an additional method to drain laser energy from the nodular defect. Finally, laser light resonates within nodular defects longer than within the adjacent nondefective multilayer structure.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Controlling light emission by engineering atomic geometries in silicon photonics

By engineering atomic geometries composed of nearly 1000 atomic segments embedded in micro-resonators, we observe Bragg resonances induced by the atomic lattice at the telecommunication wavelength. The geometrical arrangement of erbium atoms into a lattice inside a silicon nitride (SiN) microring resonator reduces the scattering loss at a wavelength commensurate with the lattice. We confirm dependency of light emission to the atomic positions and lattice spacing and also observe Fano interference between resonant modes in the system.

42 ENGINEERING↗
DOE OSTI2023

Suppression of coherent interference to electric-field-induced second-harmonic (E-FISH) signals for the measurement of electric field in mesoscale confined geometries

We present spatially enhanced electric-field-induced second-harmonic (SEEFISH) generation with a chirped femtosecond beam for measurements of electric field in mesoscale confined geometries subject to destructive spurious second-harmonic generation (SHG). Spurious SHG is shown to interfere with the measured E-FISH signal coherently, and thus simple background subtraction is not sufficient for single-beam E-FISH approaches, especially in a confined system with a large surface-to-volume ratio. The results show that a chirped femtosecond beam is effective in preventing higher-order mixing and white light generation in windows near the beam focal point which further contaminates the SEEFISH signal. The successful measurements of electric field of a nanosecond dielectric barrier discharge in a test cell demonstrated that spurious SHG detected with a congruent traditional E-FISH approach can be eliminated using the SEEFISH approach.

Vorenkamp, Madeline↗
DOE OSTI2024

Ultrabroadband two-dimensional electronic spectroscopy in the pump–probe geometry using conventional optics

We report ultrabroadband two-dimensional electronic spectroscopy (2D ES) measurements obtained in the pump–probe geometry using conventional optics. A phase-stabilized Michelson interferometer provides the pump-pulse delay interval,τ 1 , necessary to obtain the excitation-frequency dimension. Spectral resolution of the probe beam provides the detection-frequency dimension,ω 3 . The interferometer incorporates active phase stabilization via a piezo stage and feedback from interference of a continuous-wave reference laser detected in quadrature. To demonstrate the method, we measured a well-characterized laser dye sample and obtained the known peak structure. The vibronic peaks are modulated as a function of the waiting time,τ 2 , by vibrational wave packets. The interferometer simplifies ultrabroadband 2D ES measurements and analysis.

Optics↗
DOE OSTI2022

Hyporheic hydraulic geometry: Conceptualizing relationships among hyporheic exchange, storage, and water age

Hyporheic exchange is now widely acknowledged as a key driver of ecosystem processes in many streams. Yet stream ecologists have been slow to adopt nuanced hydrologic frameworks developed and applied by engineers and hydrologists to describe the relationship between water storage, water age, and water balance in finite hydrosystems such as hyporheic zones. Here, in the context of hyporheic hydrology, we summarize a well-established mathematical framework useful for describing hyporheic hydrology, while also applying the framework heuristically to visualize the relationships between water age, rates of hyporheic exchange, and water volume within hyporheic zones. Building on this heuristic application, we discuss how improved accuracy in the conceptualization of hyporheic exchange can yield a deeper understanding of the role of the hyporheic zone in stream ecosystems. Although the equations presented here have been well-described for decades, our aim is to make the mathematical basis as accessible as possible and to encourage broader understanding among aquatic ecologists of the implications of tailed age distributions commonly observed in water discharged from and stored within hyporheic zones. Our quantitative description of “hyporheic hydraulic geometry,” associated visualizations, and discussion offer a nuanced and realistic understanding of hyporheic hydrology to aid in considering hyporheic exchange in the context of river and stream ecosystem science and management.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Incorporating Full Elastodynamic Effects and Dipping Fault Geometries in Community Code Verification Exercises for Simulations of Earthquake Sequences and Aseismic Slip (SEAS)

ABSTRACT Numerical modeling of earthquake dynamics and derived insight for seismic hazard relies on credible, reproducible model results. The sequences of earthquakes and aseismic slip (SEAS) initiative has set out to facilitate community code comparisons, and verify and advance the next generation of physics-based earthquake models that reproduce all phases of the seismic cycle. With the goal of advancing SEAS models to robustly incorporate physical and geometrical complexities, here we present code comparison results from two new benchmark problems: BP1-FD considers full elastodynamic effects, and BP3-QD considers dipping fault geometries. Seven and eight modeling groups participated in BP1-FD and BP3-QD, respectively, allowing us to explore these physical ingredients across multiple codes and better understand associated numerical considerations. With new comparison metrics, we find that numerical resolution and computational domain size are critical parameters to obtain matching results. Codes for BP1-FD implement different criteria for switching between quasi-static and dynamic solvers, which require tuning to obtain matching results. In BP3-QD, proper remote boundary conditions consistent with specified rigid body translation are required to obtain matching surface displacements. With these numerical and mathematical issues resolved, we obtain excellent quantitative agreements among codes in earthquake interevent times, event moments, and coseismic slip, with reasonable agreements made in peak slip rates and rupture arrival time. We find that including full inertial effects generates events with larger slip rates and rupture speeds compared to the quasi-dynamic counterpart. For BP3-QD, both dip angle and sense of motion (thrust versus normal faulting) alter ground motion on the hanging and foot walls, and influence event patterns, with some sequences exhibiting similar-size characteristic earthquakes, and others exhibiting different-size events. These findings underscore the importance of considering full elastodynamics and nonvertical dip angles in SEAS models, as both influence short- and long-term earthquake behavior and are relevant to seismic hazard.

Geochemistry & Geophysics↗
DOE OSTI2023

Geometry Optimization for a Quadrupole Resonator at Jefferson Lab

The quadrupole resonator (QPR) is a sample characterization tool to measure the RF properties of superconducting materials using the calorime-try method at different temperatures, magnetic fields, and frequencies. Such resonators are cur-rently operating at CERN and HZB but suffer from Lorentz force detuning and modes overlap-ping, resulting in higher uncertainties in surface resistance measurement. Using the two CERN's QPR model iterations, the geometry was opti-mized via electromagnetic and mechanical simu-lations to eliminate these issues. The new QPR version was modelled for an increasing range of magnetic fields. The magnetic field is concentrat-ed at the center of the sample to reduce the un-certainty in surface resistance measurements sig-nificantly. This paper discusses the QPR geome-try optimization for the new version of QPR, which is now progressing toward fabrication.

Bira, S.↗
DOE OSTI2025

Geometry-driven modeling of electron localization in InAs/GaAs double quantum dots

The coupled electronic states in two-dimensional (2D) and three-dimensional (3D) double quantum dot (DQD) systems are investigated using a phenomenological model applied to InAs/GaAs heterostructures. The single-band k • p effective potential approach previously proposed by our group is employed to numerically calculate the energy spectrum and spatial localization of a single electron, serving as an indicator of the coupling strength within the binary system. For identical quantum dots (QDs) in a DQD, the electronic states exhibit ideal coherence. We systematically vary the DQD geometry and the strength of the confinement potential (via an applied electric field) to examine the effects of symmetry breaking and the sensitivity of electron localization in both identical and nearly identical DQDs. Our results show that coherence in DQDs is highly sensitive to these subtle variations. This sensitivity can be harnessed to detect changes in the surrounding environment, such as fluctuations in chemical or electrical properties that affect the DQD system.

electron localization↗
DOE OSTI2020

Estimation of Respirable Aerosol Release Fractions through Stress Corrosion Crack-Like Geometries

The formation of a stress corrosion crack (SCC) in the canister wall of a dry cask storage system (DCSS) has been identified as a potential issue for the long-term storage of spent nuclear fuel. The presence of an SCC in a storage system could represent a through-wall flow path from the canister interior to the environment. Modern, vertical DCSSs are of particular interest due to the significant backfill pressurization of the canister, up to approximately 800 kPa. This pressure differential offers a relatively high driving potential for blowdown of any particulates that might be present in the canister. In this study, the carrier gas flow rates and aerosol transmission properties were evaluated for an engineered microchannel with characteristic dimensions similar to those of an SCC. The microchannel was formed by mating two gage blocks with a slot orifice measuring 28.9 μm (0.0011 in.) tall by 12.7 mm (0.500 in.) wide by 8.86 mm (0.349 in.) long (flow length). Surrogate aerosols of cerium oxide, Ce0 2 , were seeded and mixed inside a pressurized tank. The aerosol characteristics were measured immediately upstream and downstream of the simulated SCC at elevated and ambient pressures, respectively. These data sets are intended to demonstrate a new capability to characterize SCCs under well-controlled boundary conditions. Separate modeling efforts are also underway that will be validated using these data. The test apparatus and procedures developed in this study can be easily modified for the evaluation of more complex SCC-like geometries including laboratory-grown SCC samples.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Improvements to CTF for Modeling of Boiling Water Reactor Geometry and Operating Conditions

This report documents work done to improve CTF modeling of boiling water reactor (BWR) designs during FY 2020. Several improvements were made to the geometry modeling capabilities, including adding support for modeling of inlet flow orifices, axially varying large water rods, channel boxes with large corner round radii, and quarter symmetry models. Work was also done to improve on the numerical stability and run-time performance of CTF in Virtual Environment for Reactor Applications (VERA) for core-scale BWR problems. An outer iteration loop was implemented that ensures mass and energy are conserved before moving on to the next time step in the solution. The ability to use a direct solver for the pressure matrix was also added, which reduces the matrix solution error and improves code stability and run-time performance. Improvements to the pressure balance iteration loop were needed to allow CTF to be used in a coupled VERA solution of a model of Hatch that was developed. Finally, various improvements have led to a 2×reduction in the total VERA walltime for the Hatch model.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Exploring new small system geometries in heavy ion collisions [Slides]

We explore various new collision geometries in the context of the publicly available hydrodynamic model SONIC. We incorporate full A-nucleon configurations for 4 He, 12 C, and 16 O obtained from quantum Monte Carlo calculations with realistic nuclear potentials, and use a Monte Carlo Glauber calculation for the initial conditions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Monte Carlo Simulation of Background and Source Measurements with CSG and CAD Geometries

Detecting radioactive sources in an urban environment is difficult due to the large magnitude and variability of the background radiation. To support the search mission of the National Nuclear Security Administration, a project was undertaken to determine if first-principles modeling and simulation could be used to accurately predict the response of radiation detectors in urban environments. This study demonstrates that the simulated detector responses to photon radiation from both an isotopic source and the background compare well to benchmark-quality measurements in a large urban-like test environment. Simulation results using a traditional combinatorial solid geometry and a CAD model generated from LiDAR measurements gave very similar detector responses. With this validation, further simulations could be used to predict detector responses for various situations in real urban settings.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Simulation of Low-Rm physics in complex geometries on GPUs with LGR

Efficient modeling of low magnetic Reynolds number (low-Rm) magnetohydrodynamics is often challenging and requires the implementation of innovative techniques to avoid key barriers experienced with prior approaches. We detail a new paradigm for first-principles simulation of the solution to the low-Rm governing equations in complex geometries. As a result of a number of innovative numerical advances, the next-generation GPU (graphics processing unit) accelerated physics code LGR has been successfully applied to the modeling of exploding wire problems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗