Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “configuration mixing”

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,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Precipitation Processes developed during ARM (1997), TOGA COARE(1992), GATE(1 974), SCSMEX(1998) and KWAJEX(1999): Consistent 2D and 3D Cloud Resolving Model Simulations

Real clouds and clouds systems are inherently three dimensional (3D). Because of the limitations in computer resources, however, most cloud-resolving models (CRMs) today are still two-dimensional (2D). A few 3D CRMs have been used to study the response of clouds to large-scale forcing. In these 3D simulations, the model domain was small, and the integration time was 6 hours. Only recently have 3D experiments been performed for multi-day periods for tropical cloud system with large horizontal domains at the National Center for Atmospheric Research. The results indicate that surface precipitation and latent heating profiles are very similar between the 2D and 3D simulations of these same cases. The reason for the strong similarity between the 2D and 3D CRM simulations is that the observed large-scale advective tendencies of potential temperature, water vapor mixing ratio, and horizontal momentum were used as the main forcing in both the 2D and 3D models. Interestingly, the 2D and 3D versions of the CRM used in CSU and U.K. Met Office showed significant differences in the rainfall and cloud statistics for three ARM cases. The major objectives of this project are to calculate and axamine: (1)the surface energy and water budgets, (2) the precipitation processes in the convective and stratiform regions, (3) the cloud upward and downward mass fluxes in the convective and stratiform regions; (4) cloud characteristics such as size, updraft intensity and lifetime, and (5) the entrainment and detrainment rates associated with clouds and cloud systems that developed in TOGA COARE, GATE, SCSMEX, ARM and KWAJEX. Of special note is that the analyzed (model generated) data sets are all produced by the same current version of the GCE model, i.e. consistent model physics and configurations. Trajectory analyse and inert tracer calculation will be conducted to identify the differences and similarities in the organization of convection between simulated 2D and 3D cloud systems.

Tao, W.-K.↗
DOE OSTI2022

First snowflake divertor experiments in MAST-U tokamak

First snowflake (SF) divertor experiments in the MAST-U tokamak demonstrated steady-state snowflake-plus divertor configurations in 450 kA ohmic L-mode plasmas. The SF divertor configuration features a second poloidal field (PF) null in the divertor region close by or overlapping with the main X-point. The resulting low PF region and two additional divertor legs (strike points) may lead to additional power and particle flux sharing via a hypothesized convective cell, and increased plasma-wetted area and radiation. The free-boundary Grad–Shafranov equilibrium code FIESTA was used to design SF configurations with several inter-null distances and orientations. In the experiment, the SF configurations with inter-null distances 0.13–0.20 m and lasting 0.2–0.3 s were obtained. Parallel connection lengths between the midplane and the outer strike point in the SF configurations (evaluated at field lines 1–2 mm from the separatrix in the midplane) were 25–30 m, higher than in the standard divertor (20–25 m) or the Super-X divertor (25 m). Diagnostic measurements highlighted salient SF features. The infra-red video bolometer diagnostic showed that the radiated power peaking in the PF null region was not as pronounces as in the standard divertor. Divertor ion fluxes measured by target Langmuir probes showed increased ion flux in the plate region where a secondary SF strike point landed, concomitantly with the SF formation. These measurements may suggest that some particle and heat redistribution was taking place in the convective SF zone. The first SF experiments provide a basis for future SF studies in MAST-U tokamak with higher input power, improved plasma control and diagnostic measurements, to be compared with the modeling predictions of plasma convective SF mixing and lower density strike point detachment threshold.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1999

Three-Dimensional Gear Crack Propagation Studied

