Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coding standards”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

A robust method for handling low density regions in hybrid simulations for collisionless plasmas

A robust method to handle vacuum and near vacuum regions in hybrid simulations for space and astrophysical plasmas is presented. The conventional hybrid simulation model dealing with kinetic ions and a massless charge-neutralizing electron fluid is known to be susceptible to numerical instability due to divergence of the whistler-mode wave dispersion, as well as division-by-density operation in regions of low density. Consequently, a pure vacuum region is not allowed to exist in the simulation domain unless some ad hoc technique is used. To resolve this difficulty, an alternative way to introduce finite electron inertia effect is proposed. Contrary to the conventional method, the proposed one introduces a correction to the electric field rather than the magnetic field. It is shown that the generalized Ohm's law correctly reduces to Laplace's equation in a vacuum which therefore does not involve any numerical problems. In addition, a variable ion-to-electron mass ratio is introduced to reduce the phase velocity of high frequency whistler waves at low density regions so that the stability condition is always satisfied. It is demonstrated that the proposed model is able to handle near vacuum regions generated as a result of nonlinear self-consistent development of the system, as well as pure vacuum regions set up at the initial condition, without losing the advantages of the standard hybrid code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Recent advances in combustion science related to hydrogen safety

Hydrogen is a key pillar in the global Net Zero strategy. Rapid scaling up of hydrogen production, transport, distribution and utilization is expected. This entails that hydrogen, which is traditionally an industrial gas, will come into proximity of populated urban areas and in some situations handled by the untrained public. To realize all their benefits, hydrogen and its technologies must be safely developed and deployed. The specific properties of hydrogen involving wide flammability range, low ignition energy and fast flame speed implies that any accidental release of hydrogen can be easily ignited. Comparing with conventional fuels, combustion systems fueled by hydrogen are also more prone to flame instability and abnormal combustion. This paper aims to provide a comprehensive review about combustion research related to hydrogen safety. It starts with a brief introduction which includes some overview about risk analysis, codes and standards. The core content covers ignition, fire, explosions and deflagration to detonation transition (DDT). Considering that DDT leads to detonation, and that detonation may also be induced directly under special circumstances, the subject of detonation is also included for completeness. The review covers laboratory, medium and large-scale experiments, as well as theoretical analysis and numerical simulation results. While highlights are provided at the end of each section, the paper closes with some concluding remarks highlighting the achievements and key knowledge gaps.

08 HYDROGEN↗
DOE OSTI2024

Ambient energy for buildings: Beyond energy efficiency

The following Key Messages comprise the salient findings of this study: 1. Ambient energy (from sun, air, ground, and sky) can heat and cool buildings; provide hot water, ventilation and daylighting; dry clothes; and cook food. These services account for about three-quarters of building energy consumption and a third of total US demand. Biophilic design (direct and indirect connections with nature) is an intrinsic adjunct to ambient energy systems, and improves wellness and human performance. 2. The current strategy of electrification and energy efficiency for buildings will not meet our climate goals, because the transition to an all-renewable electric grid is too slow. Widespread adoption of ambient energy is needed. Solar-heated buildings also flatten the seasonal demand for electricity compared to all-electric buildings, reducing required production capacity and long-term energy storage. In addition, ambient-conditioned buildings improve resilience by remaining livable during power outages. 3. National policies, incentives, and marketing should be enacted to promote ambient energy use. Federal administrative priorities should reflect the importance of ambient energy for buildings. Use of ambient energy should be encouraged through existing and new building codes and standards. 4. Ambient energy system design tools are needed for architects, engineers, builders, building scientists, realtors, appraisers, and consumers. PVWatts is used over 100 million times per year for photovoltaic system design. A similar, simple, and accessible tool for ambient design is crucial. 5. Training on ambient energy is needed throughout secondary, post-secondary, and continuing education for workforce development. Currently, only about 10% of colleges teach courses on passive heating and cooling systems. 6. Ambient-conditioned buildings should be demonstrated in all US climate zones. Performance should be monitored and reported, with quantitative case studies made widely available. 7. While current technology is sufficient to build high-performance ambient buildings now, research is needed to develop new technologies to harness ambient energy more effectively and more economically. Such advancements will facilitate adoption of ambient energy technologies in a wider range of buildings, including retrofits. Examples include windows with much lower thermal losses, use of the building shell as thermal storage, alternative light-weight thermal storage systems, sky radiation cooling systems, automated controls for solar gains and passive cooling, and ground coupling.

