Search NASA⌕ Search

SEARCH · Search NASA

Results for “Experimental and numerical techniques”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Effect of Computational Schemes on Coupled Flow and Geo-Mechanical Modeling of CO 2 Leakage through a Compromised Well

Carbon capture, utilization, and storage (CCUS) describes a set of technically viable processes to separate carbon dioxide (CO 2 ) from industrial byproduct streams and inject it into deep geologic formations for long-term storage. Legacy wells located within the spatial domain of new injection and production activities represent potential pathways for fluids (i.e., CO 2 and aqueous phase) to leak through compromised components (e.g., through fractures or micro-annulus pathways). The finite element (FE) method is a well-established numerical approach to simulate the coupling between multi-phase fluid flow and solid phase deformation interactions that occur in a compromised well system. We assumed the spatial domain consists of a three-phases system: a solid, liquid, and gas phase. For flow in the two fluids phases, we considered two sets of primary variables: the first considering capillary pressure and gas pressure (PP) scheme, and the second considering liquid pressure and gas saturation (PS) scheme. Fluid phases were coupled with the solid phase using the full coupling (i.e., monolithic coupling) and iterative coupling (i.e., sequential coupling) approaches. The challenge of achieving numerical stability in the coupled formulation in heterogeneous media was addressed using the mass lumping and the upwinding techniques. Numerical results were compared with three benchmark problems to assess the performance of coupled FE solutions: 1D Terzaghi’s consolidation, Liakopoulos experiments, and the Kueper and Frind experiments. We found good agreement between our results and the three benchmark problems. For the Kueper and Frind test, the PP scheme successfully captured the observed experimental response of the non-aqueous phase infiltration, in contrast to the PS scheme. These exercises demonstrate the importance of fluid phase primary variable selection for heterogeneous porous media. We then applied the developed model to the hypothetical case of leakage along a compromised well representing a heterogeneous media. Considering the mass lumping and the upwinding techniques, both the monotonic and the sequential coupling provided identical results, but mass lumping was needed to avoid numerical instabilities in the sequential coupling. Additionally, in the monolithic coupling, the magnitude of primary variables in the coupled solution without mass lumping and the upwinding is higher, which is essential for the risk-based analyses.

deformation flow↗
DOE OSTI2025

Audible enclaves crafted by nonlinear self-bending ultrasonic beams

Delivering audible content to a targeted listener without disturbing others is paramount in audio engineering. However, achieving this goal has long been challenging due to the diffraction of low-frequency (long-wavelength) audio waves in linear acoustics. Here, we introduce an approach for creating remote audio spots, dubbed audible enclaves, by harnessing the local nonlinear interaction of two self-bending ultrasonic beams with distinct spectra. The self-bending ultrasonic beams created by acoustic metasurfaces, though inaudible, can bypass obstacles such as human heads. At their intersection behind obstacles, highly localized audible enclaves are formed due to the local nonlinear interactions. Additionally, we demonstrate the ultrabroadband capabilities of our metasurface-based implementation both numerically and experimentally, spanning from 125 Hz to 4 kHz (6 octave bands), covering the majority of the audible frequency range. The practicality of our proposed technique is underscored by its compact implementation size (0.16 m, equivalent to 0.06 wavelengths at 125 Hz), as well as its robust performance under wideband transient audio signal excitation and in a common room with reverberations. Our proposed audible enclaves hold significant potential for various applications in advanced audio engineering, including private speech communications, immersive spatial audio reproduction, and high-resolution sound/quiet zone control.

Science & Technology - Other Topics↗
DOE OSTI2023

Parameter identification for symmetrical Prandtl-Ishlinskii hysteresis model using Gauss-Newton method

The Prandtl-Ishlinskii (PI) hysteresis model has been used to describe nonlinear hysteretic behaviors. In this paper, the PI model is first simplified to accommodate the symmetrical behavior of ferromagnetic hysteresis. A parameter identification method is then developed to extract the parameters of the symmetrical PI model from experimental data. Based on the Gauss-Newton method, regularization technique, and line search method, the model parameters can be retrieved efficiently from the measured data. Numerical examples are presented to demonstrate the performance of the proposed method.

Yan, Su↗
DOE OSTI2021

Techniques for studying materials under extreme states of high energy density compression

