Search NASA⌕ Search

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.

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

At least 1,027 records · Page 57

EFFECT OF FINITE LARMOR RADIUS CORRECTIONS ON THE THERMAL INSTABILITY OF THERMALLY CONDUCTING VISCOUS PLASMA WITH HALL CURRENT AND ELECTRON INERTIA

The thermal instability of an infinite homogeneous, thermally conducting, and rotating plasma, incorporating finite electrical resistivity, finite electron inertia, and an arbitrary radiative heat-loss function in the presence of finite Larmor radius corrections and Hall current, has been studied. Analysis has been made with the help of linearized magnetohydrodynamics (MHD) equations. A general dispersion relation is obtained using the normal mode analysis method, and the dispersion relation is discussed for longitudinal propagation and transverse propagation separately. The dispersion relation has been solved numerically to obtain the dependence of the growth rate on the various parameters involved. The conditions of modified thermal instability and stability are discussed in the different cases of interest.

79 ASTRONOMY AND ASTROPHYSICS↗

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.↗

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↗

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↗

Exact relativistic expressions for polarization of incoherent Thomson scattering

We present a derivation of the degree of polarization for incoherent Thomson scattering (TS) using Mueller matrix formalism. An exact analytic solution is obtained for spectrum-integrated matrix elements. The solution is valid for the full range of incident polarizations, scattering angles, and electron thermal motion from non-relativistic to ultra-relativistic. It is based on a newly developed theoretical model, a finite transit time (FTT) correction to previous theoretical work on TS polarization. The Mueller matrix elements are substantially different from previous calculations without the FTT correction, even to the lowest linear order in T e =m e c 2 <<1. Mathematically, the derivation is a unique example of fully analytical integration of the 3D scattering operator over a relativistic Maxwellian distribution function; experimentally, the results have application to the use of the polarization properties of Thomson scattered light as a method of electron temperature measurement. The results can also be used as a reliable tool for benchmarking and verification of numerical codes for frequency resolved properties of TS polarization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SU-F-T-490: Separating Effects Influencing Detector Response in Small MV Photon Fields

Purpose: Different detector properties influence their responses especially in field sizes below the lateral electron range. Due to the finite active volume, the detector density and electron perturbation at other structural parts, the response factor is in general field size dependent. We aimed to visualize and separate the main effects contributing to detector behavior for a variety of detector types. This was achieved in an experimental setup, shielding the field center. Thus, effects caused by scattered radiation could be examined separately. Methods: Signal ratios for field sizes down to 8 mm (SSD 90 cm, water depth 10 cm) of a 6MV beam from a Siemens Primus LINAC were recorded with several detectors: PTW microDiamond and PinPoint ionization chamber, shielded diodes (PTW P-60008, IBA PFD and SNC Edge) and unshielded diodes (PTW E-60012 and IBA SFD). Measurements were carried out in open fields and with an aluminum pole of 4 mm diameter as a central block. The geometric volume effect was calculated from profiles obtained with Gafchromic EBT3 film, evaluated using FilmQA Pro software (Ashland, USA). Results: Volume corrections were 1.7% at maximum. After correction, in small open fields, unshielded diodes showed a lower response than the diamond, i.e. diamond detector over-response seems to be higher than that for unshielded diodes. Beneath the block, this behavior was amplified by a factor of 2. For the shielded diodes, the overresponse for small open fields could be confirmed. However their lateral response behavior was strongly type dependent, e.g. the signal ratio dropped from 1.02 to 0.98 for the P-60008 diode. Conclusion: The lateral detector response was experimentally examined. Detector volume and density alone do not fully account for the field size dependence of detector response. Detector construction details play a major role, especially for shielded diodes.

61 RADIATION PROTECTION AND DOSIMETRY↗

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↗

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↗

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↗

A new approach to enforcing discrete maximum principles in continuous Galerkin methods for convection-dominated transport equations

This work introduces a set of design principles and new algorithmic tools for enforcing maximum principles and/or positivity preservation in continuous finite element approximations to convection-dominated transport problems. Enabling a linear first-order advection equation as a model problem, we address the design of first-order artificial diffusion operators and their higherorder counterparts at the element matrix level. The proposed methodology leads to a nonlinear high-resolution scheme capable of resolving moving fronts and internal/boundary layers as sharp localized nonoscillatory features. The amount of numerical dissipation depends on the difference between the solution value at a given node and a local maximum or minimum. The shockcapturing numerical diffusion coefficient is designed to vanish as the nodal values approach a mass-weighted or linearity-preserving average. The universal applicability and simplicity of the element-based limiting procedure makes it an attractive alternative to edge-based algebraic flux correction.

Kuzmin, D.↗

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↗

In Situ TEM Microcompression Pillar Size Effects in Fe-9Cr ODS

