Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic load balancing”

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.

391 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Non-linear dynamic analysis of geared systems

Under driving conditions, a typical geared system may be subjected to large dynamic loads. Also, the vibration level of the geared system is directly related to the noise radiated from the gear box. The steady state dynamic behavior of the system is examined in order to design reliable and quiet transmissions. The scope is limited to a system containing a spur gear pair with backlash and periodically time varying mesh stiffness, and rolling element bearings with clearance type nonlinearities. The internal static transmission error at the gear mesh, which is of importance from high frequency noise and vibration control view point, is considered in the formulation in sinusoidal or periodic form. A dynamic finite element model of the linear time invariant (LTI) system is developed. Effects of several system parameters, such as torsional and transverse flexibilities of the shafts and prime mover/load inertias, on free and forced vibration characteristics are investigated. Several reduced order LTI models are developed and validated by comparing their eigen solutions with the finite element model results. Using the reduced order formulations, a three degree of freedom dynamic model is developed which includes nonlinearities associated with radical clearances in the radial rolling element bearings, backlash between a spur gear pair and periodically varying gear mesh stiffness. As a limiting case, a single degree of freedom model of the spur gear pair with backlash is considered and mathematical conditions for tooth separation and back collision are defined. Both digital simulation technique and analytical models such as method of harmonic balance and the method of multiple scales were used to develop the steady state frequency response characteristics for various nonlinear and/or time varying cases.

Singh, Rajendra↗
NASA NTRSDate not supplied

Fusion of Test and Analysis: Artemis I Booster to Mobile Launcher Interface Validation

NASA is in the midst of bold and exciting next steps in human exploration and spaceflight. The designs of the new Space Launch System (SLS), the Orion spacecraft and the Exploration Ground Systems (EGS) for vehicle processing and launch are essentially complete and there has been significant progress in manufacturing and assembly of specific hardware for the Artemis I and Artemis II missions. Equally as important, the program level and integrated system level testing and analyses are also well underway to support integrated verification, validation, and Certificate of Flight Readiness (CoFR) for Artemis I. Testing and analysis are key to addressing technical challenges faced by the Artemis missions. Building block approaches are required that provide the right balance between component, element, and/or system level testing that satisfies verification and validation objectives where uncertainties are quantified and minimized. Artemis I is a system of systems that requires a fusion of test and analysis that adeptly characterizes critical interfaces between major program elements. An example of this fusion involves characterizing the interface between the SLS booster and the Mobile Launcher (ML) Vertical Support Post (VSP) interfaces. Proper characterization of this interface represents a number of challenges beginning with the fact that it is a mating of ground support structure in the form of a civil structure to flight hardware. Both sides of the interface are built to different construction standards, but are governed by interface requirements to ensure compatibility when mated. From past program experience, the flexibility at the booster to ML interface is critical in developing accurate prelaunch stacking and cryogenic preloads, squat loads, and pad separation release of preloads and squat loads. This same premise holds for Artemis I. To characterize the asymmetric characteristics at this interface, careful consideration of static forces due to gravity loading with the commensurate effects due to leveling during booster stacking (i.e., spacing and shimming) and nonlinear geometric forces are necessary for inclusion in pre-test assessments. This paper will look at these issues for the upcoming Booster Pull Test in which two boosters will be installed on the ML and one of these boosters will undergo static lateral loading followed afterwards with dynamic excitation into resonance and free-decay. This paper evaluates the booster to ML interface characteristics by characterizing the interface flexibility between the booster aft skirt and the ML VSP interfaces. Furthermore, this paper methodically evaluates the effect of the following on the test outcome: gravitational effects on the booster and ML, the effects of VSP leveling, spacing, and shimming under gravitational loading during booster stacking, the effect of geometric nonlinear follower force due to cg offset as booster is laterally displaced, and the system coupling between the booster under test, ML, and the second booster. Simulated results for a static load pull and dynamic excitation provide insight into the differences in measurement responses when boundary conditions and geometric conditions are included and not included.

Joel W Sills Jr.↗
DOE OSTI2025

Scalability analysis of heavy-duty gas turbines using data-driven machine learning