The properties of materials under extreme conditions of pressure and density are of key interest to a number of fields, including planetary geophysics, materials science, and inertial confinement fusion. In geophysics, the equations of state of planetary materials, such as hydrogen and iron, under ultrahigh pressure and density provide a better understanding of their formation and interior structure [Celliers et al., “Insulator-metal transition in dense fluid deuterium,” Science 361, 677–682 (2018) and Smith et al., “Equation of state of iron under core conditions of large rocky exoplanets,” Nat. Astron. 2, 591–682 (2018)]. The processes of interest in these fields occur under conditions of high pressure (100 GPa–100 TPa), high temperature (>3000 K), and sometimes at high strain rates (>103 s−1) depending on the process. With the advent of high energy density (HED) facilities, such as the National Ignition Facility (NIF), Linear Coherent Light Source, Omega Laser Facility, and Z, these conditions are reachable and numerous experimental platforms have been developed. To measure compression under ultrahigh pressure, stepped targets are ramp-compressed and the sound velocity, measured by the velocity interferometer system for any reflector diagnostic technique, from which the stress-density of relevant materials is deduced at pulsed power [M. D. Knudson and M. P. Desjarlais, “High-precision shock wave measurements of deuterium: Evaluation of exchange-correlation functionals at the molecular-to-atomic transition,” Phys. Rev. Lett. 118, 035501 (2017)] and laser [Smith et al., “Equation of state of iron under core conditions of large rocky exoplanets,” Nat. Astron. 2, 591–682 (2018)] facilities. To measure strength under high pressure and strain rates, experimenters measure the growth of Rayleigh–Taylor instabilities using face-on radiography [Park et al., “Grain-size-independent plastic flow at ultrahigh pressures and strain rates,” Phys. Rev. Lett. 114, 065502 (2015)]. The crystal structure of materials under high compression is measured by dynamic x-ray diffraction [Rygg et al., “X-ray diffraction at the national ignition facility,” Rev. Sci. Instrum. 91, 043902 (2020) and McBride et al., “Phase transition lowering in dynamically compressed silicon,” Nat. Phys. 15, 89–94 (2019)]. Medium range material temperatures (a few thousand degrees) can be measured by extended x-ray absorption fine structure techniques, Yaakobi et al., “Extended x-ray absorption fine structure measurements of laser-shocked V and Ti and crystal phase transformation in Ti,” Phys. Rev. Lett. 92, 095504 (2004) and Ping et al., “Solid iron compressed up to 560 GPa,” Phys. Rev. Lett. 111, 065501 (2013), whereas more extreme temperatures are measured using x-ray Thomson scattering or pyrometry. This manuscript will review the scientific motivations, experimental techniques, and the regimes that can be probed for the study of materials under extreme HED conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Time-lapse approach to correct deficiencies of 2D soil zymography

Membrane zymography is commonly used in rhizosphere ecology to map enzyme activities in intact soil samples and plant roots. The method consists of incubating a membrane saturated with an enzyme-specific fluorogenic substrate on the soil/root surface followed by measurements of fluorescence intensity of the product in the membrane. The traditional zymography is based on assumptions that fluorescence on membrane images is linearly increasing with time during zymography and an increase in product content is numerically equal to enzyme activity in soil below the membrane. These assumptions are unlikely to hold in experimental settings. In this work we introduce a new zymography technique, time-lapse zymography (TLZ); the approach that eliminates the need for assumptions of the traditional zymography and provides more realistic estimates of enzymatic activities. We assessed the new technique in a series of laboratory and modeling experiments, including quantification of the fluorescent product diffusion (e.g. MUF: 4-methylumbelliferone) from plant roots, enzyme activity measurements on the roots, and HYDRUS-2D & HP2 software calibration with obtained data. The calibrated model was used to analyze the processes governing spatial and temporal dynamics of MUF contents in the membrane. The results indicated that the enzyme diffusion within the membrane-soil system was negligible, and measured zymograms were adequately reproduced solely by accounting for substrate and product diffusions and for catalytic enzyme reaction described by the Michaelis-Menten equation. TLZ enabled identifying and using linear parts on MUF time series and accounting for MUF losses in each zymogram pixel, considerably improving the accuracy as compared to the traditional method. Results demonstrated that enzymatic activity from only a thin soil layer (~0.2 mm) is adequately represented in zymograms.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Full-Field Strain Measurement Integrated with Two Dimension Regression Analysis to Evaluate the Bi-Modulus Elastic Properties of Isotropic and Transversely Isotropic Materials

