Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Equivalence 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 181 records · Page 10

Performance Improvements for the Griffin Transport Solvers

Griffin is a Multiphysics Object-Oriented Simulation Environment based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. Griffin includes a variety of deterministic radiation transport solvers for fixed source, k-eigenvalue, adjoint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. A code assessment performed in FY-20 identified two significant issues with the transport solvers in Griffin: first, the primary heterogeneous SN (discrete ordinates) transport solver based on continuous finite element methods required significant mesh refinement and higher memory usage compared to solvers based on the method of characteristic for equivalent accuracy. Second, the homogeneous PN (spherical harmonics expansion) transport solver did not adequately support polynomial refinement, which is a feature usually required for problems with spatial homogenization and pronounced streaming, typical in fast or gas-cooled reactor systems. To address the first issue, the development effort focused on the more promising discontinuous finite element method (DFEM)-based SN transport solver in Griffin. The addition of an asynchronous parallel transport sweeper and coarse mesh finite difference (CMFD) acceleration have rendered a superior heterogeneous SN transport capability for multiphysics problems that requires far less computing resources in terms of both CPU time and memory usage. This is demonstrated with typical thermal- and fast-spectrum reactor benchmark problems, including 2D Transient Reactor Test, 3D Advanced Burner Test Reactor (ABTR), and 2D and 3D Empire microreactor. For the second issue, the development effort focused on a new transport solver based on the hybrid finite element PN method (HFEM-PN), equivalent to the variational nodal method, as well as a new diffusion solver based on HFEM-Diffusion. This solver is intended for homogenized domains with multiphysics coupling (i.e., supports mesh displacement, seamless temperature feedback, etc.). Initial calculations with the HFEM-Diffusion implementation show very good parallel efficiency for the residual evaluations with the 2D ABTR benchmark. A future development effort will be centered on further improvements to the CMFD, HFEM-PN, and DFEM diffusion solvers to ensure Griffin meets performance and software quality assurance requirements for advanced reactor design and analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Historical Perspective on Fast Coupled Loads Analysis Methods

Spacecraft structural designs are typically verified through a coupled loads analysis (CLA) process, which couples the spacecraft model with the launch vehicle (LV) model to predict low-frequency quasi-static and dynamic responses. The CLA calculations are typically the responsibility of the LV organization, but the spacecraft organization has a vested interest in being able to calculate approximate CLA results during the design of the spacecraft. Because of this, there has long been interest in a method that would allow a spacecraft organization to perform a CLA without access to the full set of LV models and forcing functions. One such method is the Norton-Thevenin Receptance Coupling (NTRC) approach, which is specifically designed to accurately transform LV free accelerations (no payload) into coupled system accelerations (LV plus payload). The purpose of this report is to provide historical context for the NTRC method and compare it with methods that have been used in the past. In particular, it is compared to a frequency-domain substitution method that had been used for a long period of time at the Jet Propulsion Laboratory, and a component-mode-based equivalent to that method.

Blelloch, Paul A.↗

Recursive Filtering And Smoothing In Robot Dynamics

Techniques developed originally for electronic systems also useful for multibody mechanical systems. Report summarizes methods developed to solve nonlinear forward-dynamics problem for robot of multiple-link arms connected by joints. Primary objective to show equivalence between recursive methods of dynamical analysis and some filtering and smoothing techniques from state-estimation theory.

Rodriguez, Guillermo↗

Methodology for Sensitivity Analysis, Approximate Analysis, and Design Optimization in CFD for Multidisciplinary Applications

An incremental iterative formulation together with the well-known spatially split approximate-factorization algorithm, is presented for solving the large, sparse systems of linear equations that are associated with aerodynamic sensitivity analysis. This formulation is also known as the 'delta' or 'correction' form. For the smaller two dimensional problems, a direct method can be applied to solve these linear equations in either the standard or the incremental form, in which case the two are equivalent. However, iterative methods are needed for larger two-dimensional and three dimensional applications because direct methods require more computer memory than is currently available. Iterative methods for solving these equations in the standard form are generally unsatisfactory due to an ill-conditioned coefficient matrix; this problem is overcome when these equations are cast in the incremental form. The methodology is successfully implemented and tested using an upwind cell-centered finite-volume formulation applied in two dimensions to the thin-layer Navier-Stokes equations for external flow over an airfoil. In three dimensions this methodology is demonstrated with a marching-solution algorithm for the Euler equations to calculate supersonic flow over the High-Speed Civil Transport configuration (HSCT 24E). The sensitivity derivatives obtained with the incremental iterative method from a marching Euler code are used in a design-improvement study of the HSCT configuration that involves thickness. camber, and planform design variables.