14 SOLAR ENERGY↗
DOE OSTI2021

Urban microclimate and its impact on building performance: A case study of San Francisco

Urban microclimate exerts an increasing influence on urban buildings, energy, and sustainability. This study uses 10-year measured hourly weather data at 27 sites in San Francisco, California, to (1) analyze and visualize the urban microclimate patterns and urban heat island effect; (2) simulate annual energy use and peak electricity demand of typical large office buildings and large hotels to investigate the influence of urban microclimate on building performance; (3) simulate indoor air temperature of a single-family house without air-conditioning during the record three-day heatwave of 2017, to quantify the divergence of climate resilience due to urban microclimate effect. Results show significant microclimate effects in San Francisco with up to 11 °C outdoor air temperature difference between the coastal and downtown areas on September 1, 2017, during the record three-day heatwave. The simulated energy results of the prototype large office and large hotel buildings using the 2017 weather data show over 100% difference in annual heating energy use and 65% difference in annual cooling energy use across different stations; as well as up to 30% difference in peak cooling electricity demand. The impacts on annual site or source energy use are minimal (less than 5%) as cooling and heating in a mild climate are a relatively small portion of overall building energy use in San Francisco. Results also show the microclimate effects influence indoor air temperature of unconditioned homes by up to 5 °C. Newer buildings and homes are much less affected by microclimate effects due to more stringent performance requirements of the building envelope and energy systems. Furthermore, these findings inform that San Francisco microclimate variations should be considered in urban energy planning, building energy codes and standards, as well as heat resilience policymaking.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

A coarse-grain force field based on quantum mechanics (CGq FF) for molecular dynamics simulation of poly(ethylene glycol)- block -poly(ε-caprolactone) (PEG- b -PCL) micelles

In order to provide the means to predict from molecular dynamics (MD) simulations the structures of copolymer-based micelles in solution, we developed coarse grain force field (CGq FF) parameters for poly(ethylene glycol) (PEG) and for poly(ε-caprolactone) (PCL). A key advance here is the use of quantum mechanics to train the parameters describing the non-bonded (NB) interactions between the CG beads. The functional forms are the same as the MARTINI CG FF so standard MD codes can be used. Our CGq FF describes well the experimentally observed properties for the polymer–air and polymer–water interfaces, indicating the accuracy of the NB interactions. Additionally, the structural properties (density, radius of gyration ( R g ), and end-to-end distance ( h )) match both experiment and all atom (AA) simulations. We illustrate the application of this CGq FF by following the formation of a spherical micelle from 250 chains of PEG 23 - b -PCL 9 diblock copolymer, each block with molecular weight of 1000 Daltons (10 500 beads, corresponding to 123 250 atoms), in a water box with 119 139 water beads (426 553 water molecules).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Thoroughly testing and integrating hundreds of Pull Requests per month: ROOT’s new Cost-efficient and Feature Rich GitHub-based CI

ROOT is an open source framework, freely available on GitHub, at the heart of data acquisition, processing and analysis of HE(N)P experiments, and beyond. It is developed collaboratively: contributions are not authored only by ROOT team members, but also by the user community at large: developers and scientists from universities, labs as well as the private sector. More than 1500 GitHub Pull Requests are merged on average per year. It is in this context that code integration acquires a primary role. The review of code contributions isn’t enough: not only they need to be thoroughly reviewed, they also need to be thoroughly tested through a powerful CI infrastructure on several different platforms to comply with the high code quality standards of the project. Since the end of 2023, ROOT moved its continuous integration system from Jenkins to GitHub Actions. In this contribution, we characterise the transition to the GitHub CI, focussing on our strategy, its implementation and the lessons learned, as well as the advantages the new system offers with respect to the previous one. Particular emphasis will be given to the evaluation of the cost-benefit ratio for Jenkins and GitHub Actions for the ROOT project. We also describe how we manage to run in less than one hour thousands of unit, integration, functional and end-to-end tests on different flavours of Windows, four versions of macOS, as well as about ten of the most used Linux distributions, taking advantage of the CERN computing infrastructure.

