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 487 records · Page 27

A meshing framework for digital twins for extrusion based additive manufacturing

Additive manufacturing (AM) allows for manufacturing of complex three-dimensional geometries not typically realizable with standard manufacturing practices. The internal microstructure of AM components has a significant impact on mechanical, vibrational, and shock properties and permits richer design space when this is controllable. Due to complex interactions of internal geometry of an extrusion-based AM component, it is common practice to assume homogeneous behavior or to perform characterization testing on specific toolpath configurations. To avoid testing or material waste, it is necessary to develop a consistently accurate numerical simulation framework with relevant boundary value problems that can handle the complicated geometry of internal material microstructure present in AM components. Herein, a framework is proposed to directly create computational meshes suitable for finite element analysis (FEA) of the fine-scale features generated from extrusion-based AM tool paths to maintain a strong process–structure–property-performance linkage. This mesh can be manually or automatically analyzed using standard FEA simulations such as quasi-static preloading or modal analysis. The framework allows an in-silico assessment of a target AM geometry where fine-scale features greatly impact quantities of design interest such as in soft elastomeric lattices where toolpath infill can greatly influence the self-contact of a structure in compression, which we use as a motivating exemplar. This approach greatly reduces both time and resource waste present in traditional build and test design cycles for non-intuitive design spaces, and acts as a tool for use in the production of a key component of a digital twin, a mesh suitable for finite element analysis. In conclusion, it also further allows for the exploration of toolpath infill to optimize component properties beyond simple linear properties such as density and stiffness.

Additive manufacturing↗

Extreme Transverse Magnetoresistance in TiZn 16

Extreme magnetoresistance (XMR) is a phenomenon characterized by an increase in resistance by factors of 10 4 –10 7 % when a magnetic field is applied. This phenomenon is found in a number of semimetals such as WTe 2 , PtSn 4 , Cd 3 As 2 , and LaSb. The origin of XMR is still hotly debated, possibly with different materials having different (or multiple) explanations. Extreme transverse magnetoresistance of up to 8000% at 14 T and 1.8 K is measured in TiZn 16 , a semimetal with a multitude of bands crossing the Fermi energy, akin to PtSn 4 . The magnetoresistance is suppressed when the magnetic field is rotated to be parallel to the applied current, similar to PtSn 4 and PdSn 4 . The resistance of TiZn 16 follows Kohler's rule, but displays different behavior under an applied transverse field and under a longitudinal magnetic field, suggesting distinct electrical phases. Also present are Shubnikov-de Haas and de Haas-van Alphen oscillations with a transverse magnetic field up to 43 T, showing that despite an insulator-like temperature-resistance curve, charge carriers are still present. This positions TiZn 16 as an interesting addition to the investigation of XMR materials as a multi-band metal with complex Fermi surface geometries.

extreme magnetoresistance↗

A phase-field model for hydraulic fracture nucleation and propagation in porous media

Many geo-engineering applications, for example, enhanced geothermal systems, rely on hydraulic fracturing to enhance the permeability of natural formations and allow for sufficient fluid circulation. Over the past few decades, the phase-field method has grown in popularity as a valid approach to modeling hydraulic fracturing because of the ease of handling complex fracture propagation geometries. However, existing phase-field methods cannot appropriately capture nucleation of hydraulic fractures because their formulations are solely energy-based and do not explicitly take into account the strength of the material. Thus, in this work, we propose a novel phase-field formulation for hydraulic fracturing with the main goal of modeling fracture nucleation in porous media, for example, rocks. Built on the variational formulation of previous phase-field methods, the proposed model incorporates the material strength envelope for hydraulic fracture nucleation through two important steps: (i) an external driving force term, included in the damage evolution equation, that accounts for the material strength; (ii) a properly designed damage function that defines the fluid pressure contribution on the crack driving force. In conclusion, the comparison of numerical results for two-dimensional test cases with existing analytical solutions demonstrates that the proposed phase-field model can accurately model both nucleation and propagation of hydraulic fractures. Additionally, we present the simulation of hydraulic fracturing in a three-dimensional domain with various stress conditions to demonstrate the applicability of the method to realistic scenarios.