Background: Measuring the physical properties of shale is critical for optimizing engineering activities such as geothermal energy generation and hydraulic fracturing. Shale is a transversely isotropic material. Furthermore, this material can also include micro and macro cracks at different locations and orientations that cause it to behave differently under tensile or compressive loading. Objective: In this work, a combined experimental–numerical approach is proposed to evaluate the bi-modulus elastic properties of isotropic and transversely isotropic materials. Methods: Full-field strain measurements for a circular disk under diametral compression are integrated with a regression analysis technique to evaluate the elastic properties of bi-modulus materials subjected to tensile and compressive loads using two loading configurations on the same specimen. Digital Image Correlation (DIC) is used to measure the full-field strains. Subsequently, in the case of an isotropic material, a linear least-squares approach is utilized to process the experimentally determined strains in conjunction with analytical expressions of the stress fields (in terms of far-field loading) to determine the elastic modulus E, the shear modulus G, and the Poisson’s ratio $v$. In the case of a transversely isotopic material, such as shale, a finite element model is implemented to determine the stress fields (again in terms of far-field loading), which is followed by repeating the previous regression analysis in an iterative process to estimate the elastic parameters. Results: The results show that the proposed technique successfully provides a complete set of elastic properties as a function of both the loading condition and the principal material directions. The technique is validated by measurements on a known isotropic material and then applied to determine the properties of shale. Conclusion: In this work, the proposed approach is successfully used to calculate the bi-modulus elastic response of poly(methyl meth- acrylate) (PMMA) and shale. As expected, PMMA exhibits an isotropic response with no bi-modulus effect, however, shale exhibits both transverse isotropy and a bi-modulus effect. Therefore, this approach holds promise for investigating the elastic properties of materials like rocks and fiber-reinforced composite laminates as functions of the principal material directions and the loading conditions.

42 ENGINEERING↗
DOE OSTI2024

Stress interference in multilayer additive friction stir deposition of AA6061 aluminum

Due to the multilayer deposition nature of metal additive manufacturing processes, each layer being printed experiences the state of thermokinetic and thermomechanical stress that in turn interfere with the state of thermokinetics and thermomechanical stress of subsequently deposited layers. Especially, this multilayer interference significantly affects the resultant properties of the component fabricated using solid-state additive friction stir deposition due to evolution of asymmetric state of planar stress. Due to the lack of comprehensive and suitable in situ diagnosis technique, the complex interference of inter- and multi-layer stresses during additive friction stir deposition was studied in an integrated approach of numerical simulation of fluidic state and experimental probing of stress influenced ultrasonic elastography. The uni-directional and bi-directional layer deposition configurations adopted during additive friction stir deposition result in the generation of constructive and destructive interference of the interlayer stress and hence, asymmetric and symmetric dynamic elasticity distribution respectively within the subsequent layers. With subsequent deposition of additional layers, the odd and even numbers of deposited layers generate asymmetric and nearly symmetric dynamic elasticity distributions.

Yang, Teng↗
DOE OSTI2022

Anion diffusion in two-dimensional halide perovskites

Commercialization of halide perovskites in the semiconductor industry is hindered by their short-term stability. The instability of perovskites is closely interlinked with ionic diffusion. Historically, attempts to study diffusion in 2D perovskites mostly utilized electrical characterizations, but these characterizations pose a challenge in deconvoluting the impact of device architecture, interlayers, and ionic species. In this Perspective, we focus our attention on simple optical characterizations employed in the literature to investigate halide diffusion in 2D perovskites using lateral and vertical heterostructure platforms. We review the various synthesis techniques used for fabrication of halide perovskite heterostructures and discuss the qualitative and quantitative diffusion studies performed using these platforms. We discuss the numerical methods used to validate and supplement the experimental halide diffusion kinetics. Finally, we highlight the need to conduct further research on the impact of device operating conditions, lattice structure, and vacancy concentration on halide diffusion. Through this Perspective, we aim to emphasize the need of developing a comprehensive understanding of halide diffusion in perovskites for their successful deployment in optoelectronics.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Using Hydrodynamic Similarity as a Verification Method for Impact Cratering Simulations in the FLAG Hydrocode