The objective of this paper is to ascertain the existence of a specimen size threshold, below which the material no longer exhibits bulk mechanical behavior when tested using in situ transmission electron microscopic (TEM) micro-compression pillars. In this work, sub-micron sized pillars of a model Fe-9Cr oxide dispersion strengthened (ODS) alloy were compressed in TEM in situ at room temperature, and attempts have been made to quantify the resulting data. Results are compared to mechanical property values in the archival literature, and a discussion on the validity of quantitative in situ TEM micro-compression data for ODS steel is included. In situ TEM mechanical testing is a relatively new technique enabling concurrent TEM imaging and mechanical testing of sub-micron-sized electron-transparent specimens. Due to limited sample volume of many nuclear materials, in situ TEM mechanical testing holds great promise for advancing our understanding of deformation phenomena. But it is well known that the dramatically reduced specimen size required for in situ testing can affect quantitative and qualitative deformation behaviors. However, Kiener, et al. used in situ TEM micro-compression pillars on irradiated and unirradiated single-crystal copper. They observed size independence in the irradiated material, while the unnirradiated material displayed a size threshold below which the mechanical properties deviated from those of bulk Cu. They assert that the more densely-spaced obstacles in irradiated Cu allow it to display bulk properties even in specimens having ∼ 100 nm dimensions. Similar to Kiener's findings, we hypothesize that unirradiated ODS will demonstrate a size threshold, but the in situ TEM micro-compression technique can be extended from single crystal materials to engineering alloys. An Fe-9Cr model ODS alloy was chosen for this study because of the broad interest in ODS for advanced nuclear reactors. ODS alloys exhibit high irradiation tolerance and high temperature creep resistance. New insight into deformation processes gained through in situ TEM mechanical testing could be valuable for the validation of ODS alloys for next generation nuclear reactors. The fact that the compressive yield stress and elastic modulus do not follow distinct trends as a function of specimen size is likely due to the high density of obstacles to dislocation motion. The ODS microstructure contains such obstacles as oxide nanoparticles, a fine grain structure, a high density of dislocation lines, and carbide precipitates. Since the deformation behavior is governed by the minimum separation distance between obstacles [1], it is plausible that the high obstacle density enables ODS alloy micropillars to exhibit bulk-like mechanical properties even prior to irradiation. Future work will focus on: (1) developing finite element models to help explain the difference in elastic moduli between bulk and micropillar tests, (2) understanding the spread in yield stress measurements at fixed micropillar volumes, and (3) extending this study to irradiated ODS. (authors)

36 MATERIALS SCIENCE↗

WE-AB-202-05: Validation of Lung Stress Maps for CT-Ventilation Imaging

Purpose: To date, lung CT-ventilation imaging has been based on quantification of local breathing-induced changes in Hounsfield Units (HU) or volume. This work investigates the use of a stress map resulting from a biomechanical deformable image registration (DIR) algorithm as a metric of the ventilation function. Method: Eight lung cancer patients presenting different kinds of ventilation defects were retrospectively analyzed. Additionally, to the 4DCT acquired for radiotherapy planning, five of them had PET and three had SPECT imaging following inhalation of Ga-68 and Tc-99m, respectively. For each patient, the inhale phase of the 4DCT was registered to the exhale phase using Morfeus, a biomechanical DIR algorithm based on the determination of boundary conditions on the lung surfaces and vessel tree. To take into account the heterogeneity of the tissue stiffness in the stress map estimation, each tetrahedral element of the finite-element model was assigned a Young’s modulus ranging from 60kPa to 12MPa, as a function of the HU in the inhale CT. The node displacements and element stresses resulting from the numerical simulation were used to generate three CT-ventilation maps based on: (i) volume changes (Jacobian determinant), (ii) changes in HU, (iii) the maximum principal stress. The voxel-wise correlation between each CT-ventilation map and the PET or SPECT V image was computed in a lung mask. Results: For patients with PET, the mean (min-max) Spearman correlation coefficients r were: 0.33 (0.19–0.45), 0.36 (0.16–0.51) and 0.42 (0.21–0.59) considering the Jacobian, changes in HU and maximum principal stress, respectively. For patients with SPECT V, the mean r were: 0.12 (−0.12–0.43), 0.29 (0.22–0.45) and 0.33 (0.25–0.39). Conclusion: The maximum principal stress maps showed a stronger correlation with the ventilation images than the previously proposed Jacobian or change in HU maps. This metric thus appears promising for CT-ventilation imaging. This work was funded in part by NIH P01CA059827.

60 APPLIED LIFE SCIENCES↗

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↗

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↗

Mitigating cutting-induced plasticity in the contour method, Part 2: Numerical analysis

Cutting-induced plasticity can have a significant effect on the measurement accuracy of the contour method. The present study examines the benefit of a double-embedded cutting configuration that relies on self-restraint of the specimen, relative to conventional edge-crack cutting configurations. A series of finite element analyses are used to simulate the planar sectioning performed during double-embedded and conventional edge-crack contour cutting configurations. The results of numerical analyses are first compared to measured results to validate the cutting simulations. The simulations are then used to compare the efficacy of different cutting configurations by predicting the deviation of the residual stress profile from an original (pre-cutting) reference stress field, and the extent of cutting-induced plasticity. Comparisons reveal that while the double-embedded cutting configuration produces the most accurate residual stress measurements, the highest levels of plastic flow are generated in this process. As a result, this cutting-induced plastic deformation is, however, largely confined to small ligaments formed as a consequence of the sample sectioning process, and as such it does not significantly affect the back-calculated residual stress field.

36 MATERIALS SCIENCE↗

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↗

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↗