Gears used in current helicopters and turboprops are designed for light weight, high margins of safety, and high reliability. However, unexpected gear failures may occur even with adequate tooth design. To design an extremely safe system, the designer must ask and address the question, "What happens when a failure occurs?" With gear-tooth bending fatigue, tooth or rim fractures may occur. A crack that propagates through a rim will be catastrophic, leading to disengagement of the rotor or propeller, loss of an aircraft, and possible fatalities. This failure mode should be avoided. A crack that propagates through a tooth may or may not be catastrophic, depending on the design and operating conditions. Also, early warning of this failure mode may be possible because of advances in modern diagnostic systems. One concept proposed to address bending fatigue fracture from a safety aspect is a splittooth gear design. The prime objective of this design would be to control crack propagation in a desired direction such that at least half of the tooth would remain operational should a bending failure occur. A study at the NASA Lewis Research Center analytically validated the crack-propagation failsafe characteristics of a split-tooth gear. It used a specially developed three-dimensional crack analysis program that was based on boundary element modeling and principles of linear elastic fracture mechanics. Crack shapes as well as the crack-propagation life were predicted on the basis of the calculated stress intensity factors, mixed-mode crack-propagation trajectory theories, and fatigue crack-growth theories. The preceding figures show the effect of the location of initial cracks on crack propagation. Initial cracks in the fillet of the teeth produced stress intensity factors of greater magnitude (and thus, greater crack growth rates) than those in the root or groove areas of the teeth. Crack growth was simulated in a case study to evaluate crack-propagation paths. Tooth fracture was predicted from the crackgrowth simulation for an initial crack in the tooth fillet region. This was the desired failure mode for an ultrasafe design. Lastly, tooth loads on the uncracked mesh of the split-tooth design were up to five times greater than those on the cracked mesh if equal deflections of the cracked and uncracked teeth were considered. This effect needs to be considered in the design of a split-tooth configuration. This work was done in-house at Lewis in support of the National Rotorcraft Technology Center project, Ultra-Safe Gear Design, with the Boeing Defense and Space Group. The crack-propagation package was developed by the Cornell Fracture Group at Cornell University. The reported results, which are the initial findings of Lewis gear-crackpropagation research for the Rotorcraft Base Program, will be further investigated to develop generalized gear design guidelines.

Lewicki, David G.↗
NASA NTRS1952

Investigation of turbines for driving supersonic compressors II : performance of first configuration with 2.2 percent reduction in nozzle flow area / Warner L. Stewart, Harold J. Schum, Robert Y. Wong

The experimental performance of a modified turbine for driving a supersonic compressor is presented and compared with the performance of the original configuration to illustrate the effect of small changes in the ratio of nozzle-throat area to rotor-throat area. Performance is based on the performance of turbines designed to operate with both blade rows close to choking. On the basis of the results of this investigation, the ratio of areas is concluded to become especially critical in the design of turbines such as those designed to drive high-speed, high-specific weight-flow compressors where the turbine nozzles and rotor are both very close to choking.

TURBINES--MIXED FLOW↗
DOE OSTI2022

Evaluating Class 6 Delivery Truck Fuel Economy and Emissions Using Vehicle System Simulations for Conventional and Hybrid Powertrains and Co-Optima Fuel Blends

The US Department of Energy’s Co-Optimization of Engine and Fuels Initiative (Co-Optima) investigated how unique properties of bio-blendstocks considered within Co-Optima help address emissions challenges with mixing controlled compression ignition (i.e., conventional diesel combustion) and enable advanced compression ignition modes suitable for implementation in a diesel engine. Additionally, the potential synergies of these Co-Optima technologies in hybrid vehicle applications in the medium- and heavy-duty sector was also investigated. In this work, vehicles system were simulated using the Autonomie software tool for quantifying the benefits of Co-Optima engine technologies for medium-duty trucks. A Class 6 delivery truck with a 6.7 L diesel engine was used for simulations over representative real-world and certification drive cycles with four different powertrains to investigate fuel economy, criteria emissions, and performance. Comparisons were made between ultralow-sulfur diesel and a blend of 25% hexyl hexanoate with diesel. Model validation data were informed by 2019 model year Cummins ISB 6.7 L diesel engine maps and transient validation data in a pre-production hybrid configuration and a direct dyno coupled configuration with diesel fuel and a blend of 25% hexyl hexanoate with diesel.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS1991

Airborne lidar measurements of ozone during the 1989 airborne Arctic stratospheric expedition