Hydrodynamic codes (hydrocodes) are common tools for modeling hypervelocity impacts to provide insight into the physical phenomenon. Hydrocodes can simulate impacts from micrometer to kilometer spatial scales and reach impact velocities difficult to achieve in experimental settings. However, numerical models are approximations, and demonstrating that a numerical method is capable of providing physical results for these models is essential. In this work, we employ a hydrocode verification technique that leverages hydrodynamic similarity, a mathematical property of the conservation equations of fluid mechanics that form the basis for hydrocode models. Using the FLAG hydrocode, we simulate aluminum (Al) and basalt projectiles and targets at spatial scales spanning 7 orders of magnitude (hundreds of micrometers to kilometers). These materials were chosen because Al-6061 is a common material in spacecraft and satellites and basalt is a useful approximation of rocky astronomical bodies. Our results show that hydrodynamic similarity holds for each material model used and across spatial scales. We show that under certain conditions hydrodynamic similarity can apply in the presence of gravity and that similarity does not hold in the presence of strength models. We conclude that the FLAG hydrocode preserves important mathematical properties of fluid dynamics in hypervelocity impacts of Al-6061 and basalt.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Experimental and Numerical Investigation of a Miniature Additively Manufactured Vortex Tube

Ranque–Hilsch vortex tubes are simple devices that can produce a cooling effect using compressed air. A key advantage of vortex tubes is the lack of moving solid parts; however, their efficiencies are relatively low. The present study focuses on the development of a miniature variable-diameter tube using additive manufacturing. A metal-based 3D printing technique was utilized to fabricate this vortex tube monolithically. Computational fluid dynamics simulations employing software STAR-CCM + with a compressible Reynolds-Averaged Navier–Stokes (RANS) approach and the elliptic-blending lag k-epsilon turbulence model have been applied to model thermofluid processes inside the vortex tube, to good agreement with the experiment. In conclusion, a temperature decrease of 13.3 °C and a cooling power of approximately 4 W were experimentally achieved with a pressure ratio of 4 in the air at normal conditions. This result shows promise for the goal of utilizing additive manufacturing to design and build complex-geometry vortex tubes intended for use with cryogenic fluids.

08 HYDROGEN↗
DOE OSTI2021

Numerical and experimental evaluation of mechanical performance of the multifunctional energy storage composites

This work presents numerical simulation methods to model the mechanical behavior of the multifunctional energy storage composites (MESCs), which consist of a stack of multiple thin battery layers reinforced with through-the-hole polymer rivets and embedded inside carbon fiber composite laminates. MESC has been demonstrated through earlier experiments on its exceptional behavior as a structural element as well as a battery. However, the inherent complex infrastructure of the MESC design has created significant challenges in simulation and modeling. A novel homogenization technique was adopted to characterize the multi-layer properties of battery material using physics-based constitutive equations combined with nonlinear deformation theories to handle the interface between the battery layers. Second, mechanical damage and failure modes among battery materials, polymer reinforcements, and carbon fiber-polymer interfaces were characterized through appropriate models and experiments. The model of MESCs has been implemented in a commercial finite element code in ABAQUS. A comparison of structural response and failure modes from numerical simulations and experimental tests are presented. The results of the study showed that the predictions of elastic and damage responses of MESCs at various loading conditions agreed well with the experimental data. © 2021

Materials Science↗
DOE OSTI2021

A simple numerical method for evaluating heat dissipation from curved wires with periodic applied heating

In many situations, the dual-purpose heater/thermometer wires used in the three-omega method – one of the most precise and sensitive techniques for measuring the thermal conductivity of thin films and interfaces – must include bends and curves to avoid obstructions on the surface of a sample. Although the three-omega analysis assumes that the heating wire is infinitely long and straight, recent experimental work has demonstrated that in some cases curved-wire geometries can be used without introducing detectable systematic error. We describe a general numerical method that can be used to calculate the temperature of three-omega heating wires with arbitrary wire geometries. Furthermore, this method provides experimentalists with a simple quantitative procedure for calculating how large the systematic error caused by a particular wire asymmetry will be. We show calculations of two useful cases: a straight wire with a single bend of arbitrary angle and a wire that forms a circle. We find that the amplitude of the in-phase temperature oscillations near a wire that forms a circle differs from the prediction using the analytic straight-line source solution by <12%, provided that the thermal penetration depth is less than ten times the radius of curvature of the wire path. The in-phase temperature amplitude 1.5 wire widths away from a 90° bend in a wire is within 11% of the straight-line source prediction for all penetration depths greater than the wire width. Our calculations indicate that the straight-line source solution breaks down significantly when the wire bend angle is less than 45°.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Glow Discharge Boronization and Real-Time Boronization Using an Impurity Powder Dropper in LHD

