Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite difference methods”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Massively parallel axisymmetric fluid model for streamer discharges

A highly parallelizable fluid plasma simulation tool based upon the first-order drift-diffusion equations is discussed. Atmospheric pressure plasmas have densities and gradients that require small element sizes in order to accurately simulate the plasm resulting in computational meshes on the order of millions to tens of millions of elements for realistic size plasma reactors. To enable simulations of this nature, parallel computing is required and must be optimized for the particular problem. Here, a finite-volume, electrostatic drift-diffusion implementation for low-temperature plasma is discussed. The implementation is built upon the Message Passing Interface (MPI) library in C++ using Object Oriented Programming. The underlying numerical method is outlined in detail and benchmarked against simple streamer formation from other streamer codes. Electron densities, electric field, and propagation speeds are compared with the reference case and show good agreement. Convergence studies are also performed showing a minimal space step of approximately 4 μm required to reduce relative error to below 1% during early streamer simulation times and even finer space steps are required for longer times. Additionally, strong and weak scaling of the implementation are studied and demonstrate the excellent performance behavior of the implementation up to 100 million elements on 1024 processors. Lastly, different advection schemes are compared for the simple streamer problem to analyze the influence of numerical diffusion on the resulting quantities of interest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Rupture Jumping and Seismic Complexity in Models of Earthquake Cycles for Fault Stepovers with Off-Fault Plasticity

ABSTRACT Fault stepovers are prime examples of geometric complexity in natural fault zones that may affect seismic hazard by determining whether an earthquake rupture continues propagating or abruptly stops. However, the long-term pattern of seismicity near-fault stepovers and underlying mechanisms of rupture jumping in the context of earthquake cycles are rarely studied. Leveraging a hybrid numerical scheme combining the finite element and the spectral boundary integral methods, FEBE, we carry out fully dynamic simulations of sequences of earthquakes and aseismic slip for both compressive and tensile stepovers with off-fault plasticity. We consider a rate-and-state friction law for the fault friction and pressure-sensitive Drucker–Prager plasticity for the off-fault bulk response. We observe that the accumulation of plastic deformation, an indication of off-fault damage, is significantly different in the two cases, with more plastic deformation projected in the overlapping region for the tensile stepover. The seismic pattern for a tensile stepover is more complex than for a compressive stepover, and incorporating plasticity also increases complexity, relative to the elastic case. A tensile stepover with off-fault plasticity shows rupture segmentation, temporal clustering, and frequent rupture jumping from one fault to another. These results shed light on possible mechanisms of rupture jumping in fault stepovers as well as the long-term evolution of the fault zone.

Geochemistry & Geophysics↗
DOE OSTI2022

EEL: A MOOSE-based Application for 3D Electro-Chemo-Thermo-Mechanical Modeling of Solid-State Batteries

This report describes a tool for three dimensional, high fidelity, coupled electro-chemo-thermo-mechanical modeling of solid-state batteries. A complete version of the tool is available as open source software at https://github.com/hugary1995/eel.git. The theoretical framework of the tool revolves around an inf-sup statement on a total potential, comprising the Helmholtz free energy, the electrical kinetic potential, the chemical potential, the Fourier potential, the chemical reaction potential, and the external power expenditure. The tool uses the finite element framework of Multiphysics Object-Oriented Simulation Environment (MOOSE) and a variational formulation to solve the boundary value problem for solid-state battery incorporating the full set of multiphysics couplings. The variational formulation also enables a modular software architecture for the tool so additional physics can be easily included by specifying the new contribution to the total potential. The report discusses several method of manufactured solutions that were used to verify the implementation of different physics and the Butler-Volmer reaction kinetics at the electrolyte-electrode interface. Finally, the report discusses results from complete charge/discharge simulations of a quasi-1D and a 3D solid-state battery.

25 ENERGY STORAGE↗
DOE OSTI2021

Striving to translate shale physics across ten orders of magnitude: What have we learned?

