Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Complex geometry”

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

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

At least 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Tensor Network Space-Time Spectral Collocation Method for Time-Dependent Convection-Diffusion-Reaction Equations

Emerging tensor network techniques for solutions of partial differential equations (PDEs), known for their ability to break the curse of dimensionality, deliver new mathematical methods for ultra-fast numerical solutions of high-dimensional problems. Here, we introduce a Tensor Train (TT) Chebyshev spectral collocation method, in both space and time, for the solution of the time-dependent convection-diffusion-reaction (CDR) equation with inhomogeneous boundary conditions, in Cartesian geometry. Previous methods for numerical solution of time-dependent PDEs often used finite difference for time, and a spectral scheme for the spatial dimensions, which led to a slow linear convergence. Spectral collocation space-time methods show exponential convergence; however, for realistic problems they need to solve large four-dimensional systems. We overcome this difficulty by using a TT approach, as its complexity only grows linearly with the number of dimensions. We show that our TT space-time Chebyshev spectral collocation method converges exponentially, when the solution of the CDR is smooth, and demonstrate that it leads to a very high compression of linear operators from terabytes to kilobytes in TT-format, and a speedup of tens of thousands of times when compared to a full-grid space-time spectral method. These advantages allow us to obtain the solutions at much higher resolutions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Predicting Wind Loading and Instability in Solar Tracking PV Arrays

Wind loading and the fluctuating pressure loads it creates on PV panel surfaces are associated with multiple degradation mechanisms and failures. Modest wind speeds create reversing loads that can initiate cell cracks and weather cracked cells. Stronger wind speeds and extreme weather events can lead to larger scale forces and the aerodynamic instability known as torsional galloping. All these effects are dependent on the complex coupling between wind speed, panel orientation, and a myriad of other hardware and site-specific factors. In this work, we present the latest developments from our work to build an open-source, high-performance computing (HPC) fluid dynamics solver to predict and mitigate these effects. This simulation package allows users to easily specify different array layouts, solar-tracking angles, panel geometries, and weather conditions before automatically generating a refined computational mesh and solving for the unsteady loading on each panel surface. Small domains (e.g., a single panel row in isolation) can be solved on a modern laptop, while larger domains or very high-fidelity studies can be solved on distributed or HPC resources with minimal modifications to the underlying problem specification. We present preliminary case studies obtained using this simulation package and highlight how increased wind speeds combined with sub-optimal tracking angles can exacerbate degradation drivers.

aerodynamics↗
DOE OSTI2023

Sensitivity analysis of a layered piezoelectric system using ZFEM

The complex variable finite element method (ZFEM) is a numerical technique which aims to find the partial derivatives of the independent variables with respect to variation in dependent parameters declared in the physics. This is done by combining the complex Taylor series expansion within the weak formulation of the governing equation in a coupled system of linear equations forming a complex valued block matrix given by the Cauchy–Riemann matrix representation. In this work, two-dimensional linear first-order elements have been implemented in ZFEM to predict the design derivatives of the mechanical displacement field and the voltage potential field for a layered piezoelectric system in a steady-state study with Dirichlet boundary condition applied at the top and bottom edges of the geometry. This approach allows the standard FEM solution to quantify the sensitivity of the mechanical displacement and voltage potential fields with respect to small variations in the material properties through the information obtained from the computation of the derivatives. The domain is formed by a layered body with PZT-4 and PZT-5 stacked together. For result verification, the numerical solution obtained with ZFEM was compared to results from a commercial FEM package and the solution from the imaginary part was compared to the exact solution of a well-known benchmark problem. In conclusion, comparison of the results showed good agreement for both the real and imaginary parts of the solution and the largest sensitivities were found in PZT-5 specifically in C 13 , C 33 , and ε 33 .

42 ENGINEERING↗
DOE OSTI2020

Dust and gas absorption in the high mass X-ray binary IGR J16318-4848

