Engineering Papers⌕ Search

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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 955 records · Page 53

Analysis and improvement of the hot disk transient plane source method for low thermal conductivity materials

The hot disk transient plane source (TPS) method is a widely used standard technique (ISO 22007-2) for the characterization of thermal properties of materials, especially the thermal conductivity, k. Despite its well-established reliability for a wide variety of common materials, the hot disk TPS method is also known to suffer from a substantial systematic errors when applied to low-k thermal insulation materials. Here, we present a combined numerical and experimental study on the influence of the geometry of hot disk sensor on measured value of low-k materials. We demonstrate that the error is strongly affected by the finite thickness and thermal mass of the sensor's insulation layer was well as the corresponding increase of the effective heater size beyond the radius of the embedded metal heater itself. We also numerically investigate the dependence of the error on the sample thermal properties, confirming that the errors are worse in low-k samples. A simple correction function is also provided, which converts the apparent (erroneous) result from a standard hot disk TPS measurement to a more accurate value. A standard polyimide sensor was also optimized using both wet and dry etching to provide more accurate measurement directly. Experimentally corrected value of k for Airloy x56 aerogel and a commercial silica aerogel using the numerical correction factor derived based on the standard TPS sensor is in excellent agreement with the directly measured value from the TPS sensor using the optimized polyimide sensor. Finally, both of these methods can reduce the errors to less than 4% as compared to around 40% error of overestimation from raw values measured with the pristine sensor. Such results show that both the numerical correction to a pristine senor or an optimized sensor are capable of providing highly accurate value of thermal conductivity for such materials.

36 MATERIALS SCIENCE↗

Multimodal sensor fusion for real-time standoff estimation in directed energy deposition

In Laser Powder-based Direct Energy Deposition (LP-DED) systems, achieving consistency, precision and quality of produced parts requires tight control over printing parameters. One of the critical parameters is the standoff distance. Maintaining an optimal standoff height is crucial for achieving correct laser power density and powder catchment efficiency, as both laser and powder streams are focused at this distance. Here, this study introduces a novel approach using multimodal sensor fusion to predict standoff height in real-time. The proposed system integrates two low-profile, cost-effective sensors: an RGB coaxial camera and a high frequency and high dynamic range microphone. By utilizing a simple fully connected neural network, trained on a limited dataset, data fusion of these sensors allowed for the real-time prediction of the standoff height. The results demonstrate high resolution and accuracy of the predictions across multiple geometries and a wide range of standoff heights. This approach offers a simple, and cost-effective solution for real-time standoff height monitoring and lays the groundwork for future integration into commercial LP-DED systems.

42 ENGINEERING↗

Exact Gaussian processes for massive datasets via non-stationary sparsity-discovering kernels

Abstract A Gaussian Process (GP) is a prominent mathematical framework for stochastic function approximation in science and engineering applications. Its success is largely attributed to the GP’s analytical tractability, robustness, and natural inclusion of uncertainty quantification. Unfortunately, the use of exact GPs is prohibitively expensive for large datasets due to their unfavorable numerical complexity of $$O(N^3)$$ O ( N 3 ) in computation and $$O(N^2)$$ O ( N 2 ) in storage. All existing methods addressing this issue utilize some form of approximation—usually considering subsets of the full dataset or finding representative pseudo-points that render the covariance matrix well-structured and sparse. These approximate methods can lead to inaccuracies in function approximations and often limit the user’s flexibility in designing expressive kernels. Instead of inducing sparsity via data-point geometry and structure, we propose to take advantage of naturally-occurring sparsity by allowing the kernel to discover—instead of induce—sparse structure. The premise of this paper is that the data sets and physical processes modeled by GPs often exhibit natural or implicit sparsities, but commonly-used kernels do not allow us to exploit such sparsity. The core concept of exact, and at the same time sparse GPs relies on kernel definitions that provide enough flexibility to learn and encode not only non-zero but also zero covariances. This principle of ultra-flexible, compactly-supported, and non-stationary kernels, combined with HPC and constrained optimization, lets us scale exact GPs well beyond 5 million data points.

