Search NASA⌕ Search

SEARCH · Search NASA

Results for “load pattern modeling”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

High-Temperature Sapphire Pressure Sensors for Harsh Environments

The project aimed to develop and characterize an optical pressure sensor capable of withstanding high temperatures, up to 1000 °C. Survivability in harsh environments was enabled using sapphire for both the substrate and the optical fibers, therefore mitigating thermal expansion mismatch. Such material cannot be machined through conventional microfabrication techniques, and unique additive and subtractive methods had to be studied and implemented. Subtractive methods included laser-machining the sapphire substrate using a picosecond laser machining tool. At first, the parameters of importance were theoretically and experimentally evaluated such that the fluence, the pulse spacing, and the number of passes were found to strongly influence the laser ablation process. A simulation tool was designed to predict the final pattern when inputting a set of parameters and a profile. This model was complemented with laser ablation material physics models and Bayesian uncertainty analysis for one dimensional and three dimensional ablation processes. An additive process called thermocompression bonding was studied in order to seal the diaphragm substrate to the back cavity. At first, bonding strength using this technique was examined using the chevron-notch method that allow to determine the fracture toughness of this bond. This parameter is independent of the geometry used and therefore can be extended to other geometries used in actual sensor bonding. A theoretical and experimental study of internal damage caused by the laser micromachining of the sapphire material was carried out. The research identified enhanced strength and fracture toughness along laser ablated surfaces. This was determined using bend bar specimens and Vicker’s microindentations. Softening of the plasticity behavior (as indicated by nanoindentation experiments) was also investigated using finite deformation plasticity modeling via a contact mechanics finite element simulation. Good correlations between experiments and the model were achieved to determine material parameters that can be used in future MEMS based sensor design and modeling to minimize sensor failure loads. The enhanced strength and fracture toughness was determined to be due to compressive residual stress near laser ablated surfaces. This was confirmed through x-ray spectra taken along different crystal structure planes in combination and Bayesian uncertainty analysis. The Bayesian analysis inferred the entire residual strain tensor (six components) and the residual stress tensor. The results illustrated compressive in-plane stress was the most likely cause of enhanced strength and toughness after laser machining. Due to major equipment failures, a parallel route for fabricating the sensor was taken, using commercially available sapphire removal method (subtractive) and ceramic adhesive bonding (additive). The fully-packaged sensor utilized laser micromachining for the sapphire substrate fabrication and was comprised only of high-temperature-compatible materials. Calibration tests of the sensor in both room-temperature and high-temperature plane-wave tubes showed good linearity and consistent frequency response, in the tested frequency range. The sensor also showed good thermal stability in tests up to 300 °C, proving the fiber-optic lever and its packaging process to be a promising approach for high-temperature measurements. The calibration from 23 °C to 300 °C gave a raw sensitivity of 3.9 × 10 -7 V/Pa and a normalized sensitivity of 9.8 × 10 -8 Pa -1 which were consistent within the temperature range. The sensor had a resonance at 16.3 kHz and a flat band (±3 dB) up to 10.3 kHz. The noise floor of the sensor at 1 kHz was 316 nV/Hz 1/2 , resulting in a MDP of 0.8 Pa. Previously, diaphragm delamination due to CTE mismatch was found at 500 °C, causing packaging failure. Unfortunately high temperature validation in hot jet environment could not be made in time. In the future, it is suggested that the two sapphire substrates be bonded by thermal compression bonding with an intermediate platinum layer instead of using an alumina-based adhesive, and therefore achieve higher temperature capacity up to 1000 °C.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
NASA NTRS1995

An Experimental Study of Incremental Surface Loading of an Elastic Plate: Application to Volcano Tectonics

Models of surface fractures due to volcanic loading an elastic plate are commonly used to constrain thickness of planetary lithospheres, but discrepancies exist in predictions of the style of initial failure and in the nature of subsequent fracture evolution. In this study, we perform an experiment to determine the mode of initial failure due to the incremental addition of a conical load to the surface of an elastic plate and compare the location of initial failure with that predicted by elastic theory. In all experiments, the mode of initial failure was tension cracking at the surface of the plate, with cracks oriented circumferential to the load. The cracks nucleated at a distance from load center that corresponds the maximum radial stress predicted by analytical solutions, so a tensile failure criterion is appropriate for predictions of initial failure. With continued loading of the plate, migration of tensional cracks was observed. In the same azimuthal direction as the initial crack, subsequent cracks formed at a smaller radial distance than the initial crack. When forming in a different azimuthal direction, the subsequent cracks formed at a distance greater than the radial distance of the initial crack. The observed fracture pattern may explain the distribution of extensional structures in annular bands around many large scale, circular volcanic features.

Williams, K. K.↗
DOE OSTI2022

Enhancing scalability of a matrix-free eigensolver for studying many-body localization

