Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Operator approximation”

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,423 records · Page 79

Projective embedding of dynamical systems: Uniform mean field equations

Herein we study embeddings of continuous dynamical systems in larger dimensions via projector operators. We call this technique PEDS, projective embedding of dynamical systems, as the stable fixed point of the original system dynamics are recovered via projection from the higher dimensional space. In this paper we provide a general definition and prove that for a particular type of rank-1 projector operator, the uniform mean field projector, the equations of motion become a mean field approximation of the dynamical system. While in general the embedding depends on a specified variable ordering, the same is not true for the uniform mean field projector. We prove a variety of results on the relationship between the spectrum of the Jacobian for fixed points in the original and in the embedded system. Direct applications of PEDS can be non-convex optimization and machine learning.

97 MATHEMATICS AND COMPUTING↗

The 4th order GISS model of the global atmosphere

The new GISS 4th order model of the global atmosphere is described. It is based on 4th order quadratically conservative differences with the periodic application of a 16th order filter on the sea level pressure and potential temperature equations, a combination which is approximately enstrophy conserving. Several short range forecasts indicate a significant improvement over 2nd order forecasts with the same resolution (approximately 400 km). However the 4th order forecasts are somewhat inferior to 2nd order forecasts with double resolution. This is probably due to the presence of short waves in the range between 1000 km and 2000 km, which are computed more accurately by the 2nd order high resolution model. An operation count of the schemes indicates that with similar code optimization, the 4th order model will require approximately the same amount of computer time as the 2nd order model with the same resolution. It is estimated that the 4th order model with a grid size of 200 km provides enough accuracy to make horizontal truncation errors negligible over a period of a week for all synoptic scales (waves longer than 1000 km).

Kalnay-Rivas, E.↗

Switching regulator emission control circuit for ion sources

An electron emission control circuit of the switching regulator type operating at 100 kHz has been developed which maintains a constant emission current within 0.1% for a cathode power demand variation of approximately 100%. The power output stage has an efficiency of 67%, and the overall efficiency is 45% when driving a thoria-coated iridium cathode having a nominal resistance at operating temperature of 2.5 ohms. Under optimum conditions, the bus power demand is 1.75 W. The circuit is useful in controlling the electron emission current of ion sources in applications which involve a substantial variation of the cathode work function, such as oxygen partial pressure measurements over a large dynamic range.

Clay, F. P., Jr.↗

On the two-dimensional modelling of horizontal CVD reactors and its limitations

In order to clarify under what conditions two-dimensional theoretical models can realistically account for the actual growth rate distribution in horizontal CVD reactors, a systematic numerical study was performed, and the results were compared with published experimental data for the MOCVD of gallium arsenide. Three models of increasing complexity were studied. It is shown that two-dimensional models can produce realistic predictions only for reactors with large width-to-height (aspect) ratios, that are operated at subcritical Rayleigh number,i.e., under dominant forced flow conditions. But even then, depending on the carrier gas, thermal (Soret) diffusion must be included in the model. Furthermore, velocity corrections for finite aspect ratios must be made, and buoyancy effects can be significant in the entrance region. On the other hand, the results are relatively insensitive to axial diffusion, temperature dependence of transport properties (versus averaged properties), and an assumed Poiseuille velocity profile (versus a rigorous solution of the velocity field). Yet, boundary layer approximations are prone to give unrealistic results at the low Reynolds numbers typical for CVD operations.

Ouazzani, Jalil↗

A 3-D Navier-Stokes CFD study of turbojet/ramjet nozzle plume interactions at Mach 3.0 and comparison with data

Advanced airbreathing propulsion systems used in Mach 4-6 mission scenarios, usually consist of a single integrated turboramjet or as in this study, a turbojet housed in an upper bay with a separate ramjet housed in a lower bay. As the engines transition from turbojet to ramjet, there is an operational envelope where both engines operate simultaneously. One nozzle concept under consideration has a common nozzle, where the plumes from the turbojet and ramjet interact with one another as they expand to ambient conditions. In this paper, the two plumes interact at the end of a common 2-D cowl, when they both reach an approximate Mach 3.0 condition and then jointly expand to Mach 3.6 at the common nozzle exit plane. At this condition, the turbojet engine operated at a higher NPR than the ramjet, where the turbojet overpowers the ramjet plume, deflecting it approximately 12 degrees downward and in turn the turbojet plume is deflected 6 degrees upward. In the process, shocks were formed at the deflections and a shear layer formed at the confluence of the two jets. This particular case was experimentally tested and the data used to compare with the PARC3D code with k-kl two equation turbulence model. The 2-D and 3-D centerline CFD solutions are in good agreement, but as the CFD solutions approach the outer sidewall, a slight variance occurs. The outer wall boundary layers are thin and do not present much of an interaction, however, where the confluence interaction shocks interact with the thin boundary layer on the outer wall, strong vortices run down each shock causing substantial disturbances in the boundary layer. These disturbances amplify somewhat as they propagate downstream axially from the confluence point. The nozzle coefficient (CFG) is reduced 1/2 percent as a result of this sidewall interaction, from 0.9850 to 0.9807. This three-dimensional reduction is in better agreement with the experimental value of 0.9790.

