Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Complex 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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

MOOSE-based Tritium Migration Analysis Program, Version 8 (TMAP8) for advanced open-source tritium transport and fuel cycle modeling

Tritium management is critical for the safety, sustainability, and economics of fusion energy systems, and advanced and reliable modeling tools help accelerate the development of tritium technologies. This paper presents the Tritium Migration Analysis Program, Version 8 (TMAP8), an open-source, MOOSE-based application developed to provide state-of-the-art tritium transport and fuel cycle modeling capabilities. TMAP8 aims to expand the capabilities of previous versions (i.e., TMAP4 and TMAP7) by leveraging modern computational techniques, ensuring high software quality assurance standards (key to building trust), and enabling multispecies, multiscale, and multiphysics simulations for integrated tritium transport modeling in complex geometries. This paper outlines TMAP8’s scope and rigorous development practices, emphasizing its transparency, accessibility, modularity, and reliability. We present the current suite of verification and validation cases based on those from TMAP4, demonstrating TMAP8’s accuracy and reliability against analytical solutions and experimental data. Additionally, the paper showcases TMAP8’s integrated fuel cycle modeling capabilities, highlighting its applicability at various scales and levels. The TMAP8 code and documentation are openly available, promoting collaborative development and widespread adoption within the fusion community. Future work will soon expand TMAP8’s verification and validation suite to include those from TMAP7 and other recent experimental studies for validation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

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↗
DOE OSTI2026

Predicting roughness effects in additively manufactured coolant channels with helical enhancements

Additive manufacturing (AM) is a promising technique for fabrication of complex geometries such as those expected to be utilized in the blanket, first wall, and divertor. In the case of cooling, metallic AM may be exploited to embed geometric enhancements (ribs, rifling, etc.) to improve cooling performance. However, due to the roughness of these unfinished internal AM surfaces, prediction of thermal hydraulic performance in such channels is difficult. In this work, we consider a methodology for predicting pressure drop and heat transfer in AM channels containing helical enhancements (e.g. helical ribs, twisted tapes) that allows the incorporation of roughness data through conventional pipe flow correlations. This methodology is tested using experimental friction factor and heat transfer coefficient data from high-pressure helium coolant flow measurements in AM stainless steel tubes fabricated by laser powder bed fusion. Both a featureless AM tube and one containing helical ribs were considered alongside a conventionally manufactured smooth tube. The AM surface roughness is obtained by profilometry and used to predict an equivalent sand-grain roughness, with this equivalent roughness confirmed through AM featureless tube measurements. Under the proposed methodology, this roughness information is incorporated into predictions of friction factor and Nusselt number for the rifled tube. Furthermore, these predictions agree well with experimental data across a large range of Reynolds numbers, encouraging the use of this methodology for thermal hydraulic analysis of similar systems and design of future coolant channel geometries.

Additive manufacturing↗
DOE OSTI2024

A robust spectral element implementation of the $k - τ$ RANS model in Nek5000/NekRS

The $k - ω$ Reynolds Averaged Navier Stokes (RANS) model is one of the industry standard approaches for modeling of turbulent flows. It performs better than the $k - ϵ$ model for low Reynolds number flows and is also more suitable for boundary layers with adverse pressure gradients. Major drawback of the model, however, is that the asymptotic value of $ω$ at the walls is singular, necessitating the use of a contrived “sufficiently” large value for $ω$ as the boundary condition for its transport equation. Here, this invariably leads to the solution being sensitive to near wall grid spacing. While an acceptable solution for low order (finite volume) methods, the excessive near wall gradients lead to persistent numerical stability issues in high order codes. To alleviate the problem, specifically in the context of the high order spectral element code Nek5000, a regularized $k - ω$ approach was formulated in our prior work (Tomboulides et al., 2018). The formulation, however, relies on the use of wall distance and its gradients for modeling the closure terms and can pose problems for simulations in complex geometries. This work presents a novel implementation of the $k - τ$ RANS model in Nek5000, where $τ = 1/ω$, eliminating the need for regularization, owing to the asymptotically bounded behavior of the source terms in the $τ$ transport equation, and also eliminating dependence on wall distance. Robustness and stability of the $k - τ$ model is ensured through implicit treatment of the source terms and their careful numerical implementation and demonstrated through several cases aimed at verification and validation. Studies include both canonical and engineering relevant problems, viz., turbulent channel flow, pipe flow, backward facing step, flow over NACA 0012 airfoil and flow in a T-junction. Results from the $k - τ$ model are shown to be consistent with regularized $k - ω$ model and also with the $k - ω$ SST model in OpenFOAM (for select studies). Comparison with experimental data is also shown, where available, to bolster validation efforts for the $k - τ$ model implementation through prediction of key turbulent quantities of interest.