Taylor, Arthur C., III↗

Current distribution in a three-dimensional IC analyzed by a perturbation method. Part 1: A simple steady state theory

The steady state current distribution in a three dimensional integrated circuit is presented. A device physics approach, based on a perturbation method rather than an equivalent lumped circuit approach, is used. The perturbation method allows the various currents to be expressed in terms of elementary solutions which are solutions to very simple boundary value problems. A Simple Steady State Theory is the subtitle because the most obvious limitation of the present version of the analysis is that all depletion region boundary surfaces are treated as equipotential surfaces. This may be an adequate approximation in some applications but it is an obvious weakness in the theory when applied to latched states. Examples that illustrate the use of these analytical methods are not given because they will be presented in detail in the future.

Edmonds, Larry D.↗

Cross Section Generation Capability in Griffin

The Griffin code is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. The code includes a variety of steady-state solvers for fixed-source, k-eigenvalue, adjoint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. The code reads multigroup cross sections in the ISOXML format generated from external deterministic or Monte Carlo cross section generation codes. The implementation of the cross section generation capability in Griffin was initiated last year by plugging in the cross section application programming interface (CSAPI) and reviewing the methodologies for treating particulate fuels. The focus this year was on improving the CSAPI integration and implementing advanced self-shielding methods for applications to advanced reactor problems with TRISO fuels. First, the process for cross section library generation was updated to accurately and rigorously produce isotopic cross section data. Second, the equivalent Dancoff factor cell method performing slowing down calculations on the fly for the resonance treatment was implemented in CSAPI to improve the accuracy of effective multigroup cross sections in the resonance energy range. Third, the iterative local spatial self-shielding method was implemented under the calculation framework of the equivalent Dancoff factor cell method to accurately deal with the double heterogeneity effect of particulate fuel. The updated CSAPI with the advanced self-shielding methods, together with the cross section libraries generated based on the improved process, were tested for pin-cell, unit-cell, and fuel assembly problems with various resonance self-shielding conditions based on very high temperature reactor, high temperature test reactor, and Empire benchmark cores, indicating that the updated CSAPI in Griffin is able to produce multigroup cross sections accurately and efficiently. We also showed that the methodology worked well for pebble bed fuel from HTR-10, but the capability still needs to be fully integrated into CSAPI. In the future, further benchmark tests will be performed for various thermal reactor core problems, including particulate fuel-based pebble bed reactors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Managing Solar Photovoltaic Integration in the Western United States: Resource Adequacy Considerations

This study examines the impact of reserve margin-based reliability assessment, as commonly used in capacity expansion models, on planning resource-adequate power systems under high penetrations of solar photovoltaics (PV). As a generation resource, PV is operationally different from the conventional dispatchable resources for which most capacity expansion models were designed. The question this study attempts to answer is whether large amounts of PV on a system (in this case, the Western Interconnection of North America) would bias the results of conventional reserve margin-based capacity expansion modeling towards an over- or under-provisioning of resource adequacy. This analysis used NREL’s Resource Planning Model (RPM) for capacity expansion modeling and NREL’s Probabilistic Resource Adequacy Suite (PRAS) for resource adequacy assessment. RPM uses a reserve margin requirement to enforce resource adequacy. PRAS, a collection of tools for studying the resource adequacy of power systems and the adequacy contributions of individual resources on a probabilistic basis, was used to compute multiple resource adequacy metrics across a number of simulated scenarios and system representations with differing regional detail. In all cases, including high PV penetrations (up to 33% annual generation from PV, interconnection-wide), RPM was able to produce resource-adequate systems as measured by normalized expected unserved energy and loss-of-load expectation results from PRAS. The accuracy of reserve margin approaches depends heavily on the underlying assumptions informing the capacity credit assigned to variable and energy-limited resources, particularly when such resources are abundant in the modeled system. RPM’s standard methodology for estimating variable and flexible resources’ capacity contributions, which is based on the top 100 hours of net load, did not appear to systematically undervalue or overvalue variable generation relative to a more rigorous equivalent firm capacity assessment using PRAS, although both over- and under-valuations were observed in specific scenarios. In the worst cases, the top 100 hour method underestimated the equivalent firm capacity of PV by two percentage points, and overestimated the equivalent firm capacity of PV by five percentage points. Calculating capacity contributions based on the top 10 hours of net load systematically underestimated equivalent firm capacities at more modest PV penetrations, but was often a better approximation of equivalent firm capacity than the existing 100-hour approach in scenarios with higher PV penetrations.

