Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mesh network”

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 127 records · Page 7

Mesh-based super-resolution of fluid flows with multiscale graph neural networks

A graph neural network (GNN) approach is introduced in this work which enables mesh-based three-dimensional super-resolution of fluid flows. In this framework, the GNN is designed to operate not on the full mesh-based field at once, but on localized meshes of elements (or cells) directly. To facilitate mesh-based GNN representations in a manner similar to spectral (or finite) element discretizations, a baseline GNN layer (termed a message passing layer, which updates local node properties) is modified to account for synchronization of coincident graph nodes, rendering compatibility with commonly used element-based mesh connectivities. Furthermore, the architecture is multiscale in nature, and is comprised of a combination of coarse-scale and fine-scale message passing layer sequences (termed processors) separated by a graph unpooling layer. The coarse-scale processor embeds a query element (alongside a set number of neighboring coarse elements) into a single latent graph representation using coarse-scale synchronized message passing over the element neighborhood, and the fine-scale processor leverages additional message passing operations on this latent graph to correct for interpolation errors. Demonstration studies are performed using hexahedral mesh-based data from Taylor–Green Vortex and backward-facing step flow simulations at Reynolds numbers of 1600 and 3200. Through analysis of both global and local errors, the results ultimately show how the GNN is able to produce accurate super-resolved fields compared to targets in both coarse-scale and multiscale model configurations. Reconstruction errors for fixed architectures were found to increase in proportion to the Reynolds number. Geometry extrapolation studies on a separate cavity flow configuration show promising cross-mesh capabilities of the super-resolution strategy.

Backward-facing step↗

Highly swollen ROMP-based gels