Piparo, Danilo [CERN]↗
DOE OSTI2014

Dispersion relations for slow and fast resistive wall modes within the Haney-Freidberg model

The dispersion relation for the resistive wall modes (RWMs) is derived by using the trial function for the magnetic perturbation proposed in S. W. Haney and J. P. Freidberg, Phys. Fluids B 1, 1637 (1989). The Haney-Freidberg (HF) approach is additionally based on the expansion in d{sub w}/s≪1, where d{sub w} is the wall thickness and s is the skin depth. Here, the task is solved without this constraint. The derivation procedure is different too, but the final result is expressed in a similar form with the use of the quantities entering the HF relation. The latter is recovered from our more general relation as an asymptote at d{sub w}≪s, which proves the equivalence of the both approaches in this case. In the opposite limit (d{sub w}≫s), we obtain the growth rate γ of the RWMs as a function of γ{sub HF} calculated by the HF prescription. It is shown that γ∝γ{sub HF}{sup 2} and γ≫γ{sub HF} in this range. The proposed relations give γ for slow and fast RWMs in terms of the integrals calculated by the standard stability codes for toroidal systems with and without a perfectly conducting wall. Also, the links between the considered and existing toroidal and cylindrical models are established with estimates explicitly showing the relevant dependencies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Ion-kinetic simulations of D- 3 He gas-filled inertial confinement fusion target implosions with moderate to large Knudsen number

Experiments designed to investigate the transition to non-collisional behavior in D 3 He-gas inertial confinement fusion target implosions display increasingly large discrepancies with respect to simulations by standard hydrodynamics codes as the expected ion mean-free-paths λ c increase with respect to the target radius R (i.e., when the Knudsen number N K = λ c /R grows). To take properly into account large N K 's, multi-ion-species Vlasov-Fokker-Planck computations of the inner gas in the capsules have been performed, for two different values of N K , one moderate and one large. The results, including nuclear yield, reactivity-weighted ion temperatures, nuclear emissivities, and surface brightness, have been compared with the experimental data and with the results of hydrodynamical simulations, some of which include an ad hocmodeling of kinetic effects. The experimental results are quite accurately rendered by the kinetic calculations in the smaller-N K case, much better than by the hydrodynamical calculations. The kinetic effects at play in this case are thus correctly understood. However, in the higher-N K case, the agreement is much worse. Furthermore, the remaining discrepancies are shown to arise from kinetic phenomena (e.g., inter-species diffusion) occurring at the gas-pusher interface, which should be investigated in the future work.

42 ENGINEERING↗
DOE OSTI2018

Speeding up simulations by slowing down particles: Speed-limited particle-in-cell simulation

Based on the particle-in-cell (PIC) plasma simulation method, the speed-limited PIC (SLPIC) method delivers faster kinetic plasma simulation in cases where the particle distributions evolve slowly compared with the maximum stable PIC timestep. SLPIC thus offers more feasible, fully kinetic simulation in regimes that historically have required fluid approaches, such as magnetohydrodynamic (MHD), two-fluid, or Boltzmann electron treatments. In particular, SLPIC allows an explicit time advance with steps much larger than the inverse plasma frequency, avoiding the instability explicit PIC faces with large timesteps. SLPIC avoids this instability by slowing down fast particles (e.g., electrons) in a way that is rigorously underpinned by an approximate Vlasov equation; unlike MHD, two-fluid, and Boltzmann electron approaches, SLPIC does not fundamentally neglect any first-principles plasma physics, although the choices of the grid cell size, timestep, and number of macroparticles per cell naturally limit the physical phenomena that can be accurately represented. SLPIC can be implemented with minor modifications of a standard PIC code and does not require an implicit time advance. It enables large timesteps in first-principles kinetic plasma simulation of appropriately slow phenomena, and it can handle many of the same complications as PIC, such as boundary conditions and collisions. Finally, in an argon plasma sheath test problem, a SLPIC simulation achieved a speed-up of a factor of 160 over the corresponding PIC simulation, without loss of accuracy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Design and Analysis of the US ITER Central Solenoid Module Lifting Fixture

