Search NASA⌕ Search

SEARCH · Search NASA

Results for “Complex Geometry”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

In Situ High Energy X-ray Diffraction Characterization of Phase Transformations and Mechanical Behaviors in Rapidly Solidified Titanium and Stainless Steel Alloys [Thesis]

Advanced manufacturing techniques like additive manufacturing (AM) have poised themselves to revolutionize metal manufacturing. A wide range of AM techniques are capable of manufacturing metal components with unique, complex geometries and hastening the scientific-engineering-development cycle. Metal AM relies on a layer-by-layer rapid manufacturing process to build components from the substrate up. Rapid solidification is a large departure from traditional metal manufacturing due to its complex physics. Characterization of rapid solidification is difficult, stemming from the small volumes used in AM and the fast dynamics of the process. High energy X-ray diffraction (HEXRD) is a solution to the characterization problems of rapidly solidified alloys and AM. HEXRD can probe small volumes at fast rates and provides a wide range of thermomechanical and kinetic information. This thesis presents the application of HEXRD to rapidly solidified titanium and stainless steel alloys through a series of case studies. In the first two studies, HEXRD is applied to rapidly solidified titanium and stainless steel welds. The materials are characterized for their temperature history, phase changes, kinetics, and microstructural evolution. In the next case study, HEXRD is applied to characterize phase changes in elastocaloric NiTi shape memory alloys (SMAs) under thermomechanical load. HEXRD, in conjunction with other tools, is used to explain the superior performance of the additively manufactured SMAs. In the final two case studies, HEXRD is used to measure the mechanical response of AM parts with complex geometries; namely, the octet truss lattice. Diffraction reveals a wide range of materials information about the AM microstructure including unexpected phases, texture, and mechanical response to loading. The mechanical results from HEXRD and then compared with theoretical predictions about the performance of octet truss lattices. Summarily, HEXRD is a diverse tool that is poised to address the complex characterization problems of many aspects of the additive manufacturing process.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A comprehensive and fair comparison of two neural operators (with practical extensions) based on $\mathrm{FAIR}$ data

Neural operators can learn nonlinear mappings between function spaces and offer a new simulation paradigm for real-time prediction of complex dynamics for realistic diverse applications as well as for system identification in science and engineering. Herein, we investigate the performance of two neural operators, which have shown promising results so far, and we develop new practical extensions that will make them more accurate and robust and importantly more suitable for industrial-complexity applications. The first neural operator, DeepONet, was published in 2019 (Lu et al., 2019), and its original architecture was based on the universal approximation theorem of Chen & Chen (1995). The second one, named Fourier Neural Operator or FNO, was published in 2020, and it is based on parameterizing the integral kernel in the Fourier space. DeepONet is represented by a summation of products of neural networks (NNs), corresponding to the branch NN for the input function and the trunk NN for the output function; both NNs are general architectures, e.g., the branch NN can be replaced with a CNN or a ResNet. According to Kovachki et al. (2021), FNO in its continuous form can be viewed conceptually as a DeepONet with a specific architecture of the branch NN and a trunk NN represented by a trigonometric basis. In order to compare FNO with DeepONet computationally for realistic setups, we develop several extensions of FNO that can deal with complex geometric domains as well as mappings where the input and output function spaces are of different dimensions. We also develop an extended DeepONet with special features that provide inductive bias and accelerate training, and we present a faster implementation of DeepONet with cost comparable to the computational cost of FNO, which is based on the Fast Fourier Transform. Here we consider 16 different benchmarks to demonstrate the relative performance of the two neural operators, including instability wave analysis in hypersonic boundary layers, prediction of the vorticity field of a flapping airfoil, porous media simulations in complex-geometry domains, etc. We follow the guiding principles of FAIR (Findability, Accessibility, Interoperability, and Reusability) for scientific data management and stewardship. The performance of DeepONet and FNO is comparable for relatively simple settings, but for complex geometries the performance of FNO deteriorates greatly. We also compare theoretically the two neural operators and obtain similar error estimates for DeepONet and FNO under the same regularity assumptions.

42 ENGINEERING↗
DOE OSTI2018

Additive Manufacturing: A New Paradigm for the Next Generation of High-Power-Density Direct-Drive Electric Generators