Nek5000↗
DOE OSTI2022

Cost-efficient copper-nickel alloy for active cooling applications

It is recently proposed that active (Peltier) and passive (conductive) cooling could be combined to design highly-efficient thermal conductors as active heat sinks. Here we study the potential of low-cost copper-nickel alloys for active cooling applications. We show that for this application, copper-nickel alloys prepared using industrial-grade powders to lower the manufacturing cost, are as good as copper-nickel alloys made out of high purity metals. Furthermore, directed energy deposition, a type of additive manufacturing, can produce copper-nickel alloys which have about 10 percent lower thermoelectric power factor compared to mechanical-alloying methods but have larger thermal conductivity values, hence, possessing similar effective thermal conductivity values, reaching 403 W m -1 K -1 at 390K. Much larger values are expected at higher temperatures. The usage of additive manufacturing enables complex geometries and large-scale production.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Direct observation of pore formation mechanisms during LPBF additive manufacturing process and high energy density laser welding

Laser powder bed fusion (LPBF) is a 3D printing technology that can print parts with complex geometries that are unachievable by conventional manufacturing technologies. However, pores formed during the printing process impair the mechanical performance of the printed parts, severely hindering their widespread application. Here, we report six pore formation mechanisms that were observed during the LPBF process. Our results in this study reconfirm three pore formation mechanisms - keyhole induced pores, pore formation from feedstock powder and pore formation along the melting boundary during laser melting from vaporization of a volatile substance or an expansion of a tiny trapped gas. We also observe three new pore formation mechanisms: (1) pore trapped by surface fluctuation, (2) pore formation due to depression zone fluctuation when the depression zone is shallow and (3) pore formation from a crack. The results presented here provide direct evidence and insight into pore formation mechanisms during the LPBF process, which may guide the development of pore elimination/mitigation approaches. Since certain laser processing conditions studied here are similar to the situations in high energy density laser welding, the results presented here also have implications for laser welding.

42 ENGINEERING↗
DOE OSTI2023

Microstructure and high temperature properties of tungsten processed via electron beam melting additive manufacturing

There is considerable interest in the adoption of additive manufacturing for processing refractory metals. The layer-wise fabrication approach enables opportunities for producing complex geometries which cannot be otherwise be achieved via powder metallurgy. However, the processing science is still in its nascent stages and structure–property relations are relatively unexplored. Fundamental research is needed to further develop the technology and enable the fabrication of refractory metals for high temperature applications. Here we focus on the processing of pure tungsten using electron beam melting additive manufacturing. Experimentally we develop a suitable processing window for achieving high density crack free material. Microstructural analysis reveals that the microstructure generally consists of a columnar structure with a (111) build direction fiber preference, although, fiber switching was observed. Process induced deformation is believed to drive the formation of subgrains whose boundaries exhibit a high dislocation density. Additionally, high temperature tensile testing reveals that the material exhibits excellent properties closer to that of annealed tungsten. Significant mechanical anisotropy was observed to be present which is likely driven by strong crystallographic texture.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Foundations for high-order, conservative cut-cell methods: Stable discretizations on degenerate meshes

Cut-cell methods for unsteady flow problems can greatly simplify the grid generation process and allow for high-fidelity simulations on complex geometries. However, cut-cell methods have been limited to low orders of accuracy. This is driven, largely, by the variety of procedures typically introduced to evaluate derivatives in a stable manner near the highly irregular embedded geometry. Here, a completely new approach, termed TEMO (truncation error matching and optimization), is taken to solve this problem. The approach is based on two simple and intuitive design principles. These principles directly allow for the construction of stable 8th To the best of the authors' knowledge, these are the highest orders ever achieved for a cut-cell discretization by a significant margin. This is done for both explicit and compact finite differences and is accomplished without any geometric transformations or artificial stabilization procedures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Bayesian sparse learning with preconditioned stochastic gradient MCMC and its applications

