SEARCH · Engineering Papers

Results for “geometry 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 667 records · Page 37

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Spatial core-edge coupling of the particle-in-cell gyrokinetic codes GEM and XGC

Two existing particle-in-cell gyrokinetic codes, GEM for the core region and XGC for the edge region, have been successfully coupled with a spatial coupling scheme at the interface in a toroidal geometry. Additionally, a mapping technique is developed for transferring data between GEM's structured and XGC's unstructured meshes. Two examples of coupled simulations are presented to demonstrate the coupling scheme. The optimization of GEM for graphics processing unit is also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

A Parallel Cut-Cell Algorithm for the Free-Boundary Grad--Shafranov Problem

A parallel cut-cell algorithm is described to solve the free-boundary problem of the Grad--Shafranov equation. The algorithm reformulates the free-boundary problem in an irregular bounded domain and its important aspects include a searching algorithm for the magnetic axis and separatrix, a surface integral along the irregular boundary to determine the boundary values, an approach to optimize the coil current based on a targeting plasma shape, Picard iterations with Aitken's acceleration for the resulting nonlinear problem, and a Cartesian grid embedded boundary method to handle the complex geometry. Here the algorithm is implemented in parallel using a standard domain-decomposition approach and a good parallel scaling is observed. Numerical results verify the accuracy and efficiency of the free-boundary Grad--Shafranov solver.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Observations and Modeling of Fiber Optic Strain on Hydraulic Fracture Height Growth in Hydraulic Fracturing Test Site 2 (HFTS-2)

Summary Understanding fracture height growth can be of great significance to optimizing field development and improving recovery. The Hydraulic Fracturing Test Site 2 (HFTS-2) has provided a unique opportunity and an advanced data set to allow us to observe and understand fracture geometries rigorously. Low-frequency distributed acoustic sensing (LF-DAS) data from a vertical well in HFTS-2 showed three key observations: (i) excessive upward height growth (>1,000 ft) and limited downward growth of the hydraulic fractures during pumping, (ii) considerable additional upward fracture height growth (~300 ft) after well shut in, and (iii) very complex LF-DAS strain rate patterns for a small fiber-to-stage offset. Advanced geomechanical modeling was performed to simulate the hydraulic fracture propagation and the resulting strain responses in the vertical direction. The modeling results demonstrated asymmetric upward and downward fracture height growths as observed in HFTS-2 with a similar upward height growth rate. Simulated waterfall plots of vertical strain rate showed distinct patterns for different fiber-to-fracture distances. The upward-growing fracture tip can be clearly identified by the interfaces between compressing and extending zones. It was also found that the complex strain rate patterns observed in HFTS-2 for small fiber-stage offsets were not caused by the mechanical layering but possibly result from the simultaneous propagation of multiple hydraulic fractures at different rates. The simulation results improved the understating of the HFTS-2 LF-DAS data, and the simulated strain rate patterns could also serve as templates for fracture height interpretation from LF-DAS data in future.

Engineering↗

DOE OSTI2021

Free Energies of Hydrated Halide Anions: High Through-Put Computations on Clusters to Treat Rough Energy-Landscapes

With a longer-term goal of addressing the comparative behavior of the aqueous halides F−, Cl−, Br−, and I− on the basis of quasi-chemical theory (QCT), here we study structures and free energies of hydration clusters for those anions. We confirm that energetically optimal (H2O)nX clusters, with X = Cl−, Br−, and I−, exhibit surface hydration structures. Computed free energies, based on optimized surface hydration structures utilizing a harmonic approximation, typically (but not always) disagree with experimental free energies. To remedy the harmonic approximation, we utilize single-point electronic structure calculations on cluster geometries sampled from an AIMD (ab initio molecular dynamics) simulation stream. This rough-landscape procedure is broadly satisfactory and suggests unfavorable ligand crowding as the physical effect addressed. Nevertheless, this procedure can break down when n≳4, with the characteristic discrepancy resulting from a relaxed definition of clustering in the identification of (H2O)nX clusters, including ramified structures natural in physical cluster theories. With ramified structures, the central equation for the present rough-landscape approach can acquire some inconsistency. Extension of these physical cluster theories in the direction of QCT should remedy that issue, and should be the next step in this research direction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Optimization using pathwise algorithmic derivatives of electromagnetic shower simulations