With the increasing integration of variable renewable energy sources into power systems, the role of flexible power generation technologies like gas turbines (GT) in rapid grid balancing remains crucial. This sustained importance underscores the need for scaled and precise modeling of GT to ensure effective integration within evolving energy frameworks. While physics-driven GT models integrate thermodynamics, fluid dynamics, and combustion principles, they often rely on approximate mathematical representations to accommodate scaling that may not capture the actual complex dynamics for GTs and inertial effects associated to GTs with different ratings. In this study, a data-driven model is proposed using machine learning (ML) techniques to conduct GT scalability analysis and performance evaluation with high accuracy. The ML model, trained on data from various operating conditions and performance parameters, aims to uncover intricate relationships and patterns, resembling GT characteristics at different scales (ratings). The model is developed to capture complex system interaction and to adapt to changing operational scenarios at different capacities, providing valuable insights of power system dynamics. In this study, the real-time digital simulator platform was employed to generate training data for the ML model and assess its dynamic characteristics. The ultimate objective was to develop a detailed modeling framework based on governing equations and data-driven ML capable of predicting key performance indicators, in thermal systems such as GTs, including power output, speed, fuel consumption, and exhaust temperature under diverse operating conditions at different scales. The developed ML framework demonstrated high accuracy, with mean relative errors for GT power prediction, reference speed, exhaust temperature, and compressor pressure ratio (CPR) parameters consistently below 0.1% across typical load fluctuation scenarios. Maximum deviations were limited to approximately 0.5 K for exhaust temperature and 0.009 for CPR, underscoring the model’s ability to replicating dynamic GT behavior with high precision. The adaptability of the ML model enables its application across diverse operational conditions and its extension to other thermal systems. By leveraging advanced ML techniques, this study presents a robust and scalable modeling framework that enhances GT simulation precision, facilitating improved integration into evolving power systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2013

The Martian Dust Cycle: Observations and Modeling

The dust cycle is critically important for Mars' current climate system. Suspended atmospheric dust affects the radiative balance of the atmosphere, and thus greatly influences the thermal and dynamical state of the atmosphere. Evidence for the presence of dust in the Martian atmosphere can be traced back to yellow clouds telescopically observed as early as the early 19th century. The Mariner 9 orbiter arrived at Mars in November of 1971 to find a planet completely enshrouded in airborne dust. Since that time, the exchange of dust between the planet's surface and atmosphere and the role of airborne dust on Mars' weather and climate has been studied using observations and numerical models. The goal of this talk is to give an overview of the observations and to discuss the successes and challenges associated with modeling the dust cycle. Dust raising events on Mars range in size from meters to hundreds of kilometers. During some years, regional storms merge to produce hemispheric or planet encircling dust clouds that obscure the surface and raise atmospheric temperatures by tens of kelvin. The interannual variability of planet encircling dust storms is poorly understood. Although the occurrence and season of large regional and global dust storms are highly variable from one year to the next, there are many features of the dust cycle that occur year after year. A low-level dust haze is maintained during northern spring and summer, while elevated levels of atmospheric dust occur during northern autumn and winter. During years without global-scale dust storms, two peaks in total dust loading are generally observed: one peak occurs before northern winter solstice and one peak occurs after northern winter solstice. Numerical modeling studies attempting to interactively simulate the Martian dust cycle with general circulation models (GCMs) include the lifting, transport, and sedimentation of radiatively active dust. Two dust lifting processes are commonly represented in these models: wind-stress lifting (i.e., saltation) and dust devil lifting. Although the predicted patterns of dust lifting and atmospheric dust loading from these simulations capture some aspects of the observed dust cycle, there are many notable differences between the simulated and observed dust cycles. For example, it is common for models to predict one peak in global dust loading near northern winter solstice due to excessive dust lifting in the Hellas basin at this season. Additionally, it is difficult for models to realistically capture the observed interannual variability in global dust storms. New avenues of dust cycle modeling research include exploring the effects of finite surface dust reservoirs and the effects of coupling the dust and water cycles on the predicted dust cycle.

Kahre, Melinda A.↗
NASA NTRS1995

NASA Tech Briefs, November 1995

