Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multiphase”

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 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Modeling supercritical CO2 injection induced rupture of a minor fault embedded in a poroelastic layered reservoir-caprock system

CO2 injection for geologic carbon sequestration involves hydromechanical processes that lead to changes in fluid pressure and stresses that can activate existing faults. This paper presents a new method and workflow of modeling fault activation considering more complex three-dimensional geometry of natural faults using the TOUGH-FLAC multiphase fluid flow and geomechanical simulator. In this method and workflow, FLAC3D mechanical interfaces and TOUGH3 finite volume elements are discretized using computer aided design and gridding software along with a tailored mesh translation routine. The method and workflow are demonstrated with a model of a curved minor fault embedded in a poro-elastic layered reservoir-caprock system. The model is used for a comprehensive sensitivity analysis of fault responses to fault length, injection mass rate, injection schedule, well-fault distance, and well locations versus fault location. Four metrics (CO2 plume, shear state of fault, pressure and stress path at fault monitoring points) are selected to assess CO2 migration, pressure change, and the reactivation of faults. The results reveal that CO2 can bypass around the tip of the minor impermeable fault, building up pressure and poro-elastic stress on both sides that tends to impede fault rupture. Our study shows the benefit of carefully designing the injection to achieve the targeted final storage volume, starting at a relatively low rate for considerable time, and then ramping up the injection rate to the full rate of injection. The initial low injection has two distinct benefits: (1) it allows for the formation of an extensive CO2 plume with a much higher mobility through a low viscosity that will result in a lower pressure for a given injection rate, and (2) it allows for gradual build-up of horizontal poro-elastic stress within the reservoir that will tend to impede activation of steeply dipping faults. The injection scenario starting at a low injection rate, denoted here as conservative injection, can significantly reduce the risk of fault activation as high fluid mobility and reservoir strengthening poro-elastic stress has been established long before reaching the peak injection rates. Moreover, simultaneous injection in two injection wells on both sides of fault can provide further reservoir strengthening through poro-elastic stress buildup acting on a fault under normal faulting stress regime. The findings presented in the paper can provide practical and effective guidance on long-term, safe, and reliable geological CO2 storage.

Cao, Meng↗
DOE OSTI2025

Factors controlling injection-induced rupture of intersecting faults during geological sequestration of CO 2

This study addresses coupled multiphase fluid flow and geomechanics effects on potential fault activation associated with subsurface CO 2 injection around intersecting faults. An enhanced fault-representation model is used to capture geomechanical responses of two intersecting faults with finite length during CO 2 injection. The faults are embedded in a strike-slip stress regime of a caprock-reservoir-basement system with the faults represented by zero-thickness interfaces with adjacent finite-thickness damage zones. A sensitivity analysis is conducted to study the effect of fault permeability, slip-weakening behavior, well location relative to the orientation of faults, and well placement (the number and location of injection wells). Five metrics (pressure, CO 2 plume, shear state on the fault, as well as shear displacement and stress path at selected fault monitoring points) are selected to assess CO 2 migration and reactivation of intersecting faults. The results show that induced ruptures are favored by low permeability faults due to high pressure buildup and by slip-weakening behavior resulting from fault strength reduction. The location of one injection well relative to fault orientation determines the magnitude of changes in effective normal stress and shear stress, affecting the location of induced ruptures. Well placement (two injection wells used in the paper) dominates pressure diffusion around the intersection and tips of faults. This redistributes changes in effective normal stress caused by each injection well, influencing the spatial distribution of ruptures along faults. A larger injection volume induces far-field ruptures that are controlled by stress transfer within the injection layer. The findings presented here can provide valuable insights into engineering operations for a long-term, safe, and reliable geologic CO 2 storage.

Fault permeability↗
DOE OSTI2014

Interobserver Variability in Target Definition for Hepatocellular Carcinoma With and Without Portal Vein Thrombus: Radiation Therapy Oncology Group Consensus Guidelines