With an absorption column density on the order of 10 24 cm -2 , IGR J16318-4848 is one of the most extreme cases of a highly obscured high mass X-ray binary. In addition to the overall continuum absorption, the source spectrum exhibits a strong iron and nickel fluorescence line complex at 6.4 keV. Previous empirical modeling of these features and comparison with radiative transfer simulations raised questions about the structure and covering fraction of the absorber and the profile of the fluorescence lines. We aim at a self-consistent description of the continuum absorption, the absorption edges, and the fluorescence lines to constrain the properties of the absorbing material, such as ionization structure and geometry. We further investigate the effects of dust absorption on the observed spectra and the possibility of fluorescence emission from dust grains. We used XMM-Newton and NuSTAR spectra to first empirically constrain the incident continuum and fluorescence lines. Next we used XSTAR to construct a customized photoionization model where we vary the ionization parameter, column density, and covering fraction. In the third step we modeled the absorption and fluorescence in a dusty olivine absorber and employed both a simple analytical model for the fluorescence line emission and a Monte Carlo simulation of radiative transfer that generates line fluxes, which are very close to the observational data. Our empirical spectral modeling is in agreement with previous works. Our second model, the single gas absorber does not describe the observational data. In particular, irrespective of the ionization state or column density of the absorber, a much higher covering fraction than previously estimated is needed to produce the strong fluorescence lines and the large continuum absorption. A dusty, spherical absorber (modeled as consisting of olivine dust, although the nature of dust cannot be constrained) is able to produce the observed continuum absorption and edges. A dense, dusty absorber in the direct vicinity of the source consisting of dust offers a consistent description of both the strong continuum absorption and the strong emission features in the X-ray spectrum of IGR J16318-4848. In particular, for low optical depth of individual grains, which is the case for typical volume densities and grain size distribution models, the dust will contribute significantly to the fluorescence emission.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

M-STAR: Magnetism second target advanced reflectometer at the Spallation Neutron Source

M-STAR is a next generation polarized neutron reflectometer with advanced capabilities. A new focusing guide concept is optimized for samples with dimensions down to a millimeter range. A proposed hybrid pulse-skipping chopper will enable experiments at constant geometry at one incident angle in a broad range of wavevector transfer Q up to 0.3 A –1 for specular, off-specular, and GISANS measurements. M-STAR will empower nanoscience and spintronics studies routinely on small samples (~2 × 2 mm 2 ) and of atomic-scale thickness using versatile experimental conditions of magnetic and/or electric fields, light, and temperature applied in situ to novel complex device-like nanosystems with multiple buried interfaces. M-STAR will enable improved grazing incidence diffraction measurements, as a surface-sensitive depth-resolved probe of, e.g., the out-of-plane component of atomic magnetic moments in ferromagnetic, antiferromagnetic, and more complex structures as well as in-plane atomic-scale structures inaccessible with contemporary diffractometry and reflectometry. New horizons will be opened by the development of an option to probe near-surface dynamics with inelastic grazing incidence scattering in the time-of-flight mode. These novel options in combination with ideally matched parameters of the second target station will place M-STAR in the world’s leading position for high resolution polarized reflectometry.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Dataset describing two reference models for full-spectral lighting and daylight simulations together with implementations for two software systems