We propose several techniques to enhance the parallel scalability of a matrix-free eigensolver designed for studying many-body localization (MBL) of quantum spin chain models with nearest-neighbor interactions and on-site disorder. This type of problem is computationally challenging because the dimension of the associated Hamiltonian matrix grows exponentially with respect to the number of spins L, and we need to average over different realizations of the random disorder to obtain relevant statistical behavior. For each disorder realization, we need to compute eigenvalues from different regions of the spectrum and their corresponding eigenvectors. In previous work, the interior eigenstates for a single eigenvalue problem are computed via the shift-and-invert Lanczos algorithm. Due to the extremely high memory footprint of the LU factorizations, this technique is not well suited for large L’s. For example, we need thousands of compute nodes on modern high performance computing infrastructures to go beyond L = 24. The matrix-free approach does not suffer from this memory bottleneck, however, its scalability is limited by a computation and communication load imbalance. To reduce this imbalance and to significantly enhance the scalability of the matrix-free eigensolver, we reorder the matrix and leverage the consistent space runtime, CSPACER. We also show its efficiency in managing irregular communication patterns at scale compared to optimized MPI non-blocking two-sided and one-sided RMA implementation variants. This effort enables us to study MBL for spin chains with a larger number of spins. The efficiency and effectiveness of the proposed algorithm is demonstrated by computing eigenstates on a massively parallel many-core high performance computer.

METIS↗
DOE OSTI2021

The impact of simultaneous inoculation of Pseudomonas aeruginosa, Staphylococcus aureus , and Candida albicans on rodent burn wounds

Burn wound infection often involves a diverse combination of bacterial and fungal pathogens. In this study, we characterize the mixed species burn wound infection by inoculating the burn surface with 1 × 10 3/4/5 CFU of Pseudomonas aeruginosa, Staphylococcus aureus, and Candida albicans in a 1:1:1 ratio. Using the revised Walker–Mason scald burn rat model, 168 male Sprague-Dawley rats (350–450 g) subject to ~10% TBSA burn injury, with or without inoculation, were evaluated for 11 days after burn. In the wound, P. aeruginosa and S. aureus formed robust biofilms as determined by the bacterial tissue load, ~1 × 10 9 CFU/g, and expression of key biofilm genes. Interestingly, within 3 days C. albicans achieved tissue loads of ~1 × 10 6 CFU/g, but its numbers were significantly reduced beyond the limit of detection in the burn wound by day 7 in partial-thickness injuries and by day 11 in full-thickness injuries. The pathogenic biofilms contributed to burn depth progression, increased release of HMGB-1 into circulation from injured tissue, and significantly elevated the numbers of circulating innate immune cells (Neutrophils, Monocytes, and Basophils). This robust model of multi-species burn wound infection will serve as the basis for the development of new antimicrobials for combating biofilm-based wound infections.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2015

Numerical Experiments of Counterflowiing Jet Effects on Supersonic Slender-Body Configurations

Previous studies have demonstrated that the use of counterflowing jets can greatly reduce the drag and heat loads on blunt-body geometries, especially when the long penetration mode jet condition can be established. Previously, the authors had done some preliminary numerical studies to determine the ability to establish long penetration mode jets on a typical Mach 1.6 slender configuration, and study its impact on the boom signature. The results indicated that a jet with a longer penetration length was required to achieve any impact on the boom signature of a typical Mach 1.6 slender configuration. This paper focuses on an in-depth parametric study, done using the space-time conservation element solution element Navier-Stokes flow solver, for investigating the effect of various counterflowing jet conditions/configurations on two supersonic slender-body models (cone-cylinder and quartic body of revolution). The study is aimed at gaining a better understanding of the relationship between the shock penetration length and reduction of drag and boom signature for these two supersonic slender-body configurations. Different jet flow rates, Mach numbers, nozzle jet exit diameters and jet-to-base diameter ratios were examined. The results show the characteristics of a short-to-long-to-short penetration-mode pattern with the increase of jet mass flow rates, observed across various counterflowing jet nozzle configurations. Though the optimal shock penetration length for potential boom-signature mitigation is tied to the long penetration mode, it often results in a very unsteady flow and leads to large oscillations of surface pressure and drag. Furthermore, depending on the geometry of the slender body, longer jet penetration did not always result in maximum drag reduction. For the quartic geometry, the maximum drag reduction corresponds well to the longest shock penetration length, while this was not the case for the cone-cylinder-as the geometry was already optimized for drag. Numerical results and assessments obtained from this parametric study along with the recommendation for future implementation of counterflowing jets as a means for drag and noise reduction are detailed in this paper.

Venkatachari, Balaji Shankar↗
DOE OSTI2025

Deconvoluting thermomechanical effects in X-ray diffraction data using machine learning