Purpose: Defining hepatocellular carcinoma (HCC) gross tumor volume (GTV) requires multimodal imaging, acquired in different perfusion phases. The purposes of this study were to evaluate the variability in contouring and to establish guidelines and educational recommendations for reproducible HCC contouring for treatment planning. Methods and Materials: Anonymous, multiphasic planning computed tomography scans obtained from 3 patients with HCC were identified and distributed to a panel of 11 gastrointestinal radiation oncologists. Panelists were asked the number of HCC cases they treated in the past year. Case 1 had no vascular involvement, case 2 had extensive portal vein involvement, and case 3 had minor branched portal vein involvement. The agreement between the contoured total GTVs (primary + vascular GTV) was assessed using the generalized kappa statistic. Agreement interpretation was evaluated using Landis and Koch's interpretation of strength of agreement. The S95 contour, defined using the simultaneous truth and performance level estimation (STAPLE) algorithm consensus at the 95% confidence level, was created for each case. Results: Of the 11 panelists, 3 had treated >25 cases in the past year, 2 had treated 10 to 25 cases, 2 had treated 5 to 10 cases, 2 had treated 1 to 5 cases, 1 had treated 0 cases, and 1 did not respond. Near perfect agreement was seen for case 1, and substantial agreement was seen for cases 2 and 3. For case 2, there was significant heterogeneity in the volume identified as tumor thrombus (range 0.58-40.45 cc). For case 3, 2 panelists did not include the branched portal vein thrombus, and 7 panelists contoured thrombus separately from the primary tumor, also showing significant heterogeneity in volume of tumor thrombus (range 4.52-34.27 cc). Conclusions: In a group of experts, excellent agreement was seen in contouring total GTV. Heterogeneity exists in the definition of portal vein thrombus that may impact treatment planning, especially if differential dosing is contemplated. Guidelines for HCC GTV contouring are recommended.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2020

Online ion mobility spectrometry of nanoparticle formation by non-thermal plasma conversion of metal salts in liquid aerosol droplets

The synthesis of nanoparticles by reaction of liquid aerosol droplets containing precursors in a flow-through, atmospheric-pressure, non-thermal plasma offers a continuous, scalable, substrate- and stabilizer-free approach for direct deposition into liquids or onto soft substrates. However, the combination of multiphase and non-equilibrium chemistry makes the process complicated and poorly understood. Here, we present ion mobility spectrometry measurements of liquid water droplets containing silver nitrate passing through an atmospheric-pressure dielectric barrier discharge reactor that allows us to monitor silver nanoparticle formation online for the first time. Mobility diameter distributions were obtained with the plasma on and off, and exhibited a shift, which was related to the degree of conversion of silver nitrate. The silver nanoparticles were also collected and characterized by UV–visible absorbance spectroscopy and transmission electron microscopy to support the online measurements. Importantly, negligible conversion was found when the water was removed by a diffusion dryer, suggesting that the key reducing species are in the liquid phase, such as solvated electrons. Finally, the study demonstrates how ion mobility spectrometry measurements can be applied to provide insight into this approach to nanoparticle synthesis.

42 ENGINEERING↗
DOE OSTI2013

Single phase melt processed powellite (Ba,Ca)MoO4 for the immobilization of Mo-rich nuclear waste

Cs is one of the more problematic fission product radionuclides to immobilize due to its high volatility at elevated temperatures, ability to form water soluble compounds, and its mobility in many host materials. The hollandite structure is a promising crystalline host for Cs immobilization and has been traditionally fabricated by solid state sintering methods. This study presents the structure and performance of Ba{sub 1.0}Cs{sub 0.3}A{sub 2.3}Ti{sub 5.7}O{sub 16}; A = Cr, Fe, Al hollandite fabricated by melt processing. Melt processing is considered advantageous given that melters are currently in use for High Level Waste (HLW) vitrification in several countries. This work details the impact of Cr additions that were demonstrated to i) promote the formation of a Cs containing hollandite phase and ii) maintain the stability of the hollandite phase in reducing conditions anticipated for multiphase waste form processing.

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

Surface enhanced performance of La 0.6 Sr 0.4 Co 0.2 Fe 0.8 O 3-δ cathodes by infiltration Pr-Ni-Mn-O progress