A dataset of two spectral lighting simulation reference models - one office and one factory hall - is presented. It aims to demonstrate and support full-spectral daylight and electric lighting simulations and facilitate evaluation of non-visual effects of light. The dataset includes Rhino CAD geometry, comprehensive spectral material and light source data and window system BSDF data. Example implementations in the two software tools, Radiance and OWL, enable reproducible workflows and support adoption in other software. The dataset is openly available on Zenodo. The office model reproduces Room 518 at the University of Innsbruck, including a west-facing façade and interior furnishings. The factory hall model follows the proposed geometry in the European standard 15193 for building energy performance. Interior reflectances in the office were measured in-situ using a handheld spectrometer. Exterior spectra and factory hall materials matching specified reflectances were obtained from an online spectral materials database. Glazing transmittance was derived from IGDB data using LBNL Optics/WINDOW. BSDFs for venetian blinds at various tilt angles, and for a diffusing pane adapted from the Complex Glazing Database, were generated in WINDOW. Luminaires in both models are specified with photometric files (Eulumdat/IES) and lamp spectra (Fluorescent 840, 4000 K LED). The provided example implementations (Radiance, OWL) include prepared input data and scripts to run first spectral simulations; example results are also included. The dataset is prepared to support reuse by researchers, designers and software developers for method validation, software engineering and comparison, and development of spectral metrics and controls.

Geisler-Moroder, David↗
DOE OSTI2021

Fracture Description of the HFTS-2 Slant Core, Delaware Basin, West Texas

Previous work in the Midland Basin Hydraulic Fracture Test Site project (HFTS1) has shown the value of core recovered from a stimulated volume in providing information on the geometry and extent of hydraulic fractures. This project (HFTS2) is in the Wolfcamp Formation of the Delaware Basin. A slant core through the stimulated volume was acquired above and below a stimulated well and the results of the fracture description are presented here. Here, nine, 3-inch-diameter whole cores from the slant well were examined for fractures prior to slabbing, together with a CT scan of the core. We used criteria developed during the HFTS1 project to identify hydraulic, natural, drilling-induced and core-handling fractures. Material collected from a few fracture faces was cleaned, sieved, and examined under a microscope to look for proppant. A total of 1261 fractures, including 500 hydraulic fractures, were described in 948 feet of core. Hydraulic fractures are not evenly distributed through the cores; they tend to occur in clusters. Many are remarkably planar, smooth, and featureless, while others have twist hackles, steps, plumose, and kinks. Natural, subvertical fractures, mostly less than 1 mm in width, sealed with blocky calcite, are oriented NE-SW (Set 1) and WNW-ESE (Set 2). Two Set 1 fractures have drusy cement on parted faces and were possibly open prior to drilling. A few subvertical sealed fractures with different orientations are also present. Bed-parallel natural fractures, most having fibrous cement (beef), and some showing top-to-north bed-parallel shear are present. Faults with normal/oblique movement were found in one part of the core. Drilling induced and core-handling breaks were identified in all cores, and were also noted to have reactivated some natural fractures. The results from this work will allow hydraulic fracture distribution relative to the stimulated well to be assessed. The number of fractures is generally greater than the number of perforations, but these can be compared for different stages and perforation clusters. Surface features provide information about the direction of hydraulic fracture propagation and about segmentation and bifurcation that might increase the number of fracture strands. Reactivation of natural fractures, including bedding-parallel fractures, occurs although the mechanism of parting is mostly unknown. The wider HFTS2 project will make use of this fracture description to improve hydraulic fracture modeling (geomechanics and engineering), to help explain pressure depletion and fiber optics observations (reservoir depletion and frac geometry), to help verify indirect fracture diagnostic techniques such as microseismic monitoring (geophysics), and to help calibrate horizontal well image logs.

58 GEOSCIENCES↗
DOE OSTI2024

Point containment algorithms for constructive solid geometry with unbounded primitives

Here, we present several algorithms for evaluating point containment in constructive solid geometry (CSG) trees with unbounded primitives. Three algorithms are presented based on postfix, prefix, and infix notations of the CSG binary expression tree. We show that prefix and infix notations enable short-circuiting logic, which reduces the number of primitives that must be checked during point containment. To evaluate the performance of the algorithms, each algorithm was implemented in the OpenMC Monte Carlo particle transport code, which relies on CSG to represent solid bodies through which subatomic particles travel. Two sets of tests were carried out. First, the execution time to generate a rasterized image of a 2D slice of three CSG models of varying complexity was measured. Use of both prefix and infix notations offered significant speedup over the postfix notation that has traditionally been used in particle transport codes, with infix resulting in a 6 x reduction in execution time relative to postfix for a model of a tokamak fusion device. We then measured the execution time of neutron transport simulations of the same three models using each of the algorithms. The results and performance improvements reveal the same trends as for the rasterization test, with a 5.52 x overall speedup using the infix notation relative to the original postfix notation in OpenMC for the tokamak model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Hydrogen-Bonded Organic Frameworks: A Rising Class of Porous Molecular Materials

