Influence of bimodal particle distribution on material properties of BaTiO3 fabricated by paste extrusion 3D printing
Not provided.
Materials Science↗
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.
Physicochemical properties of digital light processing 3D-Printed alumina and mullite ceramics
Alumina ceramics fabricated using conventional techniques such as uniaxial pressing or injection molding are popular due to their low density, excellent insulation, and mechanical properties. However, these methods often limit the fabrication of complex geometries with high dimensional accuracy due to tooling constraints and limited design freedom. To overcome these limitations, this study employed Digital Light Processing (DLP) additive manufacturing (AM), which enables the production of precise structures. In addition, the optimization of sintering parameters to enhance the densification and performance of alumina and mullite ceramics was investigated, with a specific focus on how varying sintering temperatures and hold times affect part shrinkage, geometric accuracy, and material integrity. Flexural strength of both alumina and mullite specimens was clearly influenced by the way the layers was stacked. When layers were arranged across the direction of the applied load (Z = 4; XZ), the strength was higher than when they were stacked along the same direction as the load (Z = 3; XY). Weibull analysis based on these results showed high modulus values across all samples, indicating good reliability in the flexural strength measurements. This reliability, combined with the clear influence of sintering parameters, highlights how processing conditions effect the material properties and mechanical performance of parts produced through DLP additive manufacturing. In conclusion, these findings open new avenues for ceramic AM across various applications, enhancing the potential for innovation in fields such as aerospace, biomedical engineering, and energy.
36 MATERIALS SCIENCE↗
Enriched immersed boundary method (EIBM) for interface-coupled multi-physics and applications to convective conjugate heat transfer
It is challenging to develop numerical methods that simultaneously maintain accuracy of resolving boundary conditions and mesh flexibility to handle the interface in interface-coupled multi-physics problems involving large property discontinuity. On the one hand, boundary-fitted methods possess high accuracy in capturing the interfacial phenomenon but involve complicated volumetric mesh generation, mesh-motion, and even re-meshing procedures. On the other hand, immersed boundary methods (IBM) provide mesh flexibility but sometimes suffer from inferior interface representations, leading to poor enforcement of boundary conditions. This paper presents an enriched immersed boundary method (EIBM) to overcome this challenge and demonstrates its efficacy in conjugate heat transfer, a representative example in interface-coupled multi-physics systems that have implications for many industrial processes. Further, the core technique of the method is to enhance IBM’s accuracy of the fluid–solid interface by enriching the degrees of freedom of the cut elements to enforce temperature and flux compatibilities and resolve all the physical unknowns on the background mesh to simplify the volumetric mesh generation. We implement the EIBM under the framework of a variational multiscale formulation for coupled Navier–Stokes and thermodynamics equations. The enriched DoFs enable better enforcement of temperature and flux compatibilities with large conductivity ratios across the fluid–solid interface. At the same time, the immersed nature of the proposed method still attains mesh flexibility. The EIBM’s accuracy is thoroughly evaluated through a set of examples, ranging from benchmark problems with analytical solutions to real-world cooling processes of a moving metallic structure with complex geometry.
42 ENGINEERING↗
Defects and anomalies in powder bed fusion metal additive manufacturing
Metal additive manufacturing is a disruptive technology that is revolutionizing the manufacturing industry. Despite its unrivaled capability for directly fabricating metal parts with complex geometries, the wide realization of the technology is currently limited by microstructural defects and anomalies, which could significantly degrade the structural integrity and service performance of the product. Accurate detection, characterization, and prediction of these defects and anomalies have an important and immediate impact in manufacturing fully-dense and defect-free builds. As such, this review seeks to elucidate common defects/anomalies and their formation mechanisms in powder bed fusion additive manufacturing processes. They could arise from raw materials, processing conditions, and post-processing. While defects/anomalies in laser welding have been studied extensively, their formation and evolution remain unclear. Additionally, the existence of powder in powder bed fusion techniques may generate new types of defects, e.g., porosity transferring from powder to builds. Practical strategies to mitigate defects are also addressed through fundamental understanding of their formation. Such explorations enable the validation and calibration of models and ease the process qualification without costly trial-and-error experimentation.
42 ENGINEERING↗
OpenSn: A massively parallel, open-source simulation environment for discrete ordinates radiation transport
OpenSn is an open-source, massively parallel deterministic radiation transport code for solving the discrete-ordinates ( S N ) form of the Boltzmann transport equation on unstructured, arbitrary polyhedral meshes. It supports high-fidelity simulations involving steady-state, eigenvalue, and adjoint problems for neutral particles (e.g., neutrons, photons, multi-particles), using the multigroup approximation in energy. OpenSn combines angular discretization via discrete ordinates with a discontinuous Galerkin finite element method (DGFEM) in space, enabling accurate resolution of transport physics on arbitrary polyhedral cells, included locally refined spatial grids. It includes multiple angular quadrature types, including locally refined angular quadratures. Written in modern C++ with a Python API, OpenSn runs efficiently on platforms ranging from laptops to supercomputers. The transport sweep algorithm is implemented using a task-based, directed-acyclic-graph (DAG) approach for each angle and supports asynchronous parallelism across thousands of MPI ranks. Group-set aggregation improves compute intensity, and synthetic acceleration techniques (e.g., diffusion synthetic acceleration, second-moment method) enhance solver convergence. OpenSn has been verified on reactor physics problems and demonstrated excellent weak and strong scaling performance on more than 32,768 processes, making it a versatile and robust platform for large-scale transport simulations in complex geometries.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Bridging adsorption behavior of confined CH 4 -CO 2 binary mixtures across scales
An accurate understanding of the competitive adsorption of CH 4 -CO 2 binary mixtures in nano-confined systems is critical for engineering CO 2 storage in shale gas reservoirs. Due to difficulties in making reliable experimental observations in nano-scale, atomistic simulations (ASs), such as the Grand Canonical Monte Carlo (GCMC) method, provide a viable approach to studying the adsorption behavior of confined fluids. ASs are, however, limited in the size of the compositional domain due to the high computational cost. This work proposes a framework that combines AS and the lattice Boltzmann (LB) method to bridge the physics of confined fluids across scales. The Peng–Robinson equation of state (PR-EoS) produces fugacity coefficients, which serve as input for conducting multi-component GCMC simulations. These GCMC simulations explore the competitive adsorption behavior of CH 4 -CO 2 in nano-slits at various composition, pressure, and channel-width conditions. Both components generate adsorption layers with high densities near the walls with CO 2 preferentially adsorbing compared to CH 4 on the organic walls of carbon sheets. At the mesoscale, a pseudopotential model represents the intermolecular forces in multi-component, multiple-relaxation-time LB simulations. The LB simulations are in good agreement with the GCMC results, allowing us to obtain values for tunable LB parameters. We then extend the use of LB to simulate adsorption behavior in complex networks with nano-sized channels. The phase behavior and fluid properties in the complex geometries of nano-channels differ from nano-slits and bulk systems. Furthermore, the bridging of physics from GCMC (microscale) to LB (mesoscale) via the macroscale PR-EoS connects the adsorption behavior of binary systems across scales.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
3D high-fidelity automated neutronics guided optimization of fusion blanket designs
The compact Fusion Pilot Plant (FPP) is defined in the recent National Academies of Sciences, Engineering, and Medicine report as the next step of fusion energy demonstration with a $50$ MWe peak net electricity production, $Q_e$ greater than $1$, and at least $3$ hours of continuous operation. This fusion pilot plant will be a test bed enabling materials, designs, and fuel management assessment, and it will represent an engineering challenge because of its high-fusion power and compact design targets. Previous reactor data is limited to experiments operating in different design space ranges. Therefore, design iterations and assessments should rely on high-fidelity first-principle theoretical and computational models. The high-fidelity integrated modeling of the plasma is a fundamental part of fusion energy research. However, the whole device modeling is often neglected, utilizing low-fidelity, system-level analysis. Recently, the need for high-fidelity multi-physics modeling was recognized, resulting in a selection of integrated tools. Further, autonomous design optimization requires a streamlined framework that perturbs the design point, reruns the analysis, and examines the outputs. However, high-fidelity analysis requires complex geometry specification that is difficult to perturb. This work presents the parametric CAD generation tool TRACER and a new neutronic workflow. TRACER allows the perturbation of the geometry representation, creating geometry files ready for further analysis. The streamlined neutronic workflow allows efficient and accurate calculations. The two new tools coupled together were used to perform a 3D high-fidelity multi-objective, multi-input optimization of an "ARC Class" compact tokamak design. The workflow was driven by an optimization driver for full automation.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
First operation and validation of simulations for the divertor cryo-vacuum pump in Wendelstein 7-X
Ten cryo-vacuum pumps (CVPs) were installed in the subdivertor region of each island divertor in the stellarator Wendelstein 7-X (W7-X) and operated for the first time during the recently completed plasma campaign OP2.1. A pumping speed of 70 ± 1 $\frac{m^3}{s}$ was measured during dedicated tests with known hydrogen gas injection. Based on a conductance model, the estimated pumping speed ranges from 86-93 $\frac{m^3}{s}$ for different sticking coefficients between 0.6 and 0.8. After completion of the initial tests the CVPs were operated successfully throughout the campaign, with regeneration performed once a week. Neutral gas pressures in the subdivertor in the range of 10 −4 mbar are well within the molecular flow regime and limit the particle exhaust capabilities of the CVPs. Simulations of the neutral gas pressure in the three-dimensional complex geometry of the subdivertor were performed using the DIVGAS code based on the direct simulation Monte Carlo method and a model implemented in the steady-state thermal package in ANSYS, which are in agreement with the measured values during plasma operation.
Cryo-vacuum pumping↗
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↗
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↗
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↗
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↗
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↗
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↗
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↗
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↗
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↗
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↗