97 MATHEMATICS AND COMPUTING↗

NSRD-20, Functional Testing of Novel MTC HEPA Filtration Media

To improve the safety of Department of Energy (DOE) nuclear facilities, this project aimed to develop ceramic high efficiency particulate air (HEPA) filters that also create a low pressure-drop (dP) in use. The filters under consideration use mini-tubular ceramic (MTC) HEPA filtration media made with different architectures. Once successful, this project will reduce life-cycle costs, including safety basis, operational, and waste-disposal costs. The technology’s mini-tubular geometry is known to reduce dP when compared to flow through membranes of equivalent mass and surface area. Low dP facilitates retrofitting advanced filters into existing DOE facilities to make them safer while simultaneously reducing operational costs. Operations of ventilation systems are a key cost driver for DOE nuclear facilities; reducing dP can significantly reduce DOE’s operational costs. However, the filtration efficiency (FE) of this type of filtration media depends on the architecture of the filter media in an unknown way and requires testing. Functional testing of filtration media with different architectures will establish the performance of the filtration media and will enable development of guidance on design requirements to optimize filter performance.

36 MATERIALS SCIENCE↗

Material Characterization for Large Scale Additive Manufacturing (AM)

The objective of this project was to enable prediction of residual stresses and warpage for extrusion deposition additive manufacturing (EDAM) fabricated carbon fiber thermoplastic matrix parts with three different materials systems produced by Techmer. The materials selected include Polyphenylene sulfide reinforced with 50% by weight of carbon fiber (CF-PPS), Polyethersulfone reinforced with 25% by weight of carbon fiber (CF-PESU), Polysulfone reinforced with 25% by weight of carbon fiber (CF-PSU). ADDITIVE3D © , a physics-based simulation workflow for EDAM, provided the simulation capabilities required for this project. Simulation predictions were validated against measurements carried out during and after the printing process carried out in the CAMRI and LSAM systems . Predictions for temperature, degree of crystallinity and deformation were carried out for the three material systems and two different geometries. Predictions for temperature were correlated very well with the experimental measurements for the two geometries printed using the three-material systems. The crystallinity level was verified for CF-PPS. The predictions for part deformation were in good agreement with the experimental measurements. In the best-case, predictions were within 8% of the maximum displacement observed in the 3-direction whereas predictions were within 14% for the worst-case. The adoption of this technology in commercial large-scale EDAM production processes would dramatically reduce the costs associated with producing articles by that process. Many thousands of dollars in materials, energy, and machine time could be saved by utilizing this simulation technology to develop articles rather than iterative printings to arrive at the optimal or correct design. The avoidance of a single failed print has the possibility of saving tens of thousands of dollars involved in the cost of material, machine and operators’ time.

36 MATERIALS SCIENCE↗

Methane Adsorption on Heteroatom-Modified Maquettes of Porous Carbon Surfaces

Experimental and theoretical studies disagree on the energetics of methane adsorption on carbon materials. However, this information is critical for the rational design and optimization of the structure and composition of adsorbents for natural gas storage. The delicate nature of dispersion interactions, polarization of both the adsorbent and the adsorbate, interplay between H-bonding and tetrel bonding, and induced dipole/Coulomb interactions inherent to methane physisorption require computational treatment at the highest possible level of theory. In this study, we employed the smallest reasonable computational model, a maquette of porous carbon surfaces with a central site for substitution and methane binding. The most accurate predictions of methane adsorption energetics were achieved by electron-correlated molecular orbital theory CCSD(T) and hybrid density functional theory MN15 calculations employing a saturated, all-electron basis set. Here, the characteristic geometry of methane adsorption on a carbon surface (“lander approach”) arises due to bonding interactions of the adsorbent π-system with the proximal H–C bonds of methane, in addition to tetrel bonding between the antibonding orbital of the distal C–H bond and the central atom of the maquette (C, B, or N). The polarization of the electron density, structural deformations, and the comprehensive energetic analysis clearly indicate a ∼3 kJ mol –1 preference for methane binding on the N-substituted maquette. The B-substituted maquette showed a comparable or lower binding energy than the unsubstituted, pure C model, depending on the level of theory employed. The calculated thermodynamic results indicate a strategy for incorporating electron-enriched substitutions (e.g., N) into carbon materials as a way to increase methane storage capacity over electron-deficient (e.g., B) modifications. The thermochemical analysis was revised for establishing a conceptual agreement between the experimental isosteric heat of adsorption and the binding enthalpies from statistical thermodynamics principles.