15 GEOTHERMAL ENERGY↗

Towards a classification of rank r $\mathscr{N}$ = 2 SCFTs. Part II. Special Kahler stratification of the Coulomb branch

We study the stratification of the singular locus of four dimensional $\mathscr{N}$ = 2 Coulomb branches. We present a set of self-consistency conditions on this stratification which can be used to extend the classification of scale-invariant rank 1 Coulomb branch geometries to two complex dimensions, and beyond. The calculational simplicity of the arguments presented here stems from the fact that the main ingredients needed - the rank 1 deformation patterns and the pattern of inclusions of rank 2 strata - are discrete topological data which satisfy strong self-consistency conditions through their relationship to the central charges of the SCFT. This relationship of the stratification data to the central charges is used here, but is derived and explained in a companion paper by one of the authors. We illustrate the use of these conditions by re-analyzing many previously-known examples of rank 2 SCFTs, and also by finding examples of new theories. The power of these conditions stems from the fact that for Coulomb branch stratifications a conjecturally complete list of physically allowed “elementary slices” is known. By contrast, constraining the possible elementary slices of symplectic singularities relevant for Higgs branch stratifications remains an open problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Simulations of Tidal Disruption Events

Numerical simulations have historically played a major role in understanding the hydrodynamics of the tidal disruption process. Given the complexity of the geometry of the system, the challenges posed by the problem have indeed stimulated much work on the numerical side. Smoothed Particles Hydrodynamics methods, for example, have seen their very first applications in the context of tidal disruption and still play a major role to this day. Likewise, initial attempts at simulating the evolution of the disrupted star with the so-called affine method have been historically very useful. Here, we provide an overview of the numerical techniques used in the field and of their limitations, and summarize the work that has been done to simulate numerically the tidal disruption process.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Defect and satellite characteristics of additive manufacturing metal powders

Metal additive manufacturing (AM) requires high-quality metal powders to three-dimensionally (3D) print metallic components with complex and customizable geometries. The lack of quantification of AM metal powders creates quality control challenges for 3D printed components, increases the uncertainty of printing reliability and net cost of inspected and certified printed components, and reduces the recyclability of used powders. However, critical characteristics of AM metal powders that are decisive factors for the 3D printing process, such as internal porosity, contamination, and satellite feature, remain ambiguous. In this work, we developed a novel approach to 3D quantify key characteristics of AM metal powders down to individual particles by using high-resolution synchrotron x-ray computed tomography. Empowered by the penetrative capability of high-energy x-ray, internal porosity and contamination within as-atomized metal powders from high-entropy alloys to nickel-based superalloys were evaluated. Additionally, the newly-developed dispersion method enables the homogeneous separation of individual particles, and consequently, results in the implementation of 3D particle shape analysis. To resolve a major challenge of identification and quantification of satellite-feature particles in as-atomized AM metal powders, the satellite features were quantitated by modeling and analyzing the shape parameter of local thickness variance. Furthermore, the 3D analytical methods of particle assessment in this study can be applied to other materials systems like rock, food, and pharmaceutical particles, and provide insights for process optimization across powder metallurgy, concrete, food and pharmaceutical manufacturing, and AM industries.

36 MATERIALS SCIENCE↗

Prediction of local concentration fields in porous media with chemical reaction using a multi scale convolutional neural network

