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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

All-order splits and multi-soft limits for particle and string amplitudes

The most important aspects of scattering amplitudes have long been thought to be associated with their poles. But recently a very different sort of “split” factorizations for a wide range of particle and string tree amplitudes have been discovered away from poles. In this paper, we give a simple, conceptual origin for these splits arising from natural properties of the binary geometry of the curve integral formulation for scattering amplitudes for Tr(Φ 3 ) theory. The most natural way of “joining” smaller surfaces to build larger ones directly produces a choice of kinematics for which higher amplitudes factor into lower ones. This gives a generalization of splits to all orders in the topological expansion. These splits allow us to access and compute loop-integrated multi-soft limits for particle and string amplitudes, at all loop orders. This includes split factorizations and multisoft limits for pion and gluon amplitudes, that are related to Tr(Φ 3 ) theory by a simple kinematical shift.

Bosonic Strings↗
DOE OSTI2021

A fast reduced model for a shell-and-tube based latent heat thermal energy storage heat exchanger and its application for cost optimal design by nonlinear programming

Numerical simulation of latent heat thermal energy storage (LHTES) systems plays a fundamental role in studying the physical process and guiding the engineering design. Discretization of the PDEs describing the nonlinear solidification/melting process of phase change materials (PCMs) leads to a large-scale complex dynamics system, where the system behavior depends on a set of parameters. In a design setting, repeated model evaluations are required over the set of parameters results in significant computational burden. In this paper, an explicit analytic solution was built for the propagation of the solidification front in a cylindrical coordinate. The analytic solution approach is further employed to develop a low computational reduced model (RM) as a module for a shell-and-tube based LHTES heat exchanger. The levelized Cost of Energy (LCOE) is used as a design metric and the RM model is used to apply system-level constraints in the nonlinear programming formulation that facilitates efficient global optimal design of the PCM properties, flow conditions and tube geometries. The use of LCOE as the design metric prevents over design of the heat transfer rate and also establishes a fair ground for evaluation of different thermal storage technologies and their integrated applications with other systems. Optimal results showed that a higher effectiveness results in a higher LCOE; the velocity of the HTF and the length of the channel are highly correlated with each other; both larger PCM latent energy and conductivity result in lower LCOE.

25 ENERGY STORAGE↗
DOE OSTI2023

Single-stage stellarator optimization: combining coils with fixed boundary equilibria

We introduce a novel approach for the simultaneous optimization of plasma physics and coil engineering objectives using fixed-boundary equilibria that is computationally efficient and applicable to a broad range of vacuum and finite plasma pressure scenarios. Our approach treats the plasma boundary and coil shapes as independently optimized variables, penalizing the mismatch between the two using a quadratic flux term in the objective function. Four use cases are presented to demonstrate the effectiveness of the approach, including simple and complex stellarator geometries. As shown here, this method outperforms previous two-stage approaches, achieving smaller plasma objective function values when coils are taken into account.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Ice rule breakdown and frustrated antiferrotoroidicity in an artificial colloidal Cairo ice

Abstract We combine experiments and numerical simulations to investigate the low energy states and the emergence of topological defects in an artificial colloidal ice in the Cairo geometry. This type of geometry is characterized by a mixed coordination ( z ), with coexistence of both z = 3 and z = 4 vertices. We realize this particle ice by confining field tunable paramagnetic colloidal particles within a lattice of topographic double wells at a one to one filling using optical tweezers. By raising the interaction strength via an applied magnetic field, we find that the ice rule breaks down, and positive monopoles with charge q = + 2 accumulate in the z = 4 vertices and are screened by negative ones ( q = − 1 ) in the z = 3. The resulting, strongly coupled state remains disordered. Further, via analysis of the mean chirality associated to each pentagonal plaquette, we find that the disordered ensemble for this geometry is massively degenerate and it corresponds to a frustrated antiferrotoroid.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Protein drift-diffusion dynamics and phase separation in curved cell membranes and dendritic spines: Hybrid discrete-continuum methods

We develop methods for investigating protein drift-diffusion dynamics in heterogeneous cell membranes and the roles played by geometry, diffusion, chemical kinetics, and phase separation. Our hybrid stochastic numerical methods combine discrete particle descriptions with continuum-level models for tracking the individual protein drift-diffusion dynamics when coupled to continuum fields. We show how our approaches can be used to investigate phenomena motivated by protein kinetics within dendritic spines. The spine geometry is hypothesized to play an important biological role regulating synaptic strength, protein kinetics, and self-assembly of clusters. We perform simulation studies for model spine geometries varying the neck size to investigate how phase-separation and protein organization is influenced by different shapes. We also show how our methods can be used to study the roles of geometry in reaction-diffusion systems including Turing instabilities. Furthermore, our methods provide general approaches for investigating protein kinetics and drift-diffusion dynamics within curved membrane structures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Transport in Stochastic Media with Random Chord Length Distributions

