Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “expansion cost”

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 721 records · Page 40

Energy-Storage Modeling: State-of-the-Art and Future Research Directions

Given its physical characteristics and the range of services that it can provide, energy storage raises unique modeling challenges. This paper summarizes capabilities that operational, planning, and resource-adequacy models that include energy storage should have and surveys gaps in extant models. Existing models that represent energy storage differ in fidelity of representing the balance of the power system and energy-storage applications. Modeling results are sensitive to these differences. The importance of capturing chronology can raise challenges in energy-storage modeling. Some models ‘decouple’ individual operating periods from one another, allowing for natural decomposition and rendering the models relatively computationally tractable. Energy storage complicates such a modeling approach. Improving the representation of the balance of the system can have major effects in capturing energy-storage costs and benefits.

25 ENERGY STORAGE↗

Viability Assessment of Wind and Solar Renewable Energy Generation in Support of Nationwide Vehicle Electrification

In 2022, the U.S. transportation sector was the largest source of greenhouse gas emissions in the country, with the combination of passenger and commercial vehicles contributing 80% of these emissions. As adoption of passenger electric vehicles continues to climb, sights are being set on the electrification of heavy-duty commercial vehicle (HDCV) fleets. The sustainability of these shifts relies in part on the addition of significant renewable energy generation resources to both bolster the grid in the face of increased demand, and to prevent a shift in the source of greenhouse gas (GHG) emissions to the grid, as opposed to a true net reduction. Additionally, it is necessary to quantify the variations in economic viability across the country for these technologies as it pertains to their productive capabilities. Doing so will encourage investment and ensure that the transition to electrified HDCV fleets is commercially viable, as well as sustainable. In an effort to meet these goals, multiple computational frameworks are used to locate suitable land for renewable infrastructure development, and to quantify spatiotemporal variations in the potential energy generation and financial viability of development sites across the Unites States. First, the Oak Ridge Siting Analysis for power Generation Expansion tool (OR-SAGE) is used to assess the suitability of land for potential wind and solar energy development across the contiguous U.S. From there, resource data from the National Solar Radiation Database (NSRDB) and the Wind Integration National Dataset (WIND) are used in concert with the National Renewable Energy Laboratory (NREL) Renewable Energy Potential (ReV) model to calculate the variation in potential generation capacity for each resource. Additionally, the capital and operational expenditures are calculated for an example configuration of each renewable technology. These measures are then used to calculate the levelized cost of energy (LCOE) of potential sites. All of these results are then processed and analyzed to determine where in the U.S. solar and wind energy are most viable. This viability is based on available generation potential, consistency and stability of energy generation over time, and economic viability with respect to LCOE.

Miller, Brandon [ORNL] (ORCID:0009000300169201)↗

Life after 30 years - a PV system in Colorado

The PV system consisted of 155 monocrystalline Si Arco M53 panels, installed in 1987. 80 were installed on a fixed tilt rack, and 75 on 5 separate dual-axis passive trackers. The system was decommissioned in 2017 due a campus re-design and building expansion at the Platte River Power Authority (PPRA) site. Prior to decommissioning a site visit established the following failure type observations: EVA browning (100%), cracked cells (60%), and cracked wires or grommets (9%, or more if replacements were made selectively on cracked parts). EVA browning is typical in modules of this vintage, whereas the occurrence of so many cracked cells is surprising and is likely weather-related. The fairly frequent occurrence of cracked wires and grommets serves as a reminder that accelerated testing of these components is important for goals of low-cost 50-year energy generation. In smaller numbers, thermal effects (7%), bent frames (1%), and cell corrosion (<1%) were observed. A sample of 20 modules was selected and brought to NREL for in-depth characterization. Because baseline measurements were not available using the nameplate rating, which can be inaccurate, revealed a mean module sample degradation of 0.7 %/year. Sporadic production data during the 30 years determined a system degradation rate closer to 1.2 %/year. Subsequently, four of these modules were re-mounted outdoors at NREL and are now approaching a total lifetime of 34 years.

durability↗

Potential of porous nodal-line semi-metallic carbon for sodium-ion battery anode