The study of solute transport in porous media is of interest in many chemical engineering systems. Some example applications include packed bed catalytic reactors, filtration devices, and batteries. The pore scale modeling of these systems is time consuming and may require large computing resources, for this reason computational fluid dynamics (CFD) simulations are not practical if a large number of simulations is required, like in multiscale modeling, where a model at a large scale calls for pore scale simulations. It has been shown that neural networks can be trained with a dataset of flow simulations and then predict fields orders of magnitude faster, and with less computational resources, in new domains. However, it is crucial to provide the neural network with an effective description of the domain and the undergoing operating conditions to be able to train models that generalize accurately in unseen samples. Therefore, research is needed to employ neural networks in new complex systems. The appropriate training of a network for predicting coupled flow and solute transport processes is an outstanding problem due to the complex interplay between geometry and operating conditions. In this work, we train a multi scale convolutional neural network (MSNet) with a diverse dataset of simulations of transport and chemical reaction in porous media to predict the local concentration fields in images of porous media. Our dataset contains a wide diversity of sphere pack arrangements under different operating conditions (Péclet and Reynolds numbers). Further, we train a robust model by employing different input descriptors that represent the medium and the different operating conditions of each system. Our trained model is able to provide nearly instantaneous predictions, compared to around twenty hours of the CFD workflow, with less than 3.5% error on new geometries and transport conditions. Thus the model could be easily integrated in a multiscale workflow where fast response is needed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Calculations versus measurements for residual dose rates from SNS spent structures

Here, residual dose rate measurements were conducted on target vessel #13 and proton beam window #5 after extraction from their service locations. These measurements are used to verify calculation methods of radionuclide inventory assessment that are typically performed for nuclear waste characterization and transportation of these structures. Neutronics analyses for predicting residual dose rates are carried out using the transport code MCNPX and the transmutation code CINDER90. For transport analyses a complex and rigorous geometry model of the structures and their surroundings are applied. The neutronics analyses are carried out using the Bertini and CEM high energy physics models for simulating particles interactions above the table-based cross section range. Finally, obtained calculational results are analyzed and compared to the measured dose rates and overall show good agreement within 25%, which shows applicability of the methods used in analyses.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Modeling the Transmission of Optical Lightning Signals Through Complex 3‐D Cloud Scenes

Abstract Space‐based lightning imagers have shown that complex cloud scenes consisting of multiple tall convective features, anvil clouds, and warm boundary cloud layers are illuminated by lightning in many different ways depending on where the lightning occurs and how energetic it is. Modifications to the optical lightning signals from radiative transfer in the cloud medium can lead to reductions in detection efficiency and location accuracy for these instruments and can also cause some of the optical signals that are detected to have unexpected spatial energy distributions. In this study, we perform Monte Carlo radiative transfer simulations of optical lightning emissions in clouds with complex 3‐D geometries to shed some light on the origins of certain irregular radiance patterns that have been recorded from orbit. We show that diffuse reflections off nearby cloud faces can explain lightning signals in nonelectrified clouds, tall clouds can result in poor optical transmission and suppressed radiances that could lead to missed events and that particularly favorable viewing conditions can cause otherwise normal lightning to produce a superbolt that is orders of magnitude brighter than the same flash seen from a different angle.

54 ENVIRONMENTAL SCIENCES↗

Microscale Mechanical‐Chemical Modeling of Granular Salt: Insights for Creep

Abstract Microscale numerical modeling is potentially an effective approach for understanding salt creep, but it is quite challenging because natural salt grain boundaries can have complex and evolving geometry, and because contacts between these deformable salt mineral grains are dynamic as a result of compaction, chemical reaction, fluid flow, and heat transfer. In this study, we have overcome these challenges and developed a new microscale mechanical‐chemical (MC) model to analyze creep of salt at the microscale, accounting for coupled deformation, dynamic contacts, and chemical reaction in granular systems with realistic geometric representations. The MC model was realized by linking a new microscale mechanical code based on the numerical manifold method (NMM) to a reactive transport code named Crunch. We simulated the processes of reorganization of the salt grains, microfracturing, and pressure solution that contribute to creep at larger scales. Based on this first quantitative microscale model, we found that sharp corners of mineral grains can dominate the contact dynamics, microfracturing, and pressure solution when salt is compacted, thus governing the structural changes and porosity loss of the system. We found that pressure solution, which preferentially dissolves sharp corners and edges, can lead to relatively high porosity loss in the system, thus playing an important role in the creep of salt. Our analysis shows that the dynamic changes of salt granular systems involving grain relocation and pressure solution can occur repeatedly and continuously, thus contributing to longer‐term creep of salt at larger scales.