Herein the present study reports the enhancement of electrochemical oxygen reduction activity of porous La 0.6 Sr 0.4 Co 0.2 Fe 0.8 O 3–δ (LSCF) cathodes by coating a thin film of Pr-Ni-Mn oxide (PNM5) using a multi-step infiltration process. XRD examination reveals that PNM5 mainly contains a multiphase mixture of Pr 6 O 11 , PrNiO 3 , MnO and NiO. SEM morphology shows a thin PNM5 film and small particles are formed on the surface of LSCF backbone particles. Impedance spectrum analysis indicates that PNM5 infiltrated LSCF exhibits dramatically reduced polarization resistance (R p ), reaching R p of 0.244 Ω cm 2 at 973 K, which is one-half of the baseline LSCF cathode. The activation energy of LSCF cathodes infiltrated with PNM5 is 1.45 eV, slightly lower than the baseline LSCF cathode (1.77 eV). Distribution of relaxation time (DRT) function analysis shows PNM5 infiltration layer significantly promotes the oxygen reduction reaction (ORR) of cathode surface. With the increase in firing temperature, the total resistance increases and R p changes from ion transport to oxygen reduction reaction. Degradation rate of the PNM5-infiltrated LSCF is also lower, 0.02168% vs 0.07093% for the baseline LSCF over a 200-h period. A single cell testing indicated that the peak power density of the PNM5-infiltrated cell was increased by 140.66%. Overall, PNM5 could be a potential catalyst for boosting the performance of a commercial LSCF cathode for solid oxide fuel cells (SOFC).

36 MATERIALS SCIENCE↗
DOE OSTI2025

Structural evolution and magnetic hardness of (Sm,Zr)(Fe,Co,Ti) 12 alloy particles via reduction-diffusion

The quest for achieving high coercivity in Sm(Fe,Co,Ti) 12 alloys, despite their inherent strong magnetocrystalline anisotropy, has posed significant challenges. Recently, (Sm,Zr)(Fe,Co,Ti) 12 monocrystalline particles have exhibited coercivity μ 0 H c > 1.2 T, showcasing promising prospects and significant potential for both manufacturing and research endeavors. This study delves into the structural evolution of (Sm,Zr)(Fe,Co,Ti) 12 (1:12) alloy particles made via the calciothermic reduction-diffusion synthesis process as influenced by the molar ratios of Ca atoms to O 2- ions (Ca/O), annealing time and annealing temperature. Critical insight that informs conditions to optimize the magnetic response is gained via systematic experimentation and advanced electron microscopy. Complex structural features, including core-shell morphologies and intricate multiphase compositions within individual particles, are unveiled. An optimal Ca/O ratio of 1.30 produces particles with a coercivity up to μ 0 H c = 1.63 T, while higher Ca/O ratios induce the formation of a Sm-rich TbCu 7 -type (1: 7 ) phase, which only partially transforms into the desired 1:12 phase during annealing. Persistent remnants of the 1:7 phase locally impact atomic structure, particle morphology, and coercivity. Furthermore, these findings underscore the complex interplay between synthesis parameters, resulting structures, and magnetic properties, informing the design and optimization of high-performance permanent magnets comprised of the (1:12) compound.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Higher-order compositional modeling of three-phase flow in 3D fractured porous media based on cross-flow equilibrium

Numerical simulation of multiphase compositional flow in fractured porous media, when all the species can transfer between the phases, is a real challenge. Despite the broad applications in hydrocarbon reservoir engineering and hydrology, a compositional numerical simulator for three-phase flow in fractured media has not appeared in the literature, to the best of our knowledge. In this work, we present a three-phase fully compositional simulator for fractured media, based on higher-order finite element methods. To achieve computational efficiency, we invoke the cross-flow equilibrium (CFE) concept between discrete fractures and a small neighborhood in the matrix blocks. We adopt the mixed hybrid finite element (MHFE) method to approximate convective Darcy fluxes and the pressure equation. This approach is the most natural choice for flow in fractured media. The mass balance equations are discretized by the discontinuous Galerkin (DG) method, which is perhaps the most efficient approach to capture physical discontinuities in phase properties at the matrix-fracture interfaces and at phase boundaries. In this work, we account for gravity and Fickian diffusion. The modeling of capillary effects is discussed in a separate paper. We present the mathematical framework, using the implicit-pressure-explicit-composition (IMPEC) scheme, which facilitates rigorous thermodynamic stability analyses and the computation of phase behavior effects to account for transfer of species between the phases. A deceptively simple CFL condition is implemented to improve numerical stability and accuracy. We provide six numerical examples at both small and larger scales and in two and three dimensions, to demonstrate powerful features of the formulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

A numerical method for the quasi-incompressible Cahn–Hilliard–Navier–Stokes equations for variable density flows with a discrete energy law