Hydrogen-bonded organic frameworks (HOFs) are a class of porous molecular materials that rely on the assembly of organic building blocks by means of hydrogen-bonding interactions to form two-dimensional (2D) and three-dimensional (3D) crystalline networks. The reversible nature of the hydrogen-bond formation endows HOFs with the attributes of solution processability and simple regeneration. High-quality single crystals of HOFs can be grown easily for unambiguous superstructure determination by single-crystal X-ray diffraction, which is crucial for the elucidation of superstructure–property relationships. During the past decade, considerable progress has been achieved in realizing stable HOFs with permanent porosities by focusing on the design of molecular building blocks in order to introduce rigidity, auxiliary [π···π] interactions, and interpenetration of their frameworks to sustain the extended networks. The applications of HOFs are far-reaching, spanning catalysis, energy, and biomedical products as well as the storage and separation of fine chemicals. In this paper, we, first of all, provide an overview of the chronological development of HOFs, starting from the seminal work by Marsh and Duchamp in 1969 on the crystal superstructure of the hydrogen-bonded networks of trimesic acid. We identify the development of novel hydrogen-bonding motifs such as diaminotriazine (DTA), the introduction of the concept of molecular tectonics, and the establishment of permanent porosity in HOFs as being some of the milestones, which incentivized the current burgeoning research endeavors on developing HOFs as multifunctional materials. This Account is focused primarily on surveying the strategies for constructing porous 3D HOFs based on organic building blocks with peripheral carboxyl groups. These strategies are presented in the following categories: (1) the polycatenation of 2D networks by trigonal building blocks to form global 3D frameworks, (2) the utilization of building blocks with 3D geometries—tetrahedral and trigonal prismatic—that are predisposed to form 3D networks, and (3) the docking by shape-fitting of geometrically labile building blocks. We emphasize how the molecular geometry of the building blocks plays an important role in modulating the superstructures of extended frameworks so as to address specific applications. Recognizing that the in silico design of HOFs is the ultimate goal of researchers in this field, we also discuss the recent advances in superstructure prediction that lead to the formation of porous supramolecular crystals and assess the complications in implementing computational methods for HOFs with complex superstructures. We hope this Account will inspire the development of new supramolecular designs and creative approaches to crystal engineering that aid and abet the assembly of multifunctional HOFs with customizable properties.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Importance of van der Waals interactions for ab initio studies of topological insulators

We investigate the lattice and electronic structures of the bulk and surface of the prototypical layered topological insulators Bi 2 Se 3 and Bi 2 Te 3 using ab initio density functional methods, and systematically compare the results of different methods of including van der Waals (vdW) interactions. We show that the methods utilizing semi-empirical energy corrections yield accurate descriptions of these materials, with the most precise results obtained by properly accounting for the long-range tail of the vdW interactions. The bulk lattice constants, distances between quintuple layers and the Dirac velocity of the topological surface states (TSS) are all in excellent agreement with experiment. In Bi 2 Te 3 , hexagonal warping of the energy dispersion leads to complex spin textures of the TSS at moderate energies, while in Bi 2 Se 3 these states remain almost perfectly helical away from the Dirac point, showing appreciable signs of hexagonal warping at much higher energies, above the minimum of the bulk conduction band. Furthermore, our results establish a framework for unified and systematic self-consistent first principles calculations of topological insulators in bulk, slab and interface geometries, and provides the necessary first step toward ab initio modeling of topological heterostructures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Plane-Strain Tension Testing of ATF FeCrAl Cladding (Progress Report)