The contents include: 1) Mission Accomplished; 2) Resource Report: Marshall Space Flight Center; 3) NASA 1995 Software of the Year Award; 4) Microbolometers Based on Epitaxial YBa2Cu3O(sub 7-x) Thin Films; 5) Garnet Random-Access Memory; 6) Fabrication of SNS Weak Links on SOS Substrates; 7) High-Voltage MOSFET Switching Circuit; 8) Asymmetric Switching for a PWM H-Bridge Power Circuit; 9) Better Ohmic Contacts for InP Semiconductor Devices; 10) Low-Bandgap Thermovoltaic Materials and Devices; 11) Digital Frequency-Differencing Circuit; 12) Imaging Magnetometer; 13) Computer-Assisted Monitoring of a Complex System; 14) Buffered Telemetry Demodulator; 15) Compact Multifunction Inspection Head; 16) Optical Detection of Fractures in Ceramic Diaphragms; 17) Eddy-Current Detection of Cracks in Reinforced Carbon/Carbon; 18) Apparent Thermal Conductivity of Multilayer Insulation; 19) Optimizing Misch-Metal Compositions in Metal Hydride Anodes; 20) Device for Sampling Surface Contamination; 21) Probabilistic Failure Assessment for Fatigue; 22) Probabilistic Fatigue and Flaw-Propagation Analysis; 23) Windows Program for Driving the TDU-850 Printer; 24) Subband/Transform MATLAB Functions for Processing Images; 25) Computing Equilibrium Chemical Compositions; 26) Program Processes Thermocouple Readings; 27) ICAN-Second-Generation Integrated Composite Analyzer; 28) Integrated Composite Analyzer with Damping Capabilities; 29) Computing Efficiency of Transfer of Microwave Power; 30) Program Calculates Power Demands of Electronic Designs; 31) Cost-Estimation Program; 32) Program Estimates Areas Required by Electronic Designs; 33) Program to Balance Mapped Turbopump Assemblies; 34) BiblioTech; 35) Controlling Mirror Tilt With a Bimorph Actuator; 36) Burst-Disk Device Simulates Effect of Pyrotechnic Device; 37) Bearing-Mounting Concept Accommodates Thermal Expansion; 38) Parallel-Plate Acoustic Absorbers for Hot Environments; 39) Adjustable-Length Strut Withstands Large Cyclic Loads; 40) Tool Indicates Contact Angles in Bearing Raceways; 41) Gravity Slides With Magnetic Braking; 42) High-Torque, Lightweight, Pneumatically Driven Wrench for Small Spaces; 43) Device for Testing Compatibility of an O-Ring; 44) Magnetic Heat Pump Containing Flow Diverters; 45) Variable-Tilt Helicopter Rotor Mast; 46) "Beach-Ball" Robotic Rovers; 47) Apparatus Would Measure Temperatures of Ball Bearings; 48) Flexible Borescope for Inspecting Ducts; 49) Texturing Copper To Reduce Secondary Emission of Electrons; 50) Automated Laser Cutting in Three Dimensions; 51) Algorithm Helps Monitor Engine Operation; 52) Flexible Revision of Data-Processing Communications; 53) Software for Managing the Use of Land; 54) Thermal Strap Increases Cryocooling Efficiency; 55) Reversible Nut With Engagement Indication; 56) Control Algorithms for Kinematically Redundant Manipulators; 57) Computed Hydrogen-Flow Splits in a Rocket Engine; 58) Pressure and Thermal Modeling of Rocket Launches; 59) Field of View of a Spacecraft Antenna: Analysis and Software; 60) Digital Controller for Laser-Beam-Steering Subsystem; 61) More About Beam-Steering Subsystem for Laser Communication; 62) Digital Controller for Laser-Beam-Steering Subsystem: Part 2; 63) Interface Circuit Board for Space-Shuttle Communications; 64) Automated Planning of Spacecraft Telecommunications; 65) Artifacts of Spectral Analysis of Instrument Readings; 66) Neural-Network Controller for Vibration Suppression; 67) Adaptive Finite-Element Computation in Fracture Mechanics; 68) Attitude Control for the Cassini Spacecraft; 69) Analytical Model for Fluid Dynamics in a Microgravity Environment; 70) Study of Rocket-Engine Joints Bonded by NVCU/NARloy-Z; 71) Improved Silicon Nitride for Advanced Heat Engines; 72) Parameters for Welding Aluminum/Lithium Alloys; 73) Lightweight Composite Intertank Structure; 74) Foil Patches Seal Small Vacuum Leaks; 75) Data Base on Cables and Connectors; 76) Effect of Clock Mode on Radiation Hardnessf an ADC; and 77) Fault-Tolerant Control for a Robotic Inspection System.

Source record↗
DOE OSTI2020

Comparison of aircraft measurements during GoAmazon2014/5 and ACRIDICON-CHUVA