As an alternative to lithium-ion batteries (LIBs), sodium-ion batteries (SIBs) have attracted much attention because sodium is abundant, cost-effective, and environmentally benign. However, due to the larger size of Na-ion compared to the Li-ion, graphite does not work well as an anode for SIBs, as it does for LIBs. In this work, based on first-principles calculations, we show that the recently proposed topological semi-metallic carbon, HZGM-42, has great potential as an anode material for SIBs. With a large capacity of 318.5 mAhg -1 , low diffusion barrier of 0.08–0.21 eV, low average open-circuit voltage of 0.43 V, and a small volume expansion of 1.85%, the all-carbon topological quantum HZGM-42 has high energy density, good rate capability, and excellent cycling stability and can serve as a universal anode material for both LIBs and SIBs.

25 ENERGY STORAGE↗

Implications of a regional resource adequacy program on utility integrated resource planning: Study for the Western United States

LBNL collaborated with the Western Interstate Energy Board and the University of Texas, Austin to investigate the implications of a regional resource adequacy (RA) program on utility integrated resource planning. This report is focused on an active policy discussion to develop a novel voluntary program to share capacity resource and improve RA in the Western Interconnection. This report identifies two RA components of IRP that will be highly impacted by a regional RA program: resource capacity accreditation and RA targets. There are at least four resources that will require specific attention for their capacity credit calculation: (1) variable renewable resources, (2) demand-side resources, (3) hydropower, and (4) contracts. It will then be necessary to decide on a RA target reliability metric (e.g., a planning reserve margin) that is at least the minimum requirement in IRPs to ensure consistency in RA requirement calculations. The report finds that load forecasting and transmission expansion analyses will be moderately impacted, but that most of IRP components will not be significantly impacted. The report includes (i) a review of traditional resource adequacy practices in IRP; (ii) a case study of an existing regional RA program that interacts with IRP (Southwest Power Pool) and (iii) an overview of the NWPP regional resource adequacy proposal that is the object of this work. This paper does not (1) advocate for or against a regional RA program for the NWPP, (2) make detailed design recommendations for this program, or (3) assess its benefits and costs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Titanium Carbide MXene Hole Contacts for CdTe Photovoltaics

CdTe is a high‐efficiency thin‐film photovoltaic technology that has seen tremendous commercial success over the past decade. Yet despite the improvement of other device characteristics, the fabrication of an ohmic hole back contact layer has remained a challenge due to the high ionization potential of CdTe, which limits the external potential that can be feasibly reached even as other characteristics of the device improve. MXenes, a family of 2D materials with rapidly growing scientific and commercial interest, offer a promising route to forming low‐cost, low‐barrier contacts due to their demonstrated high work function, metallic conductivity, and facile solution processing from benign solvents. Here, it is shown that Ti 3 C 2 T x MXene films processed from an aqueous colloidal dispersion can perform as a highly efficient hole contact material for CdTe solar cells, resulting in high power‐conversion efficiencies. The role of the Schottky barrier formation in Ti 3 C 2 T x ‐contacted CdTe devices is probed, and potential pathways for the future development of this potent combination of materials are elucidated. The modularity of the expansive MXene family of materials presents a promising strategy for developing next‐generation hole contacts for CdTe solar cells.

14 SOLAR ENERGY↗

Increasing the power of weak lensing data with multipole-based intrinsic alignment estimators

ABSTRACT It has long been known that galaxy shapes align coherently with the large-scale density field. Characterizing this effect is essential to interpreting measurements of weak gravitational lensing, the deflection of light from distant galaxies by matter overdensities along the line of sight, as it also produces coherent galaxy alignments that we wish to interpret in terms of a cosmological model. Existing direct measurements of intrinsic alignments using galaxy samples with high-quality shape and redshift measurements typically use well-understood but sub-optimal projected estimators, which do not make good use of the information in the data when comparing those estimators to theoretical models. We demonstrate a more optimal estimator, based on a multipole expansion of the correlation functions or power spectra, for direct measurements of galaxy intrinsic alignments. We show that even using the lowest order multipole alone increases the significance of inferred model parameters using simulated and real data, without any additional modelling complexity. We apply this estimator to measurements of parameters of the non-linear alignment model using data from the Sloan Digital Sky survey, demonstrating consistent results with a factor of ∼2 greater precision in parameter fits to intrinsic alignments models. This result is functionally equivalent to quadrupling the survey area, but without the attendant costs – thereby demonstrating the value in using this new estimator in current and future intrinsic alignments measurements using spectroscopic galaxy samples.