The NASA/NOAA Airborne Arctic Stratospheric Expedition (AASE) was conducted during the winter to study the conditions leading to possible ozone (O3) destruction in the wintertime Arctic stratosphere. As part of this experiment, the NASA-Langley airborne differential absorption lidar (DIAL) system was configured for operation on the NASA-Ames DS-8 aircraft to make measurements of O3 profiles from about 1 km above the aircraft to altitudes of 22 to 26 km. The airborne DIAL system remotely sensed O3 above the DC-8 by transmitting two laser beams at 10 Hz using wavelengths of 301.5 and 311 nm. Large scale distributions of O3 were obtained on 15 long range flights into the polar vortex during the AASE. Selected data samples are presented of O3 observed during these flights, general trends observed in O3 distributions, and correlations between these measurements and meteorological and chemical parameters. The O3 distribution observed on the first flight of the DC-8 into the polar vortex on Jan. 6 reflected the result of diabatic cooling of the air inside the vortex during the winter compared to the warmer air outside the vortex. On a potential temperature surface, the O3 mixing ratio generally increases when going from outside to inside the vortex.

Browell, Edward V.↗
NASA NTRS2012

Simulations of Magnetic Reconnection - Kinetic Mechanisms Underlying the Fluid Description of Ions

Because of its ability to transfer the energy stored in magnetic field together with the breaking of the flux freezing constraint, magnetic reconnection is considered as one of the most important phenomena in plasma physics. When it happens in a collision less environment such as the terrestrial magnetosphere, it should a priori be modelled with in the framework of kinetic physics. The evidence of kinetic features has incidentally for a long time, been shown by researchers with the help of both numerical simulations and satellite observations. However, most of our understanding of the process comes from the more intuitive fluid interpretation with simple closure hypothesis which do not include kinetic effects. To what extent are these two separate descriptions of the same phenomenon related? What is the role of kinetic effects in the averaged/fluid dynamics of reconnection? This thesis addresses these questions for the proton population in the particular case of anti parallel merging with the help of 2D Hybrid simulations. We show that one can not assume, as is usually done, that the acceleration of the proton flow is only due to the Laplace force. Our results show, for symmetric and asymmetric connection, the importance of the pressure force, opposed to the electric one on the separatrices, in the decoupling region. In the symmetric case, we emphasize the kinetic origin of this force by analyzing the proton distribution functions and explain their structure by studying the underlying particle dynamics. Protons, as individual particles, are shown to bounce in the electric potential well created by the Hall effect. The spatial divergence of this well results in a mixing in phase space responsible for the observed structure of the pressure tensor. A detailed energy budget analysis confirms the role of the pressure force for the acceleration; but, contrary to what is sometimes assumed, it also reveals that the major part of the incoming Poynting flux is transferred to the thermal energy flux rather than to the convective kinetic energy flux, although the latter is generally supposed dominant. In the symmetric case, we propose the pressure tensor to be an additional proxy of the ion decoupling region in satellite data and verify this suggestion by studying a reconnection event encountered by the Cluster spacecrafts. Finally, the last part of this thesis is devoted to the study of the kinetic structure of asymmetric tangential current sheets where connection can develop. This theoretical part consists in finding a steady state solution to the Vlasov-Maxwell system for the protons in such a configuration. We present the theory and its first confrontation to numerical tests.

Aunai, icolas↗
NASA NTRS2009

SHJAR Jet Noise Data and Power Spectral Laws

High quality jet noise spectral data measured at the Aeroacoustic Propulsion Laboratory at the NASA Glenn Research Center is used to examine a number of jet noise scaling laws. Configurations considered in the present study consist of convergent and convergent-divergent axisymmetric nozzles. The measured spectral data are shown in narrow band and cover 8193 equally spaced points in a typical Strouhal number range of 0.0 to 10.0. The measured data are reported as lossless (i.e., atmospheric attenuation is added to measurements), and at 24 equally spaced angles (50deg to 165deg) on a 100-diameter (200-in.) arc. Following the work of Viswanathan, velocity power factors are evaluated using a least squares fit on spectral power density as a function of jet temperature and observer angle. The goodness of the fit and the confidence margins for the two regression parameters are studied at each angle, and alternative relationships are proposed to improve the spectral collapse when certain conditions are met. As an immediate application of the velocity power laws, spectral density in shockcontaining jets are decomposed into components attributed to jet mixing noise and shock noise. From this analysis, jet noise prediction tools can be developed with different spectral components derived from different physics.

Khavaran, Abbas↗
NASA NTRS2001

Simple, Scalable, Script-Based Science Processor (S4P)