03 NATURAL GAS↗

Effect of sintering temperature on feature resolution and flexural strength of ceramics fabricated through vat photopolymerization additive manufacturing

Although ceramic additive manufacturing (AM) could be used to fabricate complex, high-resolution parts for diverse, functional applications, one ongoing challenge is optimizing the post-process, particularly sintering, conditions to consistently produce geometrically accurate and mechanically robust parts. This study aims to investigate how sintering temperature affects feature resolution and flexural properties of silica-based parts formed by vat photopolymerization (VPP) AM. Test artifacts were designed to evaluate features of different sizes, shapes and orientations, and three-point bend specimens printed in multiple orientations were used to evaluate mechanical properties. Sintering temperatures were varied between 1000°C and 1300°C. Deviations from designed dimensions often increased with higher sintering temperatures and/or larger features. Higher sintering temperatures yielded parts with higher strength and lower strain at break. Many features exhibited defects, often dependent on geometry and sintering temperature, highlighting the need for further analysis of debinding and sintering parameters. To the best of the authors’ knowledge, this is the first time test artifacts have been designed for ceramic VPP. This work also offers insights into the effect of sintering temperature and print orientation on flexural properties. These results provide design guidelines for a particular material, while the methodology outlined for assessing feature resolution and flexural strength is broadly applicable to other ceramics, enabling more predictable part performance when considering the future design and manufacture of complex ceramic parts.

36 MATERIALS SCIENCE↗

Breakup dynamics in a pressure-swirl injector for urea-water solution applications: A computational study

The co-optimization of in-cylinder combustion and after-treatment technology has become a major aspect in engine design and development, with the goal of meeting the increasingly restrictive emission regulations in the transportation industry. Selective Catalytic Reduction is a robust technology to control the emission of NO x , and the injection of urea in water solution is the exhaust tailpipe is a key aspect of its operation. The proposed work uses high-fidelity Computational Fluid Dynamics to characterize the atomization dynamics of the liquid jet in relevant cross-flow conditions. The study focuses on a commercial low-pressure (9 bar) pressure-swirl injector which is characterized in its internal geometry through high-resolution X-ray micro-computational tomography. The internal two-phase flow has been modeled according to the volume-of-fluid approach in a large eddy simulation framework and validated against near-nozzle X-ray radiography measurement. Moreover, characterizing the breakup dynamics for the swirling hollow cone formation, and assessing the influence of the cross-flow in the breakup dynamics was completed. The results have been reported proposing Re-Oh maps and probability density functions of the spray kinematics. Higher cross-flow momentum generates an increase in the jet intact length and a reduction of the liquid droplet diameters. The axial momentum of the jet is affected by the cross-flow already in the near-nozzle region, determining a relevant deviation of the spray velocities. In conclusion, this work aims to inform the initialization of Eulerian-Lagrangian spray models through the assignment of droplet kinematics and static one-way coupling between volume-of-fluid results and Lagrangian spray parcels, to be used for system-size domain simulations.

33 ADVANCED PROPULSION SYSTEMS↗

MADE3D: Enabling the Next-Generation High-Torque- Density Wind Generators by Additive Design and 3D Printing