In recent years, there has been a growing demand for high-power-density direct-drive generators in the wind industry owing to their high reliability, torque per unit volume, and conversion efficiencies. However, direct-drive wind turbine generators are very large, low-speed electric machines, which pose remarkable design and manufacturing issues that challenge their upscaling potential and cost of implementation. With air-gap tolerance as the main design driver, the need for high stiffness shifts the focus toward support-structure design that forms a significant portion of the generator's total mass. Existing manufacturing processes allow the use of segmented-steel-weldment disk or spoke-arm assemblies that yield stiffer structures per unit mass but tend to be heavier and more expensive to build. As a result, there is a need for a transformative approach to realize lightweight designs that can also facilitate series production at competitive costs. Inspired by recent developments in metal additive manufacturing (AM), we explore a new freedom in the structural design space with a high potential for weight savings in direct-drive generators. This includes the feasibility of using nonconventional complex geometries, such as lattice-based structures as structurally efficient options. Powder-binder jetting of a sand-cast mold was identified as the most feasible AM technology to produce large-scale generator rotor structures with complex geometry. A parametric optimization study was performed and optimized results within deformation and mass constraints were found for each design. The response to the maximum Maxwell stress due to unbalanced magnetic pull was also explored for each design. Further, a topology optimization was applied for each parameter-optimized design to validate results and provide insights into further mass reduction. These novel designs catered for AM are compared in both deflection and mass to conventional rotor designs using NREL's systems engineering design tool, GeneratorSE. The optimized lattice design with a U-beam truss resulted in a 24% reduction in structural mass of the rotor and 60% reduction in radial deflection. It is demonstrated that additive manufacturing shifts the focus from manufacturability constraints toward lower mass.

17 WIND ENERGY↗
NASA NTRS1994

Asymptotic behavior of solutions of the renormalization group K-epsilon turbulence model

Presently, the only efficient way to calculate turbulent flows in complex geometries of engineering interest is to use Reynolds-average Navier-Stokes (RANS) equations. As compared to the original Navier-Stokes problem, these RANS equations posses much more complicated nonlinear structure and may exhibit far more complex nonlinear behavior. In certain cases, the asymptotic behavior of such models can be studied analytically which, aside from being an interesting fundamental problem, is important for better understanding of the internal structure of the models as well as to improve their performances. The renormalization group (RNG) K-epsilon turbulence model, derived directly from the incompresible Navier-Stokes equations, is analyzed. It has already been used to calculate a variety of turbulent and transitional flows in complex geometries. For large values of the RNG viscosity parameter, the model may exhibit singular behavior. In the form of the RNG K-epsilon model that avoids the use of explicit wall functions, a = 1, so the RNG viscosity parameter must be smaller than 23.62 to avoid singularities.

Yakhot, A.↗
DOE OSTI2018

Additive Manufacturing: A New Paradigm for the Next Generation of High-Power-Density Direct-Drive Electric Generators: Preprint

In recent years, there has been a growing demand for high-power-density direct-drive generators in the wind industry owing to their high reliability, torque per-unit volume, and conversion efficiencies. However, direct-drive wind turbine generators are very large, low-speed electric machines, which pose remarkable design and manufacturing issues that challenge their upscaling potential and cost of implementation. With air-gap tolerance as the main design driver, the need for high stiffness shifts the focus toward support-structure design that forms a significant portion of the generator's total mass. Existing manufacturing processes allow the use of segmented-steel-weldment disk or spoke-arm assemblies that yield stiffer structures per unit mass but tend to be heavier and more expensive to build. As a result, there is a need for a transformative approach to realize lightweight designs that can also facilitate series production at competitive costs. Inspired by recent developments in metal additive manufacturing (AM), we explore a new freedom in the structural design space with a high potential for weight savings in direct-drive generators. This includes the feasibility of using nonconventional complex geometries, such as lattice-based structures as structurally efficient options. Powder-binder jetting of a sand-cast mold was identified as the most feasible AM technology to produce large-scale generator rotor structures with complex geometry. A parametric optimization study was performed and optimized results within deformation and mass constraints were found for each design. The response to the maximum Maxwell stress due to unbalanced magnetic pull was also explored for each design. Further, a topology optimization was applied for each parameter-optimized design to validate results and provide insights into further mass reduction. These novel designs catered for AM are compared in both deflection and mass to conventional rotor designs using NREL's systems engineering design tool, GeneratorSE. The optimized lattice design with a Ubeam truss resulted in a 24% reduction in structural mass of the rotor and 60% reduction in radial deflection. It is demonstrated that additive manufacturing shifts the focus from manufacturability constraints toward lower mass.