The development and deployment of data processing systems to process Earth Observing System (EOS) data has proven to be costly and prone to technical and schedule risk. Integration of science algorithms into a robust operational system has been difficult. The core processing system, based on commercial tools, has demonstrated limitations at the rates needed to produce the several terabytes per day for EOS, primarily due to job management overhead. This has motivated an evolution in the EOS Data Information System toward a more distributed one incorporating Science Investigator-led Processing Systems (SIPS). As part of this evolution, the Goddard Earth Sciences Distributed Active Archive Center (GES DAAC) has developed a simplified processing system to accommodate the increased load expected with the advent of reprocessing and launch of a second satellite. This system, the Simple, Scalable, Script-based Science Processor (S42) may also serve as a resource for future SIPS. The current EOSDIS Core System was designed to be general, resulting in a large, complex mix of commercial and custom software. In contrast, many simpler systems, such as the EROS Data Center AVHRR IKM system, rely on a simple directory structure to drive processing, with directories representing different stages of production. The system passes input data to a directory, and the output data is placed in a "downstream" directory. The GES DAAC's Simple Scalable Script-based Science Processing System is based on the latter concept, but with modifications to allow varied science algorithms and improve portability. It uses a factory assembly-line paradigm: when work orders arrive at a station, an executable is run, and output work orders are sent to downstream stations. The stations are implemented as UNIX directories, while work orders are simple ASCII files. The core S4P infrastructure consists of a Perl program called stationmaster, which detects newly arrived work orders and forks a job to run the appropriate executable (registered in a configuration file for that station). Although S4P is written in Perl, the executables associated with a station can be any program that can be run from the command line, i.e., non-interactively. An S4P instance is typically monitored using a simple Graphical User Interface. However, the reliance of S4P on UNIX files and directories also allows visibility into the state of stations and jobs using standard operating system commands, permitting remote monitor/control over low-bandwidth connections. S4P is being used as the foundation for several small- to medium-size systems for data mining, on-demand subsetting, processing of direct broadcast Moderate Resolution Imaging Spectroradiometer (MODIS) data, and Quick-Response MODIS processing. It has also been used to implement a large-scale system to process MODIS Level 1 and Level 2 Standard Products, which will ultimately process close to 2 TB/day.

Lynnes, Christopher↗
NASA NTRS1991

Advances In Li-TiS2 Cell Techno1ogy

JPL is involved in a NASA sponsored program to develop ambient temperature secondary Lithium - Titanium Disulfide Cells for future space missions. After several years of research on various lithium systems, the Li-TiS2 system was selected for development in view of its practically realizable high specific energy. In the last two years, the efforts were focused on improving the cycle life of the system and optimizing the cell design. A number of approaches, such as the use of mixed solvent electrolytes, the use of alternate anode materials, the operation of cells at low temperature, and the cycling of cells under optimized voltage limits, are examined to improve cycle life performance of this system. Cycling studies in small 150 mAh cell have identified 1.5 M LiAsF6/2-MeTHF, 1.5 M LiAsF6/EC+2-eTHF, and 1.5 M LiAsF6/THF+2-MeTHF+2-MeF as promising electrolytes for Li-TiS2 cells. Li-Al and Li-C were selected for further assessment as candidate anode materials after a detailed theoretical and experimental evaluation. Cycling of the cells at low temperature did not result in improving the cycle life of the cells. Charge and discharge voltage limits were found to have significant influence on the cycle life of the cells activated with 1.5 M LiAsF6/THF+2-MeTHF+2-MeF electrolyte. The influence of design variables, such as ratio of electrode capacity, quantity of electrolyte, pack tightness, cell configuration, etc., on the cycle life performance are being examined as a part of the design optimization study. Spiral-wound 1 Ah cells fabricated for the design studies have completed more than 500 cycles at 50% DOD. This paper summarizes the advances made in the Li-TiS2 technology at JPL since 1989.

Surampudi, S.↗
DOE OSTI2022

Real-time dispatch optimization for concentrating solar power with thermal energy storage