Chang, Ing↗

A parametric analysis of transport aircraft system weights and costs

In determining unit and operating costs for advanced aircraft, it has been found that by having first-order weight and performance approximations for the aircraft systems and structural components, a step increase in cost prediction accuracy results. This paper presents first-order approximation equations for these systems and components. These equations were developed from data for most current jet transports, and they have been ordered to use a minimum number of performance parameters such as aircraft style, number of passengers, empty and gross weight, cargo load, and operating range. A NASA Ames Research Center aircraft cost program has been used to compare calculated and actual weights for the same aircraft. Good aircraft cost correlation is shown to exist between calculated first-order and actual aircraft weight data.

Anderson, J. L.↗

Operation of helium sub-atmospheric multistage cryogenic centrifugal compressor trains: Part 1 – Steady state modeling and speed ratio selection

Helium cryogenic systems which can provide cooling below the normal boiling point of helium (approximately 4.2 K) are often required by superconducting radio-frequency niobium resonators utilized in modern high-energy particle accelerators. Achieving temperatures below 4.2 K generally involves operating a cryogenic vessel with liquid helium under sub-atmospheric conditions, thereby lowering the saturation pressure and corresponding saturation temperature. Over the last several decades, multi-stage cryogenic centrifugal compressor trains (CC’s) have been operated efficiently and reliably within large-scale cryogenic systems to continuously evacuate helium vapor generated by a device within the vessel, maintaining sub-atmospheric conditions in the vessel while pressurizing the return vapor to above atmospheric conditions. Traditionally, these CC systems have been operated using empirically derived control philosophies and insight gathered from previous operational experience. Recent efforts at the Facility for Rare Isotope Beams (FRIB) have been aimed at the development of a theoretical basis to characterize the operation of multi-stage cryogenic centrifugal compressor train and utilizing predictive model results to generate control parameters. The objective of this research was identifying operational points which adequately balance cryogenic system efficiency, stability, and overall ease of operation. Furthermore, this manuscript provides an overview of the predictive model development, characterization of the FRIB cryogenic centrifugal compressors and implementation of the predicted performance results during steady-state system operation.

Compressor train control↗

Resilient Operating Constraints for Power Distribution Systems under Setpoint Attacks

Integration and operation of distributed generation (DG) and energy storage (ES) in power distribution systems are enabled by communication networks and embedded sensor and control devices that increase the vulnerability of the systems to cyber-threats, broadening the attack surface and making adversary actions more unpredictable. This paper proposes a methodology that uses ellipsoidal approximations to quantify the potential damage caused by successful attacks that affect, directly or indirectly, the desired operation setpoints and may drive the power distribution operation to unsafe states by violating the limits of voltage or line flows. More specifically, a new methodology is introduced to find the optimal non-symmetric operating constraints that can be imposed to each DG and ES in order to guarantee that the power distribution system is resilient to any malicious setpoints. The proposed method takes as inputs the system topology, DG and ES capabilities, and load limits to solve a convex optimization problem formulated using linear matrix inequalities (LMIs) and the power flow equations. The proposed solution is agnostic to the attacker's action or load profile and it does not require any assumption about the location or means of the attack. The numerical results on a test distribution feeder with several DG and ES illustrate how the proposed resilient operating constraints guarantee the security of the power distribution system under setpoint attacks.

Giraldo, Jairo↗

Fermi’s Golden Rule Rate Expression for Transitions Due to Nonadiabatic Derivative Couplings in the Adiabatic Basis