Direct-drive wind turbine generators are increasing in popularity, thanks to recent project developments - especially offshore, where reliability and efficiency are major cost drivers. Yet, high capital costs are forcing many original equipment manufacturers to consider lightweight, high-torque density generators for next-generation multi-megawatt turbines that may be difficult to realize by traditional design or manufacturing methods. In this study, we present a new design framework enabled by advanced machine learning and multimaterial additive manufacturing to perform a magnetic topology optimization that maximizes the torque per rotor active mass for a 15-megawatt direct-drive permanent magnet wind generator. A comparison of the proposed approach against conventional topology optimization demonstrated a significant increase in computational efficiency and accuracy in performance predictions. Results using single and multimaterial compositions for rotor core and magnets identify a wider choice of 3D printable designs for a given specification. A hybrid combination of sintered and dysprosium-free polymer-bonded magnets shows good potential for torque performance by saving material costs up to 8.75%. More than 30% improvement in rotor torque densities is identified which can marginally improve the overall generator torque density. With the rapid evolution of multipowder deposition technologies, this study can greatly inspire a new paradigm for design-driven manufacturing with novel material compositions and lightweight, low-cost, high-strength multimaterial geometries that were previously unexplored for direct-drive generators.

3D printing↗

Programming approaches for scalability, performance, and portability of combustion physics codes

Here, this paper presents the process, strategy, and results associated with porting a typical combustion physics flow solver to current state-of-the-art and future massively-parallel computer architectures. Major focus is placed on the distinct algorithmic structure of these types of codes and how it can be integrated with modern programming paradigms for heterogeneous platforms (i.e., distributed many-core systems with accelerators). An end-to-end case study is presented that exemplifies the process in a generic manner, which then serves as a clear guide with respect to the strategy and best practices leading to a robust and adaptable framework that performs well, is durable over time, is portable, and requires minimal human-effort. This end is accomplished beginning with the use of a mature, validated, structured, multiblock code framework optimized for application of both Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS). This code has been ported to a variety of platforms over the past decade, including most recently the Oak Ridge Leadership Computing Facility’s “Summit” Platform. The experience gained on these multiple platforms provides general insights and thus the results presented are not specific to any one code or platform other than the overarching trend toward distributed many-core systems with accelerators in order to move toward exascale performance. The resultant performance and scalability of the ported code is demonstrated on a real-world application; a state-of-the-art rotating detonation rocket engine simulation that matches the complex geometry and boundary conditions imposed as part of a companion experimental campaign.

97 MATHEMATICS AND COMPUTING↗

Internal 2-D Surface Temperature Measurements for Large Complex Geometries

For 2D-temperature monitoring applications, a variant of EIT (Electrical Impedance Tomography) is evaluated computationally in this work. Literature examples of poor sensor performance in the center of the 2D domains away from the side electrodes motivated this study which seeks to overcome some of the previously noted shortcomings. In particular, the use of ‘sensing skins’ with novel tailored baseline conductivities were examined using the EIDORS package for EIT. It was found that the best approach for detecting a hot spot depends on several factors such as the current injection (stimulation) patterns, the measurement patterns, and the reconstruction algorithms. For a well-performing combination of these factors, tailored baseline conductivities were assessed and compared to the baseline uniform conductivity. It was discovered that for some EIT applications, a tailored distribution needs to be smooth and that sudden changes in the conductivity gradients should be avoided. Still, the benefits in terms of improved EIT performance were small for conditions for which the EIT measurements had been ‘optimized’ for the uniform baseline case. Within the limited scope of this study, only two specific cases showed benefits from tailored distributions. For one case, a smooth tailored distribution with increased baseline conductivity in the center provided a better separation of two centrally located hot spots. For another case, a smooth tailored distribution with reduced conductivity in the center provided better estimates of the magnitudes of two hot spots near the center of the sensing skin.

47 OTHER INSTRUMENTATION↗

Understanding extraction limits of plasma cathodes with experiment and simulation

The project is focused on computations for enhanced ionization near the exit orifice of a plasma cathode and how the cathode extraction electrode geometry and stray magnetic field structure influence available extractable current. The computational effort will take place at the Princeton Collaborative Research Laboratory (PCRF). The specific modeling platform to be utilized from PCRF is the electrostatic direct implicit particle in cell (EDIPIC) code. EDIPIC will be used to study ionization processes in the electron extraction region of the cathode. The control of electron energy is a key aspect to minimizing power losses in the plasma. Exiting electrons ionize gas leaving the orifice and in the process produce the dense plasma from which large amounts of current can be extracted. Project activities include studying processes that determine the energy distribution of the exiting electrons thus providing insight into how to optimize the source. The model will be validated and developed using experimental measurements as warranted. The familiarity and expertise of the team at PPPL with this type of plasma and the challenges posed greatly enhances the likelihood of gaining insight into the extraction sheath under the duration of the proposed activity. The effort is expected to shed a great deal of insight into electron extraction processes and pathways for optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Geometric concepts for stellarator permanent magnet arrays