Deep neural networks have been successfully employed in an extensive variety of research areas, including solving partial differential equations. Despite its significant success, there are some challenges in effectively training DNN, such as avoiding overfitting in over-parameterized DNNs and accelerating the optimization in DNNs with pathological curvature. Here, we propose a Bayesian type sparse deep learning algorithm. The algorithm utilizes a set of spike-and-slab priors for the parameters in the deep neural network. The hierarchical Bayesian mixture will be trained using an adaptive empirical method. That is, one will alternatively sample from the posterior using preconditioned stochastic gradient Langevin Dynamics (PSGLD), and optimize the latent variables via stochastic approximation. The sparsity of the network is achieved while optimizing the hyperparameters with adaptive searching and penalizing. A popular SG-MCMC approach is Stochastic gradient Langevin dynamics (SGLD). However, considering the complex geometry in the model parameter space in nonconvex learning, updating parameters using a universal step size in each component as in SGLD may cause slow mixing. To address this issue, we apply a computationally manageable preconditioner in the updating rule, which provides a step-size parameter to adapt to local geometric properties. Moreover, by smoothly optimizing the hyperparameter in the preconditioning matrix, our proposed algorithm ensures a decreasing bias, which is introduced by ignoring the correction term in the preconditioned SGLD. According to the existing theoretical framework, we show that the proposed algorithm can asymptotically converge to the correct distribution with a controllable bias under mild conditions. Numerical tests are performed on both synthetic regression problems and learning solutions of elliptic PDE, which demonstrate the accuracy and efficiency of the present work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A sharp immersed method for 2D flow-body interactions using the vorticity-velocity Navier-Stokes equations

Immersed methods discretize boundary conditions for complex geometries on background Cartesian grids. Here, this makes such methods especially suitable for two-way coupled flow-body problems, where the body mechanics are partially driven by hydrodynamic forces. However, for the vorticity-velocity form of the Navier-Stokes equations, existing immersed geometry discretizations for two-way coupled problems only achieve first order spatial accuracy near solid boundaries. Here we introduce a sharp-interface approach based on the immersed interface method to handle the one- and two-way coupling between an incompressible flow and one or more rigid bodies using the 2D vorticity-velocity Navier-Stokes equations. Our main contributions are three-fold. First, we develop and analyze a moving boundary treatment for sharp immersed methods that can be applied to PDEs with implicitly defined boundary conditions, such as those commonly imposed on the vorticity field. Second, we develop a two-way coupling methodology for the vorticity-velocity Navier-Stokes equations based on control-volume momentum balance that does not require the pressure field. Third, we show through extensive testing and validation that our resulting flow-body solver reaches second-order accuracy for most practical scenarios, and provides significant efficiency benefits compared to a representative first-order approach.

42 ENGINEERING↗
DOE OSTI2021

Review of additive manufacturing and densification techniques for the net- and near net-shaping of geometrically complex silicon nitride components

Silicon nitride (Si 3 N 4 ) is an advantageous material due its unique combination of mechanical, thermal, chemical, and electrical properties both at ambient and elevated temperatures. Because of these properties there are a wide range of applications for Si 3 N 4 components. Applications include heat exchangers, environmental barrier coatings, osteointegration scaffolds, radomes, and integrated circuitry. Furthermore, such applications often require geometric complexity for efficient and/or effective operation. However, traditional ceramics processing methods such as hot-pressing or die extrusion are typically limited to simple axis-symmetric shapes. With the advent of additive manufacturing, there has been significant advancement into the forming of geometrically complex Si 3 N 4 components. This review documents additive manufacturing advancements that have demonstrated, or are capable of, fabricating Si 3 N 4 components with complex geometry.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Fabrication of yttrium hydride for high-temperature moderator application

We report given the superior thermal stability and highly attainable hydrogen density, yttrium hydride is an excellent high-temperature moderator material in advanced thermal neutron spectrum reactors that require small core volumes. Large-scale, crack-free, bulk yttrium hydride is in high demand; however, fabrication of yttrium hydride is challenging and has not been demonstrated for nearly half century. The associated challenges are discussed herein. In response to these challenges, a hydriding system was designed and constructed at Oak Ridge National Laboratory and was used to successfully fabricate crack-free yttrium hydride in complex geometries at large scales. This was accomplished by precisely controlling the hydrogen’s partial pressure and the retort temperature, which was informed by the well-established thermodynamic properties of the binary H–Y system. Hydrogen content in as-fabricated hydride was determined by weight change, an approach which was considered reliable based on the use of ultra–high-purity yttrium, the absence of oxide phases up to levels detectable using x-ray diffraction (XRD), and the significant weight gain. Hydrogen distribution along one yttrium hydride rod was evaluated with XRD analysis on materials extracted from different locations on the rod. The results indicated a relatively homogeneous hydrogen distribution along the hydride rod, with <3% uncertainty in the fraction of the δ-phase hydride. In addition, significant efforts are being dedicated to establish a complete database summarizing the thermomechanical and physical properties of as-fabricated yttrium hydride and the irradiation response to facilitate its deployment as a high-temperature moderator in advanced nuclear reactors.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Affinity propagation clustering of full-field, high-spatial-dimensional measurements for robust output-only modal identification: A proof-of-concept study