14 SOLAR ENERGY↗

Determination of the Fracture Parameters in a Stiffened Composite Panel

A modified J-integral, namely the equivalent domain integral, is derived for a three-dimensional anisotropic cracked solid to evaluate the stress intensity factor along the crack front using the finite element method. Based on the equivalent domain integral method with auxiliary fields, an interaction integral is also derived to extract the second fracture parameter, the T-stress, from the finite element results. The auxiliary fields are the two-dimensional plane strain solutions of monoclinic materials with the plane of symmetry at x(sub 3) = 0 under point loads applied at the crack tip. These solutions are expressed in a compact form based on the Stroh formalism. Both integrals can be implemented into a single numerical procedure to determine the distributions of stress intensity factor and T-stress components, T11, T13, and thus T33, along a three-dimensional crack front. The effects of plate thickness and crack length on the variation of the stress intensity factor and T-stresses through the thickness are investigated in detail for through-thickness center-cracked plates (isotropic and orthotropic) and orthotropic stiffened panels under pure mode-I loading conditions. For all the cases studied, T11 remains negative. For plates with the same dimensions, a larger size of crack yields larger magnitude of the normalized stress intensity factor and normalized T-stresses. The results in orthotropic stiffened panels exhibit an opposite trend in general. As expected, for the thicker panels, the fracture parameters evaluated through the thickness, except the region near the free surfaces, approach two-dimensional plane strain solutions. In summary, the numerical methods presented in this research demonstrate their high computational effectiveness and good numerical accuracy in extracting these fracture parameters from the finite element results in three-dimensional cracked solids.

Lin, Chung-Yi↗

Equalizing the Pixel Response of the Imaging Photoelectric Polarimeter Onboard the IXPE Mission

The Gas Pixel Detector is a gas detector, sensitive to the polarization of X-rays, currently flying onboard the Imaging X-ray Polarimetry Explorer (IXPE)—the first observatory dedicated to X-ray polarimetry. It detects X-rays and their polarization by imaging the ionization tracks generated by photoelectrons absorbed in the sensitive volume, and then reconstructing the initial direction of the photoelectrons. The primary ionization charge is multiplied and ultimately collected on a finely pixellated ASIC specifically developed for X-ray polarimetry. The signal of individual pixels is processed independently and gain variations can be substantial, of the order of 20%. Such variations need to be equalized to correctly reconstruct the track shape, and therefore its polarization direction. The method to do such equalization is presented here and is based on the comparison between the mean charge of a pixel with respect to the other pixels for equivalent events. The method is shown to finely equalize the response of the detectors onboard IXPE, allowing a better track reconstruction and energy resolution, and can in principle be applied to any imaging detector based on tracks.

79 ASTRONOMY AND ASTROPHYSICS↗

Dark matter and the equivalence principle

A survey is presented of the current understanding of dark matter invoked by astrophysical theory and cosmology. Einstein's equivalence principle asserts that local measurements cannot distinguish a system at rest in a gravitational field from one that is in uniform acceleration in empty space. Recent test-methods for the equivalence principle are presently discussed as bases for testing of dark matter scenarios involving the long-range forces between either baryonic or nonbaryonic dark matter and ordinary matter.

Frieman, Joshua A.↗

Combined, nonlinear aerodynamic and structural method for the aeroelastic design of a three-dimensional wing in supersonic flow

An iterative procedure for the static aeroelastic design of a flexible wing at supersonic speeds has been developed. The procedure combines a nonlinear, full-potential solver (NCOREL) with an equivalent plate structural analysis method. The NCOREL method yields significantly improved aerodynamic estimates compared to linear theory. The equivalent plate structural analysis method demonstrates an order of magnitude reduction in computer memory and execution time compared to finite-element methods. A highly swept wing is analyzed at high lift using this aeroelastic procedure. The results indicate that the wing deforms favorably due to aerodynamic loading and, consequently, that the inviscid drag levels do not vary at the required lift coefficient although the angle of attack varies significantly. A sensitivity analysis of the type required for optimization studies was also performed with the aeroelastic design procedure.