In this paper, we investigate numerically a diffuse interface model for the Navier–Stokes equation with fluid–fluid interface when the fluids have different densities [48]. Under minor reformulation of the system, we show that there is a continuous energy law underlying the system, assuming that all variables have reasonable regularities. It is shown in the literature that an energy law preserving method will perform better for multiphase problems. Thus for the reformulated system, we design a C{sup 0} finite element method and a special temporal scheme where the energy law is preserved at the discrete level. Such a discrete energy law (almost the same as the continuous energy law) for this variable density two-phase flow model has never been established before with C{sup 0} finite element. A Newton method is introduced to linearise the highly non-linear system of our discretization scheme. Some numerical experiments are carried out using the adaptive mesh to investigate the scenario of coalescing and rising drops with differing density ratio. The snapshots for the evolution of the interface together with the adaptive mesh at different times are presented to show that the evolution, including the break-up/pinch-off of the drop, can be handled smoothly by our numerical scheme. The discrete energy functional for the system is examined to show that the energy law at the discrete level is preserved by our scheme.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Unconditionally energy stable time stepping scheme for Cahn–Morral equation: Application to multi-component spinodal decomposition and optimal space tiling

An unconditionally energy stable time stepping scheme is introduced to solve Cahn–Morral-like equations in the present study. It is constructed based on the combination of David Eyre's time stepping scheme and Schur complement approach. Although the presented method is general and independent of the choice of homogeneous free energy density function term, logarithmic and polynomial energy functions are specifically considered in this paper. The method is applied to study the spinodal decomposition in multi-component systems and optimal space tiling problems. A penalization strategy is developed, in the case of later problem, to avoid trivial solutions. Extensive numerical experiments demonstrate the success and performance of the presented method. According to the numerical results, the method is convergent and energy stable, independent of the choice of time stepsize. Its MATLAB implementation is included in the appendix for the numerical evaluation of algorithm and reproduction of the presented results. -- Highlights: •Extension of Eyre's convex–concave splitting scheme to multiphase systems. •Efficient solution of spinodal decomposition in multi-component systems. •Efficient solution of least perimeter periodic space partitioning problem. •Developing a penalization strategy to avoid trivial solutions. •Presentation of MATLAB implementation of the introduced algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Multiscale modelling of evolving foams

We present a set of multi-scale interlinked algorithms to model the dynamics of evolving foams. These algorithms couple the key effects of macroscopic bubble rearrangement, thin film drainage, and membrane rupture. For each of the mechanisms, we construct consistent and accurate algorithms, and couple them together to work across the wide range of space and time scales that occur in foam dynamics. These algorithms include second order finite difference projection methods for computing incompressible fluid flow on the macroscale, second order finite element methods to solve thin film drainage equations in the lamellae and Plateau borders, multiphase Voronoi Implicit Interface Methods to track interconnected membrane boundaries and capture topological changes, and Lagrangian particle methods for conservative liquid redistribution during rearrangement and rupture. We derive a full set of numerical approximations that are coupled via interface jump conditions and flux boundary conditions, and show convergence for the individual mechanisms. We demonstrate our approach by computing a variety of foam dynamics, including coupled evolution of three-dimensional bubble clusters attached to an anchored membrane and collapse of a foam cluster.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Lagrangian particle method for compressible fluid dynamics

A new Lagrangian particle method for solving Euler equations for compressible inviscid fluid or gas flows is proposed. Similar to smoothed particle hydrodynamics (SPH), the method represents fluid cells with Lagrangian particles and is suitable for the simulation of complex free surface / multi-phase flows. The main contributions of our method, which is different from SPH in all other aspects, are (a) significant improvement of approximation of differential operators based on a polynomial fit via weighted least squares approximation and the convergence of prescribed order, (b) a second-order particle-based algorithm that reduces to the first-order upwind method at local extremal points, providing accuracy and long term stability, and (c) more accurate resolution of entropy discontinuities and states at free inter-faces. While the method is consistent and convergent to a prescribed order, the conservation of momentum and energy is not exact and depends on the convergence order . The method is generalizable to coupled hyperbolic-elliptic systems. As a result, numerical verification tests demonstrating the convergence order are presented as well as examples of complex multiphase flows.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

A robust and efficient finite volume method for compressible inviscid and viscous two-phase flows [A robust and efficient finite volume method for compressible viscous two-phase flows]