The central solenoid (CS) consists of six large high field superconducting magnets (also known as modules) approximately 4 m in diameter and 2 m tall that weigh approximately 120 tonnes each. These large and complex modules create challenges during assembly of the CS that require the development of custom assembly tooling such as the CS lifting fixture. The CS module lifting fixture is designed to lift and stack the six CS modules in the assembly building on the ITER site. Because of its unique design, fabrication, and assembly features, no lifting attachments could be incorporated within or under the CS modules. This limitation motivated the development of a friction-based lift fixture. The design and evaluation of the CS module lifting fixture considered both worker safety and investment protection, and the assessments were performed to international codes and standards. The CS module lifting fixture consists of two principal subassemblies: spider assembly and ring weldment. These subassemblies allow the frictional force to be augmented by the mechanical advantage of shallow-angle wedges. Large radial preloads created by both screw jack assemblies and the weight of a CS module develop frictional forces capable of performing a lift with a safety factor of at least 2. The design effort resulted in the use of low friction linear bearings on angle surfaces to ensure constant pressure, integrated jacks for pretensioning the fixture prior to lifting, and load pin strain gauges for monitoring the normal force. Testing of various materials and surface treatments led to the selection of laminated aluminum and rubber pads as the friction interface on the CS lifting fixture side and a grit-blasted Nitronic 50 stainless steel band on the CS module side. A redundant lifting method using the module slings between the spider and module is also utilized after the initial friction lift. The CS lifting fixture provides a safe and reliable solution for lifting and stacking the CS module during assembly.

Morris, Eric↗
DOE OSTI2021

Structural Analysis of the US ITER ECH Transmission Line System

The electron cyclotron heating (ECH) and current drive system is one of the main plasma heating systems for ITER. It uses high-power microwave beams with the power deposition location steerable across the plasma cross section. Microwave power is conveyed via transmission lines (TLs) that run from the gyrotrons in the radio frequency building through the assembly hall and tokamak building to the ECH launchers within the tokamak vacuum vessel. The ECH system includes a vast array of interconnected TL waveguides, in-line components, and support structures.Finite element (FE) modeling provides an essential means of simulating the system, applying loads and determining deflections, rotations, forces, moments, and stresses in order to evaluate various structural and microwave transmission performance metrics. A representative FE model of the overall ECH TL system is developed in ANSYS®. This top-level model defines the centerline of the waveguide system. Waveguide segments are represented by line elements (beams and pipes) with equivalent section properties, and support structures are represented by boundary conditions. A systematic approach is used to model each ECH component with lumped masses and structurally equivalent stiffness matrices or ANSYS superelements.The top-level TL FE model is used to evaluate the various loads (thermal, vacuum, seismic, etc.) and operating scenarios. The top-level model directly calculates stresses in the straight aluminum waveguide segments. The model provides the forces and moments acting on the in-line components for detailed submodel assessments. Displacement results from the top-level analysis feed into a separate microwave performance model to help determine operational efficiency. All TL performance and thermal-structural requirements are met, as specified by the applicable codes and standards, and successfully documented in numerous technical reports and demonstrated at the final design review.

Shanmugasundaram, Aravind↗
DOE OSTI2015

ON VERTICALLY GLOBAL, HORIZONTALLY LOCAL MODELS FOR ASTROPHYSICAL DISKS