Concentrating solar power (CSP) plants present a promising path towards utility-scale renewable energy. The power tower, or central receiver, configuration can achieve higher operating temperatures than other forms of CSP, and, like all forms of CSP, naturally pairs with comparatively inexpensive thermal energy storage, which allows CSP plants to dispatch electricity according to market price incentives and outside the hours of solar resource availability. Currently, CSP plants commonly include a steam Rankine power cycle and several heat exchange components to generate high-pressure steam using stored thermal energy. The efficiency of the steam Rankine cycle depends on the temperature of the plant's operating fluid, and so is a main concern of plant operators. However, the variable nature of the solar resource and the conservatism with which the receiver is operated prevent perfect control over the receiver outlet temperature. Therefore, during periods of solar variability, collection occurs at lower-than-design temperature. To support operator decisions in a real-time setting, we develop a revenue-maximizing non-convex mixed-integer, quadradically-constrained program which determines a dispatch schedule with sub-hourly time fidelity and considers temperature-dependent power cycle efficiency. The exact nonlinear formulation proves intractable for real-time decision support. Here we present exact and inexact techniques to improve problem tractability that include a hybrid nonlinear and linear formulation. Our approach admits solutions within approximately 3% of optimality, on average, within a five-minute time limit, demonstrating its usability for decision support in a real-time setting.

14 SOLAR ENERGY↗
DOE OSTI2022

Local cation order and ferrimagnetism in compositionally complex spinel ferrites

We present an exploration of a family of compositionally complex cubic spinel ferrites featuring combinations of Mg, Fe, Co, Ni, Cu, Mn, and Zn cations, systematically investigating the average and local atomic structures, chemical short-range order, magnetic spin configurations, and magnetic properties. All compositions result in ferrimagnetic average structures with extremely similar local bonding environments; however, the samples display varying degrees of cation inversion and, therefore, differing apparent bulk magnetization. Additionally, first-order reversal curve analysis of the magnetic reversal behavior indicates varying degrees of magnetic ordering and interactions, including potentially local frustration. Finally, reverse Monte Carlo modeling of the spin orientation demonstrates a relationship between the degree of cation inversion and the spin collinearity. Collectively, these observations correlate with differences in synthesis procedures. This work provides a framework for understanding magnetic behavior reported for “high-entropy spinels,” revealing many are likely compositionally complex oxides with differing degrees of chemical short-range order—not meeting the community established criteria for high or medium entropy compounds. Moreover, this work highlights the importance of reporting complete sample processing histories and investigating local to long-range atomic arrangements when evaluating potential entropic mixing effects and assumed property correlations in high entropy materials.

36 MATERIALS SCIENCE↗
NASA NTRS1996

Experimental Investigation of the Performance of a Mach-2.7 Two-dimensional Bifurcated Duct Inlet with 30 Percent Internal Contraction

An experimental study was conducted to determine the performance of a two-dimensional, mixed-compression bifurcated duct inlet system designed for a free-stream Mach number of 2.7. Thirty percent of the supersonic area contraction occurred internally. A movable ramp was used to vary the contraction ratio for off-design operation. Boundary layer bleed regions were located on the cowl, centerbody, and sidewall surfaces. There were also provisions for vortex generators on the cowl and centerbody of the subsonic diffuser. Data were obtained over the Mach number range of 2.0 to 2.8 and at angles of yaw from 0 deg. to the maximum value prior to inlet un-start. The test at Mach 2.8 was to obtain data for an over- speed condition. The Reynolds number varied from 2.5 to 2.3 million/ft for Mach numbers above 2.5. At Mach numbers of 2.5 and lower, the Reynolds number was set at 2.5 million/ft. Bleed patterns, vortex generator patterns, and ramp position were varied, and three inlet configurations were selected for more extensive study. Two of these configurations had self-starting capability. The self-starting configuration that was developed produced 89 percent total pressure recovery at the compressor face station with 6.8 percent total bleed. The compressor face distortion was about 16 percent. Vortex generators were extremely effective in re-distributing flow but were not as effective in reducing distortion. Excellent flow symmetry was achieved between the separated halves of the inlet, and twin-duct instability was not observed. The ramp tip shock was steeper than expected. This caused the cowl lip shock to be reflected from the ramp instead of being cancelled at the shoulder. However, peak recovery at the throat was still obtained with the ramp near the design position.

Wasserbauer, J. F.↗
NASA NTRS1999