Determination of the model order is a challenging problem in system identification, especially in output-only or operational modal identification where some modes are weakly excited. While existing methods such as the stabilization diagram method (spectral information) are effective, they do not scale to high-dimensional data, which is usually needed for high-fidelity characterization of structural dynamics and has been made available in the emerging full-field measurement techniques using optical methods such as photogrammetry and laser vibrometers. In this proof-of-concept study we present a new non-parametric, data-driven approach for robust output-only identification of high-spatial-dimensional modal parameters of basic structures by efficiently processing and interactively exploiting the full-field measurement (i.e., very dense spatial measurement points). Specifically, we first over-estimate the system model once, producing a pool of candidate modes associated with their modal frequencies and full-field, high-spatial-dimensional mode shapes. This is accomplished by a data-driven method termed affinity propagation clustering (APC), where the active clusters, which are the active modes in our formulations, emerge from the "message-passing" procedure and does not require a pre-determination of the cluster number (mode or model order). Next, rather than using the spectral information to distinguish the physical and spurious modes in the stabilization diagram method, we exploit and visualize the spatial, full-field mode shape associated with each candidate mode to do so. We conduct extensive experiments on basic structural models with comparisons to a few existing methods. The results indicate that the new method is computationally efficient for identifying high-spatial-dimensional modal parameters, and robust to identify weak modes by exploiting the full-field measurement. We also discuss its applicability and limitations for structures with complex geometry (shapes).

42 ENGINEERING↗
DOE OSTI2021

A Poisson equation method for prescribing fully developed non-Newtonian inlet conditions for computational fluid dynamics simulations in models of arbitrary cross-section

Prescribing inlet boundary conditions for computational fluid dynamics (CFD) simulations of internal flow in complex geometries such as anatomical vascular models is challenging. In the absence of patient-specific inlet velocity data, a common approach for long blood vessels is to assume that the inlet flow is fully developed. In vessels of irregular cross section, however, prescribing fully developed conditions is complicated due to the lack of a general closed-form analytical solution. In this study, we develop a simple Poisson equation method for prescribing fully developed inlet conditions for the flow of either Newtonian or non-Newtonian fluids in CFD models of arbitrary cross-section. We first derive the generalized Poisson equation for fully developed flow of a non-Newtonian fluid and we then develop and verify a methodology for numerically computing the solution on any planar boundary domain. In addition, we develop a simple extension of the method for prescribing a non-orthogonal inlet velocity that represents fully developed flow from an upstream tube that is connected to the CFD inlet at a non-orthogonal angle. This may be used to investigate a common source of uncertainty in CFD simulations of internal flow that is due to a lack of information concerning the exact streamwise flow direction at the inlets. Comparison to several Newtonian and non-Newtonian benchmark verification solutions shows the method to be extremely accurate. As a practical demonstration case, we use the method to prescribe fully developed conditions on multiple non-circular inlets for the non-Newtonian flow of blood in a patient-specific model of the inferior vena cava (IVC). Finally, we further demonstrate the utility of the method by performing a sensitivity study using the patient-specific IVC model, wherein we investigate the influence of inlet velocity flow direction on the non-Newtonian IVC hemodynamics. Given its simplicity and computational efficiency, the method is shown to be far superior to alternative approaches for prescribing fully developed inlet conditions in such complicated geometries. In conclusion, to facilitate the adoption of our Poisson equation method, we have distributed our OpenFOAM source code and the associated test cases from this study as open-source software.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Texture development in magnetostrictive Fe-Ga alloys processed by laser powder bed fusion

Iron-gallium (Fe-Ga, Galfenol) alloys are promising magnetostrictive materials for actuators, sensors, and energy harvesting, but their performance is highly sensitive to microstructure and texture. Additive manufacturing by laser powder bed fusion (LPBF) offers a pathway to engineer texture and integrate functional materials into complex geometries. Here, we fabricate Fe-Ga alloys (Fe 82.2 Ga 17.8 ) by LPBF of gas-atomized powders and systematically optimize laser power and scan speed to maximize density and control texture. Nearly full-density parts (up to 99.6 %) are achieved within a narrow processing window. Electron backscatter diffraction (EBSD) reveals a strong <100> fiber texture aligned with the build direction and columnar grains up to 1 mm long. Magnetostriction measurements show saturation magnetostriction of 190 ppm in the build direction. Correlating texture data with macroscopic magnetostriction, we estimate intrinsic magnetostriction constants (λ 100 = 228 ppm, λ 111 = 12 ppm), closely matching single crystal-derived values. These results demonstrate the critical interplay between processing, texture, and functional performance in additively manufactured Fe-Ga alloys and establish LPBF as a viable route for high-performance magnetostrictive materials.

