Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reduced Order Models”

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,459 records · Page 81

Raytracing the long trace profiler

The optical design of the Long Trace Profiler optical system is explored with a commercial raytrace program, Zemax OpticStudio™ (ZOS)1, with the intent of finding and correcting sources of systematic error. ZOS provides both geometric raytracing tools and physical optics Gaussian beam propagation and diffraction image calculation tools, and two design modes, sequential component (SC) and non-sequential component (NSC) that are optimized for different aspects of the design process. The original LTP-II system employs a singlet lens with a 1250 mm focal length. It is optimized to provide minimum distortion over a surface slope angle range of ±5 mrad. Using the ZOS tools, we are able to simulate ghost rays that produce distortion in the beam spot image and can minimize the distortion by deliberate misalignment of the beamsplitter (BS) components. Unfortunately, the reference beam is compromised because of the component tilts. The most recent LTP500 system design simplifies the optical system and makes the reference beam usable again, even with misalignment of the polarizing beamsplitter (PBS). Two lenses are designed for the LTP500 - a cemented doublet that has been fabricated, and a singlet with one aspheric surface. Both have focal distances of 500 mm with an expanded angular measurement range of ±10 mrad. The aspheric singlet provides superior performance. ZOS allows the import of wavefront measurement data produced by commercial interferometer software. We apply the wavefront error measurement from the cemented doublet to the model to show that the 19nm RMS wavefront error needs to be improved by at least a factor of ten in order to reduce the systematic error to a level that will allow the LTP to approach its design limit of a few tens of nanoradians.

Takacs, Peter Z.↗

Design Considerations and Analysis of Multi-Level Erasure Coding in Large-Scale Data Centers

Multi-level erasure coding (MLEC) has seen large deployments in the field, but there is no in-depth study of design considerations for MLEC at scale. In this paper, we provide comprehensive design considerations and analysis of MLEC at scale. We introduce the design space of MLEC in multiple dimensions, including various code parameter selections, chunk placement schemes, and various repair methods. We quantify their performance and durability, and show which MLEC schemes and repair methods can provide the best tolerance against independent/correlated failures and reduce repair network traffic by orders of magnitude. To achieve this, we use various evaluation strategies including simulation, splitting, dynamic programming, and mathematical modeling. We also compare the performance and durability of MLEC with other EC schemes such as SLEC and LRC and show that MLEC can provide high durability with higher encoding throughput and less repair network traffic over both SLEC and LRC.

Wang, Meng↗

A method for high order linear system reduction and nonlinear system simplification

Least-squares-type algorithms for reducing the order of linear systems in the frequency domain and simplifying nonlinear systems in time domain are developed and demonstrated. The possible model structures are represented as nodes in a tree, and costs along the branches are assigned using the repeated-Gram-Schmidt orthogonalization procedure of Desrochers and Saridis (1980), permitting identification of the optimal n-term model by searching the tree to depth n, with no need for parameter identification. The efficiency and flexibility of the algorithms is shown in applications to the eighth-order linear system studied by Hsia (1972), a three-state eight-nonlinear-term aircraft-dynamics problem, and the related linear-controller problem (Garrard and Jordan, 1977).

Desrochers, A. A.↗

A survey of the state of the art and focused research in range systems, task 2

Many communication, control, and information processing subsystems are modeled by linear systems incorporating tapped delay lines (TDL). Such optimized subsystems result in full precision multiplications in the TDL. In order to reduce complexity and cost in a microprocessor implementation, these multiplications can be replaced by single-shift instructions which are equivalent to powers of two multiplications. Since, in general, the obvious operation of rounding the infinite precision TDL coefficients to the nearest powers of two usually yield quite poor system performance, the optimum powers of two coefficient solution was considered. Detailed explanations on the use of branch-and-bound algorithms for finding the optimum powers of two solutions are given. Specific demonstration of this methodology to the design of a linear data equalizer and its implementation in assembly language on a 8080 microprocessor with a 12 bit A/D converter are reported. This simple microprocessor implementation with optimized TDL coefficients achieves a system performance comparable to the optimum linear equalization with full precision multiplications for an input data rate of 300 baud. The philosophy demonstrated in this implementation is dully applicable to many other microprocessor controlled information processing systems.

Yao, K.↗

ORBIT: Offshore Renewables Balance-of-System and Installation Tool