Thermal radiation transport computations in binary Markovian random mixtures rely almost exclusively on the Levermore-Pomraning (LP) model which is obtained by applying a heuristic closure to the ensemble averaged random medium transport equation. The validity of this model has been extensively tested by comparing numerical results over a broad parameter range (material types and mixing parameters) against benchmark solutions in planar geometry. The conditions under which the LP-model provides useful results and when it breaks down are now well established, but work to date has been largely restricted to homogeneous mixing statistics, i.e., the mean chord lengths of both materials are taken to be spatially constant. In recent work, this limitation was relaxed by allowing the mean chord lengths and, in a consistent fashion, the volume fractions in the LP-model to vary spatially and in a follow-up investigation benchmark solutions were obtained by ensemble averaging results over material realizations sampled from a nonhomogeneous Poisson process (NHPP). Numerical experiments in rod geometry with specifically linear and quadratic spatial dependence of chord lengths showed that the material averaged radiation intensities vary nonmonotonically with depth into the medium, in stark contrast to solutions obtained assuming uniform chord lengths. Moreover, depending on the local optical thickness and strength of scattering, the LP-model results showed locally more nuanced deviations from the benchmark solutions than was the case with constant chord lengths. These limited numerical investigations highlight the nontrivial qualitative and quantitative consequences of nonhomogeneous mixing statistics, in particular that closure approximations may not be uniformly valid or invalid over the problem domain.

42 ENGINEERING↗
DOE OSTI2023

Athena-I Modeling with MCNP6.3 [Poster]

This objective of this project is to develop a model of the Athena-I experiment by creating unstructured mesh (UM) geometries in CUBIT and to perform calculations of the models in Monte Carlo N-Particle (MCNP) transport code version 6.3 (MCNP6.3).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Expanded verification and validation studies of hypersonic aerodynamics with multiple physics-fidelity models

Hypersonic aerothermodynamics is an important domain of modern multiphysics simulation. The Multi-Fidelity Toolkit is a simulation tool being developed at Sandia National Laboratories to predict aerodynamic properties for compressible flows from a range of physics fidelities and computational speeds. These models include the Reynolds-averaged Navier–Stokes (RANS) equations, the Euler equations with momentum-energy integral technique (MEIT), and modified Newtonian aerodynamics with flat-plate boundary layer (MNA+FPBL) equations, and they can be invoked independently or coupled with hierarchical Kriging to interpolate between high-fidelity simulations using lower-fidelity data. However, as with any new simulation capability, verification and validation are necessary to gather credibility evidence. This work describes formal code- and solution-verification activities, as well as model validation with uncertainty considerations. Code verification activities on the MNA+FPBL model build on previous work by focusing on the viscous portion of the model. Viscous quantities of interest are compared against those from an analytical solution for flat-plate, inclined-plate, and cone geometries. The code verification methodology for the MEIT model is also presented. Test setup and results of code verification tests on the laminar and turbulent models within MEIT are shown. Solution-verification activities include grid-refinement studies on simulations that model the HIFiRE-1 wind tunnel experiments. These experiments are used for validation of all model fidelities. A thorough validation comparison with prediction error and uncertainty is also presented. Three additional HIFiRE-1 experimental runs are simulated in this study, and the solution verification and validation work examines the effects of the associated parameter changes on model performance. Finally, a study is presented that compares the computational costs and fidelities from each of the different models.

42 ENGINEERING↗
DOE OSTI2025

Microcarbonation of Naphthalene: An Experimental and Computational Study of Photoionization in Naphthalene-Carbon Dioxide Clusters

The photoionization of naphthalene (N)-carbon dioxide (CO 2 ) clusters was studied using tunable vacuum ultraviolet (VUV) radiation from a synchrotron in the photon range of 8.0 to 13.7 eV, in combination with time-of-flight mass spectrometry. Clusters of monomer, dimer, and trimer naphthalene with CO 2 (N­(CO 2 ) 0–6 , N 2 ­(CO 2 ) 0–3 , N3) were observed. The lowest-energy conformers were obtained via a conformer search, followed by geometry optimizations at the ωB97X-V2/aug-cc-pVTZ (monomer) and ωB97X-V2/aug-cc-pVDZ (dimer) levels of theory. Carbon dioxide was found to preferentially cluster on top of the naphthalene molecule (in an out-of-plane configuration). From the mass spectra, photoionization intensity curves (PICs) were constructed, and appearance energies (AEs) were determined. No substantial trend in AE was observed with increasing size of the naphthalene-carbon dioxide clusters; rather, AE oscillations around the value for pure naphthalene were observed. These AE oscillations are also observed in a recently studied naphthalene-water cluster system, though in this system a slight downward trend in AE for pure naphthalene clusters (N1/2/3/4) was observed. The differences between the two systems are attributed to differing interaction strengths. In conclusion, understanding these differences may aid in determining how photoprocessing can proceed differently depending on the dominant matrix component in interstellar ices.