Shales will play an important role in the successful transition of energy from fossil-based resources to renewables in the coming decades. Aside from being a significant source of low-carbon intensity fuels, like natural gas, they also serve as geologic seals of subsurface formations that may be used to isolate nuclear waste, sequester CO 2 , or store intermittent energy (e.g., solar hydrogen). Despite their importance, shales pose significant engineering and environmental challenges due to their nanoporous structure and extreme heterogeneity that spans at least ~10 orders of magnitude in spatial scale. Two challenges inhibit a system-level understanding: (1) the physics of fluid flow and phase behavior in shales are poorly understood due to the dominant molecular interactions between minerals and fluids under confinement, and (2) the apparent lack of scale separation that prevents a reliable (closed) description of the physics at any single scale of observation. In this review, we focus on the latter issue and discuss scale translation, which in its broadest sense is transforming data or simulations from one spatiotemporal scale to another. While effective scale translation is not exclusive to shales, but all geologic porous media, the need for it is especially acute in shales given their high degree of heterogeneity. Classical theories like homogenization, while indispensable, fail when scales are not separated. Other methods, like numerical upscaling, scale-translate in only one direction: small to large, but not the reverse, called downscaling. However, the confluence of advances in three areas are bringing challenging problems such as shales within reach: increased computational power and scalable algorithms; high-resolution imaging and multi-modal data acquisition; and machine learning to process massive amounts of data. While these advances equip geoscientists with a wide array of experimental and computational tools, no individual tool can probe the entire gamut of heterogeneity in shales. Their effective use, therefore, requires an ability to bridge between various data types obtained at different scales. The aim of this review is to present a coherent account of computational and experimental methods that may be used to achieve just that, i.e., to perform scale translation. We provide a broader definition of scale translation, one that transcends classical homogenization and upscaling methods, but is consistent with them and accommodates notions like downscaling and data translation. After a brief introduction to homogenization, we review hybrid methods, numerical upscaling and its recent extensions, multiscale computing, high-resolution imaging, and machine learning. We place particular emphasis on multiscale computing and propose an algorithmic framework to bridge between the pore (micro) and Darcy (macro) scales. Throughout the paper, we draw comparisons between the various methods and highlight their (often hidden) similarities, differences, benefits, and pitfalls. We finally conclude with two case studies on shales that exemplify some of the methods presented.

58 GEOSCIENCES↗
DOE OSTI2022

Development and testing of a solar flux heating freeze recovery system for molten salt parabolic troughs

A major obstacle preventing commercialization of molten salt parabolic trough plants is freeze protection and recovery of the solar field. One of the most common nitrate salt compositions under consideration allows for solar field outlet temperatures up to 565°C and begins solidifying above 240°C. Using molten salt directly in the solar field can significantly reduce thermal energy storage costs and the higher temperature compared to oil heat transfer fluids can power a more efficient power cycle, but the very high freezing point is cause for concern. This work builds on a previous modeling effort investigating different methods for melting Solar Salt frozen in the solar field, which suggested the viability of a novel solar heating method to lower costs of the freeze protection system. In this work higher fidelity models are developed to better investigate the melting process with a more detailed 3D geometry which includes the insulated bellows. The thermal-fluid model of the melting process is updated to include temperature dependent density and data for solid phase properties. In addition, a finite element model is developed to resolve thermal stresses at the point with highest thermal gradients. Results with the new model confirm the viability of using controllable solar flux heating to thaw salt frozen in the solar field without damaging the receiver; however, the presence of non-illuminated sections significantly slows down the melting process. Solar heating simulations with pauses off-sun of 120 and 60 s require 18.5 and 11.9 h to melt salt from a night-time temperature of 10°C with 1000 W m2 DNI, which results in maximum thermal stresses of 48.2 and 53.5 MPa, respectively. Adding 150 W m−1 to the solar heating simulation with a 60 s pause off-sun decreases both the melt time and thermal stresses due to the more uniform heating, suggesting a combined heating method may be the best option.

14 SOLAR ENERGY↗
DOE OSTI2022

Two-body weak currents in heavy nuclei