In this study, the leakage of organic solvent during the transportation and storage is one of the major threats to the environment. Networks/gels acting as the absorbing materials for organic solvents are the most promising solution for the organic solvent pollution. Here, we report the synthesis of highly swollen networks/gels in high yield (≥95%) via ring opening metathesis polymerization (ROMP) using Grubbs’ first generation catalyst/initiator (G1), a crosslinker (1,1'-(methylenedi-4,1-phenylene)bis-exo-7-oxanorbornenecarboximide) and monofunctional monomer/diluent (exo-N-phenyl-7-oxanorbornenecarboximide). The mesh size of the network was varied by keeping the concentration of monomer constant while changing the concentrations of crosslinker and G1 (initiator). The swelling capacity of the network in N,N-dimethylformamide (DMF) solvent increased from 41 to 53 and to 61 with the increase of theoretical mesh size (molar ratio of monomer to crosslinker) from 20 to 40 and to 60, respectively. By further increasing the mesh size from 60 to 80 there was no significant change in swelling capacities (63). Thus, the gel with the mesh size of 60 was selected for swelling capacity analysis in other polar organic solvents such as dimethyl sulfoxide (DMSO), dichloromethane (DCM), chloroform and lithium-ion battery electrolyte solution (1.0 M LiPF 6 in ethylene carbonate (EC)/diethyl carbonate (DEC) = 50/50 (v/v)). They showed the decreased swelling capacity value from 61 to 57, 40, 31 and 5, respectively. For a controlled mesh size ROMP-based gel, we observed the highest swelling capacity values to date.

36 MATERIALS SCIENCE↗

Toward more mechanistic representations of biogeochemical processes in river networks: Implementation and demonstration of a multiscale model

A recently introduced multiscale model for representing the effects of small-scale hyporheic-zone biogeochemical processes is extended from the reach scale to river network scales. The model uses advection-dispersion-reaction equations for the channel network and one-dimensional advection-reaction subgrid models for the hyporheic zone. We summarize implementation in the integrated surface/subsurface hydrology modeling system Amanzi-ATS. The extension uses topologically defined meshes to represent stream/river networks and associates a hyporheic subgrid model with each channel grid cell. Biogeochemical reaction modeling capability residing in community software is accessed through an application programming interface. The implementation is verified against independent numerical solutions on a single reach. Mesh convergence studies show that commonly used semi-distributed representations can introduce significant spatial discretization error. Denitrification of farm runoff in a subbasin of the Portage River Basin in Ohio, USA is used to demonstrate the general-purpose reactive transport capability.

54 ENVIRONMENTAL SCIENCES↗

dfnWorks: A discrete fracture network framework for modeling subsurface flow and transport

DFNWORKS is a parallelized computational suite to generate three-dimensional discrete fracture networks (DFN) and simulate flow and transport. Developed at Los Alamos National Laboratory over the past five years, it has been used to study flow and transport in fractured media at scales ranging from millimeters to kilometers. The networks are created and meshed using DFNGEN, which combines FRAM (the feature rejection algorithm for meshing) methodology to stochastically generate three-dimensional DFNs with the LaGriT meshing toolbox to create a high-quality computational mesh representation. The representation produces a conforming Delaunay triangulation suitable for high performance computing finite volume solvers in an intrinsically parallel fashion. Flow through the network is simulated in dfnFlow, which utilizes the massively parallel subsurface flow and reactive transport finite volume code PFLOTRAN. A Lagrangian approach to simulating transport through the DFN is adopted within DFNTRANS to determine pathlines and solute transport through the DFN. Example applications of this suite in the areas of nuclear waste repository science, hydraulic fracturing and CO 2 sequestration are also included.

42 ENGINEERING↗

LES of Vaporizing Gasoline Sprays Considering Multi-Injection Averaging and Grid-Convergent Mesh Resolution

A state-of-the-art spray modeling methodology, recently presented by Senecal et al. [ , , ], is applied to Large Eddy Simulations (LES) of vaporizing gasoline sprays. Simulations of non-combusting Spray G (gasoline fuel) from the Engine Combustion Network are performed. Adaptive mesh refinement (AMR) with cell sizes from 0.09 mm to 0.5 mm are utilized to further demonstrate grid convergence of the dynamic structure LES model for the gasoline sprays. Grid settings are recommended to optimize the accuracy/runtime tradeoff for LES-based spray simulations at different injection pressure conditions typically encountered in gasoline direct injection (GDI) applications. The influence of LES sub-grid scale (SGS) models is explored by comparing the results from dynamic structure and Smagorinsky based models against simulations without any SGS model. Twenty different realizations are simulated by changing the random number seed used in the spray sub-models. It is shown that for global quantities such as spray penetration, comparing a single LES simulation to experimental data is reasonable. Through a detailed analysis using the relevance index (RI) criteria, recommendations are made regarding the minimum number of LES realizations required for accurate prediction of the gasoline sprays.

Som, Sibendu↗

Wireless Sensor Networks - Node Localization for Various Industry Problems

Fast, effective monitoring following airborne releases of toxic substances is critical to mitigate risks to threatened population areas. Wireless sensor nodes at fixed predetermined locations may monitor such airborne releases and provide early warnings to the public. A challenging algorithmic problem is determining the locations to place these sensor nodes while meeting several criteria: 1) provide complete coverage of the domain, and 2) create a topology with problem dependent node densities, while 3) minimizing the number of sensor nodes. This manuscript presents a novel approach to determining optimal sensor placement, Advancing Front mEsh generation with Constrained dElaunay Triangulation and Smoothing (AFECETS) that addresses these criteria. A unique aspect of AFECETS is the ability to determine wireless sensor node locations for areas of high interest (hospitals, schools, high population density areas) that require higher density of nodes for monitoring environmental conditions, a feature that is difficult to find in other research work. The AFECETS algorithm was tested on several arbitrary shaped domains. AFECETS simulation results show that the algorithm 1) provides significant reduction in the number of nodes, in some cases over 40%, compared to an advancing front mesh generation algorithm, 2) maintains and improves optimal spacing between nodes, and 3) produces simulation run times suitable for real-time applications.

42 ENGINEERING↗

Topological Relationship-Based Flow Direction Modeling: Stream Burning and Depression Filling

Flow direction modeling consists of (a) an accurate representation of the river network and (b) digital elevation model (DEM) processing to preserve characteristics with hydrological significance. In part 1 of our study, we presented a mesh-independent approach to representing river networks on different types of meshes. This follow-up part 2 study presents a novel DEM processing approach for flow direction modeling. This approach consists of (a) a topological relationship-based hybrid breaching-filling method to conduct stream burning for the river network and (b) a modified depression removal method for rivers and hillslopes. Our methods reduce modifications to surface elevations and provide a robust two-step procedure to remove local depressions in DEM. They are mesh-independent and can be applied to both structured and unstructured meshes. We applied our new methods with different model configurations to the Susquehanna River Basin. The results show that topological relationship-based stream burning, and depression-filling methods can reproduce the correct river networks, providing high-quality flow direction and other characteristics for hydrologic and Earth system models.

54 ENVIRONMENTAL SCIENCES↗

Predicting Two-Dimensional Airfoil Performance Using Graph Neural Networks