A robust and efficient density-based finite volume method is developed for solving the six-equation single pressure system of two-phase flows at all speeds on hybrid unstructured grids. Unlike conventional approaches where an expensive exact Riemann solver is normally required for computing numerical fluxes at the two-phase interfaces in addition to AUSM-type fluxes for single-phase interfaces in order to maintain stability and robustness in cases involving interactions of strong pressure and void-fraction discontinuities, a volume-fraction coupling term for the AUSM + -up fluxes is introduced in this work to impart the required robustness without the need of the exact Riemann solver. The resulting method is significantly less expensive in regions where otherwise the Riemann solver would be invoked. A transformation from conservative variables to primitive variables is presented and the primitive variables are then solved in the implicit method in order for the current finite volume method to be able to solve, effectively and efficiently, low Mach number flows in traditional multiphase applications, which otherwise is a great challenge for the standard density-based algorithms. Here, a number of benchmark test cases are presented to assess the performance and robustness of the developed finite volume method for both inviscid and viscous two-phase flow problems. The numerical results indicate that the current density-based method provides an attractive and viable alternative to its pressure-based counterpart for compressible two-phase flows at all speeds.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Latent space dynamics identification for interface tracking with application to shock-induced pore collapse

Capturing sharp, evolving interfaces remains a central challenge in reduced-order modeling, especially when data is limited and the system exhibits localized nonlinearities or discontinuities. Here, we propose LaSDI-IT (Latent Space Dynamics Identification for Interface Tracking), a data-driven framework that combines low-dimensional latent dynamics learning with explicit interface-aware encoding to enable accurate and efficient modeling of physical systems involving moving material boundaries. At the core of LaSDI-IT is a revised autoencoder architecture that jointly reconstructs the physical field and an indicator function representing material regions or phases, allowing the model to track complex interface evolution without requiring detailed physical models or mesh adaptation. The latent dynamics are learned through linear regression in the encoded space and generalized across parameter regimes using Gaussian process interpolation with greedy sampling. We demonstrate LaSDI-IT on the problem of shock-induced pore collapse in high explosives, a process characterized by sharp temperature gradients and dynamically deforming pore geometries. The method achieves relative prediction errors below 9% across the parameter space, accurately recovers key quantities of interest such as pore area and hot spot formation, and matches the performance of dense training with only half the data. This latent dynamics prediction was 10 6 times faster than the conventional high-fidelity simulation, proving its utility for multi-query applications. These results highlight LaSDI-IT as a general, data-efficient framework for modeling discontinuity-rich systems in computational physics, with potential applications in multiphase flows, fracture mechanics, and phase change problems.

Gaussian process↗
DOE OSTI2021

Phase evolution during reactive flash sintering of Li 6.25 Al 0.25 La 3 Zr 2 O 12 starting from a chemically prepared powder

Here, reactive flash sintering (RFS) of a chemically prepared multiphase precursor powder was performed to fabricate Li 6.25 Al 0.25 La 3 Zr 2 O 12 (Al-LLZO) ceramics. This approach allowed for obtaining single-phase dense samples in a remarkably short processing time of 30 s, at a furnace temperature of 600 °C, with an electric field of 50 V cm -1 and a current limit of 150 mA mm -2 . The ceramics display high bulk conductivity of 0.18 mS cm -1 at room temperature. Furthermore, phase evolution is studied by in-situ X-ray diffraction during: i) conventional heating and ii) RFS under current rate mode. As expected, the intermediate phases progressively dissolved into the Al-LLZO matrix by conventional heating. On the other hand, RFS promoted the growth of the intermediate La 2 Zr 2 O 7 , an effect that was overcome by the thermally driven formation of Al-LLZO at higher temperatures. The observed different reaction pathway suggests that RFS can be used for stabilizing phases that are not thermodynamically favored upon conventional heating.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A deep learning-based workflow for fast prediction of 3D state variables in geological carbon storage: A dimension reduction approach