The development of stellarators that use permanent magnet arrays to shape their confining magnetic fields has been a topic of recent interest, but the requirements for how such magnets must be shaped, manufactured, and assembled remain to be determined. To address these open questions, we have performed a study of geometric concepts for magnet arrays with the aid of the newly developed MAGPIE code. A proposed experiment similar to the National Compact Stellarator Experiment (NCSX) is used as a test case. Two classes of magnet geometry are explored: curved bricks that conform to a regular grid in cylindrical coordinates, and hexahedra that conform to the toroidal plasma geometry. In addition, we test constraints on the magnet polarization. While magnet configurations constrained to be polarized normally to a toroidal surface around the plasma are unable to meet the required magnetic field parameters when subject to physical limitations on the strength of present-day magnets, configurations with unconstrained polarizations are shown to satisfy the physics requirements for a targeted plasma.

stellarator, permanent magnet, optimization↗

Application of Koopman operator for model-based control of fracture propagation and proppant transport in hydraulic fracturing operation

This work explores the application of the recently developed Koopman operator approach for model identification and feedback control of a hydraulic fracturing process. Controlling fracture propagation and proppant transport with precision is a challenge due in large part to the difficulty of constructing approximate models that accurately capture the characteristic moving boundary and highly-coupled dynamics exhibited by the process. Koopman operator theory is particularly attractive here as it offers a way to explicitly construct linear representations for even highly nonlinear dynamics. The method is data-driven and relies on lifting the states to an infinite-dimensional space of functions called observables where the dynamics are governed by a linear Koopman operator. Here this work considers two problems: (a) fracture geometry control, and (b) proppant concentration control. In both cases, an approximate linear model of the corresponding dynamics is constructed and used to design a model predictive controller (MPC). The manuscript shows that in the case of highly nonlinear dynamics, as observed in the proppant concentration, use of canonical functions in the observable basis fails. In such cases, a priori system knowledge can be leveraged to choose the required basis. The numerical experiments demonstrate that the Koopman linear model shows excellent agreement with the real system and successfully achieves the desired target values maximizing the oil and gas productivity. Additionally, due to its linear structure, the Koopman models allow convex MPC formulations that avoid any issues associated with nonlinear optimization.

42 ENGINEERING↗

Testing and Expertise for Marine Energy (TEAMER) Program Support - Numerical Modeling Assistance for iProTech's "PIP" WEC Device: Cooperative Research and Development Final Report, CRADA Number CRD-20-17303

Develop a time-domain, 3 degree of freedom numerical model of IProTech's PIP device, including power take-off (PTO), that can be used to develop control algorithms and establish a baseline technology performance level (TPL). This is important to IProTech because having a good numerical model allows IProTech to optimize design and control algorithms on an ongoing basis. A baseline TPL enables the PIP technology to be benchmarked against other wave energy converter (WEC) concepts, providing (i) justification (or not) for investment in wave tank model testing, (ii) identifying optimum design parameters and control algorithms for the test model, (iii) specifying test conditions and (iv) interpreting wave tank test results. An accredited WEC-Sim model of the PIP device is an essential tool for further development of the technology. 1) The development of numerical models of the PIP concept will enable iProTech to evaluate the potential performance of this concept, and to investigate how sensitive the device is to different variables within the system (e.g. hydraulic components, geometry dimensions, etc.) 2) The numerical models were developed with WEC-Sim and PTO-Sim -- the de-facto industry standard WEC numerical modeling tools. A detailed model of the system's hydraulic circuit was developed with PTO-Sim. However, WEC-Sim is based on linear hydrodynamics and has several limitations; validation against high fidelity models/physical scale models would help to build confidence in the results.