Astronomy & Astrophysics↗

Inventory for Crystalline Silicon Module Recycling: Cooperative Research and Development (Final Report)

A critical challenge for the continued expansion of photovoltaics (PV) is to develop technically feasible, inexpensive and environmentally friendly practices for handling and recycling modules at the end of their usable life. The National Renewable Energy Laboratory (NREL) is requested by the Electric Power Research Institute (EPRI) to collect primary data regarding the environmental performance of currently operational PV module recycling facilities in Europe. Very little has been published regarding crystalline silicon (C-Si) module recycling. Thus, much effort will be needed in direct industry outreach, collection of information and other business intelligence strategies similar to NREL's approaches for developing cost models for PV manufacturing. The goal of this work effort is to produce a detailed inventory that accounts for physical (e.g., energy, water, materials) flows through each step of a C-Si recycling process. The inventory (a life cycle inventory, or LCI) shall be designed so that it can be extended to include an accounting of costs for each process step, inputs, etc. This work effort shall leverage prior LCI data collection NREL performed for the United States Department of Energy, Solar Energy Technologies Office, under the auspices of the U.S. contribution to International Energy Agency's Photovoltaics Power Systems (PVPS) Task 12 (Environmental Health and Safety), which SETO nominated NREL to chair. The primary purpose of this work effort is to augment the prior data collection to increase the sample size of manufacturers' primary data in the LCI.

14 SOLAR ENERGY↗

NLSCIDNT user's guide maximum likehood parameter identification computer program with nonlinear rotorcraft model

A nonlinear, maximum likelihood, parameter identification computer program (NLSCIDNT) is described which evaluates rotorcraft stability and control coefficients from flight test data. The optimal estimates of the parameters (stability and control coefficients) are determined (identified) by minimizing the negative log likelihood cost function. The minimization technique is the Levenberg-Marquardt method, which behaves like the steepest descent method when it is far from the minimum and behaves like the modified Newton-Raphson method when it is nearer the minimum. Twenty-one states and 40 measurement variables are modeled, and any subset may be selected. States which are not integrated may be fixed at an input value, or time history data may be substituted for the state in the equations of motion. Any aerodynamic coefficient may be expressed as a nonlinear polynomial function of selected 'expansion variables'.

Source record↗

Subsonic panel method for designing wing surfaces from pressure distribution

An iterative method has been developed for designing wing section contours corresponding to a prescribed subcritical distribution of pressure. The calculations are initialized by using a surface panel method to analyze a baseline wing or wing-fuselage configuration. A first-order expansion to the baseline panel method equations is then used to calculate a matrix containing the partial derivative of potential at each control point with respect to each unknown geometry parameter. In every iteration cycle, the matrix is used both to calculate the geometry perturbation and to analyze the perturbed geometry. The distribution of potential on the perturbed geometry is established by simple linear extrapolation from the baseline solution. The extrapolated potential is converted to pressure by Bernoulli's equation. Not only is the accuracy of the approach good for very large perturbations, but the computing cost of each complete iteration cycle is substantially less than one analysis solution by a conventional panel method.

Bristow, D. R.↗

Micromachined electron tunneling infrared sensors

The development of an improved Golay cell is reported. This new sensor is constructed entirely from micromachined silicon components. A silicon oxynitride (SiO(x)N(y)) membrane is deflected by the thermal expansion of a small volume of trapped gas. To detect the motion of the membrane, an electron tunneling transducer is used. This sensor detects electrons which tunnel through the classically forbidden barrier between a tip and a surface; the electron current is exponentially dependent on the separation between the tip and the surface. The sensitivity of tunneling transducers constructed was typically better than 10(exp -3) A/square root of Hz. Through use of the electron tunneling transducer, the scaling laws which have prevented the miniaturization of the Golay cell are avoided. This detector potentially offers low cost fabrication, compatibility with silicon readout electronics, and operation without cooling. Most importantly, this detector may offer better sensitivity than any other uncooled infrared sensor, with the exception of the original Golay cell.