Performance of advanced accident-tolerant fuel (ATF) claddings must be assessed under all relevant scenarios, including normal and off-normal operating conditions, as well as during a reactivity-initiated accident (RIA). During RIA, the fuel pellet may expand and cause significant loads on cladding, a condition known as pellet cladding mechanical interaction (PCMI), which may even lead to cladding failure. This report describes the fiscal year (FY) 2020 research activities and results from the development of separate-effects testing for ATF cladding candidate iron-chromium-aluminum (FeCrAl) tube specimens under loading conditions expected during RIA/PCMI events. To assess material performance under PCMI, a high strain rate burst test, often termed the modified burst test (MBT), is commonly employed. However, execution of such a test in a hot cell requires significant research. To overcome the complexities of performing MBTs in a hot-cell environment on irradiated samples, a plane-strain tension (PST) test was proposed by several authors as an alternative approach for PCMI testing. A new fixture was developed so that PST specimens could be tested using a more accommodating design for hot-cell manipulators. This report presents the PST fixture design, evaluates the possibility of modifying the specimen geometry, and discusses experimental results with unirradiated Al 6061 and FeCrAl tubes at various lengths. Furthermore, this report discusses the relevance and compatibility of the PST results to RIA scenarios and compares the PST and MBT results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

A Review of the Research and Development of Brayton Cycle Technology in Nuclear Power Applications with a Focus on Compressor Technology

This study reviews the integration of Brayton Cycle (BC) systems in nuclear power generation, emphasizing their potential to enhance thermal efficiency and operational flexibility over traditional Rankine Cycle (RC) systems. Key working fluids, such as helium (He), supercritical carbon dioxide (sCO 2 ), nitrogen (N 2 ), and air, are evaluated for their performance, efficiency, and compatibility with nuclear systems. He is recognized for its high thermal conductivity and inertness at elevated temperatures, while sCO 2 demonstrates advantages in compactness and efficiency in midrange temperatures. This article also highlights the importance of compressor designs in optimizing BC performance and reviews, available compressor technologies. Axial and centrifugal compressor designs enable efficient gas compression while managing the thermal and mechanical stresses associated with high-pressure operations in nuclear systems. Combined with variable geometry components and advanced materials, these technologies address the challenges posed by varying load conditions. Despite the promising features of BC systems, several challenges persist, including high leakage rates and material degradation under extreme conditions, which necessitate robust sealing technologies and thorough testing. The insights gained from operational experiences at facilities, such as the Oberhausen II plant and the High-Temperature He Test Facility (HHV), underscore the complexities involved in designing high-temperature gas turbines for nuclear applications. This review concludes that as the nuclear industry evolves, BC systems hold significant promise for contributing to a sustainable energy future, particularly in the context of small modular reactors (SMRs) and microreactors. Further exploration of combined cycle configurations that combine BCs with RCs may enhance overall efficiency and flexibility in power generation.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

M3 milestone: Advanced contact 2021

This report describes advanced contact approaches and implementations in FY21 for use in emerging reactors with geometries that differ from LWR/EBR-II cylindrical fuel. Contact remains a challenging and crucially important component of our modeling toolkit. While significant improvements have been made through mortar and automatic differentiation, we explored advanced approaches (i.e., dual mortar) that yields better conditioned systems and will improve the robustness and performance of BISON contact problems. This effort is relevant to emerging reactor types, such as micro reactors and the corresponding complex and numerous contact surfaces. In this report, improvements of the mortar-based contact formulation are described. Benchmark problems and BISON assessment cases are included to demonstrate the improvements and advantages that are brought by the recent capability. Remaining issues are discussed and corresponding development plans are outlined.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Engineering metal–metal oxide surfaces for high-performance oxygen reduction on Ag–Mn electrocatalysts