16 TIDAL AND WAVE POWER↗

Heating, DPA, and He production in the Proton Beam Window Assembly for the Second Target Station

This report details the analysis of the heating rates, dpa values, and He production in the Proton Beam Window (PBW), the assembly around the PBW, and the shielding above it. Neutronics calculations use a detailed Unstructured Mesh (UM) converted from the solid CAD model of the PBW design of June 2022. The rest of the geometry is taken from the CSG master model (as of June 2022). For this analysis, a proton beam profile with a footprint of ≈90 cm 2 has been used. This proton profile was created by a Ring to Second Target magnet array configuration using octupole magnets. Maximum values in the PBW (aluminum) are 5.4 J/cc/pulse heating rate, 1.1 dpa per year, and 660 appm He-production per year. These results within ~ 10% with the previous calculations from June 2020. Average heating and dpa rates have been calculated for each component in the PBW assembly and will be used to evaluate and design the cooling needs and make decisions on removable components for the PBW. It is desirable to optimize and simplify the PBW design by only replacing components that are at the end of their lifetime. Additional analysis for the energy deposition in the PBW is performed for two super-Gaussian beam profiles with a footprint of 60 cm 2 and 30 cm 2 respectively. Also, an estimate of the PBW lifetime is performed for two realistic quadrupole beam profiles.

43 PARTICLE ACCELERATORS↗

An implementation of a high-order generalized finite difference method for solving the time-harmonic cold plasma wave equation in toroidal geometry

A high-order physics-informed meshless finite difference numerical technique is introduced for solving the time-harmonic cold plasma wave equation in toroidal geometries, presenting a novel application of the generalized finite difference (GFD) method to plasma wave simulations. The algorithm employs an irregular distribution of computational points, with local point density informed by the shortest wavelength derived from the cold plasma dispersion relation. Numerical stability and robustness are addressed using regularization techniques. The algorithm, implemented for two spatial dimensions, solves for the wave electric field and is demonstrated to achieve convergence rates of $\mathcal{O}$($\mathcal{h}$ $\mathcal{P}$ )⁠. Verification tests reproduce plane wave solutions, and example simulations of ion cyclotron resonance heating and electron cyclotron resonance heating demonstrate its capability, approaching realistic tokamak plasma scenarios. This work contributes to laying a foundation for the GFD method to be used in more sophisticated, optimized, and physically realistic full-wave simulations in time-harmonic plasma wave research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Compatibility of molten plutonium with wrought and additively manufactured metal crucibles

Understanding plutonium’s interaction with metals is crucial for optimizing pyrochemical operations, nuclear fuel containment, and various actinide processing techniques. Traditionally, tantalum crucibles are employed for plutonium processing due to their high durability, excellent temperature stability, and low solubility in plutonium. However, tantalum faces challenges such as plutonium wetting and diffusion, making surface coatings particularly important for crucibles in pyrochemical applications to enhance corrosion resistance against plutonium. Tantalum is also expensive and difficult to machine, prompting the need for advanced manufacturing techniques to address these challenges. Here, in this work, we investigate the interaction of Pu with tantalum and titanium crucibles fabricated using both traditional machining methods and laser powder bed fusion (LPBF) additive manufacturing (AM). LPBF-AM is an advanced technique that allows for the creation of complex geometries from traditionally difficult-to-machine metals by using a high-powered laser to build parts. Previous studies of conventional manufactured tantalum have utilized oxidation and carburization of the surface to mitigate plutonium wetting; however, no studies of surface modified LPBF-AM material have been undertaken. These studies are crucial, given the typical differences in the grain structure between conventional and LPBF-AM materials. All crucibles underwent differential scanning calorimetry to confirm the melting of plutonium. Subsequently, the crucibles were sectioned and mounted in epoxy for microstructural analysis using optical microscopy and scanning electron microscopy. This investigation, comparing the performance of wrought vs AM metal crucibles, provides a basis for future tooling applications in actinide processing techniques and can address the challenges associated with traditional machining, particularly in pyrochemical applications.

Actinides↗