Starting from a general molecular Hamiltonian expressed in the basis of adiabatic electronic and nuclear position states, where a compact and complete expression for the nonadiabatic derivative coupling (NDC) Hamiltonian term is obtained, we provide a general analysis of the Fermi’s golden rule (FGR) rate expression for nonadiabatic transitions between adiabatic states. We then consider a quasi-adiabatic approximation that uses crude adiabatic states and NDC couplings, both evaluated at the minimum potential energy configuration of the initial adiabatic state, for the definition of the zeroth and first-order terms of the Hamiltonian. Although the application of this approximation is rather limited, it allows deriving a general FGR rate expression without further approximation while accounting for non-Condon contribution to the FGR rate arising from momentum operators of NDC terms and its coupling with vibronic displacements. For a generic and widely used model where all nuclear degrees of freedom and environmental effects are represented as linearly coupled harmonic oscillators, we derive a closed-form FGR rate expression that requires only Fourier transform. The resulting rate expression includes quadratic contributions of NDC terms and their couplings to Franck–Condon modes, which require evaluation of two additional bath spectral densities in addition to the conventional one that appears in a typical FGR rate theory based on the Condon approximation. Model calculations for the case where nuclear vibrations consist of both a sharp high-frequency mode and an Ohmic bath spectral density illustrate new features and implications of the rate expression. We then apply our theoretical expression to the nonradiative decay from the first excited singlet state of azulene, which illustrates the utility and implications of our theoretical results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Changes in the performance characteristics of a GaAs near infrared light emitting diode when exposed to various current and thermal stresses

The changes that occurred in the optical and electrical characteristics of a near infrared, GaAs light emitting diode, when operated under various levels and combinations of current and thermal stresses are discussed. A total of forty parts were operated for two thousand hours under eight different sets of dc current and ambient temperature conditions. Degradation in the radiant optical power of these devices was thirty-four percent when operated at their rated current and an ambient temperature of 298K (25 C). Derating the current and/or the thermal stress reduced the degradation of this parameter in approximately a linear manner. All degraded devices behaved similarly, exhibiting rapid nonlinear degradation followed by a gradual linear degradation and finally a period of stable operation. An attempt was made to correlate initial device condition to degradation during stress testing, but met with little success.

Thomas, E. F., Jr.↗

An asymptotic-preserving 2D-2P relativistic Drift-Kinetic-Equation solver for runaway electron simulations in axisymmetric tokamaks

We propose an asymptotic-preserving (AP), uniformly convergent numerical scheme for the relativistic collisional Drift-Kinetic Equation (rDKE) to simulate runaway electrons in axisymmetric toroidal magnetic field geometries typical of tokamak devices. The approach is derived from an exact Green's function solution with numerical approximations of quantifiable impact, and results in a simple, two-step operator-split algorithm, consisting of a collisional Eulerian step, and a Lagrangian orbit-integration step with analytically prescribed kernels. The AP character of the approach is demonstrated by analysis of the dominant numerical errors, as well as by numerical experiments. We demonstrate the ability of the algorithm to provide accurate answers regardless of plasma collisionality on a circular axisymmetric tokamak geometry.

97 MATHEMATICS AND COMPUTING↗

Disentangling multiple gas kinematic drivers in the Perseus galaxy cluster

Galaxy clusters, the Universe’s largest halo structures, are filled with an X-ray-emitting gas with a temperature between 10 million and 100 million degrees. Their evolution is shaped by energetic processes such as feedback from supermassive black holes (SMBHs) and mergers with other cosmic structures. The imprints of these processes on gas kinematics remain largely unknown, restricting our understanding of energy conversion within clusters4. High-resolution spectral mapping with the X-Ray Imaging and Spectroscopy Mission (XRISM) observatory offers a way forward. Here we present XRISM kinematic measurements of the Perseus cluster, radially covering the extent of its cool core. We find direct evidence for at least two dominant drivers of gas motions operating on distinct physical scales: a small-scale driver in the inner approximately 60 kpc, probably associated with the SMBH feedback; and a large-scale driver in the outer core, powered by mergers. This finding suggests that, during the active phase, SMBH feedback drives gas motions, which, if fully dissipated into heat, could have a substantial role in offsetting radiative cooling losses in the Perseus core. Furthermore, our study underscores the necessity of kinematic mapping observations of extended sources to robustly characterize the velocity fields and their role in the evolution of massive halos. It further offers a kinematic diagnostic for SMBH feedback models.

Astronomy and AstroPhysics↗

Open-Loop Control Co-Design of Semisubmersible Floating Offshore Wind Turbines Using Linear Parameter-Varying Models