This report describes the Offshore Renewables Balance-of-system Installation Tool ("ORBIT"), a new model developed by the National Renewable Energy Laboratory (NREL) to evaluate the balance-of-system (BOS) costs associated with offshore wind projects. In the context of wind energy projects BOS costs encompass all expenses required to construct a project other than the capital expenditures (CapEx) of the turbines and towers, including the procurement costs for all other components (such as substructures, cables, and electrical infrastructure), offshore and land-based construction costs, port costs, site surveying fees, permitting fees, and leasing fees are all categorized as BOS costs. BOS costs significantly contribute to the levelized cost of energy (LCOE) for offshore wind, typically comprising over 50\% of the CapEx for a fixed-bottom offshore wind project and 60\% for a floating project. In addition, technology solutions and installation methods vary drastically between projects as they are impacted by factors such as vessel availability, geographic considerations and site geotechnical conditions. These effects require a model with appropriate fidelity to understand how these costs scale as turbine rating increases and the offshore wind supply chain, particularly offshore construction vessels, is expanded. For offshore wind, cost savings attributed to increased turbine rating are primarily realized through BOS procurement and project installation as fewer substructures and less cable are required. BOS costs represent both a modeling challenge, as well as an opportunity, for project developers to optimize solutions to reduce costs. It is critical to understand how these costs are affected by novel technologies, innovative installation processes, and operational constraints in order identify meaningful cost reductions for offshore wind energy.

17 WIND ENERGY↗

Modeling of a 43 GHz microplasma within a photonic crystal

The electron density (n e ) and argon metastable density (1 s 5 ) of a 43GHz microplasma are obtained using the zero-dimensional plasma kinetics solver (ZDPlasKin®) for an incident power of 160 and 1000mW and argon pressures of 10–600Torr (1.3 ×10 3 to 8.0 × 10 4 Pa). These simulations are compared with previously published experimental data. To determine the self-consistent electric field in the modeled plasma, the three-dimensional millimeter wave fields are computed as a function of electron density using ANSYS EM19.2, HFSS®. This electromagnetic field model is coupled to ZDPlasKin such that any increase in the simulated plasma density correctly attenuates the simulated electric field within the microplasma. The electron density is found to be sensitive to argon gas temperature, so a two-zone temperature model was needed to obtain agreement with experimental measurements. The temperature in the central core of the microplasma was determined by a previous experimental study. That temperature was used as an input to the model for the simulation of volume recombination losses. The outer regions of the microplasma are assumed to be in equilibrium with the walls (300K). This second temperature was used in the model to determine diffusion losses. The modeled electron and metastable densities are of the order 10 20 and 10 18 m –3 , respectively. This is in good agreement with those measured experimentally as long as the two-zone temperature model is used. In the absence of a hot gas core, the modeled three-body recombination rates are excessive and the simulation severely under-predicts the electron density and over-estimates the metastable density. We conclude that the millimeter wave microplasma has a hot core (2500K at 600 Torr) that rarifies the argon gas and effectively reduces three-body recombination. Furthermore, this allows one to achieve high electron density on the order of 10 20 m –3 with only 100mW of wave power.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electrochemical Model Assessment of Strategies for Fast Charge

The United States Department of Energy (DOE) has identified extreme fast charging (XFC) as a critical challenge that must be overcome in order to achieve widespread adoption of electric vehicles. DOE’s eXtreme Fast Charge and Cell Evaluation (XCEL) program seeks to reduce charge time to 10-15 minutes for a 200+ Wh/kg Li-ion battery. This presentation for DOE’s Vehicle Technologies Program Annual Merit Review summarizes the XCEL team’s work developing electrochemical models to interpret data and applying them to identify best strategies to achieve fast charge. Strategies include advanced electrode architectures, elevated temperature charging, novel charge protocols, advanced electrolytes, thermal management and suppression of heterogeneities that cause early onset of the damaging Li plating side reaction.

47 OTHER INSTRUMENTATION↗

MULTI-DIMENSIONAL MODELING OF THE CFR ENGINE FOR THE INVESTIGATION OF SI NATURAL GAS COMBUSTION AND CONTROLLED END-GAS AUTOIGNITION