In light and medium-mass nuclei, two-body weak currents from chiral effective field theory account for a significant portion of the phenomenological quenching of Gamow-Teller transition matrix elements. Here we examine the systematic effects of two-body axial currents on Gamow-Teller strength and β-decay rates in heavy nuclei within energy-density functional theory. Using a Skyrme functional and the charge-changing finite amplitude method, we add the contributions of two-body currents to the usual one-body linear response in the Gamow-Teller channel, both exactly and though a densitymatrix expansion. The two-body currents, as expected, usually quench both summed Gamow-Teller strength and decay rates, but by an amount that decreases as the neutron excess grows. They can enhance individual low-lying transitions, however, leading to decay rates that are quite different from those that an energy-independent quenching would produce, particularly in neutron-rich nuclei. We show that both these effects are related to changes in the total nucleon density as the number of neutrons increases.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Performance of a Hybrid HVAC-Integrated Thermal Storage Device

Thermal equipment in buildings is a primary contributor to peak loads on the electrical grid. Thermal energy storage is a cost-effective strategy to decouple electric use from thermal loads, thus reducing grid peak costs for building owners. One method for storing thermal energy in a building is to integrate a phase change material (PCM) directly into the heating, ventilation, and air conditioning system. These systems often require additional glycol loops, pumps, valves, and heat exchangers to couple the storage to the cooling system and building space, which increases the complexity and cost. This work will discuss an alternate approach where the storage is added directly into the heat pump evaporator. A detailed two-dimensional finite difference heat transfer model of a PCM-refrigerant-glycol heat exchanger was developed to simulate the performance of this component. The fluid stream was discretized along the flow direction to capture changes in the fluid properties and local heat transfer rates, and the phase change material was discretized in both the x and y directions to capture the movement of the melt front. The model was used to understand the impact of different material and geometric properties on the charge and discharge characteristics of the device. Finally, a Ragone framework analogous to that used for electrochemical batteries was used to maximize the energy density and round-trip efficiency of the device while supplying loads appropriate for space cooling in buildings.

buildings↗
DOE OSTI2024

Topanga: A kinetic ion plasma code for large-scale ionospheric simulations on magnetohydrodynamic timescales

Topanga is a kinetic ion code developed for simulating large-scale plasma phenomena in the Earth's ionosphere on magnetohydrodynamic timescales. It is a domain-decomposed parallel code that runs on high-performance computing platforms. Features of Topanga include spherical geometry for simplified boundary conditions and computational efficiency; a hybrid plasma model with inertia-less fluid electrons, kinetic ions, and an electric field specified via an Ohm's law; a Maxwell-FDTD (finite difference time domain) plasma model which retains the displacement current in Maxwell's equations and models electron currents in the ionosphere with a tensor conductivity; sponge-layer boundary conditions for absorption of electromagnetic and plasma waves incident on the domain boundaries; and a novel mixed-implicit algorithm for evolving the EM fields inside the Maxwell-FDTD region that is stable over many orders of magnitude in the electron–ion collision frequency. We verify the numerical methods used in Topanga on a pair of test problems. The first test involves modeling a three-dimensional collisionless shock using the hybrid set of equations. The second test involves modeling a spherical TEM mode in vacuum using the Maxwell-FDTD set of equations. Finally, we demonstrate how using the combined set of hybrid and Maxwell-FDTD equations to model the Starfish Prime high-altitude nuclear test recovers a “missing” EM signal on the ground that is not present when using only the hybrid set of equations. The magnitude of this signal in the simulation containing the Maxwell-FDTD region agrees well with the E3a portion of the magnetohydrodynamic electromagnetic pulse from Starfish Prime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Method and apparatus for shape-based energy analysis of solids

A computer-readable medium stores instructions including storing a model of a physical structure and defining a mesh for the model. The mesh includes finite elements, each defined by a respective set of edges. The instructions include, for each finite element, identifying a governing differential equation and a set of complementary functions that exactly satisfy the governing differential equation. The instructions include determining an applied physical stimulus for the physical structure. The instructions include generating an energy optimization model that minimizes a difference between internal energy and external energy of the finite elements in response to the applied physical stimulus. The model includes a matrix of respective scalar multipliers for the complementary functions for each finite element. The instructions include transforming the matrix and calculating a physical parameter of interest. The instructions include, in response to the physical parameter not satisfying a design parameter, updating the model of the physical structure.

Spears, Robert E.↗
DOE OSTI2020

Structural modeling and validation of laminated stacks in magnetic gearing applications