The indirect effect of atmospheric aerosol particles on the Earth's radiation balance remains one of the most uncertain components affecting climate change throughout the industrial period. The large uncertainty is partly due to the incomplete understanding of aerosol–cloud interactions. One objective of the GoAmazon2014/5 and the ACRIDICON (Aerosol, Cloud, Precipitation, and Radiation Interactions and Dynamics of Convective Cloud Systems)-CHUVA (Cloud Processes of the Main Precipitation Systems in Brazil) projects was to understand the influence of emissions from the tropical megacity of Manaus (Brazil) on the surrounding atmospheric environment of the rainforest and to investigate its role in the life cycle of convective clouds. During one of the intensive observation periods (IOPs) in the dry season from 1 September to 10 October 2014, comprehensive measurements of trace gases and aerosol properties were carried out at several ground sites. In a coordinated way, the advanced suites of sophisticated in situ instruments were deployed aboard both the US Department of Energy Gulfstream-1 (G1) aircraft and the German High Altitude and Long-Range Research Aircraft (HALO) during three coordinated flights on 9 and 21 September and 1 October. Here, we report on the comparison of measurements collected by the two aircraft during these three flights. Such comparisons are challenging but essential for assessing the data quality from the individual platforms and quantifying their uncertainty sources. Similar instruments mounted on the G1 and HALO collected vertical profile measurements of aerosol particle number concentrations and size distribution, cloud condensation nuclei concentrations, ozone and carbon monoxide mixing ratios, cloud droplet size distributions, and downward solar irradiance. We find that the above measurements from the two aircraft agreed within the measurement uncertainties. The relative fraction of the aerosol chemical composition measured by instruments on HALO agreed with the corresponding G1 data, although the total mass loadings only have a good agreement at high altitudes. Furthermore, possible causes of the discrepancies between measurements on the G1 and HALO are examined in this paper. Based on these results, criteria for meaningful aircraft measurement comparisons are discussed.

47 OTHER INSTRUMENTATION↗
NASA NTRS1998

Ma'adim Vallis Estuarine Delta in Elysium Basin and Its Relevance as a Landing Site for Exobiology Exploration on Mars