X-ray diffraction is ideal for probing the sub-surface state during complex or rapid thermomechanical loading of crystalline materials. However, challenges arise as the size of diffraction volumes increases due to spatial broadening and because of the inability to deconvolute the effects of different lattice deformation mechanisms. Here, we present a novel approach that uses combinations of physics-based modeling and machine learning to deconvolve thermal and mechanical elastic strains for diffraction data analysis. The method builds on a previous effort to extract thermal strain distribution information from diffraction data. The new approach is applied to extract the evolution of the thermomechanical state during laser melting of an Inconel 625 wall specimen which produces significant residual stress upon cooling. A combination of heat transfer and fluid flow, elasto-plasticity and X-ray diffraction simulations is used to generate training data for machine-learning (Gaussian process regression, GPR) models that map diffracted intensity distributions to underlying thermomechanical strain fields. First-principles density functional theory is used to determine accurate temperature-dependent thermal expansion and elastic stiffness used for elasto-plasticity modeling. The trained GPR models are found to be capable of deconvoluting the effects of thermal and mechanical strains, in addition to providing information about underlying strain distributions, even from complex diffraction patterns with irregularly shaped peaks.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Storage Enabled Flexibility of Conventional Generation Assets (StorFlex)

The power systems have faced progressively more demanding operational requirements over the last two decades. Several factors contribute to these challenging operating conditions, including load growth, aging infrastructure, increasing penetrations of distributed energy resources (DERs), electrification of the economy, and policy initiatives such as decarbonization. The power system and its components must provide high operational flexibility to mitigate these challenges. For example, the proliferation of intermittent DERs such as wind and solar has increased the need for conventional generation assets like hydropower plants to respond to sudden load-generation imbalances. The higher flexibility requirements for hydropower plants cause more wear and tear, potentially shortening the useful lifespan of hydropower turbines. To reduce the need for hydropower plants to follow sudden changes in the dispatch signal, we investigate their combined operation with the energy storage systems (ESSs; “ESS-based hybridization”). Our analyses focuses on improving the lifespan of hydropower plants through ESS-based hybridization. Wear and tear on hydropower turbines (particularly Francis turbines) is modeled using a loss-of-life concept that is based on damage experienced by the turbine due to various cycles of operation. Then, we show that using ESSs to offset some of the high variation increases the remaining life of the hydropower plants. To demonstrate this, a few modeling tools were developed for this work: (1) a dynamic model for various components of the turbine and its governor; (2) a control strategy that assigns a slow-varying dispatch signal to a hydropower unit versus a fastmoving signal to ESS, such that the overall power request remains the same; and (3) models for the financial analysis to quantify the economic merits of such a framework. We used the models we developed to analyze the dispatch pattern of an actual hydropower plant with a power output of 50 MW and a head height of 152 m. This work showed that ESS-based hybridization could extend the life of the hydropower plant by 5% on average. This extension in life was then used to estimate the economic benefit in terms of cost deferrals associated with hydropower plant maintenance and replacement: on average, $3.6 million. Sensitivity analysis with respect to the size of ESS and cost of turbines was performed to show the variation in benefits over the range of turbine costs and ESS sizes. Crucially, stacking damage reduction and lifetime extension with other ESS value streams such as providing ancillary services could substantially increase the financial benefits of ESS-based hybridization. The higher costs associated with ESS of appropriate size would make more financial sense when multiple value streams are stacked and co-optimized to extract the maximum benefit. This dimension will be explored in future work.

13 HYDRO ENERGY↗
NASA NTRS1992

Turbine blade tip and seal clearance excitation forces

Experimental and theoretical work done as Phase 3 of a program sponsored by MSFC to investigate the magnitude, origin, and parametric variations of destabilizing forces which arise in high power turbines due to blade-tip leakage effects are described. The two facilities which were built for this purpose are first described. The larger one is a closed, 2 atm pressurized Freon-12 flow loop in which is installed a 1:1 replica of the SSME first stage hydrogen turbine, which can be driven by the flow, and which generates about 14 KW of power into a load-absorbing DC generator. The smaller facility is used to measure the forces on labyrinth seals of the same type as those used in our turbine tests with a shrouded turbine. The seals can be kinematically whirled and spun (independently), and the inlet swirl can be set to a variety of values. Air is the working fluid (with atmospheric discharge) and the data are real-time pressure distributions in the seal glands. The five different unshrouded turbine configurations were tested with static offsets, plus one with a shroud band and a two-ridge seal. Theoretical models of various degrees of complexity were developed to help interpreting and extrapolating the data. The notion of partial work done by the fluid leaking through the tip gaps was put on a quantitative basis by examining the leakage vortex roll-up dynamics. This was used to obtain a theory of the work loss due to a uniform gap. Perturbation and multiple scale arguments were then used to extend this to the case of an eccentric turbine. This yields an unsteady, 3-D theory which can predict the distribution of the approach flow, and its effect on work defect, cross-forces, pressure patterns, and dynamic damping. The predictions agree qualitatively with the data and exhibit the correct trends, but the cross-forces are generally under-predicted.

Martinez-Sanchez, M.↗
NASA NTRS1984