17 WIND ENERGY↗
DOE OSTI2019

Evaluation of Binder Jetting and Chemical Vapor Infiltration as a Production Method for Silicon Carbide Parts for Nuclear Applications

Application of silicon carbide (SiC) is often sought within nuclear energy systems due to its high thermal conductivity, mechanical and environmental stability at elevated temperatures, corrosion resistance, and exceptional displacement damage irradiation tolerance. As such, the transformational challenge reactor (TCR) program is considering its application in its additively manufactured core. SiC may be used as the matrix for the TCR fuel blocks that comprise of TRISO fuel particles embedded inside this material or used as core non-fuel-bearing peripherical structures such as bottom support manifolds. The core peripherical structures serves to support the overall core geometry, direct and guide coolant flow, and also reflect back leaked neutrons into the core. SiC has very low neutron absorption and adequate scattering cross sections that makes it an ideal material for utilization in the TCR core. Like many ceramics, however, SiC is difficult and at times impractical to shape into complex geometries via traditional machining due to its brittleness. Further, large structures comprising nuclear-grade SiC, i.e. highly pure and crystalline, are costly. Additive manufacturing (AM), however, offers the ability to shape materials by adding layer-by-layer rather than subtracting, contributing new possibilities for previously difficult-to-shape materials. Thus, exploring AM as a method to shape SiC offers a highly cost-effective method and is strategic to developing the next generation of components for nuclear reactors. This report outlines an evaluation of binder jet additive manufacturing combined with chemical vapor infiltration (CVI) to produce nearly-dense silicon carbide artifacts in complex geometries suitable for advanced nuclear reactor designs. The focus is placed on the binder jet printing technologies at large scale for nonfuel-bearing structures relevant to the TCR core and other advanced reactors. Among all AM processes, binder jet delivers the highest resolution capability while still maintaining large print sizes, making it the ideal candidate for nuclear component production.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2012

The Pursuit of K: Reflections on the Current State-of-the-Art in Stress Intensity Factor Solutions for Practical Aerospace Applications