The debouche of Ma'adim Vallis in the Elysium Basin generated a transitional transported sediment structure, which planimetric shape is controlled by the enclosing topography of a deep reentrant gulf of the Basin into the highland. We defined it as an estuarine delta. The location and the importance of this estuarine delta is supported by the theoretical model of graded profile constructed for Ma'adim Vallis, and by two approaches: (i) the reconstruction of Ma'adim Vallis downstream course from Gusev to Elysium Basin, and (ii) the survey of the sediment deposit in the alleged estuary. The longitudinal graded profile of Ma'adim Vallis finds its base-level in the Elysium Basin, at a about 1000 m elevation, which is in agreement with the observed Basin shoreline. This model is supported by observational evidence of flow between the northern rim of Gusev crater, and the Elysium Basin shoreline. This downstream course of Ma'adim Vallis can be divided into three hydrogeologic regions. into three hydrogeologic regions. (a) The first region is a flooded plain (Zephiria Mensae), consisting in chaotic terrain formed by highland rocks, and disintegrated lava of the western flank of Apollinaris. Morphologic indicators of the flood process are: (1) the sediment deposit over the Gusev crater northern rim that reflects the overspilling of the crater-lake water through a 40-km wide gap provided by an ancient impact crater, (2) the tear-drop shaped feature on the northeastern flank of Apollinaris Patera, and (3) the chaotic terrain that suggest the emergence of ground water generated by the seepage of the crater lake through high-permeable broken rampart material. This underground water circulation sustained by the hydrostatic pressure of the crater-lake has likely generated a hydrothermal system in the volcanic environment of Apollinaris Patera. The stratigraphy of the flooded area is identified as Hesperian age, with occurrences of Noachian hilly individual features, and as Amazonian flooded plain and chaotic material, (b) The second region is located on the western flank of Apollinaris Patera. It is surrounded by relics of deep valleys that suggest a former downstream course of Ma'adim Vallis. The geologic setting of this region (Lucus Planum) is interpreted to be an Amazonian formation composed by the middle and lower members of the Medusae Formation., c) The third region corresponds to the convergence of the west and east branches of Ma'adirn Vallis into a deep re-entrant wide gulf that penetrates about 100 km into the highland. This topographic depression is delineated by the 1000 in elevation contour. This gulf has formed an estuarine configuration centered at 3S/190W within the Elysium Basin. This configuration has favored the formation of a estuarine sedimentary delta, because of topographically controlled lateral migration. This estuarine structure is strongly dominated by the incoming supply of Ma'adim Vallis fluvial sediment extracted from Zephiria Mensae and Lucus Planum. The obtuse-angle geometry of the estuary increases the sedimentation rate, which is higher than in the course of the channel. The sediment deposition process is governed by the estuarine water circulation. The inflowing loaded fluvial water enters the estuary as a bottom current, and mixes with the relatively less-loaded water of the receiving basin. When they mixed. the inflowing fluvial material, and the landward basin circulating water generate an accumulation of highly-diversified estuarine deposit stratification. This accumulation of material is mostly centered in the transitional zone of the delta. The sediment trapping efficiency of the estuary is function of the energy balance between the inflowing fluvial water, and the ingoing basin current. The submergence of the delta by the rising of the water-level increases the estuary water-depth, and consequently the sediment entrapment is favored. The locus of sediment accumulation moves landward in the zone of inflowing fluvial water. This results in the rising of the channel base-level, thus in the increase of the length of the longitudinal graded-profile. The sediment deposit facies of the zone A shows a generally smooth surface. The longitudinal deposit is bordered by alluvial terraces that reflect the variations of the channel level. The waning of the Elysium Basin caused the erosion of the Basin estuarine zone by small channels, this episode being characterized by dissected tear-drop shaped mesa-like morphologies in the delta. Our estuarine delta model predicts a lithostratigraphic depositional sequence associated with the water submergence and the transgression of Elysium Basin. The thickness of the estuarine sediment corresponds to the Elysium Basin levels changes relatively to the bed floor of the estuary, The depositional sequence of Ma'adim Vallis are described: (1) a pro-current filled region (A), where fluvial are longitudinally accumulated by the inflowing water, (2) inverse current from Elysium Basin (B), where fluvial and lacustrine sediments are accumulated, and (3) zone of current equilibrium (C), where the sediments are distributed as a shoreline at the boundary of the estuarine delta. The estuary sedimentology dynamics collects and keeps the record of the geologic unit material crossed by Ma'adim Vallis, and those of the lakebed deposit of Elysium Basin. The predicted mixed stratigraphic sequence from fluvial and lacustrine sediment makes this site an exceptional environment to concentrate potential multi-origin biologic records. We envision four possible strategies to explore this sedimentologic record: (1) longitudinal surface and subsurface traverses in region A to investigate outcrop levees, (2) exploration of the mesa walls in region B, (3) deep drilling hole lodging of the sequential deposits in the zones A and B, and (4) surface and subsurface exploration of the shoreline delta. The expected results for each of these strategies are: (1) in the deepest layers of region A are predicted frequent and abundant coarse material, sandy lenses lamination grading downward from sand to cobbles. Volcanic debris from the Noachian crustal Plateau unit material, hydrothermal altered rocks, carbonates, Hesperian and possibly Amazonian volcanic material, from Apollinaris Patera, altered rocks and carbonates from Zephiria Mensae are expected. As a favorable environment for inception of life, possible biological records are expected in transported rock, (2) At the surface, and subsurface (<=100 m), large deposits sandy to silted material from Elysium paleolake basin mixed with fine-grained sediments from Ma'adim Vallis are expected mostly in the upstream part of region B, (3) on the shoreline of the estuarine delta, abundant fine material from Elysium paleolake basin (evaporites, carbonates), mostly Amazonian in age are expected. The Ma'adini estuary is a favorable landing site for all the above mentioned science aspects, and .for its location. The site lies near the equator, which is favorable for the rover solar power supply, and at 1000m elevation, which is a favorable configuration for the descent system braking. Another advantage is the extent of the area of high scientific interest (33,000 sq km), which provides a good ellipse, and potential long study traverses.

Grin, E. A.↗
DOE OSTI2023

Component level modeling of materials degradation for insights into operational flexibility of Existing Coal Power Plants

Increasingly, coal-fired power plants are required to balance power grids by compensating for the variable electricity supply from renewable energy sources. Fossil-fueled power plants, originally designed to be base loaded, will increasingly need to operate on a load following or cyclic basis. This demanding requirement for operational flexibility needs insights into accelerated material degradation arising due to the harsh operating conditions (e.g., fatigue, early oxide exfoliation due to stresses) along with current damage mechanisms (fireside corrosion, creep and erosion) observed in service. Our research objective is to develop component level modeling toolkit for materials-based degradation for two key mechanisms that can accelerate with cyclic operations. In more detail, this includes the fireside corrosion/steam oxidation/erosion/creep/fatigue of superheaters/reheaters and steam pipework and also the water droplet erosion/ fatigue of last stage steam turbine blades degradation mechanisms, that demand routine and sometimes unplanned maintenance and repair. The innovation is in developing a computational fluid dynamics/finite element (CFD/FE) modeling toolkit for the component level models of the boilers and low-pressure steam turbines in coal power plants that can tackle multidisciplinary failure mechanisms occurring concurrently for extreme environment materials. Lifetime assessment in such environments also needs to account for the unit-specific analyses, operational history and fuel feedstock; this can only be obtained by destructive analysis of components. This, in turn, enables validation of the model toolkits utilizing service feedback data, improving the probability of time/temperature dependent life prediction.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2026