Computer simulations require the use of meshes to simulate geometries. These meshes capture important geometric features of the design and can be used in machine learning modeling. This report explores the use of graph neural networks (GNNs) to learn features from two-dimensional (2D) airfoil designs represented as a set of nodes connected using edges. This type of network is common in aerospace applications: most geometries are represented as a mesh in order to perform analysis. The objective of this work is to use GNNs to predict the performance of 2D airfoils generated using the program XFOIL. The predicted performance parameters include bulk quantities such as coefficients of lift (C L ), drag (C d , C dp ), moment (C m ), and node-specific quantities such as coefficient of pressure (C p ). In this report, a spline convolutional graph-based neural network is compared with deep learning neural networks to predict both bulk and node-specific quantities. The findings indicate the GNNs are able to predict bulk quantities quite well; however, when the number of outputs is increased, the deep neural network (DNN) proves to be better in its prediction capability. Two different normalization strategies were compared in the training of both GNNs and DNNs: minmax and standard deviation. In both types of networks, standard deviation scaling proved to be the best.

machine learning↗

Three-dimensional interconnected nickel phosphide networks with hollow microstructures and desulfurization performance

Graphical abstract: Three-dimensional interconnected nickel phosphide networks with hollow microstructures and desulfurization performance. - Highlights: • Three-dimensional Ni{sub 2}P has been prepared using foam nickel as a template. • The microstructures interconnected and formed sponge-like porous networks. • Three-dimensional Ni{sub 2}P shows superior hydrodesulfurization activity. - Abstract: Three-dimensional microstructured nickel phosphide (Ni{sub 2}P) was fabricated by the reaction between foam nickel (Ni) and phosphorus red. The as-prepared Ni{sub 2}P samples, as interconnected networks, maintained the original mesh structure of foamed nickel. The crystal structure and morphology of the as-synthesized Ni{sub 2}P were characterized by X-ray diffraction, scanning electron microscopy, automatic mercury porosimetry and X-ray photoelectron spectroscopy. The SEM study showed adjacent hollow branches were mutually interconnected to form sponge-like networks. The investigation on pore structure provided detailed information for the hollow microstructures. The growth mechanism for the three-dimensionally structured Ni{sub 2}P was postulated and discussed in detail. To investigate its catalytic properties, SiO{sub 2} supported three-dimensional Ni{sub 2}P was prepared successfully and evaluated for the hydrodesulfurization (HDS) of dibenzothiophene (DBT). DBT molecules were mostly hydrogenated and then desulfurized by Ni{sub 2}P/SiO{sub 2}.

36 MATERIALS SCIENCE↗

Ultrastructure of Critical-Gel-like Polyzwitterion–Polyoxometalate Complex Coacervates: Effects of Temperature, Salt Concentration, and Shear

The structure of unconventional complex coacervates, such as polymer-nonpolymer complex coacervates, remains less investigated than that of the conventional coacervates formed between oppositely charged polyelectrolytes with symmetric charge species. Yet, their microscopic structural organization is important to further elucidate the mechanism of liquid-liquid-phase separation processes upon complexation. In this work, we characterize the microstructural organization of complex coacervates formed between inorganic polyoxometalate (POM) and polyzwitterion by small-angle X-ray scattering (SAXS) with in situ temperature and shear control. Furthermore, despite the apparent transparent and homogeneous morphology of dense coacervates as observed by optical microscopy, our previous results (Macromolecules 2018, 51, 22, 9405-9411) suggest that dense polyzwitterion-POM coacervates exhibit critical-gel-like networks containing both complex-poor region (mesh pore) and complex-rich region (connective network). SAXS results as reported in this work support that the complex-rich region is actually in the form of loosely packed POM aggregates linked by polyzwitterion, designated as complex particles. The structure of aggregating complex particles is further examined against varied composition and salt concentrations, temperature, and shear, thanks to the high X-ray scattering contrast of POMs from that of other components in the coacervates. The complex particles in the dense coacervates appear to grow with more tightly packed POM aggregates with increasing POM-to-polyzwitterion concentration ratio, in contrast to more loosely packed POM aggregates with decreasing salt concentration. Conversely, increasing temperature could result in smaller complex particles containing more loosely packed POM aggregates, consistent with temperature-dependent viscoelasticity of dense coacervates. Furthermore, such POM-based hybrid dense coacervates exhibit intriguing strain-hardening behavior, resulting from shear-enhanced POM packing inside the complex particles. Distinct from the shear-thinning behavior of most conventional crosslinked polymeric networks, the strain-hardening property combined with thermal- and salt-responsive characteristics of polyzwitterion-POM coacervates could broaden the applications of hybrid organic-inorganic macroion coacervates as smart functional materials working under extreme environmental conditions.

36 MATERIALS SCIENCE↗

Finite deformation cohesive polygonal finite elements for modeling pervasive fracture