Wannenmacher, Anna [Lawrence Berkeley National Lab↗
DOE OSTI2025

High-Pressure Effects on an Octa-Hydrated Curium Complex: An Experimental and Theoretical Investigation

An octa-hydrated curium compound [Cm(H 2 O) 8 ](Hdtp)(dtp)·H 2 O (Cm1,H 2 dtp = 2,3-di(tetrazol-5-yl)pyrazine) along with its lanthanide analogues [Ln(H 2 O) 8 ](Hdtp)-(dtp)·H 2 O (Ln1, Ln 3+ = La 3+ –Nd 3+ , Sm 3+ –Lu 3+ ) were synthesized and characterized using single crystal X-ray diffraction and spectroscopic methods. Bond length analysis of VIII Cm(III)–OH 2 (where VIII refers to the coordination number) was compared to VIII Ln(III)–OH 2 (Ln 3+ = Nd 3+ and Sm 3+ ), indicating similar VIII M(III)–OH 2 bond lengths owing to the similar eight-coordinate ionic radii of these VIII M(III) cations. Owing to the reduced coordination number, the VIII Cm(III)–OH 2 bond lengths were shorter than previously reported IX Cm(III)–OH 2 bonds in [Cm(H 2 O) 9 ](CF 3 SO 3 ) 3 . The octa-aquo complexes were also characterized by solid-state UV–vis–NIR spectroscopy in addition to variable-temperature and variable-pressure photoluminescence. Variable-pressure absorption spectra of Cm1 were compared with Ln1 and show that the Cm(III) f → f transitions have a stronger dependence on pressure than that observed in Ln1 (Ln 3+ = Nd 3+ and Sm 3+ ). The experimental and computational analyses reveal that the monotonic decrease in the computed energy difference between the ground state and the first excited state corresponds to the observed red shift of the photoluminescence peak. This is accompanied by a gradual reduction in the average Cm(III)–OH 2 bond length and a delocalization of spin densities, alongside an intensified interaction involving the 5f orbitals under increasing pressure. Furthermore, these changes accommodate the new geometry and collectively modify the energy landscape, resulting in peak broadening and quenching.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Multidimensional modeling of non-equilibrium plasma generated by a radio-frequency corona discharge

Low-temperature plasma (LTP) ignition concepts rely on the production of radical and charged species to speed up the onset of combustion in spark-ignition engines. These features are responsible for the superior performance of LTP igniters under extremely dilute combustion operation that is not achievable by conventional spark igniters. Additionally, LTP discharges extend the lifetime of the igniters, due to the avoidance of spark processes. For these reasons, the engine research community and the automotive industry have shown growing interest in this technology in the recent years. As of today, computational fluid-dynamics (CFD) codes typically used by the multi-dimensional engine modeling community do not have reliable models to describe LTP ignition processes. One key missing piece of information is the physical and chemical properties of the plasma and their effect on combustion ignition. Most non-equilibrium plasma simulations reported in literature are based on simplified, canonical geometries, with simple discharge excitation schemes. Here we conduct multi-dimensional modeling of the non-equilibrium plasma generated by an application-relevant radio-frequency (RF) corona discharge in air. Three test cases are simulated, characterized by different environmental pressure levels and peak electrode voltage values at room temperature. Streamer penetration, electron number density, atomic oxygen production, and bulk gas temperature distribution in the first 10 sinusoidal pulses are presented and discussed. This model can be used as a key tool for an in-depth understanding of RF-corona discharge for automotive applications and provides the basis for future implementations of dedicated LTP ignition models in CFD codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A Practical Analysis Framework for Assessment of Printed Circuit Heat Exchangers in High-Temperature Nuclear Service

Abstract Printed circuit heat exchangers (PCHEs) have high thermal efficiency because of the numerous minuscule channels. These minuscule channels result in a high thermal exchange area per unit volume, making PCHE a top contender for an intermediate heat exchanger (IHX) in high-temperature reactors. Thousands of minuscule channels make finite element analysis of the PCHE computationally infeasible. A two-dimensional analysis is usually performed for the PCHE core, which cannot simulate the local channel level responses reasonably because of the absence of global constraint influence. At present, there is no analysis technique available in the ASME Code or literature that is computationally efficient and suitable for engineers to estimate PCHE local responses. A novel but practical two-step analysis framework is proposed for performing PCHE analysis. In the first step, the channeled core is replaced by orthotropic solids with similar stiffness to simulate the global thermomechanical elastic responses of the PCHE. In the second step, local submodel analysis with detailed channel geometry and loading is performed using the elastic-perfectly plastic (EPP) material model. The proposed two-step analysis technique provides a unique capability to estimate the channel corner responses to be used for PCHE performance assessment. This study first developed a methodology for calculating the elastic orthotropic properties of the PCHE core. Next, the two-step analysis is performed for a realistic size PCHE core, and different issues observed in the results are scrutinized and resolved. Finally, a practical finite element analysis framework for PCHEs in high-temperature nuclear service is recommended.

Engineering↗
DOE OSTI2025

Additive Manufacturing of Leak-Free Metal Components with Thin Walls and Sealing Surfaces

Laser-powder bed fusion (L-PBF) offers the ability to print free form design components which often do not require post-processing. However, challenges arise when printing small geometries with mating surfaces. Using an AddUp FormUp 350 L-PBF machine installed at the Manufacturing Demonstration Facility (MDF) of Oak Ridge National Laboratory (ORNL), a User Agreement project was formed with intent to manufacture metal leak-free cylindrical sealing surfaces. One print of twelve 12.7mm outer diameter (OD) cylinders was designed varying wall thickness in the computer-aided design (CAD) model. A build plate was successfully printed but included two build pauses each adding about five minutes per layer. The build plate of 12 cylinders was removed from the chamber and shipped to the partner.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Analytical and numerical modeling of optical second harmonic generation in anisotropic crystals using #SHAARP package

Electric-dipole optical second harmonic generation (SHG) is a second-order nonlinear process that is widely used as a sensitive probe to detect broken inversion symmetry and local polar order. Analytical modeling of the SHG polarimetry of a nonlinear optical material is essential to extract its point group symmetry and the absolute nonlinear susceptibilities. Current literature on SHG analysis involves numerous approximations and a wide range of (in)accuracies. We have developed an open-source package called the Second Harmonic Analysis of Anisotropic Rotational Polarimetry (#SHAARP.si) which derives analytical and numerical solutions of reflection SHG polarimetry from a single interface (.si) for bulk homogeneous crystals with arbitrary symmetry group, arbitrary crystal orientation, complex and anisotropic linear dielectric tensor with frequency dispersion, a general SHG tensor and arbitrary light polarization. #SHAARP.si enables accurate modeling of polarimetry measurements in reflection geometry from highly absorbing crystals or wedge-shaped transparent crystals. The package is extendable to multiple interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Computational optimal transport for molecular spectra: The fully discrete case

In this study, the use of computational optimal transport is investigated as a tool for comparing two molecular spectra. Unlike other techniques for comparing molecular spectra in a pattern-recognition framework, transport distances simultaneously encode information about line positions and intensities. In addition, it is shown that transport distances are a useful alternative to Euclidean distances as Euclidean distances are based on line-by-line comparisons, while transport distances reflect broader features of molecular spectra and adequately compare spectra with different resolutions. This paper includes a tutorial on the use of optimal transport and investigates several well-chosen examples to illustrate the utility of computational optimal transport for comparing molecular spectra.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Uniform convergence of an upwind discontinuous Galerkin method for solving scaled discrete-ordinate radiative transfer equations with isotropic scattering

Here we present an error analysis for the discontinuous Galerkin (DG) method applied to the discrete-ordinate discretization of the steady-state radiative transfer equation with isotropic scattering. Under some mild assumptions, we show that the DG method converges uniformly with respect to a scaling parameter which characterizes the strength of scattering in the system. However, the rate is not optimal and can be polluted by the presence of boundary layers. In one-dimensional slab geometries, we demonstrate optimal convergence when boundary layers are not present and analyze a simple strategy for balance interior and boundary layer errors. Some numerical tests are also provided in this reduced setting.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Narrow Slot Gasket Grooves And Step Changes In Width

We provide corrections to the slot capacitance and inverse inductance per unit length for slot gasket groove geometries using an approximate conformal mapping approach. We also provide corrections for abrupt step changes in slot width along with boundary discontinuity conditions for implementation in the various slot models.

42 ENGINEERING↗
DOE OSTI2023

Theory and Simulations of Emergent Geometry in Quantum Gravity (Final Technical Report)

This is the final technical report for the QuantISED Award to PIs Hartman, Ginsparg, and McMahon at Cornell University, outlining the goals of the research project and summarizing its main technical achievements. These include the discovery of new non-perturbative effects in quantum gravity known as replica wormholes which play an important role in the black hole information paradox, new methods to study cosmological spacetimes using the tools of holographic duality, and new experimental results on error mitigation in quantum computing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