Here, this paper discusses a framework to design elements of the plant and control systems for floating offshore wind turbines in an integrated manner using linear parameter-varying models. Multiple linearized models derived from aero-elastic simulation software in different operating regions characterized by the incoming wind speed are combined to construct an approximate low-fidelity model of the system. The combined model is then used to generate open-loop, optimal control trajectories as part of a nested control co-design strategy that explores the system’s power production and stability using the platform pitch tilt as a proxy in the context of crucial plant and control design decisions. The radial distance between the central and outer columns and the diameter of the outer columns of the semisubmersible platform are the plant design variables. The platform stability and power production are studied for different plant design decisions. The effect of plant decisions on subsequent power production and stability response of the floating wind turbine is quantified in terms of the levelized cost of energy. The results show that the inner-loop constraints and the plant design decisions affect the turbine’s power and, subsequently, the cost of the system.

17 WIND ENERGY↗

CEBAF Injector Model for KL Beam Conditions

The Jefferson Lab KL experiment will run at the Continuous Electron Beam Accelerator Facility with a much lower bunch repetition rate (7.80 or 15.59 MHz) than nominally used (249.5 or 499 MHz). While the proposed average current of 2.5 - 5.0 muA is relatively low compared to the maximum CEBAF current of approximately 180 muA, the corresponding bunch charge is atypically high for CEBAF injector operation. In this work, we investigated the evolution and transmission of low-rep-rate, high-bunch-charge (0.32 to 0.64 pC) beams through the CEBAF injector. Using the commercial software General Particle Tracer, we have simulated and analyzed the beam characteristics for both values of bunch charge. We performed these simulations with the existing injector using a 130 kV gun voltage. We have calculated and measured the transmission as a function of the photocathode laser spot size and pulse length. We report on the findings of these simulations and optimum parameters for operating the experiment.

Pokharel, Sunil↗

Performance Assessment for the Environmental Restoration Disposal Facility (FY 2019 Annual Status Report)

DOE O 435.1 and DOE M 435.1-1 require that a determination of continued adequacy of the performance assessment (PA), composite analysis, and disposal authorization statement (DAS) be annually, and it must consider the results of data collection and analysis from research, field studies, and monitoring as well as update any Radioactive Waste Management Basis documents. Beginning in 1996, the Environmental Restoration Disposal Facility (ERDF) has accepted low-level radioactive, hazardous, and mixed wastes that were generated during cleanup activities at the Hanford Site. ERDF is composed of a series of cells (currently 10) or disposal areas and can accommodate future design expansions as needed. During this reporting period (fiscal year 2019, October 1, 2018, to September 30, 2019). approximately 139,721 U.S. tons of waste was disposed at ERDF. From inception through September 30, 2019, ERDF has received approximately 18.5 million U.S. tons of waste, which equates to consumption of approximately 88% of the disposal volume. As a condition of the DAS, disposal operations within ERDF must be in accordance with the waste acceptance criteria (WAC) that provide specific radionuclide disposal limits, waste from restrictions, and descriptions of acceptable waste packages in compliance with DOE M 435-1.1 requirements. The ERDF WAC stipulates that waste destined for disposal at ERDF be controlled based on source, physical form, and containment concentration and activity levels

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Annual Status Report (FY 2020): Performance Assessment for the Environmental Restoration Disposal Facility

DOE O 435.1 and DOE M 435.1-1 require that a determination of continued adequacy of the performance assessment (PA) (CP-60089), composite analysis, and disposal authorization statement (DAS) be made annually, and these guidelines must be used to consider the results of data collection and analysis from research, field studies, and monitoring as well as provide the need to update any radioactive waste management basis documents. Beginning in 1996, the Environmental Restoration Disposal Facility (ERDF) started accepting low-level radioactive, hazardous, and mixed wastes generated during cleanup activities at the Hanford Site. ERDF is composed of a series of cells or disposal areas and can accommodate future design expansions as needed. Currently, there are 10 cells. During this reporting period (fiscal year 2020, which extended from October 1, 2019, through September 30, 2020), approximately 3.39E+04 U.S. tons (3.07E+04 metric tons) of waste was disposed at ERDF. From ERDF inception through September 30, 2020, approximately 18.5 million U.S. tons of waste has been disposed at ERDF, which equates to the consumption of approximately 88% of the disposal volume. As a condition of the DAS, disposal operations within ERDF must be in accordance with the waste acceptance criteria (ERDF-00011) that provide specific radionuclide disposal limits, waste form restrictions, and descriptions of acceptable waste packages in compliance with DOE M 435.1-1 requirements. The ERDF waste acceptance criteria stipulate that waste destined for disposal at ERDF be controlled based on source, physical form, and contaminant concentration and activity levels. There have been no changes to the physical configuration of ERDF or to the waste forms (source, physical form, etc.). No new Unreviewed Disposal Question Screenings or Evaluations have been generated during this reporting period. Therefore, there are no noted impacts to the PA, composite analysis, DAS, or radioactive waste management basis documents resulting from the evaluations and screenings. Sum of fraction analysis shows that the disposed inventory meets both the concentration and inventory threshold requirements. A sum of fractions value is computed for ERDF sensitive radionuclides contributing to the all pathways and air pathway inventory limits. Computed values were 8.85E-02 and 1.78E-01, respectively. The disposed waste inventory remained well under the PA imposed limits. Required monitoring was satisfactorily completed during the fiscal year reporting period (fiscal year 2020). Compliance with performance objectives were met as each of the reported values were well below the established limit. Overall, there are no substantive changes to primary PA assumptions nor changes to the PA analysis conclusion; therefore, compliance with DOE O 435.1 and the DAS is maintained.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Front-End Engineering Design (FEED) Study for a Carbon Capture Plant Retrofit to a Natural Gas-Fired Gas Turbine Combined Cycle Power Plant (2x2x1 Duct-Fired 758-MWe Facility with F Class Turbines)