In this work, we introduce a framework for modeling dynamic fracture problems using cohesive polygonal finite elements. Random polygonal meshes provide a robust, efficient method for generating an unbiased network of fracture surfaces. Further, these meshes have more facets per element than standard triangle or quadrilateral meshes, providing more possible facets per element to insert cohesive surfaces. This property of polygonal meshes is advantageous for the modeling of pervasive fracture. We use both Wachspress and maximum entropy shape functions to form a finite element basis over the polygons. Fracture surfaces are captured through dynamically inserted cohesive zone elements at facets between the polygons in the mesh. Contact is enforced through a penalty method that is applied to both closed cohesive surfaces and general interpenetration of two polygonal elements. Finally, several numerical examples are presented that illustrate the capabilities of the method and demonstrate convergence of solutions.

97 MATHEMATICS AND COMPUTING↗

Closing the Loop between In Situ Stress Complexity and EGS Fracture Complexity

We present an agent-guided approach to CAD geometry decomposition that automates hex/hybrid meshing with graph neural networks (GNNs) to accelerate next-generation ModSim workflows. Our end-to-end pipeline (i) reduces 3D boundary-representation (B-Rep) models to a 2D chordal axis skeleton (CAT) and then to a 1D bipartite graph of surface and curve nodes, (ii) assigns per node labels as Cubit® WebCut actions, (iii) trains a multi-action GNN under supervised learning, and (iv) predicts five surface-node and three curve-node actions on out-of-distribution test geometries. Each graph node carries geometric, topological, and meshing attributes drawn from the B-Rep “skin” and CAT “skeleton,” with two-way mappings across 3D↔2D↔1D representations to maintain traceability back to 3D CAD. The supervised learning model exhibits stable convergence of the binary cross-entropy loss and achieves 98.7% accuracy on unseen lattice models. To operationalize decision-making, we rank predicted commands by geometric significance and prototyped the agent-guided workflow through the Cubit® Meshing PowerTool GUI. As a stretch goal, we explore reinforcement learning (RL) to reduce or remove label requirements and to learn policies for action sequences that maximize total reward (e.g., size of hex-meshable regions and resulting hex mesh quality). When all-hex meshing is not feasible, the agent assists in producing hybrid meshes—prioritizing hex in critical regions and transitioning to tetrahedral elements (tets) elsewhere—maintaining fidelity while ensuring robustness. The overarching objective is to replace manual, heuristics-based decomposition with data-driven, reproducible automation, cutting meshing turnaround time by orders of magnitude. We anticipate direct impact on simulation workflows through intelligent, scalable decomposition of complex CAD models into hex-meshable subdomains.

42 ENGINEERING↗

Development of the FHR advanced natural circulation analysis code and application to FHR safety analysis

The University of California, Berkeley (UCB) is performing thermal hydraulics safety analysis to develop the technical basis for design and licensing of fluoride-salt-cooled, high-temperature reactors (FHRs). FHR designs investigated by UCB use natural circulation for emergency, passive decay heat removal when normal decay heat removal systems fail. The FHR advanced natural circulation analysis (FANCY) code has been developed for assessment of passive decay heat removal capability and safety analysis of these innovative system designs. The FANCY code uses a one-dimensional, semi-implicit scheme to solve for pressure-linked mass, momentum and energy conservation equations. Graph theory is used to automatically generate a staggered mesh for complicated pipe network systems. Heat structure models have been implemented for three types of boundary conditions (Dirichlet, Neumann and Robin boundary conditions). Heat structures can be composed of several layers of different materials, and are used for simulation of heat structure temperature distribution and heat transfer rate. Control models are used to simulate sequences of events or trips of safety systems. A proportional-integral controller is also used to automatically make thermal hydraulic systems reach desired steady state conditions. A point kinetics model is used to model reactor kinetics behavior with temperature reactivity feedback. The underlying large sparse linear systems in these models are efficiently solved by using direct and iterative solvers provided by the SuperLU code on high performance machines. Input interfaces are designed to increase the flexibility of simulation for complicated thermal hydraulic systems. In conclusion, this paper mainly focuses on the methodology used to develop the FANCY code, and safety analysis of the Mark 1 pebble-bed FHR under development at UCB is performed.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

nn-PINNs: Non-Newtonian physics-informed neural networks for complex fluid modeling