Real Surface Effects in Elastohydrodynamic Lubrication

Following a general discussion of real effects in EHL, consideration is given in detail to the role played by stochastic roughness superimposed on the nominal solid boundaries. Particular attention is given to the full-film EHL regime where incipient asperity contact bears a negligible fraction of the load. Flow, from which an averaged Reynolds equation can be formed, is nonetheless modified by the amplitude and texture of the roughness patterns. In describing these effects by means of certain well-defined flow factors, it is found that only two extra parameters are needed - the rms surface height and the ratio of the correlation lengths in the two principal roughness directions. In a lowest order perturbation expansion of these flow factors in powers of the ratio of rms roughness to nominal film thickness, no other properties of the roughness appear. In most cases of practical interest, the factors describing Poiseuille flow are separable into the sum of two single-surface flow factors which means that a combination of a single equivalent rough surface versus an ideal smooth surface can always be found. For flow entrained by slip velocity the factors separate instead into a difference of single-surface factors and it becomes significant which of the two surfaces carries the equivalent roughness. Results are discussed of some applications of the averaged Reynolds equation based on the flow factor method to the EHL line contact problem. Finally, the partial EHL regime is considered where comparable load fractions are carried by the hydrodynamic film and by incipient mechanical contact. An extension of the method into this regime by combining it with asperity contact models appears most encouraging.

John H Tripp↗
NASA NTRS1999

Data Structure and Parallel Decomposition Considerations on a Fibonacci Grid

The Fibonacci grid, proposed by Swinbank and Purser (see companion abstract), provides attractive properties for global numerical atmospheric prediction by offering an optimally homogeneous, geometrically regular, and approximately isotropic discretization, with only the polar regions requiring special numerical treatment. It is a mathematical idealization, applied to the sphere, of the multi-spiral patterns often found in botanical structures, such as in pine cones and sunflower heads. Computationally, it is natural to organize the domain, into zones, in each of which the same pair, or triple, of "Fibonacci spirals" dominate. But the further subdivision of such zones into "tiles" of a shape and size suitable for distribution to the processors of a massively parallel computer requires very careful consideration if the subsequent spatial computations along the respective spirals, especially those computations (such as compact differencing schemes) that involve recursion, can be implemented in an efficient "load-balanced "manner without requiring excessive amounts of inter-processor communications. In this paper we show how certain "number theoretic" properties of the Fibonacci sequence (whose numbers prescribe the multiplicity of successive spirals) may be exploited in the decomposition of grid zones into tidy arrangements of triangular grid tiles, each tile possessing one side approximately parallel to the constant-latitude zone boundary. We also describe how the spatially recursive processes may be decomposed across such a tiling, and the directionality of the recursions reversed on alternate grid lines, to ensure a very high degree of load balancing throughout the execution of the computations required for one time step of a global model.

Michalakes, John↗
DOE OSTI2024

Observations of airflow around a supertall curved building and its impact on temperature and humidity in Houston's urban center

Characterization of realistically shaped skyscrapers embedded in non-uniform neighbourhoods experiencing intricate weather patterns remains inadequately investigated. Aiming to close this gap, the Center for Multiscale Applied Sensing team deployed its mobile observatory in the street canyons around the curved Wells Fargo Plaza skyscraper in downtown Houston, TX. Three deployments allowed airflow observations under different inflow wind and thermodynamic stability conditions. Doppler lidar measurements reveal that when inflow hits the curved wall of the skyscraper, perpendicular canyons experience similar vortex configurations creating two windward and two leeward circulations. Windsond measurements support that buoyancy within the deep street canyons can generate thermal updrafts as strong as 2 m s -1 which is sufficient to overturn the mechanical downwash under gentle wind conditions. Canyons experiencing venting during the daytime were observed to be more thermodynamically stable at night while thermodynamically stable canyons during the day were observed to be more thermodynamically unstable at night owing to the accumulation of heat near street-level. Fourier decomposition of the vertical velocity measurements shows that in all cases flow exhibited high Reynolds numbers and was composed of turbulent eddies of predominantly 6 and 15 min periods. Here, this observational dataset provides insights to assess wind load, pedestrian comfort, urban air mobility, and natural ventilation and may be used as a benchmark for numerical model and wind tunnel studies attempting to represent realistically complex urban conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Determining Compression Characteristics of Honeycomb Material - 19674