Understanding fundamental material–property relationships in mixed-element catalyst systems is crucial to advancing the viability of renewable electrochemical energy technologies, an important part of creating a more sustainable future. Herein, we report our insight on the nature and dynamics of highly active silver–manganese oxide (Ag–MnOx) catalyst surfaces for the oxygen reduction reaction (ORR) via a combined experimental–theoretical approach. Experimentally, we synthesize well-mixed Ag–Mn co-deposited thin films that are measurably flat and smooth, despite Mn surface migration and oxidation upon air exposure and electrochemical measurements. Cyclic voltammetry in 0.1 M KOH demonstrates up to 10-fold specific activity enhancements over pure Ag at 0.8 V vs. RHE for Ag-rich films (70–95% Ag in bulk). To further understand the Ag–Mn system, separate samples were synthesized with small amounts of Mn sequentially deposited onto the surface of a pure Ag thin film (Mn@Ag), ranging from partial to full surface coverage (down to 0.3 nm Mn $cm^{-2}_{geo}$ ~ 0.2 μg Mn $cm^{-2}_{geo}$). These sequentially deposited Mn@Ag films show analogous performance to their co-deposited counterparts indicating similar enhanced active sites. With density functional theory (DFT), we calculate that this enhancement arises from the tuned d-band of these material surfaces owing to the optimal hybridization of the electronic structures in specific Ag and MnO x geometries. Together, electrochemical measurements, DFT calculations, X-ray absorption spectroscopy, and valence-band X-ray photoelectron spectroscopy suggest synergistic electronic interactions between Ag and MnOx yield enhanced oxygen adsorption, and thus ORR activity, with DFT highlighting the Ag–MnO x interface sites as the most enhanced. This work demonstrates how combined experimental–theoretical methods can help design electrocatalysts with enhanced electrocatalytic properties and understand the nature of complex mixed metal–metal oxide surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Heat Exchanger Design for Waste Heat Recovery using Additive Manufacturing

With anticipated increases in global energy use and their associated emissions, there are many efforts to improve the emissions and efficiency of many industrial processes through waste heat recovery (WHR). To maximize improvements, the heat exchangers in WHR power cycles must be designed and optimized for a given heat source. Meanwhile, additive manufacturing (AM) is a single part manufacturing technique that provides the opportunity to develop unique and complex heat exchanger designs. Consequently, this work centers on the modeling, design, optimization, and testing of AM shell and tube heat exchangers for WHR. Specifically, a modified one-dimensional (1+1-D) and three-dimensional (3-D) heat exchanger models were developed and validated, with the 1+1-D model serving as the basis for design optimization. After modifications to the design to allow for AM, heat exchanger geometry optimizations were completed, targeting increased effectiveness across a range of physical parameters. Of the designs considered, nine were additively manufactured using Ti5553, and then tested using an external heat source and working fluid. Experimentally, the nineteen-tube, zero-baffle heat exchanger displayed the best overall performance, with the effects of dimensional accuracy, surface roughness, and baffles in the other designs also indicated. Next, the experimental results were used as inputs for the 1+1-D model. Although the original model was able to predict fluid outlet temperatures within 5% differences, larger variations in the heat exchanger performance parameters, temperature differences, and pressure drop necessitated optimizing the pressure drop and heat transfer calculations within the model. Thus, the resulting improved model was able to predict the heat transfer rates, heat source temperature difference, and effectiveness within 15% differences, while showing improvements in the pressure drops and working fluid temperature differences. The improved model serves as the foundation for future work to improve the design and performance of AM heat exchangers.

Spickler, Bailey Ann↗
DOE OSTI2021

The response field and the saddle points of quantum mechanical path integrals