The elastic stress intensity factor (SIF, commonly denoted as K) is the foundation of practical fracture mechanics (FM) analysis for aircraft structures. This single parameter describes the first-order effects of stress magnitude and distribution as well as the geometry of both structure/component and crack. Hence, the calculation of K is often the most significant step in fatigue analysis based on FM. This presentation will provide several reflections on the current state-of-the-art in SIF solution methods used for practical aerospace applications, including a brief historical perspective, descriptions of some recent and ongoing advances, and comments on some remaining challenges. Newman and Raju made significant early contributions to practical structural analysis by developing closed-form SIF equations for surface and corner cracks in simplified geometries, often based on empirical fits of finite element (FE) solutions. Those solutions (and others like them) were sometimes revised as new analyses were conducted or limitations discovered. The foundational solutions have exhibited striking longevity, despite the relatively "coarse" FE models employed many decades ago. However, in recent years, the accumulation of different generations of solutions for the same nominal geometry has led to some confusion (which solution is correct?), and steady increases in computational capabilities have facilitated the discovery of inaccuracies in some (not all!) of the legacy solutions. Some examples of problems and solutions are presented and discussed, including the challenge of maintaining consistency with legacy design applications. As computational power has increased, the prospect of calculating large numbers of SIF solutions for specific complex geometries with advanced numerical methods has grown more attractive. Fawaz and Andersson, for example, have been generating literally millions of new SIF solutions for different combinations of multiple cracks under simplified loading schemes using p-version FE methods. These data are invaluable, but questions remain about their practical use, because the tabular databases of key results needed to support practical life analysis can occupy gigabytes of storage for only a few classes of geometries. The prospect of using such advanced numerical methods to calculate in real time only those K solutions actually needed to support a specific crack growth analysis is also tempting, but the stark reality is that the computational cost is still so high that the approach is not practical except for specific, critical application problems. Some thoughts are offered about alternative paradigms. Compounding approaches are some of the earliest building blocks of SIF development for more complex geometries. These approaches are especially attractive because of their very low computational cost and their conceptual robustness; they are, in some ways, an intriguing contrast and complement to the brute-force numerical methods. In recent years, researchers at NRC-Canada have published remarkable results showing how compounding approaches can be used to generate accurate solutions for very difficult problems. Examples are provided of some successes--and some limitations--using this approach. These closed-form, tabulated numerical, and compounding approaches have typically been used for simple remote loading with simple load paths to the crack. However, many significant cracks occur in complex stress gradient fields. This is a job for weight function (WF) methods, where the arbitrary stress distribution on the crack plane in the corresponding uncracked body (typically determined using FE methods) is used to determine K. Several significant recent advances in WF methods and solutions are highlighted here. Fueled by advanced 3D numerical methods, many new solutions have been generated for classic geometries such as surface and corner cracks with wide ranges of geometrical validity. A new WF formulation has also be developed for part-through cracks considering the arbitrary stress gradients in all directions in the crack plane (so-called bivariant solutions). Basic WF methods have recently been combined with analytical expressions for crack plane stresses to develop a large family of accurate SIF solutions for corner, surface, and through cracks at internal or external notches with very wide ranges of shapes, sizes, acuities, and offsets. Finally, WF solutions are much faster than FE or boundary element solutions, but can still be much slower than simple closed-form solutions, especially for bivariant solutions that can require 2D numerical integration. Novel pre-integration and dynamic tabular methods have been developed that substantially increase the speed of these advanced WF solutions. The practical utility of advanced SIF methods, including both WF and direct numerical methods, is greatly enhanced if the FM life analysis can be directly and efficiently linked with digital models of the actual structure or component (e.g., FE models for stress analysis). Two recent advances of this type will be described. One approach directly interfaces the FM life analysis with the FE model of the uncracked component (including stress results). Through a powerful graphical user interface, simplified FM life models can be constructed (and visualized) directly on the component model, with the computer collecting the geometry and stress gradient information needed for the life calculation. An even more powerful paradigm uses expert logic to automatically build an optimum simple fracture model at any and every desired location in the component model, perform the life calculation, and even generate fatigue crack growth life contour maps, all with minimal user intervention. This paradigm has also been extended to the automatic calculation of fracture risk, considering uncertainty or variability in key input parameters such as initial crack size or location. Another new integrated approach links the engineering life analysis, the component model, and a 3D numerical fracture analysis built with the same component model to generate a table of SIF values at a specific location that can then be employed efficiently to perform the life calculation. Some attention must be given to verification and validation (V&V) issues and challenges: how good are these SIF solutions, how good is good enough, and does anyone believe the life answer? It is important to think critically about the different sources of error or uncertainty and to perform V&V in a hierarchal, building-block manner. Some accuracy issues for SIF solutions, for example, may actually involve independent material behavior issues, such as constraint loss effects for crack fronts near component surfaces, and can be a source of confusion. Recommendations are proposed for improved V&V approaches. This presentation will briefly but critically survey the range of issues and advances mentioned above, with a particular view towards assembling an integrated approach that combines different methods to create practical tools for real-world design and analysis problems. Examples will be selectively drawn from the recent literature, from recent enhancements in the NASGRO and DARWIN computer codes, and from previously unpublished research

CraigMcClung, R.↗
NASA NTRS1989

Application of unstructured grid methods to steady and unsteady aerodynamic problems

The purpose is to describe the development of unstructured grid methods which have several advantages when compared to methods which make use of structured grids. Unstructured grids, for example, easily allow the treatment of complex geometries, allow for general mesh movement for realistic motions and structural deformations of complete aircraft configurations which is important for aeroelastic analysis, and enable adaptive mesh refinement to more accurately resolve the physics of the flow. Steady Euler calculations for a supersonic fighter configuration to demonstrate the complex geometry capability; unsteady Euler calculations for the supersonic fighter undergoing harmonic oscillations in a complete-vehicle bending mode to demonstrate the general mesh movement capability; and vortex-dominated conical-flow calculations for highly-swept delta wings to demonstrate the adaptive mesh refinement capability are discussed. The basic solution algorithm is a multi-stage Runge-Kutta time-stepping scheme with a finite-volume spatial discretization based on an unstructured grid of triangles in 2D or tetrahedra in 3D. The moving mesh capability is a general procedure which models each edge of each triangle (2D) or tetrahedra (3D) with a spring. The resulting static equilibrium equations which result from a summation of forces are then used to move the mesh to allow it to continuously conform to the instantaneous position or shape of the aircraft. The adaptive mesh refinement procedure enriches the unstructured mesh locally to more accurately resolve the vortical flow features. These capabilities are described in detail along with representative results which demonstrate several advantages of unstructured grid methods. The applicability of the unstructured grid methodology to steady and unsteady aerodynamic problems and directions for future work are discussed.