Disks with a barotropic equilibrium structure, for which the pressure is only a function of the density, rotate on cylinders in the presence of a gravitational potential, so that the angular frequency of such a disk is independent of height. Such disks with barotropic equilibria can be approximately modeled using the shearing box framework, representing a small disk volume with height-independent angular frequency. If the disk is in baroclinic equilibrium, the angular frequency does generally depend on height, and it is thus necessary to go beyond the standard shearing box approach. In this paper, we show that given a global disk model, it is possible to develop approximate models that are local in horizontal planes without an expansion in height with shearing-periodic boundary conditions. We refer to the resulting framework as the vertically global shearing box (VGSB). These models can be non-axisymmetric for globally barotropic equilibria but should be axisymmetric for globally baroclinic equilibria. We provide explicit equations for this VGSB which can be implemented in standard magnetohydrodynamic codes by generalizing the shearing-periodic boundary conditions to allow for a height-dependent angular frequency and shear rate. We also discuss the limitations that result from the radial approximations that are needed in order to impose height-dependent shearing periodic boundary conditions. We illustrate the potential of this framework by studying a vertical shear instability and examining the modes associated with the magnetorotational instability.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Simulation of site and topographic effects on ground motion in Los Alamos, NM mesas

Mesas, which consist of an elevated area of land with a flat top and steep cliffs at the sides, are one of the common geological formations present in the Los Alamos region. Previous research has shown that geological formations such as sedimentary canyons can amplify soil response during earthquakes. There have also been parametric studies to understand the response of an idealized and isolated mountain/canyon under inclined plane waves. In this study, a 2-D linear anelastic soil domain, with topography modelled after the Los Alamos region with many mesas and canyons, is considered to understand site-specific topographic effects in the presence of non-isolated topographical features. Various earthquake sources configurations ranging from point sources to finite faults with varying rupture length, dip angles and spatial slip distributions are considered. To isolate the effect of topography, three other soil domains—a homogeneous half-space, homogeneous half-space with mesas and canyons, and a computational domain with just the soil properties from the Los Alamos region on a flat soil domain is also considered. 2-D site-response analyses of these soil domains under earthquake excitation show that the free-field response of the soil can be amplified or de-amplified depending on the topography of the region and the location of the station relative to the fault. These studies also show that even relatively small mesas with height less than 100 m can significantly amplify the response (by a factor of 2 or more), which differ from the much smaller amplification factors (≤1.2) specified by standard building codes such as Eurocode-8 (2000) for topographical features with a similar slope. De-amplifications are also fairly common, especially close to the canyons. The results from this study agree qualitatively with those from the ambient vibration study conducted by Stolte et al. on a mesa from the Los Alamos region. Furthermore, such site-specific studies provide important insights into the variability of the topographic amplification factors within a region of interest. Such knowledge is important in the design of safety-related critical infrastructure located within that region.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

DSPS: Differentiable stellar population synthesis

ABSTRACT Models of stellar population synthesis (SPS) are the fundamental tool that relates the physical properties of a galaxy to its spectral energy distribution (SED). In this paper, we present DSPS: a python package for SPS. All of the functionality in DSPS is implemented natively in the JAX library for automatic differentiation, and so our predictions for galaxy photometry are fully differentiable, and directly inherit the performance benefits of JAX, including portability onto GPUs. DSPS also implements several novel features, such as i) a flexible empirical model for stellar metallicity that incorporates correlations with stellar age, ii) support for the Diffstar model that provides a physically-motivated connection between the star formation history of a galaxy (SFH) and the mass assembly of its underlying dark matter halo. We detail a set of theoretical techniques for using autodiff to calculate gradients of predictions for galaxy SEDs with respect to SPS parameters that control a range of physical effects, including SFH, stellar metallicity, nebular emission, and dust attenuation. When forward modelling the colours of a synthetic galaxy population, we find that DSPS can provide a factor of 5 speed-up over standard SPS codes on a CPU, and a factor of 300-400 on a modern GPU. When coupled with gradient-based techniques for optimization and inference, DSPS makes it practical to conduct expansive likelihood analyses of simulation-based models of the galaxy–halo connection that fully forward model galaxy spectra and photometry.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Radiation of a charge flying in a partially loaded dielectric section of a waveguide