Pittman, J. L.↗

Determination of parameters of a new method for predicting alloy properties

Recently, a semiempirical method for alloys based on equivalent crystal theory was introduced. The method successfully predicts the concentration dependence of the heat of formation and lattice parameter of binary alloys. A study of the parameters of the method is presented, along with new results for (gamma)Fe-Pd and (gamma)Fe-Ni alloys.

Bozzolo, Guillermo↗

Modification of a Macromechanical Finite-Element Based Model for Impact Analysis of Triaxially-Braided Composites

A macro level finite element-based model has been developed to simulate the mechanical and impact response of triaxially-braided polymer matrix composites. In the analytical model, the triaxial braid architecture is simulated by using four parallel shell elements, each of which is modeled as a laminated composite. For the current analytical approach, each shell element is considered to be a smeared homogeneous material. The commercial transient dynamic finite element code LS-DYNA is used to conduct the simulations, and a continuum damage mechanics model internal to LS-DYNA is used as the material constitutive model. The constitutive model requires stiffness and strength properties of an equivalent unidirectional composite. Simplified micromechanics methods are used to determine the equivalent stiffness properties, and results from coupon level tests on the braided composite are utilized to back out the required strength properties. Simulations of quasi-static coupon tests of several representative braided composites are conducted to demonstrate the correlation of the model. Impact simulations of a represented braided composites are conducted to demonstrate the capability of the model to predict the penetration velocity and damage patterns obtained experimentally.

Goldberg, Robert K.↗

Smart Ventilation Controls for Occupancy and Auxiliary Fan Use Across U.S. Climates

Smart Ventilation has been developed as a way to reduce the energy associated with ventilation by changing when ventilation happens and how much ventilation occurs at any given time. In high performance buildings with low envelope and appliance related loads, ventilation is becoming a bigger part of the total building energy use and needs to be addressed if performance targets are to be met. This work explores the development and performance of smart ventilation controls based on occupancy and auxiliary fan operation that provide annual dwelling unit ventilation equivalence to ASHRAE Standard 62.2-2016. A prototype high performance home compliant with U.S. DOE Building America Zero Energy Ready program requirements was simulated using the REGCAP tool with and without smart controls across 15 U.S. DOE climate zones. Balanced and unbalanced IAQ fans were independently simulated, and all smart controlled fans had airflows double the reference 62.2-2016 airflow. Three idealized occupancy patterns were examined: 1st shift (a typical daily work/school absence), an extended 1st shift with more time spent out of the home evenings and weekends and 3rd shift (night work). The Occupancy controller turns the IAQ fan off during unoccupied periods, and it resumes ventilation upon their return. While unoccupied, contaminants are allowed to accumulate in the space, because occupants are not exposed to these contaminants. When occupants return home, they are exposed to these higher contaminant concentrations, and the controller increases the ventilation rate sufficiently to ensure equivalence with a continuous IAQ fan. Accounting for pollutant emissions that occur during unoccupied periods (as required by ASHRAE 62.2-2016), sharply distinguishes our occupancy controls from past Demand Controlled Ventilation systems. This new accounting method results in equivalent contaminant exposure, as well as lower reductions in average ventilation rates and lower energy savings. Smart controls were demonstrated that saved HVAC energy (i.e., avoided heating/cooling load, as well as fan energy)—averaging between 6 and 46% of ventilation-related energy use depending on the control strategy and occupancy pattern assumptions. The greatest savings were in the combined Auxiliary Fan + Occupancy control. Energy savings increased with climate heating demand and longer unoccupied time periods. The 3rd shift occupancy pattern had better performance, due to the thermal benefit of reducing the ventilation rate during the cold nighttime hours. This same effect provided a cooling energy benefit in hotter locations for the 1st shift. Overall, savings from occupancy-based smart controls were low, because of the recovery period required after occupants return home, during which the airflow is double the 62.2 reference. This recovery is required to maintain equivalence with the ASHRAE standard. Occupancy-based control performance was improved when combined with sensing auxiliary fans and when providing a pre-occupancy flush out of one- or two-hours. Performance was similarly improved if the ASHRAE Standard were to recognize that pollutant emissions are lower during unoccupied periods, iii allowing a lower target ventilation rate during unoccupied periods (not currently in the standard) (see Full vs. Half AEQ in this report). Finally, over-sized unbalanced fans that are cycled on-and-off by a smart controller (or timer) were found to substantially increase annual average air exchange and energy use relative to a continuous unbalanced fan due to the effects of superposition with natural infiltration.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Knowledge-Based Hazardous Waste Determinations for Solar Photovoltaic (PV) Modules: A Foundational Study