A GPU-based compressible combustion solver for applications exhibiting disparate space and time scales

High-speed chemically active flows pose significant computational challenges due to their disparate space and time scales, with stiff chemistry often dominating simulation time. While modern scientific computing programs achieve exascale performance by leveraging graphics processing units (GPUs), existing GPU-based compressible combustion solvers face critical limitations in memory management, load balancing, and handling the highly localized nature of chemical reactions. To this end, we present a high-performance compressible reacting flow solver built on the AMReX framework and optimized for multi-GPU settings. Here, our approach addresses three GPU performance bottlenecks: memory access patterns through column-major storage optimization, computational workload variability via a bulk-sparse integration strategy for chemical kinetics, and multi-GPU load distribution for adaptive mesh refinement applications. The solver adapts existing matrix-based chemical kinetics formulations to multi-grid contexts. Using representative combustion applications, including 2D and 3D detonations and a 3D jet-in-crossflow configuration, we demonstrate 1.4–5× performance improvements over initial implementations on an in-house cluster of NVIDIA H100 GPUs, and near-ideal weak scaling on the Frontier supercomputer (Oak Ridge Leadership Computing Facility) with up to 1024 AMD Instinct MI250X GPUs. Roofline analysis reveals substantial improvements in arithmetic intensity for both convection (∼ 10 ×) and chemistry (∼ 4 ×) routines, confirming efficient utilization of GPU memory bandwidth and computational resources.

42 ENGINEERING↗
NASA NTRS2004

Design of Timer Circuit for Dynamic Data System

The Branch That I work in is in the Aero Electronic Test Branch, which is part of the Research and Testing Division. The Aero Electronic Test Branch deals with electronic control and instrumentation systems. This branch supports the research and test study of wind tunnels such as the l0x10,9x15, and 8x6. Wind tunnels are used in research to test certain parts of a jet, plane, shuttle or any other flying object in certain test conditions. My assignment is to design a programmable trigger circuit on a 19 standard rack mount that will allow the circuit to latch and hold for a predefined amount of time entered by the user when receiving a signal. It should then re-arm itself within 0.25 seconds after the time is finished. The time should be able to be seen on a display showing the time entered. The time range has to be from 0-600 seconds in 0.01 second increments (600.00). From the information given, counters will be needed to design and build this circuit. A counter, in it s simplest form, is a group of flip flops that can temporarily store bits of information put into the circuit. They can be constructed in many different ways, such as in 4 flip flops (4-bit counter) or 8 flip flops and even higher. Counters are usually cascaded with other counters to reach higher bits, such as 16 or 24 bit counters. The application in which I will use the counters will be to count down from any programmable number that I input either by a keyboard or a thumbwheel. Also, I will use counters that will be used specifically as a frequency divider to divide the pulses that enter the circuit through an input signal from a crystal clock. The pulses will need to be divided so that it will function as a 100Hz clock putting out 100 pulses per second. A switch will be used to load my inputs in and more than likely a button also so that I can stop and hold the count at any point of time. I will use 5 BCD up/down programmable counters, and a certain amount (depending on what kind of "divide by N" counter I use) of frequency dividing counters for the assignment. After the design is carefully made, a task order will be written and then given to the manufacturer to create a rack mount circuit board that will match my specifications given. The applications in which this design will be used for is in the use of the six-component balance signal conditioner for measurement and electronic system control. It can be used as a timer system for the balance signal conditioner in which it does numerous tests for the Wind tunnel research, in which a preset time can be set for how long it performs its tests. Specifically, my design should be applied to the balance signal conditioner used for the 8x6 wind tunnel research. Hopefully this design should aid in more efficient research for the 8x6 wind tunnel.

Young, Nathaniel, III↗
NASA NTRS1999

Computer Modeling of Sand Transport on Mars Using a Compart-Mentalized Fluids Algorithm (CFA)