Testing of a Spray-Bar Zero Gravity Cryogenic Vent System for Upper Stages

The capability to vent in zero gravity without resettling is a fundamental technology need that involves practically all uses of subcritical cryogenics in space. Venting without resettling would extend cryogenic orbital transfer vehicle capabilities. However, the lack of definition regarding liquid/ullage orientation coupled with the somewhat random nature of the thermal stratification and resulting pressure rise rates, lead to significant technical challenges. Typically a zero gravity vent concept, termed a thermodynamic vent system (TVS), consists of a tank mixer to destratify the propellant, combined with a Joule- Thomson (J-T) valve to extract then-nal energy from the propellant. In a cooperative effort, Marshall Space Flight Center's (MSFC's) Multipurpose Hydrogen Test Bed (N4HTB) was used to test a unique "spray bar" TVS system developed by Boeing. A schematic of this system is included in Figure 1. The system consists of a recirculation pump, a parallel flow concentric tube, heat exchanger, and a spray bar positioned close to the longitudinal axis of the tank. In the mixing mode, the recirculation pump withdraws liquid from the tank and sprays it radially into the tank liquid, ullage, and exposed tank surfaces. When energy extraction is required, a small portion of the recirculated liquid is passed sequentially through the J-T expansion valve, the spray bar heat exchanger element, and is vented overboard. The vented vapor cools the circulated bulk fluid, thereby removing thermal energy and reducing tank pressure. Figure 2 is a plot of ullage pressure (P4) and liquid vapor pressure (PSAI) versus time. The pump operates alone, cycling on and off, to destratify the tank liquid and ullage until the liquid vapor pressure reaches the lower set point. At that point, the J-T valve begins to cycle on and off with the pump. Thus, for short duration missions, only the mixer may operate, thus minimizing or even eliminating boil-off losses. The primary advantage of the spray bar configuration is that pressure reduction is achieved independent of liquid and vapor location, thereby enhancing the applicability of normal gravity test data to zero gravity conditions. The in-tank components are minimized with the proposed TVS design. Because the recirculation pump is external to the tank, no electrical power penetration of the tank is required for pump or valve operation. This is especially desirable for L02 tanks since the presence of an electrical ignition source in oxygen represents a critical failure mode. Also, since the critical components (pump, motor, valve, orifice) are external to the tank, system checkout and ground servicing/replacement are easier. For zero-g operation, component replacement external to the tank may be a significant benefit. In addition to satisfying the zero g TVS design objectives, the TVS concept tested offers additional benefits to the integrated subcritical cryogenic storage and launch system.

Lak, Tibor↗
NASA NTRS1990

A decentralized linear quadratic control design method for flexible structures

A decentralized suboptimal linear quadratic control design procedure which combines substructural synthesis, model reduction, decentralized control design, subcontroller synthesis, and controller reduction is proposed for the design of reduced-order controllers for flexible structures. The procedure starts with a definition of the continuum structure to be controlled. An evaluation model of finite dimension is obtained by the finite element method. Then, the finite element model is decomposed into several substructures by using a natural decomposition called substructuring decomposition. Each substructure, at this point, still has too large a dimension and must be reduced to a size that is Riccati-solvable. Model reduction of each substructure can be performed by using any existing model reduction method, e.g., modal truncation, balanced reduction, Krylov model reduction, or mixed-mode method. Then, based on the reduced substructure model, a subcontroller is designed by an LQ optimal control method for each substructure independently. After all subcontrollers are designed, a controller synthesis method called substructural controller synthesis is employed to synthesize all subcontrollers into a global controller. The assembling scheme used is the same as that employed for the structure matrices. Finally, a controller reduction scheme, called the equivalent impulse response energy controller (EIREC) reduction algorithm, is used to reduce the global controller to a reasonable size for implementation. The EIREC reduced controller preserves the impulse response energy of the full-order controller and has the property of matching low-frequency moments and low-frequency power moments. An advantage of the substructural controller synthesis method is that it relieves the computational burden associated with dimensionality. Besides that, the SCS design scheme is also a highly adaptable controller synthesis method for structures with varying configuration, or varying mass and stiffness properties.

Su, Tzu-Jeng↗
DOE OSTI2023