The objective is to determine compression test characteristics of stainless steel honeycomb material to be able to represent honeycomb structures accurately in analytical models used to simulate hypothetical accident scenarios of shipping packages. Honeycomb is a material used primary in the aerospace industry due to its high strength-to-weight ratio. Because it doesn't have a shelf life and can withstand high heats, it is an excellent candidate for a structural material in package designs. Honeycomb has orthotropic material properties. The honeycomb currently being evaluated is made of metal ribbons spot welded together to form hexagon pattern between metal plates. The hexagon ribbons are brazed to the metal plates. The complexity of the honeycomb significantly lengthens the simulation time, which is further complicated when the brazing, welding, and imperfections of the material is considered. This means to affectively represent the honeycomb in simulation programs, such as Abaqus, the material needs to be approximated as a uniform orthotropic material. This requires material properties in each of the three directions, T, W, and L. The T direction is defined as the 'strong' direction, perpendicular to the honeycomb sheet. The L direction is parallel the ribbon and the W direction is perpendicular to the ribbon. Honeycomb has 3 different compression stages. Stage 1 is the initial compression where the honeycomb maintains its structural integrity and does not permanently deform from loads from normal operations. Stage 2 is when the initial buckling causes irreversible damage to the honeycomb. Stage 2 is the one we are most interested in because it absorbs the most energy from a hypothetical accident scenario. During Stage 2, the honeycomb fails layer-by-layer, indicating that the more layers, the longer material crushing is sustained, and thus the more energy that is absorbed. Stage 3 is final compression, similar to compressing solid metal. During stage 3 the stress increases with a diminishing rate of elongation and the honeycomb is completely failed where the plates between the honeycomb core sandwich the crushed honeycomb ribbon. The stress-strain graph bellow shows compression test of 3 different samples in the T direction. The red lines divide the three stages. The far left is Stage 1, the middle stage 2, and the right stage 3. The data collected is displayed in the table below. The expressions in the table represent the regression formula that represents each stage of the compression on a stress-strain graph. Stage 1 is assumed to intersect with the origin; data for stage 3 for the crush in W and L directions where unable to be gathered due to the nature of failure for those directions. Stage 1 and 2 are linear regressions while stage 3 is a degree 2, to best match the curve. Linear regression is used on stage 2 because ultimately that represents the energy absorbed. With more data, a degree n x 2 regression would be a more appropriate regression for stage 2, where n is the number of layers. The table below is the approximation of each scenario and stage. Each stage starts and ends at the intersection of the next stage. There where two main objectives for these tests: firstly to understand how honeycomb performs under extreme compression, and secondly to be able to numerically represent the honeycomb structure. Both of these where completed to varying degrees. It is important understood how the honeycomb fails. If it is crushed in the T direction it retains its integrity even after being crushed, however when crushed in the L and W crush direction, if it fails, it disintegrates, and loses all integrity. Also, the more layers the more time the material spends in stage 2. Failure is started by buckling, thus if there is any imperfection, the stress will not spike but transition straight into stage 2. The data collected and aggregated can be used for initial simulation of honeycomb used in packages. The initial testing has set the ground work for more data to be collected in order to verify results and to allow more confidence in the simulation results. This is only the initial data collected. The next steps is to continue to collect more data to verify results and to test more variants of honeycomb, with different brazing, layers, and shape.

36 MATERIALS SCIENCE↗
NASA NTRS2010

Analysis of Unsteady Tip and Endwall Heat Transfer in a Highly Loaded Transonic Turbine Stage

In a previous study, vane-rotor shock interactions and heat transfer on the rotor blade of a highly loaded transonic turbine stage were simulated. The geometry consists of a high pressure turbine vane and downstream rotor blade. This study focuses on the physics of flow and heat transfer in the rotor tip, casing and hub regions. The simulation was performed using the Unsteady Reynolds-Averaged Navier-Stokes (URANS) code MSU-TURBO. A low Reynolds number k-epsilon model was utilized to model turbulence. The rotor blade in question has a tip gap height of 2.1 percent of the blade height. The Reynolds number of the flow is approximately 3x10(exp 6) per meter. Unsteadiness was observed at the tip surface that results in intermittent "hot spots". It is demonstrated that unsteadiness in the tip gap is governed by inviscid effects due to high speed flow and is not strongly dependent on pressure ratio across the tip gap contrary to published observations that have primarily dealt with subsonic tip flows. The high relative Mach numbers in the tip gap lead to a choking of the leakage flow that translates to a relative attenuation of losses at higher loading. The efficacy of new tip geometry is discussed to minimize heat flux at the tip while maintaining choked conditions. In addition, an explanation is provided that shows the mechanism behind the rise in stagnation temperature on the casing to values above the absolute total temperature at the inlet. It is concluded that even in steady mode, work transfer to the near tip fluid occurs due to relative shearing by the casing. This is believed to be the first such explanation of the work transfer phenomenon in the open literature. The difference in pattern between steady and time-averaged heat flux at the hub is also explained.

Shyam, Vikram↗
DOE OSTI2023

Stabilizing the Grid and Reducing Utility Bills Through Price-Responsive Controls for Heat Pump Water Heaters