Laminated structures are commonly employed in magnetic gears based on their ability to direct magnetic flux, while inhibiting eddy currents. Unfortunately, the stiffness of the laminated structure is much lower than a solid part, which can lead to unacceptable deflections, given the large magnetic forces and small air gaps in magnetic gears. Finite element analysis of laminations to accurately predict these deflections can be computationally expensive due to the small elements needed to model the thin lamination dimensions. These analyses also have some uncertainty due to the assignment of bonding stiffness. This paper presents a means of accurately modelling the deflection of a laminated stack with reduced complexity and computational time, which is compared to experimental results. Here, the resulting methods are then used to assess the magnetic and mechanical advantages of different designs for the modulating pieces of a magnetic gear.

42 ENGINEERING↗
DOE OSTI2025

Hypercomplex Automatic Differentiation in the Eulerian Hydrocode PAGOSA

Enabling the computation of partial derivatives or sensitivities in production hydrocodes is beneficial for design, optimization, sensitivity analysis, and uncertainty quantification. Traditional finite difference approximations of these sensitivities are inefficient since convergence studies of the step size is required for each parameter of interest. For these reasons, HYPercomplex Automatic Differentiation (HYPAD) was implemented in the Eulerian hydrocode PAGOSA. HYPAD is analogous to forward-mode automatic differentiation except hypercomplex numbers (numbers with multiple imaginary parts) are used instead of dual numbers. Accurate partial derivatives can be computed of all state variables with respect to multiple input variables in a single run. The method was implemented using operator overloading to handle hypercomplex algebra. HYPAD was demonstrated and verified on Sod’s shock tube problem to compute derivatives of the state variables with respect to a material parameter, initial conditions, and geometry.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Origin of high optical contrast in zinc-zinc oxide electrodeposits for dynamic windows

The control of solar light and heat emission through windows is an important strategy for increasing the energy efficiency of buildings. Reversible Zn electrodeposition has recently emerged as a promising method for constructing electronically tintable robust dynamic windows. In Zn electrodeposits formed from dimethyl sulfoxide (DMSO) electrolytes during device tinting, we observe extraordinary absorption that is in excess of what is predicted by the Beer-Lambert law for a uniform Zn thin film. Here, the charge required to electroplate these films is abnormally low, significantly less than previously reported dynamic windows based on reversible metal electrodeposition, which facilitates the construction of large-area devices that switch uniformly. Finite-difference time-domain (FDTD) simulations are used to investigate the origin of this enhanced absorption, which arises from plasmonic effects among the Zn nanoparticles and ZnO dendrites in the film. The dielectric ZnO dendrites promote absorption via slit-like Zn-dielectric-Zn structures that from hybrid surface plasmon resonance at metal walls. Through these investigations, we provide design principles to construct low-charge and high-contrast metal and metal oxide-based dynamic windows.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Characterization of anisotropic mechanical properties of polymer composites from a hybrid additive manufacturing-compression molding process using x-ray computer tomography

In this study, anisotropic stiffness tensors were reconstructed based on fiber orientation distributions obtained from X-ray computer tomography (xCT). A preform was manufactured via a big area additive manufacturing (BAAM) system with carbon fiber (CF) filled acrylonitrile butadiene styrene (ABS). The tailored preform from additive manufacturing (AM) was used in the compression molding (CM) process to produce a low-void high-performance thermoplastic composite panel. An xCT technique was employed to detect the fiber orientations in CF/ABS composites manufactured via three different methods: AM from BAAM, extrusion compression molding (ECM), and AM-CM. The anisotropic stiffness tensor was obtained from the composite panel manufactured via the three manufacturing methods (AM, ECM, and AMCM). A micromechanics theory was used to obtain the orthotropic stiffness tensors of the composite panels and compared with the experimental values. The predicted stiffness tensors of AM and AM-CM composite panels were used to study the deformation characteristics of a steering wheel during airbag deployment by performing finite element analysis (FEA). The approach developed in this study can be utilized for evaluating high-performance composites.

Pokkalla, Deepak↗
DOE OSTI2020

Robust block preconditioners for poroelasticity

In this paper we develop and analyze several preconditioners for the linear systems arising from discretized poroelasticity problems. The preconditioners include one block preconditioner for the two-field Biot model and several preconditioners for the classical three-field Biot model. We manage to analyze these different preconditioners under a same theoretical framework and show that all of them are uniformly optimal with respect to material and discretization parameters. Numerical tests demonstrate the robustness of these preconditioners.