It has been postulated that aeolian transport on Mars may be significantly different from that on Earth. From laboratory experiments simulating martian grain transport [2], it has been observed that (saltating) grains striking the bed can cause hundreds of secondary reptation trajectories when impact occurs at speeds postulated for Mars. Some of the ballistically induced trajectories "die ouf' and effectively join the ranks on the creep population that is merely nudged along by impact. Many of the induced reptation trajectories, however, are sufficiently high for the grains to become part of the saltation load (it is irrelevant to the boundary layer how a grain attained its initial lift force). When these grains, in turn, strike the surface, they too are capable of inducing more reptating grains. This cascading effect has been discussed in connection with terrestrial aeolian transport in an attempt to dispel the notion that sand motion is divisible only into creep and saltation loads. On Earth, only a few grains are splashed by impact. On Mars, it may be hundreds. We developed a computer model to address this phenomenon because there are some important ramifications: First, this ratio may mean that martian aeolian transport is dominated by reptation flux rather than saltation. On Earth, the flux would be a roughly balanced mixture between reptation/creep and saltation. On Venus, there would be no transport other than by saltation. In other words, an understanding of planetary aeolian processes may not be necessarily understood by extrapolating from the "Earth case", with only gravity and atmospheric density/viscosity being considered as variables. Second, the reptation flux on Mars may be self sustaining, so that little input is required by the wind once transport has been initiated. The number of grains saturating the boundary layer near the bed may mean that average grain speed on Mars might conceivably be less than that on Earth. This would say much for models of sand comminution on Mars. A multiple-grain transport model using just the equations of grain motion describing lift and drag is impossible to develop owing to stochastic effects --the very effects we wish to model. Also, unless we were to employ supercomputing techniques and extremely complex computer codes that could deal with millions of grains simultaneously, it would also be difficult to model grain transport if we attempted to consider every grain in motion. No existing computer models were found that satisfactorily used the equations of motion to arrive at transport flux numbers for the different populations of saltation and reptation. Modeling all the grains in a transport system was an intractable problem within our resources, and thus we developed what we believe to be a new modeling approach to simulating grain transport. The CFA deals with grain populations, but considers them to belong to various compartmentalized fluid units in the boundary layer. In this way, the model circumvents the multigrain problem by dealing primarily with the consequences of grain transport --momentum transfer between air and grains, which is the physical essence of a dynamic grain-fluid mixture. We thus chose to model the aeolian transport process as a superposition of fluids. These fluids include the air as well as particle populations of various properties. The prime property distinguishing these fluids is upward and downward grain motion. In a normal saltation trajectory, a grain's downwind velocity increases with time, so a rising grain will have a smaller downwind velocity than a failing grain. Because of this disparity in rising and falling grain proper-ties, it seemed appropriate to track these as two separate grain populations within the same physical space. The air itself can be considered a separate fluid superimposed within and interacting with the various grain-cloud "fluids". Additional informaiton is contained in the original.

Marshall, J.↗
DOE OSTI2024

HydraGNN_Predictive_GFM_2024 - Ensemble of predictive graph foundation models for ground state atomistic materials modeling