58 GEOSCIENCES↗

Data transfers for full core heterogeneous reactor high- fidelity multiphysics studies

Multiphysics simulations for nuclear reactor analysis are usually performed by resorting to operator splitting and fixed point iterations between single-physics solvers. This enables the separate solution of each physics, such as neutronics, fuel performance, and thermal hydraulics, on meshes tailored to the requirements of the respective numerical discretizations of the equations. As the equations are coupled, several fields must be transferred between single-physics solves. Projecting fields between meshes while preserving order of accuracy, conservation properties, and mapping non-overlapping geometries is a complex endeavor. This conference paper will present the transfers as implemented in MOOSE, which can handle arbitrary meshes, arbitrary mappings, conservation of integral quantities, and are made to scale with distributed simulations on both ends of the transfers. Their adequacy for advanced nuclear reactor multiphysics coupling is shown through examples and numerical studies.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Updated Thermofluid Performance of the Simplified Flat Variant of the HEMJ

Our group has recently developed and studied “finger”-type divertors that are a simplified version of the helium-cooled modular divertor with multiple jets (HEMJ) using coupled computational fluid dynamics and thermal stress simulations. Such a simplified geometry could reduce complexity and cost given the large number of fingers required to cover the total divertor target area. Previous experimental studies for this simplified flat design reported lower heat transfer coefficients and higher pressure drops than the HEMJ, contrary to numerical predictions. Subsequent measurements determined that the original test section had significant dimensional variations in the jet exit holes. A new test section was therefore manufactured and tested in the Georgia Tech (GT) helium loop. The experimental results presented here for this test section at maximum heat flux of 7.1 MW/m 2 are in good agreement with numerical predictions. Correlations developed from these experimental data are extrapolated to predict the maximum heat flux that can be accommodated by the flat design and the coolant pumping power requirements under prototypical conditions. Lastly, numerical simulations are used to estimate the sensitivity of the flat design to geometric variations typical of manufacturing tolerances and variations in the gap width.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Optimization of VISAR Mirrors for NIF Targets

The velocity interferometer system for any reflector (VISAR) is a key diagnostic at the National Ignition Facility that is used for measuring propagation and velocities of shocks through a sample for several high energy density campaigns. To support simultaneous measurements with VISAR and other diagnostics, optical mirrors are used to divert the VISAR signal off axis from the primary diagnostic line of sight, which can then be used for other measurements. Here, in this study, we describe our work on investigating issues with these mirrors to support complex new target geometries and signal variance due to system aging and more demanding polling rates. This paper also describes the key mitigation efforts performed to improve both the flatness and reflectivity of these commonly used mirrors.

National Ignition Facility↗

W7-X and the sawtooth instability: towards realistic simulations of current-driven magnetic reconnection

Magnetic reconnection in W7-X is studied by means of global numerical simulations in a series of models of increasing complexity. The magnetic geometry ranges from that of a cylinder to the full three-dimensional field of W7-X, and the equations solved range from ideal magnetohydrodynamics (MHD) to gyrokinetics. We simulate plasmas from the first operation phase with electron cyclotron current drive (ECCD). These are characterized by an equilibrium magnetic field featuring an ECCD-distorted 'humped' profile of the rotational transform ι, with ι = 1 in two radial locations. Such plasmas generally show sawtooth activity, hence motivating the present study. We pay particular attention to the role of equilibrium current density gradients in the destabilization of reconnecting modes. When the equilibrium temperature and density gradients are artificially suppressed (to eliminate the pressure gradient drive), the perturbed electrostatic potential is radially localized between the locations at which ι = 1. This is shown with a purely collisionless gyrokinetic model, in cylindrical geometry. In the real toroidal geometry of W7-X, for a non-ideal MHD model including a uniform resistivity, electron inertia and (numerical) viscosity, the same qualitative behaviour is observed. In particular, even if a resonant (m, n) = (1, -1) perturbation is initialized, the most unstable mode is the (m, n) = (-4, 4), where m and n are the poloidal and toroidal mode numbers, respectively. Other modes are destabilized due to geometric coupling. The growth rate of this instability scales as η 1/3 , where η is the plasma resistivity, thus suggesting that ECCD drives ideal MHD stable W7-X plasmas towards non-ideal marginality. An ideal magnetohydrodynamic analysis confirms the result. A fluid-kinetic hybrid version of the EUTERPE code shows that gyrokinetic ions have a stabilizing effect on these modes. For W7-X relevant collisionalities, the growth rate scales linearly with the electron skin depth, d e . Implications of our results for sawtoothing W7-X operation are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Behavior of thermally quenched topological defects in quasicrystal artificial spin ices