We consider the electromagnetic field of a point charged particle moving along the axis of a cylindrical waveguide from a homogeneously filled area to a dielectric loading area having an axially symmetrical channel. We are interested in studying the Cherenkov radiation excited in the bilayer area. The solution is performed by expanding the field in each area in a series of orthogonal eigenmodes. The main attention is focused on investigation of the wave field in the bilayer section. We show that, at a given observation point, the "reduced wakefield" is simplified with time (the number of modes decreases). The obtained results are generalized for the case of a bunch with Gaussian longitudinal profile. The typical numerical results for wakefield formation process are presented. These results agree with simulations done by the industry standard electromagnetic code CST Particle Studio.

Grigoreva, Aleksandra A.↗
DOE OSTI2024

Hybrid Modeling of Three-Phase Grid-Supporting Inverters for Dynamic Studies

Grid technologies connected by power electronic converter (PEC) interfaces continually implement grid support functions mandated by grid codes and standards. The transition to converter-based generation demands precise PEC models to assess system dynamics, which have been previously overlooked in conventional power systems. This study proposes a hybrid method for analyzing grid-connected three-phase PEC dynamics with the IEEE standard 1547-2018 Volt-VAr mode that combines physics and data-driven techniques. The physics model reflects the PEC’s internal behavior, whereas the data-driven modeling technique evaluates the grid-supporting capabilities of the smart PEC. The system identification approach is used to generate dynamic PEC models based on changing grid voltage and measured current injected into the grid by the PEC. In the Volt-VAr support mode, a detailed topological model including switches is utilized to compare the goodness-of-fit of the extracted hybrid dynamic model. The results demonstrate that the hybrid PEC model in the Volt-VAr mode accurately matches the dynamics with the topological model.

Subedi, Sunil↗
DOE OSTI2020

Evaluating the Incident Energy of Arcs in Photovoltaic DC Systems: Comparison Between Calculated and Experimental Data

Solar Photovoltaic (PV) systems have permeated the energy generation world at a very high rate, some of the safety codes and standards are still lagging in accurately assessing the hazards and risks associated with PV array arcing energies. Safety professionals and maintenance workers using NFPA 70E have utilized the Doan, Stokes & Oppenlander or Enrique models, meant to determine arc energies in DC power systems using the maximum power method. These methods may lead to an overestimate of energy available in PV systems during a fault. Since PV modules/arrays are non-linear, current limited DC devices, some of these calculation methods may not accurately predict fault energy. This paper will validate current arc energy models for PV systems by comparing experimental and calculated data. Additionally, this data will help modify the current NFPA 70E models related to smaller solar arrays. Understanding where the real safety threshold for DC arc flash in PV systems exists will help maintenance and safety professionals better prepare for a variety of work related activities. This paper will analyze real arc data taken for PV systems <1000VDC and <60amps and compare this to the calculated incident energy models, to include 70E. Thus, using these comparisons, it may be possible to reduce the safety hazard severity and thus relax the PPE requirements for installation and maintenance crews.

14 SOLAR ENERGY↗
DOE OSTI2012

Approach for Selection of Rayleigh Damping Parameters Used for Time History Analysis

Nonlinearities, whether geometric or material, need to be addressed in seismic analysis. One good analysis method that can address these nonlinearities is direct time integration with Rayleigh damping. Modal damping is the damping typically specified in seismic analysis Codes and Standards. Modal damping is constant for all frequencies where Rayleigh damping varies with frequency. An approach is proposed here for selection of Rayleigh damping coefficients to be used in seismic analyses that is consistent with given Modal damping. The approach uses the difference between the modal damping response and the Rayleigh damping response along with effective mass properties of the model being evaluated to match overall system response levels. This paper provides a simple example problem to demonstrate the approach. It also provides results for a finite element model representing an existing piping system. Displacement, acceleration, and stress results are compared from model runs using modal damping and model runs using Rayleigh damping with coefficients selected using the proposed method.

42 ENGINEERING↗