Batina, John T.↗
DOE OSTI2026

Enabling On-Demand Aerospace Component Manufacturing: Topology Optimization of GE Engine Bracket and Fabrication Using Metal FFF

Additive Manufacturing (AM) offers advantages over conventional manufacturing processes, particularly by reducing the number of parts produced through multistage combined technologies, but these often result in low manufacturing yields or require post-processing. AM facilitates the production of complex geometries with fine features, overhangs, and lattice structures. For instance, Laser Powder Bed Fusion (LPBF) technology enables the fabrication of intricate parts that can be easily post-processed by removing residual powder. Laser powder bed AM technologies are widely discussed in the literature due to their design freedom in creating complex geometries, with and without the need for support generation. However, rapid solidification due to a thermal gradient in the build direction, which leads to the formation of columnar grains and warpage, is one of the challenges. To address this challenge, we propose layer-by-layer metal FFF technology, followed by the debinding and sintering process, as an alternative to powder- and laser-based approaches. Furthermore, design for additive manufacturing (DfAM) principles are discussed to minimize the need for support generation and enable easy post-processing, thereby improving surface finish to meet high tolerances in fabricating components for aerospace and healthcare applications.

Singh, Abhishek [University of Michigan, Ann Arbor↗
NASA NTRS2014

Parallel Adaptive High-Order CFD Simulations Characterizing Cavity Acoustics for the Complete SOFIA Aircraft

This paper presents one-of-a-kind MPI-parallel computational fluid dynamics simulations for the Stratospheric Observatory for Infrared Astronomy (SOFIA). SOFIA is an airborne, 2.5-meter infrared telescope mounted in an open cavity in the aft of a Boeing 747SP. These simulations focus on how the unsteady flow field inside and over the cavity interferes with the optical path and mounting of the telescope. A temporally fourth-order Runge-Kutta, and spatially fifth-order WENO-5Z scheme was used to perform implicit large eddy simulations. An immersed boundary method provides automated gridding for complex geometries and natural coupling to a block-structured Cartesian adaptive mesh refinement framework. Strong scaling studies using NASA's Pleiades supercomputer with up to 32,000 cores and 4 billion cells shows excellent scaling. Dynamic load balancing based on execution time on individual AMR blocks addresses irregularities caused by the highly complex geometry. Limits to scaling beyond 32K cores are identified, and targeted code optimizations are discussed.

Acoustics↗
NASA NTRS1996

Progress on LES of Flow Past a Circular Cylinder

The objective of the present research is to assess the usefulness of large-eddy simulation (LES) methodology for flows in complex geometries. Flow past a circular cylinder has been calculated using a central-difference based solver, and the results have been compared to those obtained by a solver that employs higher-order upwind biased schemes (Beaudan & Moin, 1994). This comparison allows us to assess the suitability of these schemes for LES in complex geometry flows.

Mittal, R.↗
DOE OSTI2021

Complex Fault Geometry of the 2020 M ww 6.5 Monte Cristo Range, Nevada, Earthquake Sequence

On 15 May 2020 an M ww 6.5 earthquake occurred beneath the Monte Cristo Range in the Mina Deflection region of western Nevada. Rapid deployment of eight temporary seismic stations enabled detailed analysis of its productive and slowly decaying aftershock sequence (p=0.8), which included ~18,000 autodetected events in 3.5 months. Double-difference, waveform-based relative relocation of 16,714 earthquakes reveals a complex network of faults, many of which cross the inferred 35-km-long east–northeast-striking, left-lateral mainshock rupture. Seismicity aligns with left-lateral, right-lateral, and normal mechanism moment tensors of 128 of the largest earthquakes. The mainshock occurred near the middle of the aftershock zone at the intersection of two distinct zones of seismicity. In the western section, numerous subparallel, shallow, north-northeast-striking faults form a broad flower-structure-like fault mesh that coalesces at depth into a near-vertical, left-lateral fault. We infer the near-vertical fault to be a region of significant slip in the mainshock and an eastward extension of the left-lateral Candelaria fault. Near the mainshock hypocenter, seismicity occurs on a northeast-striking, west-dipping structure that extends north from the eastern Columbus Salt Marsh normal fault. Together, these two intersecting structures bound the Columbus Salt Marsh tectonic basin. East of this intersection and the mainshock hypocenter, seismicity occurs in a narrow, near-vertical, east-northeast-striking fault zone through to its eastern terminus. At the eastern end, the aftershock zone broadens and extends northwest toward the southern extension of the northwest-striking, right-lateral Petrified Springs fault system. The eastern section hosts significantly fewer aftershocks than the western section, but has more moment release. We infer that shallow aftershocks throughout the system highlight fault-fracture meshes that connect mapped fault systems at depth. Here, comparing earthquake data with surface ruptures and a simple geodetic fault model sheds light on the complexity of this recent M 6.5 Walker Lane earthquake.

58 GEOSCIENCES↗
DOE OSTI2018

Negative Additive Manufacturing of Complex Shaped Boron Carbides

Boron carbide (B 4 C) is one of the hardest materials in existence. However, this attractive property also limits its machineability into complex shapes for high wear, high hardness, and lightweight material applications such as armors. To overcome this challenge, negative additive manufacturing (AM) is employed to produce complex geometries of boron carbides at various length scales. Negative AM first involves gelcasting a suspension into a 3D-printed plastic mold. The mold is then dissolved away, leaving behind a green body as a negative copy. Resorcinol-formaldehyde (RF) is used as a novel gelling agent because unlike traditional hydrogels, there is little to no shrinkage, which allows for extremely complex molds to be used. Furthermore, this gelling agent can be pyrolyzed to leave behind ~50 wt% carbon, which is a highly effective sintering aid for B 4 C. Due to this highly homogenous distribution of in situ carbon within the B 4 C matrix, less than 2% porosity can be achieved after sintering. This protocol highlights in detail the methodology for creating near fully dense boron carbide parts with highly complex geometries.

36 MATERIALS SCIENCE↗
NASA NTRS2003

Numerical Simulation of Regional Circulation in the Monterey Bay Region

The objective of this study is to produce a high-resolution numerical model of Mon- terey Bay area in which the dynamics are determined by the complex geometry of the coastline, steep bathymetry, and the in uence of the water masses that constitute the CCS. Our goal is to simulate the regional-scale ocean response with realistic dynamics (annual cycle), forcing, and domain. In particular, we focus on non-hydrostatic e ects (by comparing the results of hydrostatic and non-hydrostatic models) and the role of complex geometry, i.e. the bay and submarine canyon, on the nearshore circulation. To the best of our knowledge, the current study is the rst to simulate the regional circulation in the vicinity of Monterey Bay using a non-hydrostatic model. Section 2 introduces the high resolution Monterey Bay area regional model (MBARM). Section 3 provides the results and veri cation with mooring and satellite data. Section 4 compares the results of hydrostatic and non-hydrostatic models.

Tseng, Y. H.↗
DOE OSTI2024

Enforcing global constraints for the dispersion closure problem: τ 2 -SIMPLE algorithm

Permeability and effective dispersion tensors are critical parameters to characterize flow and transport in porous media at the continuum scale. Homogenization theory defines a framework in which such effective properties are first computed from solving a closure problem in a repeating unit cell of the periodic microstructure and then used in a macroscopic formulation for efficient computation. The closure problem is formulated as a local boundary value problem subjected to global constraints, which guarantee the uniqueness of the solution and can be difficult to satisfy for complex geometries and at high flow conditions. These constraints also ensure that pore-scale pressure, velocity, and concentration fields can be accurately reconstructed from the closure variable. Building on a previous work, here we present a framework that allows to satisfy global constraints associated to both the permeability and the dispersion closure problems by introducing two artificial time scales. The algorithm, called τ 2 -SIMPLE, computes both permeability and effective dispersion given an arbitrarily complex geometry and flow condition. Furthermore, this algorithm is demonstrated to be accurate for both 2D and 3D geometries across varying flow conditions, and thus it can be used to quickly characterize effective properties from porous media images in many applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Packing fraction control during additive manufacturing of powder green bodies

The demand for high performance ceramic and metal components with complex geometry necessitates developments in powder handling and green body shaping. Here, in this study, vibrational powder deposition is adapted to three-dimensional printing with the ability to modulate packing fraction during printing. Boron carbide powder with sub-micron primary particle size is printed over a packing fraction range of 25.0 % to 46.5 % (below the poured density to above the tapped density). Solid metal powders are printed at 67.5 %, which enables novel freestanding, vertical-walled features without binder. These capabilities introduce opportunities for more complex geometries in high-performance materials, e.g. large-scale uniaxial hot-pressing of ceramic ballistic armor components with graded thickness and three-dimensional curved or stepped features (by compensating for displacement differences during compaction). This technology also enables multi-material patterning of additive manufacturing powder beds with reduced feedstock quantity requirements and wastage.

36 MATERIALS SCIENCE↗
DOE OSTI2017

A finite element approach to self-consistent field theory calculations of multiblock polymers

Self-consistent field theory (SCFT) has proven to be a powerful tool for modeling equilibrium microstructures of soft materials, particularly for multiblock polymers. A very successful approach to numerically solving the SCFT set of equations is based on using a spectral approach. While widely successful, this approach has limitations especially in the context of current technologically relevant applications. These limitations include non-trivial approaches for modeling complex geometries, difficulties in extending to non-periodic domains, as well as non-trivial extensions for spatial adaptivity. As a viable alternative to spectral schemes, we develop a finite element formulation of the SCFT paradigm for calculating equilibrium polymer morphologies. We discuss the formulation and address implementation challenges that ensure accuracy and efficiency. We explore higher order chain contour steppers that are efficiently implemented with Richardson Extrapolation. This approach is highly scalable and suitable for systems with arbitrary shapes. We show spatial and temporal convergence and illustrate scaling on up to 2048 cores. Finally, we illustrate confinement effects for selected complex geometries. This has implications for materials design for nanoscale applications where dimensions are such that equilibrium morphologies dramatically differ from the bulk phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Non-conformal interface-cohesive modeling with the shifted boundary method

The accurate simulation of boundary- and interface-dominated problems on complex geometries remains challenging when boundary- or interface-fitted meshes are difficult to generate, particularly for curved boundaries, polycrystalline microstructures, and dense interface networks. The Shifted Boundary Method (SBM) alleviates this meshing burden by shifting the enforcement of boundary conditions from the true boundary to a nearby surrogate boundary and recovering the effect of the true boundary through geometric correction terms, thereby enabling standard finite element spaces on non-boundary-fitted meshes. In this report, we develop a general shiftedboundary and shifted-interface framework within the open-source MOOSE framework. We first present a general SBM implementation for complex geometries on non-boundary-fitted meshes. We then adopt the Shifted Interface Method (SIM) for internal interfaces and develop a unified shifted-interface treatment in which the interface law is enforced on a surrogate interface and the effect of the true interface is recovered through shifted jumps, fluxes, and tractions. This perspective brings scalar thermal-contact and vector-valued cohesive-zone mechanics into a single framework, the latter realized as the Shifted Cohesive Zone Method (SCZM) and coupled with history-dependent constitutive models from NEML2. We further extend the MOOSE mesh infrastructure to support cohesive-zone calculations on distributed meshes. The framework is verified and demonstrated through three progressive studies: Poisson’s equation on a smoothed starshaped domain, a manufactured thermal-contact problem on a non-interface-fitted mesh, and a two-dimensional polycrystalline representative volume element combining crystal plasticity with cohesive grain-boundary interfaces. Across these studies, the shifted formulations reproduce boundary- and interface-fitted reference solutions with high fidelity, indicating that the proposed framework provides an accurate and efficient route to boundary- and interface-dominated simulations on arbitrary geometries without requiring fitted meshes.

Yang, Cheng-Hau↗