We report that understanding the formation of topological defects in connected magnetic nanowires such as in artificial spin ices (ASIs) is critical to controlling their behavior. Moreover, emergent frustration behavior of these defects can be expected due to complex underlying lattice geometry of the ASIs. We have explored in situ thermalization and quenching of quasicrystalline ASIs, which exhibit additional frustration due to aperiodicity from near the Curie temperature, and we observe the formation of two sets of topological vortex defects: within the magnetic bars and within the vertices at which the bars meet. By varying the laser fluence, we show that the number of defects in the vertices varies as the quenching rate changes following a power law relation, which we interpret as being related to the ferromagnetic-to-paramagnetic phase transition. Conversely, the defects formed within the bars are metastable states related to the spontaneous magnetization reversal of the magnetic bars and are therefore strongly dependent on underlying lattice geometry. These results provide insight into the emergence and control of topological defects in confined frustrated magnetic systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Multifidelity validation of digital surrogates using variable-density turbulent mixing models

In this study, ensembles of experimental data are presented and utilized to compare and validate two models used in the simulation of variable-density (Atwood = 0.22), compressible turbulent mixing. Though models of this kind (Reynolds averaged NavierStokes and large-eddy simulations) have been validated extensively with more canonical flows in previous studies, here the present approach offers novelty in the complexity of the geometry, the ensemble-based validation, and the uniformity of the computational framework on which the models are tested. Moreover, all experimental and computational tasks were completed by the authors which has led to a tightly coupled experimental configuration with its “digital twin.” The experimental divergent-shock-tube facility and its data acquisition methods are described and replicated in simulation space. A 2D Euler model which neglects the turbulent mixing at the interface is optimized to experimental data using a Gaussian process. This model then serves as the basis for both the 2D RANS and 3D LES studies that make comparisons to the mixing-layer data from the experiment. A relatively simple RANS model is shown to produce good agreement with experimental data only at late flow development times. The LES ensembles generally show good agreement with experimental data but display sensitivity to the characterization of initial conditions. Resolution-dependent behavior is also observed for certain higher-order statistics of interest. Overall, the LES model successfully captures the effects of divergent geometry, compressibility, and combined nonlinear instabilities inherent to the problem. The successful prediction of mixing width and its growth rate highlight the existence of three distinct regimes in the development of the instability, each with similarities to previously studied instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Quench Protection of Direct Wind Magnets

The direct wind magnets have several unique features. The conductor layout is characterized by thin single conductors wound over long lengths and separated from the other turns with thick layer of insulation. Moreover, there are several layers of winding. The strong non-uniformity in the magnetic field distribution results in different current and thermal margins for quenching in different regions. The quench protection modeling involves multi-physics coupling between electrical, thermal and magnetic transient. The material properties vary over time and location. This coupled with the complexity in the geometry adds to long computation times. The objective of this study is to understand the quench propagation in an EIC high inductance direct wind magnet which enables one to design appropriate quench protection strategies to protect the magnet. An advanced modeling program has been developed to simulate direct wind magnet quench. In conclusion, we validate theoretical simulations with experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