Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mesh data structures”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Enabling the Broader Use of MOOSE for Nuclear Energy and Other Simulation

This Final Scientific and Technical Report summarizes work performed under the Phase IIA SBIR project “Enabling the Broader Use of MOOSE for Nuclear Energy and Other Simulation” (DE-SC0020906) from August 2023 through August 2025. The objective of the Phase IIA effort was to mature and harden capabilities developed during Phase II, with the goal of enabling practical interoperability between Coreform’s isogeometric analysis (IGA) technologies and the Multiphysics Object-Oriented Simulation Environment (MOOSE), while improving robustness, performance, and scalability for complex, nuclear-relevant geometries. Over the course of Phase IIA, the project established and validated an extraction-based interoperability pathway between Coreform tools and MOOSE. A combined mesh and matrix format was defined collaboratively with MOOSE developers and integrated into the solver, enabling standard MOOSE workflows to operate on data exported from Coreform’s IGA and Flex Representation Method (FRM) pipelines. Early demonstrations validated architectural compatibility using linear solid mechanics problems, while later efforts focused on benchmark testing and external use. By the end of the project period, engineers at BWXT were able to independently set up and execute a simulation using the Coreform–MOOSE workflow and provide direct feedback that informed further refinement. In parallel, substantial effort was devoted to improving the robustness of trimmed U-spline construction for complex CAD geometries. A growing test suite of nuclear-relevant models was compiled through collaboration with multiple stakeholders and used to drive extensive bug fixing and reliability improvements. These efforts resulted in improved robustness and performance, including the addition of fallback capabilities that enhance reliability when the underlying commercial CAD kernel fails. Performance-oriented work progressed later in the project, with the development and demonstration of methods to decompose complex geometries into structured subregions and updated data representations to support more efficient solver processing. Additionally, extensive enhancements to threadsafe parallel data structures and trimming operations established a foundation for scalable processing of large assemblies. Collaboration with Sandia National Laboratories on the SGM geometric modeling kernel advanced to a functioning interface test case, positioning the workflow for future kernel integration. Overall, the Phase IIA effort successfully transitioned the project from architectural proof-of-concept to externally exercised, solver-integrated capability, while clarifying remaining technical challenges related to standardization, performance optimization, and kernel integration.

42 ENGINEERING↗
DOE OSTI2015

Predictive Engineering Tools for Injection-Molded Long-Carbon-Fiber Thermoplastic Composites (FY 2015 Third Quarterly Report)

During the third quarter of FY 2015, the following technical progress has been made toward project milestones: 1) Magna oversaw the tool build and prepared the molding plan for the complex part of Phase II. 2) PlastiComp hosted a visit by Magna and Toyota on April 23rd to finalize the molding scope and schedule. The plan for molding trials including selection of molding parameters for both LFT and D-LFT for the U-shape complex part was established. 3) Toyota shipped the U-shape complex part tool to Magna on May 28th, 2015. 4) Plasticomp provided 30wt% LCF/PP and 30wt% LCF/PA66 compounded pellets to Magna for molding the complex part. 5) Magna performed preliminary molding trials on June 2nd, 2015 to validate wall thickness, fill profile, tool temperature and shot size requirements for the complex part. 6) Magna performed the first complex part run on June 16th and 17th, 2015 at Magna’s Composite Centre of Excellence in Concord, ON, Canada. Dale Roland of Plasticomp, and Umesh Gandhi of Toyota also attended the molding. 7) Magna discussed and finalized the plan with PNNL and the team for cutting samples from molded parts at selected locations for fiber orientation and length measurements. 8) Magna provided the computer-aided design (CAD) files of the complex parts with and without ribs to PNNL and Autodesk to build the corresponding ASMI models for injection molding simulations. Magna also provided the actual parameters used. 9) Plasticomp’s provided knowledge and experience of molding LCF materials essential to the successful molding of the parts including optimization of fill speed, tool temperatures, and plasticizing conditions for the 30wt% LCF/PP and 30wt% LCF/PA66 materials in both rib and non-rib versions. 10) Magna molded additional parts for evaluation of mechanical property testing including torsional stiffness on June 29th and 30th, 2015 at Magna’s Composite Center of Excellence. 11) Toyota began preparation for the torsion test of the specimens. Preparation of a computer-aided engineering (CAE) model to predict the performance is in progress. 12) Autodesk fixed an error in the implementation of the proper orthogonal decomposition (POD) calculation of fiber length that had caused the ASMI solution to crash and provided an updated build of ASMI containing the fix. 13)Autodesk reviewed and provided feedback for the complex part molding and measurement locations. 14) Autodesk provided support to set up the workflow for ASMI-ABAQUS® analysis, and provided a fix and workaround for a bug in the ASMI-ABAQUS® output command. 15) Autodesk helped build ASMI analysis models for the complex parts with and without ribs. 16) Autodesk worked on improving the orientation prediction accuracy in the shearing layer for 3D meshes based on comparison to measured data of the plaque moldings. 17) PNNL installed a new ASMI version received from Autodesk and performed comparative analyses to assess mid-plane versus 3D fiber length predictions using the full fiber length model and the reduced-order model (ROM) using POD. 18) PNNL presented the project scope, accomplishments, significant results and future plans to DOE and the USCAR Materials Tech Team on June 3rd, 2015. 19) PNNL discussed the cutting of samples from molded parts and finalized a plan with Magna and the team suggesting the sample size, locations and number of samples per location. 20) PNNL and Autodesk built ASMI models for the complex parts with and without ribs, and preliminary analyses of the part with ribs were conducted using the actual molding parameters received from Magna. 21) PNNL worked on a procedure to extract fiber orientation and length results from a 3D ASMI analysis to a 3D ABAQUS model. This procedure is essential to import ASMI fiber orientation and length to a 3D ABAQUS model of the part allowing future part structural analysis for weight reduction study.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Parallel Infrastructure Modeling and Inversion Module for E4D