In the Large Helical Device (LHD), boronization is utilized for wall conditioning, employing two distinct methods: glow discharge using diborane gas and the injection of boron powder into the plasma via an impurity powder dropper (IPD). Glow discharge boronization is conducted after maintenance when the vacuum vessel is exposed to air and just before the start of an experimental campaign. In contrast, real-time boronization using the IPD is performed during discharges. The distribution of boron on plasma-facing components was assessed through long-term material probes in the case of glow discharge and through numerical simulations for the IPD method. Results indicated non-uniform boron distributions in both cases. The impact of these boronization techniques on impurity levels and the lifetime of their effects were also examined. Both methods effectively reduced oxygen and iron in the plasma, with the iron impurities originating from the first wall. For glow discharge boronization, the reduction in oxygen levels persisted until the end of the experimental campaign, whereas the reduction in iron levels only lasted a few days. In the case of boron injection during long pulse discharges exceeding 200 s, the suppression of oxygen and iron was observed. The effect on iron quickly diminished within several seconds after the injection ended, whereas the effect on oxygen lasted more than 100 s. The difference in the lifetimes of the boronization effects is further discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Validation of Spent Fuel Dose Rate Calculations Using NDA Measurement Campaign Data

The predicted radiation gamma dose rate is a key parameter to determine if a spent nuclear fuel assembly is self-protected from misuse as well as an indicator for worker/inspector safety. Recent studies showed that using state-of-the-art techniques for dose rate calculations, including the one adopted for this study, the numerical calculations of the gamma dose rates are three times lower than those from the most commonly-cited reference values for spent LWR fuel used by domestic and international agencies. With only a few publically available experimental data points worldwide for this type of benchmark study, the gamma flux measurements from the Spent Fuel Non Destructive Assay (NDA) measurement campaign can be used to help validate the physics approach used in calculating gamma dose rates. Specifically, this data would validate the gamma flux calculation from spent fuel, which is a key component of these multi-physics techniques for dose rate assessments; dose rate measurements were not part of the NDA campaign. These data are primarily the passive gamma measurements of the spent nuclear fuel assemblies from the Swedish Nuclear Fuel and Waste Management Company (SKB).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Subvoxel Control of Fiber Orientation via Multidirectional Shearing in 3D Printing

Anisotropy, the characteristic of materials exhibiting different properties based on their direction, is widespread in nature. Conventional manufacturing techniques often fall short in recreating such complex anisotropy. While 3D printing allows for precise fiber deposition in anisotropic composites, previous studies have only achieved bulk reorientation of fibers at the voxel level in two dimensions. This limits the replication of localized 3D anisotropies found in natural materials. To address this, a novel 3D printing technique is presented that enables subvoxel control of fiber orientation in all three directions using multidirectional shearing via nozzle rotation and inclination. The fiber orientation control in these experimental tests is driven by a numerical model, enabling a fully digital approach to program microstructures within a strand. This programmability is demonstrated through mechanical and thermal tests, illustrating a localized and controllable response to external stimuli. Achieving such complex anisotropy holds potential across several fields, including wearables and biomedical implants, lightweight composite structures, and energy storage technologies such as batteries and supercapacitors.

3D printing↗
DOE OSTI2025

How are Heterogeneous Nucleation Rate Observations Influenced by Instrument Resolution?