This report explores using "generator knowledge" to determine whether a solar photovoltaic (PV) module must be managed as hazardous waste prior to recycling or landfilling. Generator knowledge is a legal term and existing regulatory pathway for making a hazardous waste determination that has been used by other industries but is a relatively unknown option to the PV industry. In the United States, a hazardous waste determination often acts as a pre-requisite to recycle or landfill a PV module. The results of the hazardous waste determination dictate whether the PV module must be managed as hazardous waste or nonhazardous solid waste. Managing a PV module as hazardous waste requires compliance with stringent U.S. federal and state hazardous waste law. In addition to increased management costs, which can be ten times higher, legal liability for PV modules regulated as hazardous waste is also heightened with both civil and criminal penalties for noncompliance which includes making an inaccurate or faulty hazardous waste determination. The most common reason a PV module would be regulated as hazardous is if it contains a regulated metal in an amount that equals or exceeds the toxicity characteristic limits. To determine whether a PV module exhibits a hazardous characteristic, the regulated person/entity must "apply knowledge...in light of the materials and processes used." In the absence of adequate knowledge to determine whether the PV module is hazardous, it must be tested using Test Method 1311 Toxicity Characteristic Leaching Procedure (TCLP) or an equivalent EPA-approved method. Although TCLP is the predominant method used to today to make a hazardous waste determination for PV modules, evidence from this study concludes it is not practical to TCLP test every PV module even in a single utility-scale installation, and a scalable solution is needed. This study finds that knowledge-based hazardous waste determinations may allow a regulated person/entity to make a hazardous waste determination for more than one PV module at one time - making this regulatory pathway a potential scalable solution. Through legal analysis and interviews with 44 experts, the authors explore what it means to make an accurate knowledge-based hazardous determination for PV modules considering sources and forms of information as well as potential limitations. The work aims to provide a foundation for building consensus on whether knowledge-based hazardous determinations are a viable, scalable industry approach for solar.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Radar cross section of an open-ended circular waveguide Calculation of second-order diffraction terms

The radar cross section of an open-ended, semi-infinite circular waveguide is calculated by the geometrical theory of diffraction and the equivalent-current (EC) method. Both single- and double-diffraction terms are included. It is found that the double-diffraction term is stronger than the single-diffraction term for the horizontal polarization at wide-angle incidence. Its inclusion is necessary in order to check with experimental data and the asymptotic expansion of the exact Wiener-Hopf solution. The EC method is used for the single diffraction. On the other hand the double-diffraction term is obtained by using the EC method for axial incidence and geometrical theory of diffraction for wide-angle incidence. These two solutions are matched with a proper matching function. This technique is computationally more efficient than the matching technique of integration over the equivalent currents.

Lee, C. S.↗

Determining blood and plasma volumes using bioelectrical response spectroscopy

We hypothesized that an electric field (inductance) produced by charged blood components passing through the many branches of arteries and veins could assess total blood volume (TBV) or plasma volume (PV). Individual (N = 29) electrical circuits (inductors, two resistors, and a capacitor) were determined from bioelectrical response spectroscopy (BERS) using a Hewlett Packard 4284A Precision LCR Meter. Inductance, capacitance, and resistance from the circuits of 19 subjects modeled TBV (sum of PV and computed red cell volume) and PV (based on 125I-albumin). Each model (N = 10, cross validation group) had good validity based on 1) mean differences (-2.3 to 1.5%) between the methods that were not significant and less than the propagated errors (+/- 5.2% for TBV and PV), 2) high correlations (r > 0.92) with low SEE (< 7.7%) between dilution and BERS assessments, and 3) Bland-Altman pairwise comparisons that indicated "clinical equivalency" between the methods. Given the limitation of this study (10 validity subjects), we concluded that BERS models accurately assessed TBV and PV. Further evaluations of the models' validities are needed before they are used in clinical or research settings.

NASA Discipline Number 00-00↗