The electricity grid is facing increasing challenges in cost-effectively balancing supply and demand. These challenges are exacerbated by increased penetration of photovoltaics, which causes mid-day overproduction, and electrification of gas appliances, which increases peak-period electricity demand. Decarbonization requires shifting building loads from fossil-intensive high-cost times to renewable-intensive low-cost times while maintaining quality of service to occupants. Utilities and ISO’s are investigating new ways of incentivizing this load shifting. One promising method is the use of Highly Dynamic Prices (HDPs). HDPs feature continuously changing prices that reflect real-time grid generation and distribution costs and capacity constraints, and thus incentivize consumers to shift their loads. California’s CPUC CalFUSE proposal and Hawaii’s recent changes demonstrate that electricity tariffs are moving towards this model. For this to work however, loads must have the capability to respond to these prices. Heat pump water heaters (HPWHs) are an ideal device for this purpose because their storage tanks decouple delivery of domestic hot water from electricity consumption. The storage enables control strategies that consume midday solar power to increase the energy stored in the tank, then provide evening peak domestic hot water services using the stored energy. Berkeley Lab's CalFlexHub project is pioneering price-driven load flexibility by developing and deploying cost-minimizing controls for many flexible loads - including HPWHs - in response to HDP. Control development is based on simulations using the Flexible Heat Pump Water Heater Performance Predictor which captures the control decisions of the on-board controller in a residential, integrated HPWH . The price-responsive controls a) shift load in ways that consume additional midday solar power to help stabilize the grid and reduce overall emissions, b) ensure that occupants receive equal or better hot water delivery service, and c) minimize the operating cost for each home in the fleet. On the grid level, the resulting shift will reduce utility operating costs and emissions, and can avoid expensive system capacity expansions. The control approach is customized to each home based on typical hot water consumption patterns. HPWH controllers, whether on the device or remotely, will receive a schedule of CTA-2045-B signals or set temperature adjustments customized to the current HDP price schedule and home. Simulation results for a fleet of 148 HPWHs on a summer day in Berkeley, California show cost savings of 29% and high price electricity consumption reductions of 80%, while maintaining full quality of service.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2001

Experimental and Numerical Investigation of Buoyancy Driven Convection During PDAMNA Thin Film Growth

This paper presents results from numerical simulations as well as laboratory experiments of buoyancy driven convection in an ampoule under varying heating and gravitational acceleration loadings. The modeling effort in this work resolves the large scale natural convective motion that occurs in the fluid during photodeposition of polydiacetelene films which is due to energy absorbed by the growth solution from a UV source. Consequently, the growth kinetics of the film are ignored in the model discussed here, and also a much simplified ampoule geometry is considered. The objective of this work is to validate the numerical prediction on the strength and structure of buoyancy driven convection that could occur under terrestrial conditions during nonlinear optical film growth. The validation is used to enable a reliable predictive capability on the nature and strength of the convective motion under low gravity conditions. The ampoule geometry is in the form of a parallelepiped with rectangular faces. The numerical results obtained from the solution to the Boussinesq equations show that natural convection will occur regardless of the orientation of the UV source with respect to the gravity vector. The least strong convective motion occurred with the UV beam directed at the top face of the parallelepiped. The strength of the convective motion was found to be almost linearly proportional to the total power of the UV source. Also, it was found that the strength of the convective motion decreased linearly with the gravity due to acceleration. The pattern of the convective flow on the other hand, depended on the source location.

Antar, Basil N.↗
NASA NTRS2007

Zero-Gravity Locomotion Simulators: New Ground-Based Analogs for Microgravity Exercise Simulation