Among the well-known methods to approximate derivatives of expectancies computed by Monte-Carlo simulations, averages of pathwise derivatives are often the easiest one to apply. Computing them via algorithmic differentiation typically does not require major manual analysis and rewriting of the code, even for very complex programs like simulations of particle-detector interactions in high-energy physics. However, the pathwise derivative estimator can be biased if there are discontinuities in the program, which may diminish its value for applications. This work integrates algorithmic differentiation into the electromagnetic shower simulation code HepEmShow based on G4HepEm, allowing us to study how well pathwise derivatives approximate derivatives of energy depositions in a sampling calorimeter with respect to parameters of the beam and geometry. We found that when multiple scattering is disabled in the simulation, means of pathwise derivatives converge quickly to their expected values, and these are close to the actual derivatives of the energy deposition. Additionally, we demonstrate the applicability of this novel gradient estimator for stochastic gradient-based optimization in a model example.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Molecular concentration field design using closed-form steady-state solutions

Control over spatial concentration fields represents a fundamental challenge in designing synthetic biological systems and programmable soft materials. While nature creates morphogen gradients that orchestrate complex developmental processes, synthetic approaches have largely relied on empirical optimization and computationally intensive simulations. Here, we present an analytical framework for steady-state concentration fields generated by finite-sized localized sources in diffusion–degradation systems and derive closed-form solutions for one-, two-, and three-dimensional geometries. By expressing these solutions in dimensionless form, we show that gradient steepness and spatial structure are organized by the Thiele modulus, which captures the competition between diffusion and degradation length scales. The analysis reveals distinct design regimes: in degradation-dominated systems, gradient shape is governed by exponential decay and becomes dimension-independent, whereas in diffusion-dominated systems, gradient magnitude and extent follow dimension-dependent power-law scaling. Building on these results, we introduce a quantitative design strategy that uses threshold-based criteria to program concentration ranges by tuning physically accessible parameters, most directly the production rate, while holding transport and degradation properties fixed. Comparisons with numerical solutions and reported experimental systems demonstrate consistency with the predicted scaling behavior. Together, this work provides a generalizable and physically transparent framework for designing steady-state concentration fields in synthetic biological and soft matter systems, enabling predictive control of gradient-mediated organization without reliance on extensive numerical optimization.