Highlights: • Moyal quantum mechanics and Marinov’s path integral. • Classical and semiclassical limits of Marinov’s path integral. • Oscillating functional integrals. • Instantons of the Marinov’s path integral. In quantum statistical mechanics, Moyal’s equation governs the time evolution of Wigner functions and of more general Weyl symbols that represent the density matrix of arbitrary mixed states. A formal solution to Moyal’s equation is given by Marinov’s path integral. In this paper we demonstrate that this path integral can be regarded as the natural link between several conceptual, geometric, and dynamical issues in quantum mechanics. A unifying perspective is achieved by highlighting the pivotal role which the response field, one of the integration variables in Marinov’s integral, plays for pure states even. The discussion focuses on how the integral’s semiclassical approximation relates to its strictly classical limit; unlike for Feynman type path integrals, the latter is well defined in the Marinov case. The topics covered include a random force representation of Marinov’s integral based upon the concept of “Airy averaging”, a related discussion of positivity-violating Wigner functions describing tunneling processes, and the role of the response field in maintaining quantum coherence and enabling interference phenomena. The double slit experiment for electrons and the Bohm–Aharonov effect are analyzed as illustrative examples. Furthermore, a surprising relationship between the instantons of the Marinov path integral over an analytically continued (“Wick rotated”) response field, and the complex instantons of Feynman-type integrals is found. The latter play a prominent role in recent work towards a Picard–Lefschetz theory applicable to oscillatory path integrals and the resurgence program.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Enriched immersed finite element and isogeometric analysis: algorithms and data structures

Immersed finite element methods provide a convenient analysis framework for problems involving geometrically complex domains, such as those found in topology optimization and microstructures for engineered materials. However, their implementation remains a major challenge due to, among other things, the need to apply nontrivial stabilization schemes and generate custom quadrature rules. This article introduces the robust and computationally efficient algorithms and data structures comprising an immersed finite element preprocessing framework. The input to the preprocessor consists of a background mesh and one or more geometries defined on its domain. The output is structured into groups of elements with custom quadrature rules formatted such that common finite element assembly routines may be used without or with only minimal modifications. The key to the preprocessing framework is the construction of material topology information, concurrently with the generation of a quadrature rule, which is then used to perform enrichment and generate stabilization rules. While the algorithmic framework applies to a wide range of immersed finite element methods using different types of meshes, integration, and stabilization schemes, the preprocessor is presented within the context of the extended isogeometric analysis. This method utilizes a structured B-spline mesh, a generalized Heaviside enrichment strategy considering the material layout within individual basis functions’ supports, and face-oriented ghost stabilization. Using a set of examples, the effectiveness of the enrichment and stabilization strategies is demonstrated alongside the preprocessor’s robustness in geometric edge cases. Additionally, the performance and parallel scalability of the implementation are evaluated.

Computer implementation↗
DOE OSTI2020

Current challenges and potential directions towards precision microscale additive manufacturing – Part III: Energy induced deposition and hybrid electrochemical processes

The Part III of the four-part series of articles discusses the challenges and opportunities in microscale additive manufacturing processes, specifically focusing on energy-induced deposition and electrochemical processes. Compared to the direct ink write (DIW) and laser-based processes, the energy-induced deposition methods can fabricate high-resolution, high aspect ratio and complex parts, while the hybrid electrochemical process can be used to fabricate complex parts using a wide range of conductive and photoactive materials. However, the volumetric throughput of these processes is lower than their DIW and laser-based counterparts. The processes that have been explored in this process are Focused-ion Beam Induced Deposition (FIBID), Laser Chemical Vapor Deposition (LCVD), Menicus-confined Electrodeposition (MCED) and Laser-Enabled Electrochemical Printing (LECP). The range of processable materials, feature-size resolution, geometry and volumetric throughput are used as factors to evaluate the current state-of-the-art for these processes. Finally, novel approaches have been proposed in the article to address these challenges associated with microscale AM processes.

42 ENGINEERING↗