Maintaining health and fitness in crewmembers during space missions is essential for preserving performance for mission-critical tasks. NASA's Exercise Countermeasures Project (ECP) provides space exploration exercise hardware and monitoring requirements that lead to devices that are reliable, meet medical, vehicle, and habitat constraints, and use minimal vehicle and crew resources. ECP will also develop and validate efficient exercise prescriptions that minimize daily time needed for completion of exercise yet maximize performance for mission activities. In meeting these mission goals, NASA Glenn Research Center (Cleveland, OH, USA), in collaboration with the Cleveland Clinic (Cleveland, Ohio, USA), has developed a suite of zero-gravity locomotion simulators and associated technologies to address the need for ground-based test analog capability for simulating in-flight (microgravity) and surface (partial-gravity) exercise to advance the health and safety of astronaut crews and the next generation of space explorers. Various research areas can be explored. These include improving crew comfort during exercise, and understanding joint kinematics and muscle activation pattern differences relative to external loading mechanisms. In addition, exercise protocol and hardware optimization can be investigated, along with characterizing system dynamic response and the physiological demand associated with advanced exercise device concepts and performance of critical mission tasks for Exploration class missions. Three zero-gravity locomotion simulators are currently in use and the research focus for each will be presented. All of the devices are based on a supine subject suspension system, which simulates a reduced gravity environment by completely or partially offloading the weight of the exercising test subject s body. A platform for mounting treadmill is positioned perpendicularly to the test subject. The Cleveland Clinic Zero-g Locomotion Simulator (ZLS) utilizes a pneumatic subject load device to apply a near constant gravity-replacement load to the test subject during exercise, and is currently used in conjunction with the General Clinical Research Center for evaluating exercise protocols using a bedrest analog. The enhanced ZLS (eZLS) at NASA Glenn Research Center features an offloaded treadmill that floats on a thin film of air and interfaces to a force reaction frame via variably-compliant isolators, or vibration isolation system. The isolators can be configured to simulate compliant interfaces to the vehicle, which affects mechanical loading to crewmembers during exercise, and has been used to validate system dynamic models for new countermeasures equipment designs, such as the second International Space Station treadmill slated for use in 2010. In the eZLS, the test subject and exercise device can be pitched at the appropriate angle for partial gravity simulations, such as lunar gravity (1/6th earth gravity). On both the eZLS and the NASA-Johnson Space Center standalone ZLS installed at the University of Texas Medical Branch in Galveston, Texas, USA, the subject's body weight relative to the treadmill is controlled via a linear motor subject load device (LM-SLD). The LM-SLD employs a force-feedback closed-loop control system to provide a relatively constant force to the test subject during locomotion, and is set and verified for subject safety prior to each session. Locomotion data were collected during parabolic flight and on the eZLS. The purpose was to determine the similarities and differences between locomotion in actual and simulated microgravity. Subjects attained greater amounts of hip flexion during walking and running during parabolic flight. During running, subjects had greater hip range of motion. Trunk motion was significantly less on the eZLS than during parabolic flight. Peak impact forces, loading rate, and impulse were greater on the eZLS than during parabolic while walking with a low external load (EL) and rning with a high EL. Activation timing differences existed between locations in all muscles except for the rectus femoris. The tibialis anterior and gluteus maximus were active for longer durations on the eZLS than in parabolic flight during walking. Ground reaction forces were greater with the LM-SLD than with bungees during eZLS locomotion. While the eZLS serves as a ground-based analog, researchers should be aware that subtle, but measurable, differences in kinematics and leg musculature activities exist between the environments. Aside from space applications, zero-gravity locomotion simulators may help medical researchers in the future with development of rehabilitative or therapeutic protocols for injured or ill patients. Zero-gravity locomotion simulators may be used as a ground-based test bed to support future missions for space exploration, and eventually may be used to simulate planetary locomotion in partial gravity environments, including the Moon and Mars. Figure: Zero-gravity Locomotion Simulator at the Cleveland Clinic, Cleveland, Ohio, USA

Perusek, Gail P.↗
DOE OSTI2017

Efficient Machine Learning Approach for Optimizing Scientific Computing Applications on Emerging HPC Architectures

Efficient parallel implementations of scientific applications on multi-core CPUs with accelerators such as GPUs and Xeon Phis is challenging. This requires - exploiting the data parallel architecture of the accelerator along with the vector pipelines of modern x86 CPU architectures, load balancing, and efficient memory transfer between different devices. It is relatively easy to meet these requirements for highly structured scientific applications. In contrast, a number of scientific and engineering applications are unstructured. Getting performance on accelerators for these applications is extremely challenging because many of these applications employ irregular algorithms which exhibit data-dependent control-ow and irregular memory accesses. Furthermore, these applications are often iterative with dependency between steps, and thus making it hard to parallelize across steps. As a result, parallelism in these applications is often limited to a single step. Numerical simulation of charged particles beam dynamics is one such application where the distribution of work and memory access pattern at each time step is irregular. Applications with these properties tend to present significant branch and memory divergence, load imbalance between different processor cores, and poor compute and memory utilization. Prior research on parallelizing such irregular applications have been focused around optimizing the irregular, data-dependent memory accesses and control-ow during a single step of the application independent of the other steps, with the assumption that these patterns are completely unpredictable. We observed that the structure of computation leading to control-ow divergence and irregular memory accesses in one step is similar to that in the next step. It is possible to predict this structure in the current step by observing the computation structure of previous steps. In this dissertation, we present novel machine learning based optimization techniques to address the parallel implementation challenges of such irregular applications on different HPC architectures. In particular, we use supervised learning to predict the computation structure and use it to address the control-ow and memory access irregularities in the parallel implementation of such applications on GPUs, Xeon Phis, and heterogeneous architectures composed of multi-core CPUs with GPUs or Xeon Phis. We use numerical simulation of charged particles beam dynamics simulation as a motivating example throughout the dissertation to present our new approach, though they should be equally applicable to a wide range of irregular applications. The machine learning approach presented here use predictive analytics and forecasting techniques to adaptively model and track the irregular memory access pattern at each time step of the simulation to anticipate the future memory access pattern. Access pattern forecasts can then be used to formulate optimization decisions during application execution which improves the performance of the application at a future time step based on the observations from earlier time steps. In heterogeneous architectures, forecasts can also be used to improve the memory performance and resource utilization of all the processing units to deliver a good aggregate performance. We used these optimization techniques and anticipation strategy to design a cache-aware, memory efficient parallel algorithm to address the irregularities in the parallel implementation of charged particles beam dynamics simulation on different HPC architectures. Experimental result using a diverse mix of HPC architectures shows that our approach in using anticipation strategy is effective in maximizing data reuse, ensuring workload balance, minimizing branch and memory divergence, and in improving resource utilization.