Design optimization of an ethanol heavy-duty engine using design of experiments and bayesian optimization

Diesel-fueled engines still hold a large market share in the medium and heavy-duty transportation sector. However, the increase in fossil fuel prices and the strict emission regulations are leading engine manufacturers to seek cleaner alternatives without a compromise in performance. Alcohol-based fuels, such as ethanol, offer a promising alternative to diesel fuel in meeting regulatory demands. Ethanol provides cleaner combustion and lower levels of soot due to its chemical properties, in particular its lower level of carbon content. In addition, the stoichiometric operating conditions of alcohol fueled engines enable the mitigation of NOx emissions in aftertreatment stage. With the promise of retrofitting diesel engines to run on ethanol to reduce emissions, the thermal efficiency of these engines remains the primary optimization target. In order to find the optimal ethanol-fueled engine design that maximizes the thermal efficiency, a large design space needs to be investigated using engineering tools. In this study, previous research by the authors on optimizing the design of a single-cylinder ethanol-fueled engine was extended to explore the design space for a heavy-duty multi-cylinder engine configuration. A heavy-duty engine setup with multiple operating conditions at different engine speeds and loads were considered. A design optimization analysis was performed to identify the potential designs that maximize the indicated thermal efficiency in an ethanol-fueled compression ignition engine. First, a computational fluid dynamics (CFD) model of the engine was validated using experimental data for four drive cycle points. Using a design of experiments (DoE) approach and a parameterized piston bowl geometry, the model was then exercised to explore the relationship among geometric features of the piston bowl and spray targeting angle and indicated thermal efficiency across all tested operating conditions. After evaluating 165~candidate designs, a piston bowl geometry was identified that yielded an increase between 1.3 to 2.2 percentage points in indicated thermal efficiency for all tested conditions, while satisfying the operational design constraints for peak pressure and maximum pressure rise rate. The increased performance was attributed to enhanced mixing that led to the formation of a more homogeneous distribution of in-cylinder temperature and equivalence ratio, higher combustion temperatures, and shorter combustion duration. Finally, a Bayesian optimization (BOpt) analysis was employed to find the optimal piston bowl geometry with a fixed spray injector angle for one of the operating conditions. Using BOpt, a piston candidate was identified that resulted in a 1.9~percentage point increase in thermal efficiency from the baseline design, yet only required 65\% of the design samples investigated using the DoE approach.

Tekgul, Bulut↗
DOE OSTI2024

Harnessing complexity: Nonlinear optical phenomena in L-shapes, nanocrescents, and split-ring resonators

Here, we conduct systematic studies of the optical characteristics of plasmonic nanoparticles that exhibit C 2v symmetry. In particular, we analyze three distinct geometric configurations: an L-type shape, a crescent, and a split-ring resonator shaped like the Greek letter π. Optical properties are examined using the finite-difference time-domain method. It is demonstrated that all three shapes exhibit two prominent plasmon modes associated with the two axes of symmetry. This is in addition to a wide range of resonances observed at high frequencies corresponding to quadrupole modes and peaks due to sharp corners. Next, to facilitate nonlinear analysis, we employ a semiclassical hydrodynamic model, where the electron pressure term is explicitly accounted for. This model goes beyond the standard Drude description and enables capturing nonlocal and nonlinear effects. Employing this model enables us to rigorously examine the second-order angular resolved nonlinear optical response of these nanoparticles in each of the three configurations. Two pumping regimes are considered, namely, continuous wave (CW) and pulsed excitations. For CW pumping, we explore the properties of the second harmonic generation (SHG). Polarization and angle-resolved SHG spectra are obtained, revealing strong dependence on the nanoparticle geometry and incident wave polarization. The C 2v symmetry is shown to play a key role in determining the polarization states and selection rules of the SHG signal. For pulsed excitations, we discuss the phenomenon of broadband terahertz (THz) generation induced by the difference-frequency generation . It is shown that the THz emission spectra exhibit unique features attributed to the plasmonic resonances and symmetry of the nanoparticles. The polarization of the generated THz waves is also examined, revealing interesting patterns tied to the nanoparticle geometry. To gain deeper insight, we propose an analytical theory that agrees very well with the numerical experiments. The theory shows that the physical origin of the THz radiation is the mixing of various frequency components of the fundamental pulse by the second-order nonlinear susceptibility. An expression for the far-field THz intensity is derived in terms of the incident pulse parameters and the nonlinear response tensor of the nanoparticle. The results presented in this work offer new insights into the linear and nonlinear optical properties of nanoparticles with C 2v symmetry. The demonstrated strong SHG response and efficient broadband THz generation hold great promise for applications in nonlinear spectroscopy, nanophotonics, and optoelectronics. The proposed theoretical framework also provides a valuable tool for understanding and predicting the nonlinear behavior of other related nanostructures.

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

