SEARCH · Search NASA

Results for “finite difference method”

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.

SourceResearch connectionFrom yearThrough year

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

At least 829 records · Page 46

All sources NASA NTRS DOE OSTI

DOE OSTI2019

Combined crystal plasticity and grain boundary modeling of creep in ferritic-martensitic steels: II. The effect of stress and temperature on engineering and microstructural properties

This paper describes a series of physically-based crystal plasticity finite element method (CPFEM) simulations of long-term creep and creep rupture of Grade 91 steel. It is Part 2 of a two part series of papers. Part 1 describes the simulation framework; this part focuses on specific simulations and on how the predicted long-term creep properties of Grade 91 compare to the assumptions used in current high temperature design practices. This work extends the model developed in Part 1 to look at creep properties at different temperatures, principal stresses, and multiaxial stress states. The simulations show that empirically extrapolating creep rupture stresses from short-term experimental data may substantially over predict the actual long-term creep properties of Grade 91. Additionally, the CPFEM calculations predict a transition from notch strengthening creep behavior for high values of maximum principal stress and moderate notch severity to notch weakening behavior for low principal stresses and more severe notches. The latter regime better categorizes conditions in engineering components designed for long term elevated temperature use, which implies Grade 91 may be a notch weakening material in actual service. This would have a significant impact on high temperature design practices, though confirmatory test data on long-life, low stress notched specimens is difficult to obtain. Finally, one advantage of the physically-based modeling approach adopted here is that the simulation results also elucidate the microstructural mechanisms causing the macroscopic trends in engineering properties predicted by the simulations. This paper shows that the detailed micromechanical mechanisms predicted by the CPFEM simulations can be abstracted with a simple micromechanical model that can be used to both explain the detailed results and make improved predictions of engineering properties from experimental data.

36 MATERIALS SCIENCE↗

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 OSTI2019

Reconfigurable metamaterial absorber as an optical switch based on organic-graphene control

In this paper, we have suggested an effective optical switch with the cross fractal shape. It contains the main cross element that connected to four parasitic elements by four organic’s material junctions. This absorber is modified as a tunable element by using the graphene layer which can be considered for compensating the physical distortion. Each organic and parasitic load is placed over a graphene layer. These graphene layers separated with a SiO{sub 2} thin film. We have implemented a SiO{sub 2} layer as a spacer between the graphene and metal ground layer. The organic materials are popular for their various resistances value. We have applied a bulky dielectric slab with two different conductivity values for high and low resistance (HRS and LRS). This absorber works at the range of 90–100 THz for mid-infrared spectrum and the bandwidth of 4%. In the HRS mode the reflection value is about − 35 up to − 25 dB. When the organic material switched to LRS mode, this value drastically increased to − 5 dB. The designed absorber has both the tunable and switchable attributes. We have simulated this switch by the finite integral technique method with CST microwave studio.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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 OSTI2018

Multi-objective design of microvascular panels for battery cooling applications

Building on a recently developed optimization method based on an interface-enriched generalized finite element method, multiple objective functions are considered for the optimization of 2D networks of microchannels embedded in battery-cooling panels. The objective functions considered in this study are a differentiable alternative to the maximum temperature (the p-mean of the temperature), the pressure drop and the variance of the temperature. The -constraint method and the normalized normal constraint method are used to generate the pressure-temperature Pareto optimal front of the multi-objective optimization problem. The effects of different operating constraints/conditions such as localization of heat sources, prescribed pump power and imposed flow rate on the optimal designs are investigated. In addition to the topology of the embedded network, the cross sections of the microchannels are also introduced as design parameters to further improve the pressure drop of the designs. The resulting variable-cross-section optimized design is validated with experiment.