Kenny, T. W.↗

Fast Multipole Methods for Three-Dimensional N-body Problems

We are developing computational tools for the simulations of three-dimensional flows past bodies undergoing arbitrary motions. High resolution viscous vortex methods have been developed that allow for extended simulations of two-dimensional configurations such as vortex generators. Our objective is to extend this methodology to three dimensions and develop a robust computational scheme for the simulation of such flows. A fundamental issue in the use of vortex methods is the ability of employing efficiently large numbers of computational elements to resolve the large range of scales that exist in complex flows. The traditional cost of the method scales as Omicron (N(sup 2)) as the N computational elements/particles induce velocities at each other, making the method unacceptable for simulations involving more than a few tens of thousands of particles. In the last decade fast methods have been developed that have operation counts of Omicron (N log N) or Omicron (N) (referred to as BH and GR respectively) depending on the details of the algorithm. These methods are based on the observation that the effect of a cluster of particles at a certain distance may be approximated by a finite series expansion. In order to exploit this observation we need to decompose the element population spatially into clusters of particles and build a hierarchy of clusters (a tree data structure) - smaller neighboring clusters combine to form a cluster of the next size up in the hierarchy and so on. This hierarchy of clusters allows one to determine efficiently when the approximation is valid. This algorithm is an N-body solver that appears in many fields of engineering and science. Some examples of its diverse use are in astrophysics, molecular dynamics, micro-magnetics, boundary element simulations of electromagnetic problems, and computer animation. More recently these N-body solvers have been implemented and applied in simulations involving vortex methods. Koumoutsakos and Leonard (1995) implemented the GR scheme in two dimensions for vector computer architectures allowing for simulations of bluff body flows using millions of particles. Winckelmans presented three-dimensional, viscous simulations of interacting vortex rings, using vortons and an implementation of a BH scheme for parallel computer architectures. Bhatt presented a vortex filament method to perform inviscid vortex ring interactions, with an alternative implementation of a BH scheme for a Connection Machine parallel computer architecture.

Koumoutsakos, P.↗

An Integrated Computational Materials Engineering (ICME) Approach to Design Nonlinear Transition Zones Between Dissimilar Metals

Current approaches to designing graded transition joints (GTJs) between dissimilar metals often rely on linear changes in both composition profiles and thickness of each sublayer. This increases fabrication cost and may not be optimal with respect to residual stress or the formation of undesirable phases. Here, in this study, GTJs between P91 ferritic/martensitic steel and 347H austenitic stainless steel were designed using Integrated Computational Materials Engineering (ICME) principles with nonlinear composition and length profiles. Guided by inputs from classical mechanics and CALPHAD predictions of carbon chemical potential, a novel transition zone consisting of five discrete compositions was proposed, with the thickness of each sublayer varying according to a brachistochrone-inspired distribution. In addition to carbon potential gradients, CALPHAD was used to predict coefficients of thermal expansion, which were incorporated into finite element models to evaluate stress evolution. The proposed nonlinear design resulted in a smoother carbon potential gradient, lower carbon depletion at the P91 interface, and a comparable residual stress under long-term thermal exposure, compared to a conventional linear design using ten sublayers with equal thickness. This work introduces a brachistochrone-inspired distribution for GTJ design, offering a general framework for optimizing graded interfaces between dissimilar metals.

Directed Energy Deposition↗

Oxidation performance and limitations of additively manufactured SiC/Si-O-C composites at elevated temperatures