Knowledge Oriented Graph Unified Transformer (KOGUT) v0.1

KOGUT — Knowledge Oriented Graph Unified Transformer KOGUT implements the Relational Graph Transformer (RelGT) architecture for knowledge graph link prediction in biological domains, with a primary focus on microbial growth media prediction. While the original RelGT (arXiv:2505.10960) targets relational tables, time series, and multi-table databases, KOGUT adapts this architecture for heterogeneous biological knowledge graphs, providing first-in-class AI predictive models for microbial cultivation. Key Adaptations Beyond Original RelGT: - Knowledge Graph Focus: Applied to biological KGs with semantic node types (taxa, chemicals, media, phenotypes, environments) versus generic relational database tables, trained on the KG-Microbe knowledge graph (1.3M entities, 2.9M edges, 24 relation types). - Multimodal Node Encoding: Integrates node labels, categories, descriptions, and synonyms from KG metadata through learned embedding layers—adapting relational column features to graph node attributes with textual semantics. - Extended K-Hop Subgraph Strategy: Optimized neighborhood sampling (3-hop default, configurable up to 200 nodes) tuned for sparse biological networks, building on the original local-global attention framework with biological relation preservation. - Biolink Predicate Preservation: Type-specific transformations for 24 biological edge semantics (occurs_in, consumes, produces, has_phenotype, subclass_of) beyond standard relational foreign keys, enabling multi-relation link prediction. - Inductive Learning Support: Enables zero-shot predictions for novel taxa through feature-based embeddings (temperature, oxygen requirements, gram stain, cell shape), extending the original transductive relational benchmark scope to uncultured microorganisms. CheapSOTA Performance Optimizations (This Distribution): - VQ-EMA Centroid Attention: Vector quantization with exponential moving average for improved global context modeling (+5-10% MRR improvement). - HDF5 Precomputed Data Loading: One-time preprocessing of k-hop subgraphs to eliminate redundant graph traversals (2-5× training speedup). - Distributed Data Parallel Training: Multi-GPU support for scaling to larger knowledge graphs (tested on 4× NVIDIA A100 GPUs at NERSC Perlmutter). - Mixed Precision Training: Automatic mixed precision (AMP) for memory efficiency and faster training. Advantages Over Standard Knowledge Graph Embedding Models: Combines RelGT's proven multi-element tokenization (features, type, hop, structure) with graph-native biological representations, enabling interpretable link prediction across heterogeneous entities that standard embedding models (TransE, RotatE, ComplEx) and table-based transformers cannot directly model. Achieves near-perfect performance on microbial growth media prediction (MRR: 0.9966, Precision@1: 0.9932, Hit@10: 1.0000) while maintaining explainability through attention-based reasoning over biological pathways. Training Data: - KG-Microbe merged knowledge graph: 1,379,337 nodes, 2,960,472 edges - 24 biological relation types including taxonomic hierarchies, metabolic interactions, phenotype associations, and environmental relationships - Primary prediction task: Growth media suitability for microbial taxa (biolink:occurs_in, 50K edges) - Multi-relation capability: Predicts links for any of the 24 relation types, including chemical consumption/production, phenotype associations, and taxonomic classification Citation: Original RelGT Architecture: Dwivedi et al., "Relational Graph Transformer", arXiv:2505.10960, 2025 KOGUT Implementation: Knowledge Oriented Graph Unified Transformer for Microbial Growth Media Prediction Developed at Lawrence Berkeley National Laboratory (LBNL) Trained on NERSC Perlmutter supercomputer

Joachimiak, Marcin [Lawrence Berkeley National Lab↗