Arumugam, Kamesh↗
NASA NTRS2016

Forward Skirt Structural Testing on the Space Launch System (SLS) Program

Structural testing was performed to evaluate heritage forward skirts from the Space Shuttle program for use on the Space Launch System (SLS) program. One forward skirt is located in each solid rocket booster. Heritage forward skirts are aluminum 2219 welded structures. Loads are applied at the forward skirt thrust post and ball assembly. Testing was needed because SLS ascent loads are roughly 40% higher than Space Shuttle loads. Testing objectives were to determine margins of safety, demonstrate reliability, and validate analytical models. Two forward skirts were structurally tested using the test configuration. The test stand applied loads to the thrust post. Four hydraulic actuators were used to apply axial load and two hydraulic actuators were used to apply radial and tangential loads. The first test was referred to as FSTA-1 (Forward Skirt Structural Test Article) and was performed in April/May 2014. The purpose of FSTA-1 was to verify the ultimate capability of the forward skirt subjected to ascent ultimate loads. Testing consisted of two liftoff load cases taken to 100% limit load followed by an ascent load case taken to 110% limit load. The forward skirt was unloaded to no load after each test case. Lastly, the forward skirt was tested to 140% limit and then to failure using the ascent loads. The second test was referred to as FSTA-2 and performed in July/August of 2014. The purpose of FSTA-2 was to verify the ultimate capability of the forward skirt subjected to liftoff ultimate loads. Testing consisted of six liftoff load cases taken to 100% limit load followed by the six liftoff cases taken to 140% limit load. Two ascent load cases were then tested to 100% limit load. The forward skirt was unloaded to no load after each test case. Lastly, the forward skirt was tested to 140% limit and then to failure using the ascent loads. The forward skirts on FSTA-1 and FSTA-2 successfully carried all applied liftoff and ascent load cases. Both FSTA-1 and FSTA-2 were tested to failure by increasing the ascent loads. Failure occurred in the forward skirt thrust post radius. The forward skirts on FSTA-1 and FSTA-2 had nearly identical failure modes. FSTA-1 failed at 1.72 times limit load and FSTA-2 failed at 1.62 times limit load. This difference is primarily attributed to variation in material properties in the thrust post region. Test data were obtained from strain gages, deflection gages, ARAMIS digital strain measurement, acoustic emissions, and high-speed video. Strain gage data and ARAMIS strain were compared to finite element (FE) analysis predictions. Both the forward skirt and tooling were modeled. This allows the analysis to simulate the loading as close as possible to actual test configuration. FSTA-1 and FSTA-2 were instrumented with over 200 strain gages to ensure all possible failure modes could be captured. However, it turned out that three gages provided critical strain data. One was located in the post bore and two on the post radius. More gages were not specified due to space limitations and the desire to not interfere with the use of the ARAMIS system on the post radius. Measured strains were compared to analysis results for the load cycle to failure. Note that FSTA-1 gages were lost before failure was reached. FSTA-2 gages made it to the failure load but one of the radius gages was lost before testing began. This gage was not replaced because of the time and cost associated with disassembly of the test structure. Correlation to analysis was excellent for FSTA-1. FSTA-2 was not quite as good because there was more residual strain from previous load cycles. FSTA-2 was loaded and unloaded with 12 liftoff cases and two ascent cases before taking the skirt to failure. FSTA-1 only had two liftoff cases and one ascent case before taking the skirt to failure. The ARAMIS system was used to determine strain at the post radius by processing digital images of a speckled paint pattern. Digital cameras recorded images of the speckled paint pattern. ARAMIS strain results for FSTA-2 just prior to failure. Note a high strain location develops near the left side. This high strain compares well to analysis prediction for both FSTA-1 and FSTA-2. The strain at this location was also plotted versus limit load. Both FSTA-1 and FSTA-2 had excellent correlation between ARAMIS and analysis strains. Acoustic emission (AE) sensors were used to monitor for damage formation that may occur during testing (e.g., crack formation and growth or propagation). AE was very important because after disassembly of FSTA-1, a crack was observed in the ball fitting radius. The ball fitting did not crack on FSTA-2. AE data was used to reconstruct when the crack occurred. The AE energy versus time plot for FSTA. The energy increased considerably at 850 seconds (152% limit load), indicating a crack could have formed at this point. The only visual evidence found that could have corresponded to this was the crack that initiated in the ball fitting. The cracks in the forward skirt aluminum structures would likely have been lower energy due to a lower modulus and all that were found after failure correlated to occurring after the initial crack in the post radius. This was verified by high-speed cameras used to record the failure.

Lohrer, J. D.↗