Here, SiC preforms were infiltrated with polycarbosiloxane preceramic polymer to create low-cost SiC—Si x O y C z matrix composites for possible high temperature application. The effect of specimen annealing conditions on the steam and air oxidation resistance was analyzed. Air exposures from 900–1500°C showed stability of the Si x O y C z matrix phase with moderate oxidation resistance. After 900°C exposures, oxidation of the composite is primarily limited to the matrix phase, and only at 1200 and 1500°C was SiC consumption visualized. Steam cycle exposures at 1200°C resulted in rapid linear oxidation of the composites with induced specimen swelling from volume expansion associated with internal oxidation. Specimens annealed in argon were found to be the most oxidation resistant due to increased matrix crystallinity and carbon retention. The results of this work were used to develop clear strategies for improving the high-temperature properties of SiC—Si x O y C z composites.

36 MATERIALS SCIENCE↗

Mission Options Scoping Tool for Mars Orbiters: Mass Cost Calculator (MC2)

Prior to developing the details of an advanced mission study, the mission architecture trade space is typically explored to assess the scope of feasible options. This paper describes the main features of an Excel-based tool, called the Mass-Cost-Calculator (MC2 ), which is used to perform rapid, high-level mass and cost options analyses of Mars orbiter missions. MC2 consists of a combination of databases, analytical solutions, and parametric relationships to enable quick evaluation of new mission concepts and comparison of multiple architecture options. The tool's outputs provide program management and planning teams with answers to "what if" queries, as well as an understanding of the driving mission elements, during the pre-project planning phase. These outputs have been validated against the outputs generated by the Advanced Projects Design Team (Team X) at NASA's Jet Propulsion Laboratory (JPL). The architecture of the tool allows for future expansion to other orbiters beyond Mars, and to non-orbiter missions, such as those involving fly-by spacecraft, probes, landers, rovers, or other mission elements.

orbiters↗

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗

Simulations and symmetries

ABSTRACT We investigate the range of applicability of a model for the real-space power spectrum based on N-body dynamics and a (quadratic) Lagrangian bias expansion. This combination uses the highly accurate particle displacements that can be efficiently achieved by modern N-body methods with a symmetries-based bias expansion which describes the clustering of any tracer on large scales. We show that at low redshifts, and for moderately biased tracers, the substitution of N-body-determined dynamics improves over an equivalent model using perturbation theory by more than a factor of two in scale, while at high redshifts and for highly biased tracers the gains are more modest. This hybrid approach lends itself well to emulation. By removing the need to identify haloes and subhaloes, and by not requiring any galaxy-formation-related parameters to be included, the emulation task is significantly simplified at the cost of modelling a more limited range in scale.

79 ASTRONOMY AND ASTROPHYSICS↗

Extraction of Oxygen from the Martian Atmosphere

A mechanical process was designed for direct extraction of molecular oxygen from the martian atmosphere based on liquefaction of the majority component, CO2, followed by separation of the lower-boiling components. The atmospheric gases are compressed from about 0.007 bar to 13 bar and then cooled to liquefy most of the CO2. The uncondensed gases are further compressed to 30 bar or more, and then cooled again to recover water as ice and to remove much of the remaining CO2. The final gaseous products consisting mostly of nitrogen, oxygen, and carbon monoxide are liquefied and purified by cryogenic distillation. The liquefied CO2 is expanded back to the low-pressure atmosphere with the addition of heat to recover a majority of the compression energy and to produce the needed mechanical work. Energy for the process is needed primarily as heat to drive the CO2-based expansion power system. When properly configured, the extraction process can be a net producer of electricity. The conceptual design, termed 'MARRS' for Mars Atmosphere Resource Recovery System, was based on the NASA/JSC Mars Reference Mission (MRM) requirement for oxygen. This mission requires both liquid oxygen for propellant, and gaseous oxygen as a component of air for the mission crew. With single redundancy both for propellant and crew air, the oxygen requirement for the MRM is estimated at 5.8 kg/hr. The process thermal power needed is about 120 kW, which can be provided at 300-500 C. A lower-cost nuclear reactor made largely of stainless steel could serve as the heat source. The chief development needed for MARRS is an efficient atmospheric compression technology, all other steps being derived from conventional chemical engineering separations. The conceptual design describes an exceptionally low-mass compression system that can be made from ultra-lightweight and deployable structures. This system adapts to the rapidly changing martian environment to supply the atmospheric resource to MARRS at constant conditions.

England, C.↗