Time- and rate-dependent material functions in non-Newtonian fluids in response to different deformation fields pose a challenge in integrating different constitutive models into conventional computational fluid dynamic platforms. Considering their relevance in many industrial and natural settings alike, robust data-driven frameworks that enable accurate modeling of these complex fluids are of great interest. The main goal is to solve the coupled Partial Differential Equations (PDEs) consisting of the constitutive equations that relate the shear stress to the deformation and fully capture the behavior of the fluid under various flow protocols with different boundary conditions. Here, in this work, we present non-Newtonian physics-informed neural networks (nn-PINNs) for solving systems of coupled PDEs adopted for complex fluid flow modeling. The proposed nn-PINN method is employed to solve the constitutive models in conjunction with conservation of mass and momentum by benefiting from Automatic Differentiation (AD) in neural networks, hence avoiding the mesh generation step. nn-PINNs are tested for a number of different complex fluids with different constitutive models and for several flow protocols. These include a range of Generalized Newtonian Fluid (GNF) empirical constitutive models, as well as some phenomenological models with memory effects and thixotropic timescales. nn-PINNs are found to obtain the correct solution of complex fluids in spatiotemporal domains with good accuracy compared to the ground truth solution. We also present applications of nn-PINNs for complex fluid modeling problems with unknown boundary conditions on the surface, and show that our approach can successfully recover the velocity and stress fields across the domain, including the boundaries, given some sparse velocity measurements.

42 ENGINEERING↗

Directing Colloidal Assembly and a Metal-Insulator Transition Using a Quench-Disordered Porous Rod Template

Replica and effective-medium theory methods are employed to elucidate how to massively reconfigure a colloidal assembly to achieve globally homogeneous, strongly clustered, and percolated equilibrium states of high electrical conductivity at low physical volume fractions. A key idea is to employ a quench-disordered, large-mesh rigid-rod network as a templating internal field. By exploiting bulk phase separation frustration and the tunable competing processes of colloid adsorption on the low-dimensional network and fluctuation-driven colloid clustering in the pore spaces, two distinct spatial organizations of greatly enhanced particle contacts can be achieved. As a result, a continuous, but very abrupt, transition from an insulating to metallic-like state can be realized via a small change of either the colloid-template or colloid-colloid attraction strength. Furthermore, the approach is generalizable to more complicated template or colloidal architectures.

36 MATERIALS SCIENCE↗

Gelation and Re-entrance in Mixtures of Soft Colloids and Linear Polymers of Equal Size

Liquid mixtures composed of colloidal particles and much smaller non-adsorbing linear homopolymers can undergo a gelation transition due to polymer-mediated depletion forces. We now show that the addition of linear polymers to suspensions of soft colloids having the same hydrodynamic size yields a liquid-togel-to-re-entrant liquid transition. In particular, the dynamic state diagram of 1,4-polybutadiene star-linear polymer mixtures was determined with the help of linear viscoelastic and small-angle Xray scattering experiments. While keeping the star polymers below their nominal overlap concentration, a gel was formed upon increasing the linear polymer content. Further addition of linear chains yielded a re-entrant liquid. This unexpected behavior was rationalized by the interplay of three possible phenomena: (i) depletion interactions, driven by the size disparity between the stars and the polymer length scale which is the mesh size of its entanglement network; (ii) colloidal deswelling due to the increased osmotic pressure exerted onto the stars; and (iii) a concomitant progressive suppression of the depletion efficiency on increasing the polymer concentration due to reduced mesh size, hence a smaller range of attraction. Our results unveil an exciting new way to tailor the flow of soft colloids and highlight a largely unexplored path to engineer soft colloidal mixtures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Flow and transport in three-dimensional discrete fracture matrix models using mimetic finite difference on a conforming multi-dimensional mesh

Here, we present a comprehensive workflow to simulate single-phase flow and transport in fractured porous media using the discrete fracture matrix approach. The workflow has three primary parts: (1) a method for conforming mesh generation of and around a three-dimensional fracture network, (2) the discretization of the governing equations using a second-order mimetic finite difference method, and (3) implementation of numerical methods for high-performance computing environments. A method to create a conforming Delaunay tetrahedralization of the volume surrounding the fracture network, where the triangular cells of the fracture mesh are faces in the volume mesh, that addresses pathological cases which commonly arise and degrade mesh quality is also provided. Our open-source subsurface simulator uses a hierarchy of process kernels (one kernel per physical process) that allows for both strong and weak coupling of the fracture and matrix domains. We provide verification tests based on analytic solutions for flow and transport, as well as numerical convergence. We also provide multiple expositions of the method in complex fracture networks. In the first example, we demonstrate that the method is robust by considering two scenarios where the fracture network acts as a barrier to flow, as the primary pathway, or offers the same resistance as the surrounding matrix. In the second test, flow and transport through a three-dimensional stochastically generated network containing 257 fractures is presented.

97 MATHEMATICS AND COMPUTING↗