42 ENGINEERING↗
DOE OSTI2026

Scattering Processes from Quantum Simulation Algorithms for Scalar Field Theories

We provide practical simulation methods for scalar field theories on a quantum computer that yield improved asymptotics as well as concrete gate estimates for the simulation and physical qubit estimates using the surface code. We achieve these improvements through two optimizations. First, we consider a finite volume approach for estimating the elements of the S-matrix. This approach is appropriate in general for 1+1D and for certain low-energy elastic collisions in higher dimensions. Second, we implement our approach using a series of different fault-tolerant simulation algorithms for Hamiltonians formulated both in the field occupation basis and field amplitude basis. Our algorithms are based on either second-order Trotterization or qubitization. The cost of Trotterization in occupation basis scales as O ( λ N 7 | Ω | 3 / ( M 5 / 2 ϵ 3 / 2 ) ) where λ is the coupling strength, N is the occupation cutoff, | Ω | is the volume of the spatial lattice, M is the mass of the particles and ϵ is the uncertainty in the energy calculation used for the S -matrix determination. Qubitization in the field basis scales as O ( | Ω | 2 ( k 2 Λ + k M 2 ) / ϵ ) , where k is the cutoff in the field and Λ is a scaled coupling constant. We find in both cases that the bounds suggest physically meaningful simulations can be performed using on the order of 4 × 10 6 physical qubits and 10 12 T -gates which corresponds to roughly one day on a superconducting quantum computer with surface code and a cycle time of 100 ns. This places the simulation of scalar field theory within striking distance of the gate counts for the best available chemistry simulation results.

Hardy, Andrew [Toronto U.] (ORCID:0000000235817382↗
DOE OSTI2025

A Simple, Scalable Large Deformation Solid Mechanics Implementation in the MOOSE Framework

This article describes a large deformation solid mechanics solver implemented as part of the freely available and open source MOOSE finite element simulation framework. The article documents the choices made in developing the solid mechanics framework and describes novel formulations for the gradient operator and constitutive modeling framework made to simplify implementations of different coordinate systems, stabilized gradient operators, and different constitutive model inputs and outputs. In the process, the article describes a new formulation that casts objective integration of the Cauchy stress as a linear transformation of the small stress rate. Finally, the article presents key implementation details and examines the parallel efficiency of the solid mechanics solver implemented in MOOSE. The implementation retains a good weak scaling efficiency beyond 1,000 parallel processes. The article includes a discussion of the factors limiting the parallel efficiency of implicit, large deformation solid mechanics codes on current high-performance computers, with the main current limitation being the scalability of the algebraic multigrid methods used to solve the linearized equilibrium equations.

Applied computing → Computer-aided design↗
DOE OSTI2020

Simulation scaling studies of reactor core two-phase flow using direct numerical simulation

Tremendous growth in supercomputing power in recent years has resulted in the emergence of high-resolution flow analysis methods as an advanced research tool to evaluate single and two-phase flow behavior. In particular, unstructured mesh-based methods have been applied to analyze flows in complex reactor core geometries, including those of light water reactors (LWR). The finite-element based code, PHASTA, is utilized to perform large-scale simulations of two-phase bubbly flows in LWR geometries. Given the large computational cost of direct numerical simulation (DNS) coupled with interface tracking methods (ITM), typical domains encompass a portion of a single subchannel. In the presented research, the state-of-the-art analysis of turbulent two-phase flows in complex LWR subchannel geometries are demonstrated at both prototypical reactor parameters as well as scaled low pressure conditions. Three different cases are studied, a high-pressure simulation in prototypical reactor subchannel geometry, a low-pressure case in prototypical geometry and a final low-pressure case in a geometry scaled up to conserve the ratio between the bubble size and the domain pitch. Utilizing advanced statistical processing tools, these simulation conditions are compared to shed light on the relevancy of two-phase flow characteristics given the significant differences between LWR and low-pressure conditions. These findings can lead to the generation of useful guiding principles when researchers need to scale the two-phase flow behavior captured at low pressure and temperature conditions to those at reactor operating conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