Engine knock and misfire are barriers to pathways leading to high-efficiency Spark-Ignited (SI) Natural Gas engines. The general tendency to knock is highly dependent on engine operating conditions and the fuel reactivity. The problem is further complicated by low emission limits and the wide range of chemical reactivity in pipeline quality natural gas. Depending on the region and the source of the natural gas, its reactivity, described by its methane number (analogous to the octane number for liquid SI fuels) can span from 65 - 95. In order to realize diesel-like efficiencies, SI natural gas engines must be designed to operate at high BMEP near knock limits over a wide range of fuel reactivity. This requires a deep understanding regarding the combustion-engine interactions pertaining to flame propagation and end-gas autoignition (EGAI). However, EGAI, if controlled, provides an opportunity to increase SI natural gas engine efficiency by increasing combustion rate and the total burned fuel, mitigating the effects of the slow flame speeds of natural gas fuels which generally reduce BMEP and increase unburned hydrocarbon emissions. For this reason, in order to study EGAI phenomenon, the present work highlights multi-dimensional computational fluid dynamics (CFD) models of the Cooperative Fuel Research (CFR) engine. The CFR engine models are used to investigate fuel-engine interactions that lead to EGAI with natural gas, including effects of fuel reactivity, engine operating parameters, and exhaust gas recirculation (EGR). A Three-Pressure Analysis, performed with GT-Power, was used to estimate initial and boundary conditions for the three-dimensional CFD model. CONVERGE CFD v2.4 was used for the three-dimensional CFD modeling where the level set G-Equation model and SAGE detailed chemical kinetics solver were used. An assessment of the different modeling approaches is also provided to evaluate their limitations, advantages and disadvantages, and for which situations they are most applicable. Model validation was performed with experimental data taken with a CFR engine over varying compression ratio, CA50, EGR fraction, and IMEP and shows good agreement in Peak Cylinder Pressure (PCP), PCP crank angle, and the location of the 10%, 50%, and 90% mass fraction burned (CA10, CA50, and CA90, respectively). The models can predict the onset crank angle and pressure rise rate for light, medium, and heavy EGAI under a variety of fuel reactivities and engine operating conditions.

Bestel, Diego↗

Reducing Memory Consumption in Calico with Shared Memory

This document details the work to reduce memory consumption in Calico. Calico is SimTools’ Constructive Solid Geometry (CSG) and geometry painting library. It is primarily used to paint material volume fractions in the Eulerian meshes of the physics codes. Calico provides point-in-body checks for the geometry supplied by an Oso model, which are then aggregated by the host codes. In addition, Calico can be used to build Oso models and is used by Ingen for that purpose. Oso models, and thus Calico, provide support for various CSG primitives such as spheres, cylinders, surfaces generated by rotating tabular curve data, and STL files as well as binary combinations of those primitives. Prior to refactoring Calico will run out of memory on CTS-1 machines when 36 MPI ranks are used per node when reading STL models on the order of 1.5 GB. This limitation is a bottleneck in designer workflow. This problem has been alleviated through the use of data structures to both reduce memory consumption and to leverage MPI-3 shared memory. This report details the data structures targeted for refactoring in Calico, the methods and implementation details for reducing memory consumption and leveraging shared memory, and results for one test problem. Results show a memory reduction when loading a 1 GB STL file by a factor of 27.5, from 93.4 to 3.4 GB.

97 MATHEMATICS AND COMPUTING↗

Viscoplasticity based on overstress with a differential growth law for the equilibrium stress

Two coupled, nonlinear differential equations are proposed for the modeling of the elastic and rate (time) dependent inelastic behavior of structural metals in the absence of recovery and aging. The structure of the model is close to the unified theories but contains essential differences. It is shown that the model reproduces almost elastic regions upon initial loading and in the unloading regions of the hysteresis loop. Under loading, unloading and reloading in strain control the model simulated the experimentally observed sharp transition from nearly elastic to inelastic behavior. When a formulation akin to existing unified theories is adopted the almost elastic regions reduce the points and the transition upon reloading is very gradual. For different formulations the behavior under sudden in(de)creases of the strain rate by two orders of magnitude is simulated by numerical experiments and differences are noted. The model represents cyclically neutral behavior and contains three constants and two positive, decreasing functions. The determination of constants and functions from monotonic loading with strain rate changes and relaxation periods is described.

Krempl, E.↗

Viscoplasticity based on overstress with a differential growth law for the equilibrium stress

Two coupled, nonlinear differential equations are proposed for the modeling of the elastic and rate (time) dependent inelastic behavior of structural metals in the absence of recovery and aging. The structure of the model is close to the unified theories but contains essential differences. It is shown that the model reproduces almost elastic regions upon initial loading and in the unloading regions of the hysteresis loop. Under loading, unloading and reloading in strain control the model simulated the experimentally observed sharp transition from nearly elastic to inelastic behavior. When a formulation akin to existing unified theories is adopted the almost elastic regions reduce the points and the transition upon reloading is very gradual. For different formulations the behavior under sudden in(de)creases of the strain rate by two orders of magnitude is simulated by numerical experiments and differences are noted. The model represents cyclically neutral behavior and contains three constants and two positive, decreasing functions. The determination of constants and functions from monotonic loading with strain rate changes and relaxation periods is described.