A comprehensive front-end engineering design (FEED) study has been undertaken by Bechtel National Inc. (Bechtel) for locating a post-combustion capture and compression (PCC) unit at Panda’s Sherman natural gas–combined cycle (NGCC) power plant in Sherman, Texas. This is described in the unredacted FEED Study report (Attachment 1) with all supporting documents, numbering over 150. The Study Report is publicly available. Sizing of the PCC plant is based on treating an amount of flue gas equivalent to that produced when generating 420 MW, which is approximately 68% of the total flue gas emitted by the NGCC power plant operating at guarantee condition with duct burners off. A reduced power plant capacity factor was used for sizing the PCC plant because the gas turbines at the site often operate at reduced load due to the high penetration of renewable power in the ERCOT region. The cost of carbon capture is primarily driven by capital cost (and therefore is highly sensitive to capacity factor). Sizing the capture unit so that when used it is nearly always operating at full capacity is critical to the economic viability of the proposed investment.

03 NATURAL GAS↗

In-Situ Pipeline Coatings for Methane Emissions Mitigation and Quantification from Natural Gas Pipelines

Addressing the current health of the nation’s existing 3 million miles of pipeline infrastructure is key to preventing further climate change. In 2020, natural gas production exceeded 34 trillion cubic feet (Tcf). Roughly 75% of natural gas consists of methane (CH 4 ), which is up to 25 times more powerful than carbon dioxide (CO 2 ) at trapping heat within the atmosphere over a 100-year period, and studies from the Environmental Defense Fund (EDF) estimate approximately 2% of all the natural gas produced will be lost during normal operations due to unaddressed leaks. This does not even consider the risks of major disaster due to pipeline failure, or the losses and extra fuel costs incurred due to corrosion and scale deposits in under-maintained pipelines. The objective of the proposed research is to demonstrate the protection capabilities and economic benefits of Oceanit’s internal pipe surface treatment, known as DragX™. DragX™ is a chemically resistant, water-and-oil repellent nanocomposite system that can be readily applied in-situ on natural gas transmission and distribution pipelines with a minimum of surface preparation. This makes it an ideal candidate for in-place retrofitting and refurbishment of existing pipelines without the need for expensive extraction and replacement. DragX™ is also able to significantly reduce the surface roughness, and subsequently, the frictional drag forces within a pipeline, improving throughput, decreasing energy costs of pressurization and pumping, and allowing for longer pipeline operation without interruption, reducing the methane emitted during pipe isolation and venting. As part of this project, Oceanit has utilized the Department of Energy’s support to fully develop, de-risk and prove the DragX™ core technology is both economically viable and commercially desirable to pipeline operators and energy companies alike. DragX™ material properties were optimized in this effort both for ease of applicability, to provide value in certain key parameters, and was demonstrated on pilot applications exceeding 2 miles in length. Beyond the already field demonstrated applications, this innovative nanocomposite surface treatment has the potential to be the backbone for CO 2 and Hydrogen transporting pipeline infrastructure. The learnings from this project could accelerate the deployment of surface treatment technologies related to the energy transition infrastructure, thus benefitting the clean energy initiatives in the United States and all around the world.

03 NATURAL GAS↗