Kim, Dong Woo [Johns Hopkins University, Baltimore↗

DOE OSTI2022

Understanding resistance oscillation in the CsV 3 Sb 5 superconductor

A recent demonstration of the periodic oscillation of resistance in the thin film of the CsV 3 ⁢Sb 5 superconductor with a hole in the film suggests that charge-4⁢ e and charge-6⁢⁢ e Cooper pairs may have condensed in this compound. While exciting, such interpretation calls for a precise determination of the effective area for the passage of Cooper pairs from one end of the lead to the other. Unlike the traditional Little-Parks effect where the rim around the hole is thin, the effective hole area is not obviously defined for the “thick-rim geometry” adopted in the experiment. Here, we note that the experiment was conducted in a regime where the superconductivity is strongly fluctuating, which motivates an analysis based on the spacetime formulation of the time-dependent Ginzburg-Landau theory. Further, we argue that under appropriate conditions, the optimal semiclassical path is not the geometrically shortest one, but the one that moves along the edge of the hole and takes advantage of the reduced fluctuations at the boundary of the hole. The condition for the crossover from the geometrically shortest path to the path that sticks to the wall is clarified. In such a scenario, the geometric area of the hole indeed emerges as the effective area for the flux, providing a theoretical justification to the interpretation given by J. Ge et al. (arXiv:2201.10352). The conclusion of our analysis may have implications for similar experiments in other superconductors where the geometry of the device is not obviously that of the Little-Parks experiment employing the thin wall, but of the thick-rim type such as used in the experiment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Predictive binder jet additive manufacturing enabled by clean burn-off binder design

Binder jet additive manufacturing (BJAM) presents an avenue for advanced manufacturing of various high-value materials due to high deposition rates, scalability, and geometric flexibility. However, conventional organic binders in BJAM introduce residual carbon upon pyrolysis, often leading to imprecise alloy composition in the final sintered part. The undesirable residual carbon from binder burn-off limits the application of BJAM for high-performance alloys due to their high sensitivity to carbon addition. In this study, we have designed poly(vinylpyrrolidone-co-vinyl acetate) (PVP-VAc) as a clean burn-off binder for BJAM, where excess oxygen groups in VAc enable cleaner burn-off and reduce residual carbon retention. Compared to a widely used commercial binder, the optimized PVP-VAc binder reduced residual carbon retention by over 90% in H13 tool steel. The significant reduction in residual carbon enables predictable printing and subsequent sintering of complex H13 tool steel geometries, an alloy known to have substantial challenges around distortion due to the sintering window shifting from carbon addition. The design of a clean burnout binder provides a major path forward for BJAM by enabling new AM designs and applications for composition-sensitive high-performance alloys, such as nickel-based superalloys, titanium alloys, and high alloyed steels.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Nature-inspired lotus-shaped fins combined with hybrid nanoparticles and metal foam for high-performance latent heat thermal energy storage

Latent heat thermal energy storage (LHTES) systems play a critical role in renewable energy integration by providing high energy density and nearly isothermal operation during phase transitions. However, their performance is often limited by slow melting/charging rates, which motivates the search for enhanced heat transfer designs. This study investigates the melting behavior of RT-82 phase change material (PCM) using novel lotus-shaped fins combined with copper metal foam and conductive graphene nanoparticles and carbon nanotubes. A two-dimensional enthalpy-porosity model in ANSYS Fluent was developed to simulate the charging/melting process, capturing non-thermal equilibrium between the foam and PCM/nano-PCM. In this study, effects of fin geometry, nanoparticle concentration, and foam porosity on melting dynamics and cost-performance trade-offs were investigated. Results showed that natural convection accelerated melting by ~12% compared to conduction-only scenarios. Optimized lotus-shaped fins with higher fin density (T3F4 and T3F10) achieved up to 63% faster melting relative to sparse configurations. Graphene nanoparticles improved thermal conductivity, with a 6% volume fraction, by reducing melting time by ~6.9%, while their combination with 75% porosity foam achieved a maximum reduction in the melting time of ~51% compared to pure PCM. Cost-performance analysis identified T3F4 as the most balanced design, offering rapid thermal response without excessive material costs, while moderate-density designs like T3S6 provided economical alternatives with acceptable performance. These results highlight the performance enhancement that can be achieved by integrating bio-inspired fins, nanoparticles, and foams, into compact and efficient LHTES for solar heating, building thermal management, and industrial waste-heat recovery applications.

25 ENERGY STORAGE↗

DOE OSTI2026

Metal additively manufactured wavy fin cold-plate architecture for improved thermal-hydraulic performance

Rapid growth in artificial intelligence and data center workloads demands high-performance liquid cooling to manage increasing chip power. This study presents two metal-additive-manufactured cold plates with sinusoidal fins, constant-amplitude wavy fins and linearly variable-amplitude wavy fins and compares them against metal-additive-manufactured straight fins using experiments conducted at 1 kW heat dissipation as well as high-fidelity 3D conjugate computational fluid dynamic simulations. The cold plates were printed in AlSi10Mg material and underwent design using a Python-automated workflow prior to manufacture and testing. The experiments show that wavy fins reduce the normalized thermal resistance by 35 to 45 % at water flow rates from 1 to 4 LPM. At a fixed 20 kPa pressure drop, the variable-waviness design lowered peak surface temperature by 9 °C and thermal resistance by 51 %, while edge-channel maldistribution in the constant wavy fin design limited gains. A thermal resistance breakdown revealed that 55–63 % of the total thermal resistance in wavy designs comes from base heat conduction, 27–33 % from fin heat conduction, and 9–13 % from fin heat convection, indicating the need to address conduction bottlenecks. Parametric sweeps identify a 3 mm fin pitch as optimal, and that horizontal inlet/outlet manifolds further reduce pressure drop by 30–60 % and thermal resistance by 9–16 % relative to vertical inlet-outlet manifolds. The results yield comprehensive guidelines for fin geometry, manifold alignment, material selection and additive-manufacturing constraints to realize high-performance liquid-cooled cold plates for power-dense electronics.

3d printing↗

DOE OSTI2022

Design of inertial fusion implosions reaching the burning plasma regime

Abstract In a burning plasma state 1–7 , alpha particles from deuterium–tritium fusion reactions redeposit their energy and are the dominant source of heating. This state has recently been achieved at the US National Ignition Facility 8 using indirect-drive inertial-confinement fusion. Our experiments use a laser-generated radiation-filled cavity (a hohlraum) to spherically implode capsules containing deuterium and tritium fuel in a central hot spot where the fusion reactions occur. We have developed more efficient hohlraums to implode larger fusion targets compared with previous experiments 9,10 . This delivered more energy to the hot spot, whereas other parameters were optimized to maintain the high pressures required for inertial-confinement fusion. We also report improvements in implosion symmetry control by moving energy between the laser beams 11–16 and designing advanced hohlraum geometry 17 that allows for these larger implosions to be driven at the present laser energy and power capability of the National Ignition Facility. These design changes resulted in fusion powers of 1.5 petawatts, greater than the input power of the laser, and 170 kJ of fusion energy 18,19 . Radiation hydrodynamics simulations 20,21 show energy deposition by alpha particles as the dominant term in the hot-spot energy balance, indicative of a burning plasma state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

CHESS: The future direct geometry spectrometer at the second target station

CHESS, chopper spectrometer examining small samples, is a planned direct geometry neutron chopper spectrometer designed to detect and analyze weak signals intrinsic to small cross sections (e.g., small mass, small magnetic moments, or neutron absorbing materials) in powders, liquids, and crystals. CHESS is optimized to enable transformative investigations of quantum materials, spin liquids, thermoelectrics, battery materials, and liquids. The broad dynamic range of the instrument is also well suited to study relaxation processes and excitations in soft and biological matter. The 15 Hz repetition rate of the Second Target Station at the Spallation Neutron Source enables the use of multiple incident energies within a single source pulse, greatly expanding the information gained in a single measurement. Furthermore, the high flux grants an enhanced capability for polarization analysis. This enables the separation of nuclear from magnetic scattering or coherent from incoherent scattering in hydrogenous materials over a large range of energy and momentum transfer. This paper presents optimizations and technical solutions to address the key requirements envisioned in the science case and the anticipated uses of this instrument.

47 OTHER INSTRUMENTATION↗

DOE OSTI2022

Preparation for a Measurement of Charge Asymmetry in the Bethe-Heitler Process

We have prepared a measurement of the energy asymmetry in wide- and medium-angle electron/positron pair production off protons and heavy targets. This asymmetry is caused by the interference between the first- and second-order Born diagrams and the Compton scattering diagram. It directly probes aspects of QED, as well as providing a direct measurement of the real part of the Compton amplitude. It will be conducted at the HI??S facility at Duke University, using a 60 MeV photon beam. This dissertation serves as documentation of the preparation stage of the Bethe-Heitler experiment. The major was the recommissioning of the vertical drift chambers previously used in the Q-weak experiment at the Jefferson Lab. Cosmic test runs were conducted, drift time data were collected and efficiency plateaus were measured. We made modifications to the JLAB Hall A analyzer to suit the geometry and drift characteristics of these wire chambers. The analyzer was used for the reconstruction of the trajectories of cosmic ray test runs with the results confirmed by direct measurement of trigger geometry. Spatial and angular resolution is estimated to ~300?? and 0.17° respectively. Geant 4 simulations with generated Bethe-Heitler pairs satisfying theoretical differential cross sections. It was used to check detector acceptance, optimize apparatus layout, and estimate measurable energy asymmetry. The measurable asymmetries from electron/positron pairs with polar angles around between approximately 5° and 8°, azimuthal angles differing by 180°, and energy differing by approximately 9 MeV to 15 MeV are predicted to be above 10%. The kinematics of primary vertices are reconstructed using the data from wire chambers in the simulation. The energy resolution is determined to be better than 1MeV.

Chen, Haoyu↗

DOE OSTI2022

Fourier analysis of continuous fractional diffusion synthetic acceleration schemes in slab geometry

We propose two fractional extensions of continuous diffusion synthetic acceleration (DSA) with fractional derivative order α varying over the interval 2 ≥ α ≥ 1 . We investigate the spectral properties of the corresponding continuous families of fractional preconditioners by performing Fourier analysis for a model infinite homogeneous medium problem in slab geometry. The first family results in a fractional acceleration scheme, FrDSAo, that reduces to traditional DSA for .α = 2 and scattering ratio c limiting to a unit value (c → 1) but is otherwise optimized via the Fourier analysis, to obtain the smallest possible spectral radius, for c < 1 and 2 ≥ α ≥ 1. The second family corresponds to a fractional acceleration scheme, FrDSAs, that reduces to traditional DSA for α = 2 for all values of c. The latter scheme is not optimized but has the advantage of lending itself to a more straightforward implementation. For high values of c, the results of the Fourier analysis point to the existence of an interval 2 > α > ∼1.8 where both FrDSAo and FrDSAs can achieve a lower spectral radius than DSA. For example, DSA has a spectral radius of ∼0.2246 for c = 0.9999 while FrDSAo produces a value of ∼0.1616 at α = 1.92 and FrDSAs results in ∼0.2116 at α =1.93. (author)

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Reduced Geometry for the NB5 Velocity Selector Shielding Box

This report details the reduced NB5 velocity selector shielding box that has been isolated from the complete HFIR CGH model in order to speed up detailed calculations of the velocity selector shielding. This shielding model is intended for outsourcing shielding optimization calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Predictive process mapping for laser powder bed fusion: A review of existing analytical solutions

One of the main challenges in the laser powder bed fusion (LPBF) process is making dense and defect-free components. These porosity defects are dependent upon the melt pool geometry and the processing conditions. Power-velocity (PV) processing maps can aid in visualizing the effects of LPBF processing variables and mapping different defect regimes such as lack-of-fusion, under-melting, balling, and keyholing. This work presents an assessment of existing analytical equations and models that provide an estimate of the melt pool geometry as a function of material properties. The melt pool equations are then combined with defect criteria to provide a quick approximation of the PV processing maps for a variety of materials. Finally, the predictions of these processing maps are compared with experimental data from the literature. Here, the predictive processing maps can be computed quickly and can be coupled with dimensionless numbers and high-throughput (HT) experiments for validation. The present work provides a boundary framework for designing the optimal processing parameters for new metals and alloys based on existing analytical solutions.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Microstructural Development in Inconel 718 Nickel-Based Superalloy Additively Manufactured by Laser Powder Bed Fusion

Excellent weldability and high temperature stability make Inconel 718 (IN718) one of the most popular alloys to be produced by additive manufacturing. Here in this study, we investigated the effects of laser powder bed fusion (LPBF) parameters on the microstructure and relative density of IN718. The samples were fabricated with independently varied laser power (125 - 350W), laser scan speed (200 - 2200 mm/s), and laser scan rotation (0° - 90°). Archimedes’ method, optical microscopy, and scanning electron microscopy were employed to assess the influence of LPBF parameters on the relative density and microstructure. Optimal processing windows were identified for a wide range of processing parameters and relative density greater than 99.5% was achieved using volumetric energy density between 50 and 100 J/mm3. Microstructural features including melt pool geometry, lack of fusion defect, keyhole porosity, and sub-grain cellular microstructure were examined and quantified to correlate to LPBF parameters. A simple empirical model was postulated to relate relative sample density and LPBF volumetric energy density. Melt pool dimensions were quantitatively measured and compared to estimations based on Rosenthal solution, which yielded a good agreement with the width, but underestimated the depth, particularly at high energy input, due to lack of consideration for keyhole mode. In addition, the sub-grain cellular-dendritic microstructure in the as-built samples was observed to decrease with increasing laser scan speed. Quantification of the sub-micron cellular-dendritic microstructure yielded estimated cooling rate in the order of 10 5 – 10 7 K/s.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2026

Bio-inspired alula-based winglet design for enhanced heat transfer in high temperature fin-and-tube heat exchangers

Fin-and-tube heat exchangers (FTHEs) are widely used for high-temperature flue-gas heat recovery, but their performance is often limited by wake regions and non-uniform fin-surface temperatures. This study proposes and numerically evaluates four bio-inspired longitudinal vortex generator (VG) configurations in a high-temperature FTHE with flue-gas inlet temperature ∼1230 K: double-delta, curved double-delta, alula, and a new curved-alula geometry. The reference fin is not hydraulically plain; it already incorporates leading-edge separation columns and convex protrusions, so the alula-type winglets are assessed as downstream add-ons acting on a strongly disturbed flow. In a second step, perforations (one, two and three circular holes) are introduced into the curved-alula VGs to further tailor the flow field. Three-dimensional simulations with the Shear Stress Transpor (SST) $k - ω$ model, temperature-dependent flue-gas properties and conjugate conduction are carried out for gas-side Reynolds numbers $Re_g ≈ 8.0$ x $10^2 - 3.6$ x $10^3$ (mass flow rates 0.5 – 2.5 g/s), and the designs are compared in terms of surface heat flux, Nusselt number, friction factor and hydrothermal performance factor (HTPF). For this already-promoted fin, the additional downstream winglets provide moderate, incremental hydrothermal gains. At the highest Reynolds number, the best non-perforated design (curved-alula) increases surface heat flux from 1630.9 to 1794.7 kW/m² (∼ 10 % gain) and the Nusselt number from 227.6 to 242.6 (∼ 7 % gain), while the friction factor rises from 0.26 to about 0.30, yielding HTPF values close to unity (∼ 0.9 – 1.0). Introducing circular perforations into the curved-alula winglets acts mainly as a wake-bleeding refinement: the three-hole configuration provides a heat flux of 1824.7 kW/m² and a pressure drop of 127.9 Pa, with HTPF in the range ∼ 1.03 – 1.14 and a small (∼ 1 – 3 %) improvement over the solid curved-alula design. Flow-field analysis shows that the perforated curved-alula VGs shrink tube-wake regions, thin the thermal boundary layer and homogenize the fin-surface temperature (outlet-gas temperature ∼ 510 – 520 K and fin-surface temperature ∼ 420 – 421 K for the three-hole case). An optimal flue-gas mass flow rate of ∼ 1 g/s ($Re_g ≈ 1.5$ x $10^3$) is identified, beyond which additional heat-transfer gains are offset by rapidly increasing pressure losses. Overall, the results highlight that initial fin geometry and VG placement are as important as VG shape: alula-based winglets are expected to yield larger relative gains on simpler flat-fin layouts or when positioned closer to the fin leading edge and tube

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