Experimental measurements of the heterogeneous nucleation rate rely on counting the number of nuclei with time. However, the size of a thermodynamically stable nucleus is often a few nanometers in diameter and is below the resolution of most (in situ) measurement techniques that provide a statistically valid sample. Due to the finite resolution of the instruments and analysis methods, it is challenging to capture the incipient nuclei and the subsequent evolution of nuclei density over time. In this work, we demonstrate the impact of instrument resolution on observed nuclei densities by comparing numerical modeling with experimental results. Further, to achieve this, we implemented heterogeneous nucleation within the pore-scale reactive transport modeling framework using classical nucleation theory (CNT). We compared the modeling results with nucleation rates measured using X-ray nanotomography (XnT) and evaluated how these impact the apparent values of the prefactor and interfacial energy based on CNT and the crystal growth rate. Specifically, we applied a resolution threshold (artificial resolution limit) in the model during nuclei counting to resemble an experimental resolution, ranging from 15 to 500 nm. The findings reveal that the instrument resolution significantly impacts the apparent prefactor and interfacial energy. Both apparent prefactor and interfacial energy decrease with a decrease in the instrument resolution. While deviation in the prefactor due to resolution is anticipated, those in the interfacial energy are unexpected. The approach described here allows one to correct apparent nucleation rates that depend on the instrument’s resolution to derive “intrinsic” CNT parameters for the prefactor and interfacial energy.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Compaction of crushed salt for safe containment – overview of the KOMPASS project

Abstract. In Germany, rock salt formations are possible host rock candidates for a repository for heat-emitting radioactive waste. The safety concept of a repository in salt bases on a multibarrier system consisting mainly of the geological barrier salt and geotechnical seals ensuring safe containment. Crushed salt will be used for backfilling of cavities and sealing measures in drifts and shafts due to its favourable properties and its easy availability (mined-off material). The creep of the rock salt leads to crushed salt compaction with time. Thereby, the crushed salts' porosity is reduced from the initial porosity of 30 %–40 % to a value comparable to the porosity of undisturbed rock salt (≤1 %). In such low porosity ranges, technical impermeability is assumed. The compaction behaviour of crushed salt is rather complex and involves several coupled THM processes (Kröhn et al., 2017; Hansen et al., 2014). It is influenced by internal properties like humidity and grain size distribution, as well as boundary conditions such as temperature, compaction rate or stress state. However, the current process understanding has some important gaps referring to the material behaviour, experimental database and numerical modelling. It needs to be extended and validated, especially in the low porosity range. The objective of the KOMPASS project was development of methods and strategies for the reduction of deficits in the prediction of crushed salt compaction leading to an improvement of the prognosis quality. Key results are as follows (KOMPASS Phase 1, 2020): selection of an easily available and permanently producible synthetic crushed salt mixture, acting as a reference material for generic investigations; development and proof of different techniques for producing pre-compacted samples for further investigations; establishment of a tool of microstructure investigation methods to demonstrate the comparability of grain structures of pre-compacted samples with in-situ compacted material for future investigations; execution of various laboratory experiments using pre-compacted samples, e.g. long-term creep tests which deliver reliable information about time- and stress-dependent compaction behaviour; development of a complex experimental investigation strategy to derive necessary model parameters considering individual functional dependencies. Its technical feasibility was successfully verified; benchmarking with various existing numerical models using datasets from three different triaxial long-term tests. The result was not entirely satisfactory; however, the number of influencing factors is small and further validation work has to be done. Overall, the KOMPASS project has made significant progress in the approaches to solving the outstanding question, building the basis for further investigations.

Friedenberg, Larissa↗
DOE OSTI2024

Adamantine 1.0: A Thermomechanical Simulator for Additive Manufacturing

Adamantine is a thermomechanical simulation code that is written in C++ and built on top of deal.II (Arndt et al., 2023), p4est (Burstedde et al., 2011), ArborX (Lebrun-Grandié et al., 2020), Trilinos (The Trilinos Project Team, 2020), and Kokkos (Trott et al., 2022). Adamantine was developed with additive manufacturing in mind and it is particularly well adapted to simulate fused filament fabrication, directed energy deposition, and powder bed fusion. Adamantine employs the finite element method with adaptive mesh refinement to solve a nonlinear anisotropic heat equation, enabling support for various additive manufacturing processes. It can also perform elastoplastic and thermoelastoplastic simulations. It can handle materials in three distinct phases (solid, liquid, and powder) to accurately reflect the physical state during different stages of the manufacturing process. To enhance simulation accuracy, adamantine incorporates data assimilation techniques (Asch et al., 2016). This allows it to integrate experimental data from sensors like thermocouples and infrared (IR) cameras. This combined approach helps account for errors arising from input parameters, material properties, models, and numerical calculations, leading to more realistic simulations that reflect what occurs in a particular print.

36 MATERIALS SCIENCE↗