Krempl, E.↗

Transient Thermal Analyses of Passive Systems on SCEPTOR X-57

As efficiency, emissions, and noise become increasingly prominent considerations in aircraft design, turning to an electric propulsion system is a desirable solution. Achieving the intended benefits of distributed electric propulsion (DEP) requires thermally demanding high power systems, presenting a different set of challenges compared to traditional aircraft propulsion. The embedded nature of these heat sources often preclude the use of traditional thermal management systems in order to maximize performance, with less opportunity to exhaust waste heat to the surrounding environment. This paper summarizes the thermal analyses of X-57 vehicle subsystems that don't employ externally air-cooled heat sinks. The high-power battery, wires, high-lift motors, and aircraft outer surface are subjected to heat loads with stringent thermal constraints. The temperature of these components are tracked transiently, since they never reach a steady-state equilibrium. Through analysis and testing, this report demonstrates that properly characterizing the material properties is key to accurately modeling peak temperature of these systems, with less concern for spatial thermal gradients. Experimentally validated results show the thermal profile of these systems can be sufficiently estimated using reduced order approximations.

Thermal Analyses↗

A scoping study of far-SOL main-wall protection limiters for steady-state operation of compact pilot plant tokamaks

We present a novel method for handling steady-state heat fluxes incident on the main wall of pilot plant-scale magnetic fusion devices, based on the utilization of protection limiters in the far scrape-off layer (SOL). This method helps avoid large plasma-wall gaps, without excessively compromising blanket performance. We present an optimization algorithm for determining the appropriate size and scale of these protection limiters given (1) probability distributions of SOL plasma parameters and (2) assumed risk tolerance. As part of this optimization, we have developed an analytic description of parallel heat fluxes across limiter shadows, and an objective cost function (the ‘Far-SOL Marginal Cost’) to quantify the impact that different main-wall thermal management design choices have on reactor capital cost. Applying the model to a midscale fusion pilot plant concept shows that making use of far-SOL protection limiters can reduce capital costs on the order of $500 M, relative to naively increasing the plasma-wall gap. Our analysis demonstrates that the far-SOL power decay length is the highest-leverage plasma assumption for thermal loading of the first wall, and the primary cost driver for main wall thermal management. The relative cost efficiency of protection limiters increases as assumptions on the far-SOL heat flux become more pessimistic. The concepts described in this paper motivate the further development of far-SOL protection limiters as part of larger efforts to design economical core-edge-wall compatible solutions for a fusion pilot plant.

Design under uncertainty↗

Modeling active nematics via the nematic locking principle

Active nematic systems consist of rod-like internally driven subunits that interact with one another to form large-scale coherent flows. They are important examples of far-from-equilibrium fluids, which exhibit a wealth of nonlinear behavior. This includes active turbulence, in which topological defects in the nematic order braid around one another in a chaotic fashion. One of the most studied examples of active nematics consists of a dense two-dimensional layer of microtubules, crosslinked by kinesin molecular motors that inject extensile deformations into the fluid. Though numerous theoretical studies have modeled microtubule-based active nematics, no consensus has emerged on how to fully and quantitatively capture the features of the experimental system. Here, to better understand the theoretical foundations for modeling this system, we propose a fundamental principle we call the nematic locking principle—individual microtubules cannot rotate without all neighboring microtubules also rotating. Physically, this is justified by the high density of the microtubules, their elongated nature, and their corresponding steric interactions. We assert that the nematic locking principle holds throughout the majority of the material, but breaks down in the neighborhood of topological defects and other regions of low density. We derive the most general nematic transport equation consistent with this principle and also derive the most general term that violates it, introducing fracturing into the material. We then examine the standard Beris–Edwards approach, commonly used to model this system, and show that it violates the nematic locking principle throughout the majority of the material due to fracturing. We then propose a modification to the Beris–Edwards model that enforces nematic locking nearly everywhere. This modification shuts off fracturing except in regions where the order parameter (a proxy for density) is reduced. In these regions fracturing is turned on. The resulting simulations in turn show strong nematic locking throughout the bulk of the material, with narrow bands of fracturing, consistent with experimental observation. One additional advantage of enforcing nematic locking is that nontrivial stationary state solutions, common in Beris–Edwards simulations but not seen in experiments, are eliminated.