Deep learning (DL) models are extensively used as surrogate models for high-fidelity simulations of multiphase fluid flow in porous media at large scales, enabling fast forecasts of the spatial–temporal evolution of three-dimensional (3D) state variables in geological carbon storage (GCS). However, training these models in high-dimensional space remains computationally demanding and prone to overfitting because of limited training data. This paper presents a novel workflow to address these challenges by integrating dimension reduction (DR) methods. Here, the proposed workflow employed pre-trained DR models to extract the latent variables of geological models and state variables and utilized the multi-layer perceptron (MLP) for constructing mapping functions between the input and output variables in latent spaces. Subsequently, the pre-trained reconstruction models converted the MLP-predicted latent state variables to their original high-dimensional form. Furthermore, we proposed a novel strategy for the DR and reconstruction of 3D saturation fields to account for the unique data characteristics of sparsity, nonuniformity, and discontinuity. The proposed strategy applied PCA and inverse PCA for 2D average saturation fields and developed a DL-based 3D reconstruction model, leveraging three 2D average saturation fields as input to produce a 3D saturation field as output. The pre-training of DR and reconstruction models and training of MLP models were conducted on 84 Gulf of Mexico (GoM) simulations and evaluated on 12 testing simulations. Each simulation contained 720 monthly time steps, with the first 360 months as the injection period and the rest as the post-injection period. The proposed workflow, incorporating DR and DL models, accurately predicts the normalized 3D pressure fields, achieving mean square error (MSE) of 2.92 × 10 -7 compared to the ground truth obtained from a full-physics simulator. Furthermore, the proposed strategy outperformed PCA and convolutional autoencoder (CAE) models on 3D saturation fields, resulting in minor workflow prediction errors with an MSE of 2.93 × 10 -5 . The results suggest the proposed workflow provides sufficient predictive fidelity across temporal and spatial scales, and enables a speedup of 160 times compared to the full-physics simulator, facilitating improved decision-making and risk assessment for large-scale GCS management in real-time scenarios.

3D reconstruction model↗
DOE OSTI2025

Multi-physics melt pool modeling and process optimization for laser direct energy deposition of Nb-based refractory C103: Defect formation, geometric precision, and process mapping

Recent developments in additive manufacturing (AM) technology have reignited interest in the fabrication of the Nb-based refractory C103 alloy offering solutions to the challenges posed by traditional manufacturing methods. However, the limited numerical and experimental studies on laser direct energy deposition (DED) of C103 have hindered the understanding of the relationships between process parameters and build quality. This has made it challenging to consistently produce parts with the desired quality and microstructure suitable for critical applications. In this study, we focus on optimizing the laser DED process for C103 by employing a hybrid approach that combines experimental techniques and computational fluid dynamics (CFD). This approach facilitates the development of process maps for defect detection and geometric precision. To achieve this, multi-layer C103 samples were fabricated using laser DED under various process parameters, enabling the creation of a process map for defect detection. Additionally, a multi-physics, multiphase simulation framework was developed within a high-performance computing (HPC) environment to establish process maps for geometric precision. Using these process maps, printability windows were identified for achieving both the desired geometric accuracy and defect-free prints. It was observed that prints with a power-to-velocity (P/V) ratio close to unity resulted in defect-free outcomes. This study provides a foundation for reducing design lead time and rejected parts, ultimately optimizing the laser DED process for C103.

Defect formation and geometric precision↗
DOE OSTI2017

Formulation and validation of a reduced order model of 2D materials exhibiting a two-phase microstructure as applied to graphene oxide [Formulation and validation of a reduced order model of 2D materials exhibiting a two-phase microstructure]

Novel 2D materials, e.g., graphene oxide (GO), are attractive building blocks in the design of advanced materials due to their reactive chemistry, which can enhance interfacial interactions while providing good in-plane mechanical properties. Recent studies have hypothesized that the randomly distributed two-phase microstructure of GO, which arises due to its oxidized chemistry, leads to differences in nano- vs meso–scale mechanical responses. However, this effect has not been carefully studied using molecular dynamics due to computational limitations. Herein, a continuum mechanics model, formulated based on density functional based tight binding (DFTB) constitutive results for GO nano-flakes, is establish for capturing the effect of oxidation patterns on the material mechanical properties. GO is idealized as a continuum heterogeneous two-phase material, where the mechanical response of each phase, graphitic and oxidized, is informed from DFTB simulations. A finite element implementation of the model is validated via MD simulations and then used to investigate the existence of GO representative volume elements (RVE). We find that for the studied GO, an RVE behavior arises for monolayer sizes in excess to 40 nm. Moreover, we reveal that the response of monolayers with two main different functional chemistries, epoxide-rich and hydroxyl–rich, present distinct differences in mechanical behavior. In addition, we explored the role of defect density in GO, and validate the applicability of the model to larger length scales by predicting membrane deflection behavior, in close agreement with previous experimental and theoretical observations. As such the work presents a reduced order modeling framework applicable in the study of mechanical properties and deformation mechanisms in 2D multiphase materials.

36 MATERIALS SCIENCE↗