We provide the ensemble of fifteen pre-trained graph foundation models (GFMs) for atomistic materials modeling applications. Each one of the fifteen GFMs has been trained on five open-source datasets that (once aggregated) amount to over 154 million atomistic structures, which cover over two-thirds of the natural elements of the periodic table and that comprises a broad set of organic and inorganic compounds. This vast set of atomistic structures comprises ground state configurations that are dynamically stable (i.e., equilibrated structures with atomic forces approximately close to zero values) as well as dynamically unstable structures (i.e., non-equilibrium structures with non-negligible non-zero values of atomic forces). The ensemble of datasets aggregated does NOT include excited states. The datasets have been curated to remove atomistic structures with spectral norm of the force tensor above 100 eV/angstrom. Moreover, a linear term of the energy was computed for each dataset using a linear regression model that uses the chemical concentration of each natural element as regressor. The linear term predicted by the linear regression model has been subtracted from each original energy value to perform a re-alignment of the energy values across different electronic structures approximation theories performed to generate the diverse multi-source, multi-fidelity datasets. The folder "ADIOS_files" contains the set of pre-processed datasets in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used for the development and training of GFMs in this work. The "ADIOS_files" directory contains 6 sub-directories named as follows: - ANI1x-v3.bp - MPTrj-v3.bp - OC2020-20M-v3.bp - OC2020-v3.bp - OC2022-v3.bp - qm7x-v3.bp Each sub-directory contains the pre-processed datasets converted in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used to the development, training, and performance testing of the ensemble go predictive graph foundation models. Each GFM was developed using HydraGNN (https://github.com/ORNL/HydraGNN) as underlying graph neural network (GNN) architecture. The multi-task learning (MTL) capability of HydraGNN was used to simultaneously train the GFMs on labeled values for direct predictions of energy (a total system property of an atomistic structure that measures the chemical stability) and atomic forces (an atomic level property of an atomistic structure that measures the dynamical stability). The hyper parameters of the GFM have been tuned using scalable hyperparameter optimization (HPO) algorithms implemented in the software DeepHyper (https://github.com/deephyper/deephyper). The pre-training of each HPO trial was performed using distributed data parallelism (DDP) to scale the training across 128 compute nodes of the exascale OLCF supercomputer Frontier. Each HPO trial was trained only for 10 epochs and an early stopping was performed to avoid wasting significant computational resources on GNN architectures that were clearly underperforming. For each HPO trial, the 'omnistat' tool developed by (AMD Research - Advanced Micro Device) was used to measure the total energy consumption in kWh. The ensemble of GFMs was obtained by selecting the fifteen best performing HPO trials. Four models have been selected for their clear advantage in accuracy, and these are the GFMs with IDs 229, 156, 147, 260. Additional eleven models have been selected based on judicious balance between accuracy and energy consumption needed for training, and these are the GFMs with IDs 165, 78, 137, 1, 175, 171, 181, 67, 179, 167, 351. Each selected GFM of the ensemble was continued to cumulate a total of at most 30 epochs. In some cases, the total number of epochs actually performed was les than 30 due to two combined factors: (1) the size of the GFM (i.e., the number of model parameters to train) and (2) the total wall-clock time for which the computational resources could be allocated on OLCF-Frontier. The "Ensemble_of_models" directory contains 15 sub-directories named as follows: - gfm_0.229 - gfm_0.156 - gfm_0.147 - gfm_0.260 - gfm_0.165 - gfm_0.78 - gfm_0.137 - gfm_0.1 - gfm_0.175 - gfm_0.171 - gfm_0.181 - gfm_0.67 - gfm_0.179 - gfm_0.167 - gfm_0.351 Each one of these sub-directories refers to one of the fifteen HPO trials that have been selected to continue the pre-training with at most 30 epochs. With each sub-directory associated with a specific HPO trial, the following files can be found: - config.json: file for argument parsing to develop and train an HydraGNN architecture - gfm_0.ID_epoch_N.pk: file with model parameters for HPO ID trial after N epochs of training The ensemble of fifteen GFM architectures was used for (1) ensemble averaging to stabilize the predictions of energy and atomic forces after pre-training for post-processing analysis and (2) ensemble uncertainty quantification (UQ). The code used to develop, pre-train, and load the pre-trained models for post-processing analysis is available on the ORNL-GitHub at the following link: https://github.com/ORNL/HydraGNN/tree/Predictive_GFM_2024

36 MATERIALS SCIENCE↗
NASA NTRS1994

Accelerated panel methods using the fast multipole method

Panel methods are commonly used in computational fluid dynamics for the solution of potential flow problems. The methods are a numerical technique based on the surface distribution of singularity elements. The solution is the process of finding the strength of the singularity elements distributed over the body's surface. This process involves the solution of the matrix problem Pq = p' for a set of unknowns q. The Fast Multipole Method is used to directly compute q without using matrix solvers. The algorithm works in O(N) time for N points, a great improvement over standard matrix solvers. In panel methods, the surface of a body is divided into a series of quadrilateral panels. The methods involve the computation of the influence of all other panels on each individual panel. The influence is based on the surface distribution, though this can be approximated by the area for distant panels. An alternative approximation, though with arbitrary accuracy, is to develop a multipole expansion about the center of the panel to describe the effect of a given panel on distant points in space. The expansion is based on the moments of the panel, thus allow the use of various surface distributions without changing the basic algorithm, just the computation of the various moments. The expansions are then manipulated in a tree walk to develop Taylor series expansions about a point in space which describe the effect of all distant panels on any point within a volume of convergence. The effect of near panels then needs to be computed directly, but the effect of all distant panels can be computed by simply evaluating the resulting expansion. The Fast Multipole Method has been applied to panel methods for the solution of source and doublet distributions. A major feature of the algorithm is that the algorithm does not change to derive the potential and velocity for sources and doublets. The same expansions can be used for both sources and doublets. Since the velocity is related to the potential, and the doublet potential is related to the z-component of the source velocity, all values can be derived from the same expansion by taking a series of partial derivatives. This requires more expansion terms to be kept since terms are lost in the process of taking partial derivatives. Thus to maintain accuracy for the doublet computation, more terms are required than if just evaluating for sources. The resulting Fast Multipole code should then parallelize better than classical panel methods due to the locality of data dependencies found in the Fast Multipole Method. Theoretically the parallelized code should execute in O(log N) time with O(N) processors, though this is not practical. Ongoing work includes implementing the parallel accelerated panel method, including methods to improve the load balancing of the problem by taking advantage of the known geometry of panels, and to encorporate sensitivity analysis into the algorithm.

Leathrum, James F., Jr.↗