Additive manufacturing↗
DOE OSTI2022

Volumetric energy density impact on mechanical properties of additively manufactured 718 Ni alloy

The 718 Ni-based superalloy has gained enormous attention in the additive manufacturing community for its great weldability, allowing for complex geometries to be formed, and its superb mechanical strength. Here we explore the influence of a wide range of major build parameters on microstructures and mechanical properties of 718 Ni alloys. Volumetric energy density appears to be one of the major parameters that captures the relationship between energy input and measured mechanical properties and microstructures of the alloys. A threshold energy density was identified below which the properties of the additively manufactured specimen degrade substantially. Finally, the influence of energy density on the quality of the built specimen simulated using Flow3D AM correlates well with the experimental results.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A lightweight Fe–Mn–Al–C austenitic steel with ultra-high strength and ductility fabricated via laser powder bed fusion

Lightweight Fe–Mn–Al–C steels have become a topic of significant interest for the defense and automotive industries. These alloys can maintain high strength and ductility while also reducing weight in structural applications. Conventionally processed Fe–Mn–Al–C austenitic steels with high Al content (~9 wt%) demonstrate greater than 1.5 GPa strength with 35% elongation. Several recent studies have demonstrated success in fabricating steel parts using laser powder bed fusion (L-PBF) additive manufacturing (AM), which can generate near-net-shape components with complex geometries and is capable of local microstructural control. However, studies on L-PBF processing of Fe–Mn–Al–C alloys have focused on low Al content (<5 wt%) compositional regimes representing alloys that undergo transformation-induced plasticity (TRIP) and twinning-induced plasticity (TWIP). Here, in this study, we present the effects of L-PBF processing on the microstructure and mechanical properties of an Fe–30Mn–9Al–1Si-0.5Mo-0.9C austenitic steel. A process optimization framework is employed to determine an ideal L-PBF processing space that will result in >99% density parts. Implementing this framework resulted in near-fully dense specimens fabricated over a broad range of process parameters. Additionally, two bi-directional scan rotation strategies (90° and 67°) were applied to understand their effects on texture and anisotropy in this material. As-printed specimens displayed considerable work-hardening characteristics with average strengths of up to 1.3 GPa and 36% elongation in the build direction. However, solidification microcracks oriented in the build direction resulted in anisotropy in tensile strength and ductility resulting in average strengths of 1.1 GPa and 20% elongation perpendicular to the build direction. The successful L-PBF fabrication of Fe–30Mn–9Al–1Si-0.5Mo-0.9C presented here is expected to open new avenues for weight reduction in structural applications with a high degree of control over part topology.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Influence of microstructure and temperature on impact toughness of H13 steel produced by binder jet additive manufacturing

Binder Jet Additive Manufacturing (BJAM) is a promising manufacturing pathway to produce H13 steel dies and tooling with complex geometries for applications in high pressure aluminum die casting, hot stamping, and injection molding. While fully dense H13 coupons produced using BJAM have been subjected to detailed microstructure characterization, properties which are critical and relevant to the aforementioned applications, such as impact toughness, have not been reported. Here, this work evaluated the influence of microstructural characteristics and test temperature on the impact toughness of H13 produced by BJAM. Coupons were produced from three different powder size distributions (PSDs), with nominal powder size ranges of −22 μm, 10-32 μm, and 15-53 μm. Coupons were printed, sintered, hot isostatically pressed (HIPed), and heat treated by quenching and tempering. After HIPing, the measured porosities of the different PSDs were all less than 0.015 vol%. The −22 μm PSD BJAM material exhibited the best impact toughness of all three PSDs across the entire test temperature range from 25 to 400 °C, and also exhibited an ambient temperature impact toughness of 10.9 J at a hardness of 46 HRC, equivalent to minimum threshold requirements for premium grade H13. However, all three BJAM PSDs exhibited moderately lower impact toughness than premium grade wrought H13 from RT to 400 °C. We attribute this in part due to the very large prior-austenite grain (PAG) sizes near 400 μm and segregation from former liquid channels leading to greater amounts of carbide precipitation on PAG boundaries. Technical pathways for optimizing processing and chemistry of BJAM H13 to further improve impact toughness are discussed.

Binder jet additive manufacturing↗