Electrical resistivity tomography ERT is a method of imaging the electrical conductivity of the subsurface. Electrical conductivity is a useful metric for understanding the subsurface because it is governed by geomechanical and geochemical properties that drive subsurface systems. ERT works by injecting current into the subsurface across a pair of electrodes, and measuring the corresponding electrical potential response across another pair of electrodes. Many such measurements are strategically taken across an array of electrodes to produce an ERT data set. These data are then processed through a computationally demanding process known as inversion to produce an image of the subsurface conductivity structure that gave rise to the measurements. Data can be inverted to provide 2D images, 3D images, or in the case of time-lapse 3D imaging, 4D images. ERT is generally not well suited for environments with buried electrically conductive infrastructure such as pipes, tanks, or well casings, because these features tend to dominate and degrade ERT images. This reduces or eliminates the utility of ERT imaging where it would otherwise be highly useful for, for example, imaging fluid migration from leaking pipes, imaging soil contamination beneath leaking subsurface tanks, and monitoring contaminant migration in locations with dense network of metal cased monitoring wells. The location and dimension of buried metallic infrastructure is often known. If so, then the effects of the infrastructure can be explicitly modeled within the ERT imaging algorithm, and thereby removed from the corresponding ERT image. However, there are a number of obstacles limiting this application. 1) Metallic infrastructure cannot be accurately modeled with standard codes because of the large contrast in conductivity between the metal and host material. 2) Modeling infrastructure in true dimension requires the computational mesh to be highly refined near the metal inclusions, which increases computational demands. 3) The ERT imaging algorithm requires specialized modifications to accommodate high conductivity inclusions within the computational mesh. The solution to each of these challenges was implemented within E4D (formerly FERM3D), which is a parallel ERT imaging code developed at PNNL (IPID #30249). The infrastructure modeling module implement in E4D uses a method of decoupling the model at the metallic interface(s) boundaries, into several well posed sub-problems (one for each distinct metallic inclusion) that are subsequently solved and recombined to form the global solution. The approach is based on the immersed interface method, with has been applied for similar problems in other fields (e.g. semiconductor industry). Comparisons to analytic solutions have shown the results to be very accurate, addressing item 1 above. The solution is implemented about an unstructured mesh, which enables arbitrary shapes to be efficiently modelled, thereby addressing item 2 above. In addition, the algorithm is written in parallel and shows excellent scalability, which also addresses equation 2 above. Finally, because only the boundaries of metallic inclusions are modeled, there are no high conductivity cells within the modeling mesh, and the problem described by item 3 above is no longer applicable.

Johnson, TimC.↗
DOE OSTI2019

Numerical Analysis of Detonability Assessment in a Natural Gas-Air Fueled Rotating Detonation Engine

Rotating Detonation Engines (RDE) offer an alternative combustion strategy to replace conventional constant pressure combustion with a process that could produce a pressure gain without the use of a mechanical compressor. Recent numerical and experimental publications that consider air as the oxidizer have primarily focused on the ability of these annular combustors to sustain a stable continuous detonation wave when fueled by hydrogen. However, for this to be a viable consideration for the land-based power generation it is necessary to explore the ability to detonate natural gas and air within the confines of the annular geometry of an RDE. Previous studies on confined detonations have expressed the importance of permitting detonation cells to fully form within the combustor in order to achieve stability. This poses a challenge for natural gas–air fueled processes as their detonation cell size can be quite large even at moderate pressures. Despite the practical importance, only a few studies are available on natural gas detonations for air-breathing RDE applications. Moreover, the extreme thermodynamic condition (high temperature inside the combustor) allows limited accessibility inside the combustor for detailed experimental instrumentations, providing mostly single-point data. Recent experimental studies at the National Energy Technology Laboratory (NETL) have reported detonation failure at higher methane concentration in an air-breathing RDE fueled by natural gas-hydrogen fuel blends. This encourages to perform a detailed numerical investigation on the wave characteristics of detonation in a natural gas-air fueled RDE to understand the various aspects of instability associated with the natural gas-air detonation.This study is a numerical consideration of a methane-air fueled RDE with varying operating conditions to ascertain the ability to achieve a stable, continuous detonation wave. The simulations have been performed in a 2D unwrapped RDE geometry using the open-source CFD library “OpenFOAM” employing an unsteady pressure-based compressible reactive flow solver with a <em>k</em>–<em>ε</em> turbulence model in a structured rectangular grid system. Both reduced and detailed chemical kinetic models have been used to assess the effect of the chemistry on the detonation wave characteristics and the underlying flow features. A systematic grid sensitivity study has been conducted with various grid sizes to quantify the weakly stable overdriven detonation on a coarse mesh and oscillating features at fine mesh resolutions.The main focus of the current study is to investigate the effects of operating injection pressure on detonation wave characteristics of an air-breathing Rotating Detonation Engine (RDE) fueled with natural gas-hydrogen fuel blends. Wave speeds, peak pressures and temperatures, and dominant frequencies have been computed from the time histories. The flow structures were then visualized using 2D contours of temperature and species concentration.

Saha, Pankaj↗
DOE OSTI2012

NEAMS Update: Quarterly report for January - March 2012

Quarterly highlights are: (1) The integration of Denovo and AMP was demonstrated in an AMP simulation of the thermo-mechanics of a complete fuel assembly; (2) Bison was enhanced with a mechanistic fuel cracking model; (3) Mechanistic algorithms were incorporated into various lower-length-scale models to represent fission gases and dislocations in UO2 fuels; (4) Marmot was improved to allow faster testing of mesoscale models using larger problem domains; (5) Component models of reactor piping were developed for use in Relap-7; (6) The mesh generator of Proteus was updated to accept a mesh specification from Moose and equations were formulated for the intermediate-fidelity Proteus-2D1D module; (7) A new pressure solver was implemented in Nek5000 and demonstrated to work 2.5 times faster than the previous solver; (8) Work continued on volume-holdup models for two fuel reprocessing operations: voloxidation and dissolution; (9) Progress was made on a pyroprocessing model and the characterization of pyroprocessing emission signatures; (10) A new 1D groundwater waste transport code was delivered to the used fuel disposition (UFD) campaign; (11) Efforts on waste form modeling included empirical simulation of sodium-borosilicate glass compositions; (12) The Waste team developed three prototypes for modeling hydride reorientation in fuel cladding during very long-term fuel storage; (13) A benchmark demonstration problem (fission gas bubble growth) was modeled to evaluate the capabilities of different meso-scale numerical methods; (14) Work continued on a hierarchical up-scaling framework to model structural materials by directly coupling dislocation dynamics and crystal plasticity; (15) New 'importance sampling' methods were developed and demonstrated to reduce the computational cost of rare-event inference; (16) The survey and evaluation of existing data and knowledge bases was updated for NE-KAMS; (17) The NEAMS Early User Program was launched; (18) The Nuclear Regulatory Commission (NRC) Office of Regulatory Research was introduced to the NEAMS program; (19) The NEAMS overall software quality assurance plan (SQAP) was revised to version 1.5; and (20) Work continued on NiCE and its plug-ins and other utilities, such as Cubit and VisIt.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2015

GPU-based simulation of the two-dimensional unstable structure of gaseous oblique detonations

In this paper, the two-dimensional structure of unstable oblique detonations induced by the wedge from a supersonic combustible gas flow is simulated using the reactive Euler equations with a one-step Arrhenius chemistry model. A wide range of activation energy of the combustible mixture is considered. Computations are performed on the Graphical Processing Unit (GPU) to reduce the simulation runtimes. A large computational domain covered by a uniform mesh with high grid resolution is used to properly capture the development of instabilities and the formation of different transverse wave structures. After the initiation point, where the oblique shock transits into a detonation, an instability begins to manifest and in all cases, the left-running transverse waves first appear, followed by the subsequent emergence of right-running transverse waves forming the dual-head triple point structure. This study shows that for low activation energies, a long computational length must be carefully considered to reveal the unstable surface due to the slow growth rate of the instability. For high activation energies, the flow behind the unstable oblique detonation features the formation of unburnt gas pockets and strong vortex-pressure wave interaction resulting in a chaotic-like vortical structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Orbit-averaging and deposition accuracy for runaway electron beams in hybrid kinetic-MHD simulations of the runaway plateau

We develop a new procedure that combines the kinetic orbit runaway electron code (KORC) and the NIMROD extended-magnetohydrodynamic code to simulate runaway electrons (REs) in the post-disruption plateau. KORC integrates guiding-center orbits, with a barycentric-based binary search strategy providing initial guesses for the Newton–Raphson logical-to-physical coordinate inversion, ensuring reliable particle-to-mesh mapping in NIMROD, whose fields remain static for the present study. Samples are drawn in accord with experimental parallel current profiles of RE beams during the plateau phase. Deposition in NIMROD is verified through comparison with a Python-based finite-element code that ensures periodicity in the poloidal direction and continuity at the magnetic axis. Accurate representation of near-axis fields requires finer mesh resolution to prevent under- and overshoots in current density from orbit inaccuracies. Yet, at a fixed particle count, increasing mesh resolution amplifies statistical noise in the deposited fields. An orbit-averaging method accumulates partial current deposits over multiple kinetic steps and reduces the statistical noise with little added computational cost. By coupling kinetic routines from KORC directly into the NIMROD codebase, these developments lay essential groundwork for future self-consistent KORC–NIMROD coupling.

Algorithms and data structure↗
DOE OSTI2017

An immersed boundary method for fluid–structure interaction with compressible multiphase flows

This paper presents a two-dimensional immersed boundary method for fluid–structure interaction with compressible multiphase flows involving large structure deformations. This method involves three important parts: flow solver, structure solver and fluid–structure interaction coupling. In the flow solver, the compressible multiphase Navier–Stokes equations for ideal gases are solved by a finite difference method based on a staggered Cartesian mesh, where a fifth-order accuracy Weighted Essentially Non-Oscillation (WENO) scheme is used to handle spatial discretization of the convective term, a fourth-order central difference scheme is employed to discretize the viscous term, the third-order TVD Runge–Kutta scheme is used to discretize the temporal term, and the level-set method is adopted to capture the multi-material interface. In this work, the structure considered is a geometrically non-linear beam which is solved by using a finite element method based on the absolute nodal coordinate formulation (ANCF). The fluid dynamics and the structure motion are coupled in a partitioned iterative manner with a feedback penalty immersed boundary method where the flow dynamics is defined on a fixed Lagrangian grid and the structure dynamics is described on a global coordinate. We perform several validation cases (including fluid over a cylinder, structure dynamics, flow induced vibration of a flexible plate, deformation of a flexible panel induced by shock waves in a shock tube, an inclined flexible plate in a hypersonic flow, and shock-induced collapse of a cylindrical helium cavity in the air), and compare the results with experimental and other numerical data. The present results agree well with the published data and the current experiment. Finally, we further demonstrate the versatility of the present method by applying it to a flexible plate interacting with multiphase flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

A mass-conserving lattice Boltzmann method with dynamic grid refinement for immiscible two-phase flows

A mass-conserving lattice Boltzmann method (LBM) for multiphase flows is presented in this paper. The proposed LBM improves a previous model (Lee and Liu, 2010 [21]) in terms of mass conservation, speed-up, and efficiency, and also extends its capabilities for implementation on non-uniform grids. The presented model consists of a phase-field lattice Boltzmann equation (LBE) for tracking the interface between different fluids and a pressure-evolution LBM for recovering the hydrodynamic properties. In addition to the mass conservation property and the simplicity of the algorithm, the advantages of the current phase-field LBE are that it is an order of magnitude faster than the previous interface tracking LBE proposed by Lee and Liu (2010) [21] and it requires less memory resources for data storage. Meanwhile, the pressure-evolution LBM is equipped with a multi-relaxation-time (MRT) collision operator to facilitate attainability of small relaxation rates thereby allowing simulation of multiphase flows at higher Reynolds numbers. Additionally, we reformulate the presented MRT-LBM on nonuniform grids within an adaptive mesh refinement (AMR) framework. Various benchmark studies such as a rising bubble and a falling drop under buoyancy, droplet splashing on a wet surface, and droplet coalescence onto a fluid interface are conducted to examine the accuracy and versatility of the proposed AMR-LBM. The proposed model is further validated by comparing the results with other LB models on uniform grids. A factor of about 20 in savings of computational resources is achieved by using the proposed AMR-LBM. As a more demanding application, the Kelvin–Helmholtz instability (KHI) of a shear-layer flow is investigated for both density-matched and density-stratified binary fluids. The KHI results of the density-matched fluids are shown to be in good agreement with the benchmark AMR results based on the sharp-interface approach. When a density contrast between the two fluids exists, a typical chaotic structure in the flow field is observed at a Reynolds number of 10000, which indicates that the proposed model is a promising tool for direct numerical simulation of two-phase flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2015

Visualization and Quantification of Rotor Tip Vortices in Helicopter Flows

Helicopter aeromechanics encompasses a highly vortical flow field. The vortices generated at each blade tip contain unsteady, complex, three-dimensional structures, which interact with each other, other blades, the fuselage and various components of the helicopter. It is crucial to understand vortex kinematics and their subsequent dynamic evolution. Much research has been devoted to the understanding of helicopter vortex dynamics, including a number of experimental studies.1-6 In May 2010 Particle Image Velocimetry (PIV) measurements of a full-scale UH-60A rotor were acquired in the National Full-Scale Aerodynamics Complex (NFAC) 40- by 80-Foot Wind Tunnel.1 These measurements were taken at a plane just downstream of the advancing blade in the vicinity of the blade tip—the so-called PIV plane. The resulting PIV data were then processed using an ensemble-average approach to create graphical representations of the vortical wake velocity and vorticity fields, which, in turn, have enhanced the understanding of rotorcraft vortical wake flow field physics and have provided a more detailed validation of vortical wake computer simulations.7 A common approach used to analyze flow field features is to compute and plot color contour maps of various scalar quantities such as pressure, velocity magnitude and vorticity magnitude. For example, the color map of the vorticity magnitude is typically used to determine vortical flow structure. With this approach the vortex core may appear larger or smaller, depending on the contour levels that are selected. Thus, the resulting visualization is sensitive to user-specified contour levels. For vortex core radius measurements, it is more accurate to calculate the vortex core radius using the cross-flow velocity profile across the vortex core. The task of extracting the cross-flow velocity profile can be time consuming with existing tools since the user needs to manually select the core center then specify sampling points along the profile axis. The task becomes even more challenging when the associated grid system uses AMR (Adaptive Mesh Refinement) where the profile axis could span multiple grid blocks. There are a number of existing techniques for profiling of vortex core attributes;8-9 however, these techniques are not fully automatic in that the user still needs to select the vortex core center to compute the cross-flow velocity profile. The present study introduces a new color map scheme that is based on the vortex core radius, which is fully automatic and does not require user intervention. Analysis and visualization of blade tip vortices on the PIV plane using the proposed new color map scheme are described in Section II. The new approach is evaluated using two case studies, which are described in Section III. The paper ends with a summary in Section IV.

Quantification↗
DOE OSTI2016

Predicting Fueled Experiment Swelling and Thermal Expansion Using STAR-CCM+

An experiment being designed at Idaho National Laboratory (INL) has an unusual temperature requirement. The surface temperature must remain above 300 C over the irradiation lifetime. The fuel in this experiment is composed of uranium-nitride (UN) tri-structural-isotropic (TRISO) particles encased in a cylindrical, SiC matrix. Increased swelling of SiC at lower temperatures is the basis for this lower temperature limit. For a typical thermal analysis, the focus is on some upper limit of the center-line temperature. This tends to be a more straight-forward process where the analyst makes a series of worst-case assumptions and demonstrates that the experiment - even under those conditions - cannot exceed the maximum temperature. The model used to ensure the experiment will always have a surface temperature in excess of 300 C has been developed using STAR-CCM+. Notable properties of the simulation include - Full 3D fluid and heat transfer modeling; - Geometry input taken directly from design engineer's CAD files; - Fuel swelling and thermal conductivity as functions of displacements per atom (DPA) and temperature; - Radial thermal expansion of fuel, rodlet, and capsule; - Jump gap and surface roughness implemented using equations from BISON. These things together have yielded a model with sufficient detail to predict minimum surface temperatures with the greatest precision attainable with available data. The only notable approximation is the zirconia insulator, which will actually be a sphere. This particular shape was chosen to expedite meshing, while maintaining a realistic degree of physical contact with adjacent components. Given temperature and irradiation dependent data for thermal conductivity and fuel swelling of a SiC matrix, it was possible to perform an in-depth thermal analysis to find an optimal configuration of radial gas gaps. Of greatest importance was that the fuel surface temperature remain above 300 C. This criteria was met with a healthy margin, even when considering machining tolerances of ±1 mil on the gas gaps.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Structure- and Physics-Preserving Reductions of Power Grid Models

The large size of multiscale, distribution and transmission, power grids hinders fast systemwide estimation and real-time control and optimization of operations. This paper studies graph reduction methods of power grids that are favorable for fast simulations and follow-up applications. While the classical Kron reduction has been successful in reduced order modeling of power grids with traditional, hierarchical design, the selection of reference nodes for the reduced model in a multiscale, distribution and transmission, network becomes ambiguous. In this work we extend the use of the iterative Kron reduction by utilizing the electric grid's graph topology for the selection of reference nodes, consistent with the design features of multiscale networks. Additionally, we propose further reductions by aggregation of coherent subnetworks of triangular meshes, based on the graph topology and network characteristics, in order to preserve currents and build another power-flow equivalent network. Our reductions are achieved through the use of iterative aggregation of subgraphs that include general tree structures, lines, and triangles. Important features of our reduction algorithms include that (i) the reductions are either equivalent to the Kron reduction or otherwise produce a power-flow equivalent network; (ii) due to the mentioned power-flow equivalence, the reduced network can model the dynamic of the swing equations for a lossless, inductive, steady state network; (iii) the algorithms efficiently utilize hash tables to store the sequential reduction steps. The third feature allows for easy reintroduction of detailed models into the reduced, conceptual network and makes the final reduced order model backward compatible with a sequence of intermediate, partially reduced networks with varying resolution---the ordered sequence of iterative reductions corresponds to a sequence of reduced order models. The performance of our graph reduction algorithms and features of the reduced grids are discussed on a real-word transmission and distribution grid. We produce visualizations of the reduced models through open source libraries and release our reduction algorithms with example code and toy data.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

The Thermal Plumbing System of Stromboli Volcano, Aeolian Islands (Italy) Inferred From Electrical Conductivity and Induced Polarization Tomography

Abstract We performed the first 3D island‐scale tomography of the electrical conductivity of Stromboli volcano (Aeolian Islands, Italy) using 2D acquisition lines (37.2 km) and a total of 18,880 measurements and 2,402 unique electrode locations. This 3D data set was inverted using a Gauss‐Newton algorithm, parallel‐processing on an unstructured tetrahedral mesh containing 678,420 finite‐element nodes and 3,580,145 elements to account for the topography of the volcanic island. The tomogram exhibits a conductive body (10 −2 –1.0 S m −1 ) consistent with the location of CO 2 and temperature anomalies observed at the ground surface. It corresponds to the hydrothermal system with high electrical conductivity associated with alteration. In order to confirm this interpretation, a 2.5D large‐scale induced polarization tomography was performed crossing the volcano. The joint interpretation of the conductivity and normalized chargeability is done with a petrophysical model previously tested and verified at both shield‐ and strato‐volcanoes. This model implies that alteration (through the effect of the cation exchange capacity associated with clay minerals and zeolites) plays a strong role in both controlling the electrical conductivity and normalized chargeability at Stromboli volcano. A temperature tomogram, derived from the geoelectrical measurements, is consistent with surface temperature anomalies and the Very Long Period (VLP) seismicity related to the mild‐explosive activity. This survey displays at 600 m a.s.l. a lateral shift in the highest temperature location, also corresponding to the source of VLP seismicity. Structural boundaries have a major role in the hottest hydrothermal fluids rising below the active crater terrace of Stromboli volcano.

58 GEOSCIENCES↗
DOE OSTI2015

Mass transport and turbulence in gravitationally unstable disk galaxies. I. The case of pure self-gravity

The role of gravitational instability-driven turbulence in determining the structure and evolution of disk galaxies, and the extent to which gravity rather than feedback can explain galaxy properties, remains an open question. To address it, we present high-resolution adaptive mesh refinement simulations of Milky Way-like isolated disk galaxies, including realistic heating and cooling rates and a physically motivated prescription for star formation, but no form of star formation feedback. After an initial transient, our galaxies reach a state of fully nonlinear gravitational instability. In this state, gravity drives turbulence and radial inflow. Despite the lack of feedback, the gas in our galaxy models shows substantial turbulent velocity dispersions, indicating that gravitational instability alone may be able to power the velocity dispersions observed in nearby disk galaxies on 100 pc scales. Moreover, the rate of mass transport produced by this turbulence approaches ∼1 M{sub ⊙} yr{sup −1} for Milky Way-like conditions, sufficient to fully fuel star formation in the inner disks of galaxies. In a companion paper, we add feedback to our models, and use the comparison between the two cases to understand which galaxy properties depend sensitively on feedback and which can be understood as the product of gravity alone. All of the code, initial conditions, and simulation data for our model are publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRSDate not supplied

Proof Testing Aluminum and Titanium Thin-Walled Tubes for Aerospace Applications

Proof testing aerospace components aids in verifying the structural integrity of flight hardware. This process is often an important step in the certification process for critical parts. Though proof testing is a common practice at NASA and in industry, there is often a need to understand better the effectiveness of proof testing pressure systems discovering cracklike flaws or defects. Often proof tests are used in requirements and in practice as technical rationale for screening critical defects in a part in lieu of performing a non-destructive inspection. The intent of this work is to define proof test failure envelopes for thin-walled tubes composed of two materials: aluminum and titanium. Both materials are common in aerospace thin-walled tubing applications which is why these were the materials chosen for this study. The fracture mechanics software NASGRO was used to perform a parametric study on the effect of defects in thin-walled tubes subjected to proof testing. For both materials, the parameters varied in the NASGRO analyses were aspect ratio, proof test pressure, thickness, diameter, and flaw multiplier. With each case defined as a through crack or failure, the results were then pulled into MATLAB in order to parse through the data and produce failure envelope plots. Failure envelopes were created by meshing a surface in between the minimum proof pressure failing data points. As expected, based on material properties, the plots indicated the aluminum material would fail at much lower proof tests than the titanium. The parametric data and failure envelopes may be used as general qualitative insight on the effectiveness of proof testing thin-walled tubes of similar alloys.

Leah B Davis↗
NASA NTRS1995

Large HAWT wake measurement and analysis

From the theoretical fluid dynamics point of view, the wake region of a large horizontal-axis wind turbine has been defined and described, and numerical models of wake behavior have been developed. Wind tunnel studies of single turbine wakes and turbine array wakes have been used to verify the theory and further refine the numerical models. However, the effects of scaling, rotor solidity, and topography on wake behavior are questions that remain unanswered. In the wind tunnel studies, turbines were represented by anything from scaled models to tea strainers or wire mesh disks whose solidity was equivalent to that of a typical wind turbine. The scale factor compensation for the difference in Reynolds number between the scale model and an actual turbine is complex, and not typically accounted for. Though it is wise to study the simpler case of wakes in flat topography, which can be easily duplicated in the wind tunnel, current indications are that wind turbine farm development is actually occurring in somewhat more complex terrain. Empirical wake studies using large horizontal-axis wind turbines have not been thoroughly composited, and, therefore, the results have not been applied to the well-developed theory of wake structure. The measurement programs have made use of both in situ sensor systems, such as instrumented towers, and remote sensors, such as kites and tethered, balloonborne anemometers. We present a concise overview of the work that has been performed, including our own, which is based on the philosophy that the MOD-2 turbines are probably their own best detector of both the momentum deficit and the induced turbulence effect downwind. Only the momentum deficit aspects of the wake/machine interactions have been addressed. Both turbine power output deficits and wind energy deficits as measured by the onsite meteorological towers have been analyzed from a composite data set. The analysis has also evidenced certain topographic influences on the operation of spatially diverse wind turbines.

Miller, A. H.↗
NASA NTRS1995

Navier-Stokes Aerodynamic Simulation of the V-22 Osprey on the Intel Paragon MPP

The paper will describe the Development of a general three-dimensional multiple grid zone Navier-Stokes flowfield simulation program (ENS3D-MPP) designed for efficient execution on the Intel Paragon Massively Parallel Processor (MPP) supercomputer, and the subsequent application of this method to the prediction of the viscous flowfield about the V-22 Osprey tiltrotor vehicle. The flowfield simulation code solves the thin Layer or full Navier-Stoke's equation - for viscous flow modeling, or the Euler equations for inviscid flow modeling on a structured multi-zone mesh. In the present paper only viscous simulations will be shown. The governing difference equations are solved using a time marching implicit approximate factorization method with either TVD upwind or central differencing used for the convective terms and central differencing used for the viscous diffusion terms. Steady state or Lime accurate solutions can be calculated. The present paper will focus on steady state applications, although time accurate solution analysis is the ultimate goal of this effort. Laminar viscosity is calculated using Sutherland's law and the Baldwin-Lomax two layer algebraic turbulence model is used to compute the eddy viscosity. The Simulation method uses an arbitrary block, curvilinear grid topology. An automatic grid adaption scheme is incorporated which concentrates grid points in high density gradient regions. A variety of user-specified boundary conditions are available. This paper will present the application of the scalable and superscalable versions to the steady state viscous flow analysis of the V-22 Osprey using a multiple zone global mesh. The mesh consists of a series of sheared cartesian grid blocks with polar grids embedded within to better simulate the wing tip mounted nacelle. MPP solutions will be shown in comparison to equivalent Cray C-90 results and also in comparison to experimental data. Discussions on meshing considerations, wall clock execution time, load balancing, and scalability will be provided.

Vadyak, Joseph↗
NASA NTRS2011

Numerical Implementation of a Multiple-ISV Thermodynamically-Based Work Potential Theory for Modeling Progressive Damage and Failure in Fiber-Reinforced Laminates

A thermodynamically-based work potential theory for modeling progressive damage and failure in fiber-reinforced laminates is presented. The current, multiple-internal state variable (ISV) formulation, enhanced Schapery theory (EST), utilizes separate ISVs for modeling the effects of damage and failure. Damage is considered to be the effect of any structural changes in a material that manifest as pre-peak non-linearity in the stress versus strain response. Conversely, failure is taken to be the effect of the evolution of any mechanisms that results in post-peak strain softening. It is assumed that matrix microdamage is the dominant damage mechanism in continuous fiber-reinforced polymer matrix laminates, and its evolution is controlled with a single ISV. Three additional ISVs are introduced to account for failure due to mode I transverse cracking, mode II transverse cracking, and mode I axial failure. Typically, failure evolution (i.e., post-peak strain softening) results in pathologically mesh dependent solutions within a finite element method (FEM) setting. Therefore, consistent character element lengths are introduced into the formulation of the evolution of the three failure ISVs. Using the stationarity of the total work potential with respect to each ISV, a set of thermodynamically consistent evolution equations for the ISVs is derived. The theory is implemented into commercial FEM software. Objectivity of total energy dissipated during the failure process, with regards to refinements in the FEM mesh, is demonstrated. The model is also verified against experimental results from two laminated, T800/3900-2 panels containing a central notch and different fiber-orientation stacking sequences. Global load versus displacement, global load versus local strain gage data, and macroscopic failure paths obtained from the models are compared to the experiments.

Pineda, Evan J.↗