Mitchell, Kevin A. [Univ. of California, Merced, C↗

Breaking the barrier of human-annotated training data for machine learning-aided plant research using aerial imagery

Machine learning (ML) can accelerate biological research. However, the adoption of such tools to facilitate phenotyping based on sensor data has been limited by (i) the need for a large amount of human-annotated training data for each context in which the tool is used and (ii) phenotypes varying across contexts defined in terms of genetics and environment. This is a major bottleneck because acquiring training data is generally costly and time-consuming. This study demonstrates how a ML approach can address these challenges by minimizing the amount of human supervision needed for tool building. A case study was performed to compare ML approaches that examine images collected by an uncrewed aerial vehicle to determine the presence/absence of panicles (i.e. “heading”) across thousands of field plots containing genetically diverse breeding populations of 2 Miscanthus species. Automated analysis of aerial imagery enabled the identification of heading approximately 9 times faster than in-field visual inspection by humans. Leveraging an Efficiently Supervised Generative Adversarial Network (ESGAN) learning strategy reduced the requirement for human-annotated data by 1 to 2 orders of magnitude compared to traditional, fully supervised learning approaches. The ESGAN model learned the salient features of the data set by using thousands of unlabeled images to inform the discriminative ability of a classifier so that it required minimal human-labeled training data. This method can accelerate the phenotyping of heading date as a measure of flowering time in Miscanthus across diverse contexts (e.g. in multistate trials) and opens avenues to promote the broad adoption of ML tools.

59 BASIC BIOLOGICAL SCIENCES↗

Backpropagated frame transformation theory: A reformulation

In this work, the energy-dependent frame transformation theory of Gao and Greene [H. Gao and C. H. Greene, Phys. Rev. A 42, 6946 (1990)] is extended to yield a quantitatively accurate description of the dissociative recombination process. Evidence is presented to show that direct application of the original theory leads to inaccurate cross sections. A major revision, based on an interaction-free backpropagation of the Born-Oppenheimer solutions, markedly improves the frame transformation theory, reducing its average error by orders of magnitude. The original theory and its extension are tested on the previously explored two-dimensional (2D) model that is tailored to describe the singlet ungerade states of molecular hydrogen. The 2D model can be solved exactly (within the numerical accuracy) without implementing the Born-Oppenheimer approximation. These exact results then serve as a benchmark for the frame transformation theory developed in this paper.

74 ATOMIC AND MOLECULAR PHYSICS↗

Direct link between disorder and magnetoresistance in topological semimetals

The extent to which disorder influences the properties of topological semimetals is relevant to the understanding of topological states and their use in practical applications. Using molecular beam epitaxy, we achieve systematic control of point defect concentrations in the prototypical Dirac semimetal Cd 3 As 2 to gain insight into the role of disorder on electron transport behavior. Using the guiding center diffusion model for linear magnetoresistance, we extract point defect densities as a function of deposition conditions. We find that reducing cadmium defect concentrations by an order of magnitude results in an 2x increase in the magnetoresistance from 450% to 900%. This finding yields important information in the quest to identify the origin of linear magnetoresistance in a wider range of materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

pixelvar79/ESGAN-Flowering-Detection-paper

Machine learning (ML) can accelerate biological research. However, the adoption of such tools to facilitate phenotyping based on sensor data has been limited by (i) the need for a large amount of human-annotated training data for each context in which the tool is used and (ii) phenotypes varying across contexts defined in terms of genetics and environment. This is a major bottleneck because acquiring training data is generally costly and time-consuming. This study demonstrates how a ML approach can address these challenges by minimizing the amount of human supervision needed for tool building. A case study was performed to compare ML approaches that examine images collected by an uncrewed aerial vehicle to determine the presence/absence of panicles (i.e. “heading”) across thousands of field plots containing genetically diverse breeding populations of 2 Miscanthus species. Automated analysis of aerial imagery enabled the identification of heading approximately 9 times faster than in-field visual inspection by humans. Leveraging an Efficiently Supervised Generative Adversarial Network (ESGAN) learning strategy reduced the requirement for human-annotated data by 1 to 2 orders of magnitude compared to traditional, fully supervised learning approaches. The ESGAN model learned the salient features of the data set by using thousands of unlabeled images to inform the discriminative ability of a classifier so that it required minimal human-labeled training data. This method can accelerate the phenotyping of heading date as a measure of flowering time in Miscanthus across diverse contexts (e.g. in multistate trials) and opens avenues to promote the broad adoption of ML tools.

Varela, Sebastian↗