energy storage (including batteries and capacitors↗

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 OSTI2017

Weak-lensing Power Spectrum Reconstruction by Counting Galaxies. I. The ABS Method

We propose an analytical method of blind separation (ABS) of cosmic magnification from the intrinsic fluctuations of galaxy number density in the observed galaxy number density distribution. The ABS method utilizes the different dependences of the signal (cosmic magnification) and contamination (galaxy intrinsic clustering) on galaxy flux to separate the two. It works directly on the measured cross-galaxy angular power spectra between different flux bins. It determines/reconstructs the lensing power spectrum analytically, without assumptions of galaxy intrinsic clustering and cosmology. It is unbiased in the limit of an infinite number of galaxies. In reality, the lensing reconstruction accuracy depends on survey configurations, galaxy biases, and other complexities due to a finite number of galaxies and the resulting shot noise fluctuations in the cross-galaxy power spectra. We estimate its performance (systematic and statistical errors) in various cases. We find that stage IV dark energy surveys such as Square Kilometre Array and Large Synoptic Survey Telescope are capable of reconstructing the lensing power spectrum at z≃1 and ℓ≲5000 accurately. This lensing reconstruction only requires counting galaxies and is therefore highly complementary to cosmic shear measurement by the same surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

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↗

DOE OSTI2020

Simulation of the Fast Reactor Fuel Assembly Duct-Bowing Reactivity Effect Using Monte Carlo Neutron Transport and Finite Element Analysis

This paper discusses a new method of simulating the fuel assembly duct-bowing reactivity coefficient for EBR-II run 138B. Quantification of the fuel assembly duct-bowing reactivity effect in liquid metal–cooled fast reactors has been a persistent problem since they were first designed and operated. Simulation of the duct-bowing reactivity effect is difficult because the level of detail required to simulate the effect has exceeded most modeling capabilities. The new method outlined in this paper utilizes the finite element analysis code ANSYS to analyze the thermal and structural components. Here, the displacement of the fuel assembly duct due to thermal expansion and mechanical interaction was calculated by ANSYS using recorded EBR-II run 138B temperature and power boundary value data. The displacement values were incorporated into to a Monte Carlo model of EBR-II run 138B and keff was calculated. Multiple Monte Carlo calculations were performed with duct displacement values corresponding to different reactor temperatures. Using the calculated keff values associated with the different duct displacement results allowed calculation of the duct-bowing reactivity coefficient. The duct-bowing reactivity coefficient was calculated to be –14.5 × 10 –4 $/°C/ ± 4.4%.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2019

Numerical Investigation of the Effect of Rolling on the Localized Stress and Strain Induction for Wire + Arc Additive Manufactured Structures

Cold rolling can be used in-process or post-process to improve microstructure, mechanical properties and residual stress in directed-energy-deposition techniques, such as the high deposition rate wire + arc additive manufacturing (WAAM) process. Finite element simulations of the rolling process are employed to investigate the effect of rolling parameters, in particular rolling load and roller profile radius on the residual stress field as well as plastic strain distribution for the profiled roller. The results show the response to rolling of commonly used structural metals in WAAM, i.e., AA2319, S335JR steel and Ti-6Al-4V, taking into account the presence of residual stresses. The rolling load leads to changes in the location and the maximum value of the compressive residual stresses, as well as the depth of the compressive residual stresses. However, the roller profile radius only changes the maximum value of these compressive residual stresses. Changing the rolling load influences the equivalent plastic strain close to the top surface of the wall as well as in deeper areas, whereas the influence of the roller profile radius is negligible. The plastic strain distribution is virtually unaffected by the initial residual stresses prior to rolling. Finally, design curves were generated from the simulations for different materials, suggesting ideal rolling load and roller profile combinations for microstructural improvement requiring certain plastic strains at a specific depth of the additive structure.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Measurement-driven, model-based estimation of residual stress and its effects on fatigue crack growth. Part 1: Validation of an eigenstrain model

The objective of this paper is to validate a measurement-driven, model-based approach to estimate residual stress (RS) in samples machined from quenched aluminum stock. Model input is derived from measurement of RS in the parent stock. Validation is performed for prismatic T-sections removed from bars at different locations. We find RS predicted agrees with RS measured, by contour and neutron diffraction methods, with root-mean-square model-measurement difference of 22 MPa. Follow-on work (in Part 2) applies the RS estimation to samples representative of aircraft structures and examines the effects of RS on fatigue crack growth in the RS-bearing samples.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Digital image correlation and infrared thermography data for seven unique geometries of 304L stainless steel

Material Testing 2.0 (MT2.0) is a paradigm that advocates for the use of rich, full-field data, such as from digital image correlation and infrared thermography, for material identification. By employing heterogeneous, multi-axial data in conjunction with sophisticated inverse calibration techniques such as finite element model updating and the virtual fields method, MT2.0 aims to reduce the number of specimens needed for material identification and to increase confidence in the calibration results. To support continued development, improvement, and validation of such inverse methods—specifically for rate-dependent, temperature-dependent, and anisotropic metal plasticity models—we provide here a thorough experimental data set for 304L stainless steel sheet metal. The data set includes full-field displacement, strain, and temperature data for seven unique specimen geometries tested at different strain rates and in different material orientations. Commensurate extensometer strain data from tensile dog bones is provided as well for comparison. We believe this complete data set will be a valuable contribution to the experimental and computational mechanics communities, supporting continued advances in material identification methods.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Composites and additive manufacturing for high-field coil supports for stellarators

The utilization of composites layered on an additive manufactured substrate, for the production of coil supports for modular coils in small or middle size experimental stellarators is assessed. The focus of the study is a monolithic coil support comprising the coils of a half-period of a stellarator, somewhat similar to the ones in UST_2 and ARIES-CS stellarators. However, the concept may be applicable to quasi-monolithic coil supports (coil forms of NCSX type) or individual coil casings (W7-X type).Coil supports for stellarators require high precision, stiffness and strength for large contorted parts. Traditionally, monolithic coil supports are produced by steel casting/forging and final machining. This production method and material gives accurate, stiff and strong coil supports, but the method may be expensive for monolithic supports due to the geometrical complexity and required accuracy of the structure. In relation to those matters, this work investigates whether a monolithic coil support comprising an additive manufactured resin substrate, which is externally (outward from the coils) surrounded with a thick layer of fibre-reinforced resin, may achieve enough strength and stiffness under middle/high magnetic fields. Finite element calculations are produced to obtain the direction of the principal stresses and their values in compression and tension at different areas of the monolithic support, which is relevant for anisotropic materials. The feasibility of directional application of (carbon) fibres on the winding surface of the stellarator outward from the coils is experimentally tested on a scaled-down additively manufactured prototype of a monolithic support. The strength and stiffness of the composite structure appears sufficient for common magnetic fields in experimental stellarators, and the 3D-composite design and manufacturing was technically feasible.

Queral, Vicente↗