Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiple solutions”

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 415 records · Page 23

A MultiGPU Performance-Portable Solution for Array Programming Based on Kokkos

Today, multiGPU nodes are widely used in high-performance computing and data centers. However, current programming models do not provide simple, transparent, and portable support for automatically targeting multiple GPUs within a node on application areas of array programming. In this paper, we describe a new application programming interface based on the Kokkos programming model to enable array computation on multiple GPUs in a transparent and portable way across both NVIDIA and AMD GPUs. We implement different variations of this technique to accommodate the exchange of stencils (array boundaries) among different GPU memory spaces, and we provide autotuning to select the proper number of GPUs, depending on the computational cost of the operations to be computed on arrays, that is completely transparent to the programmer. We evaluate our multiGPU extension on Summit (#5 TOP500), with six NVIDIA V100 Volta GPUs per node, and Crusher that contains identical hardware/software as Frontier (#1 TOP500), with four AMD MI250X GPUs, each with 2 Graphics Compute Dies (GCDs)for a total of 8 GCDs per node. We also compare the performance of this solution against the use of MPI + Kokkos, which is the cur-rent de facto solution for multiple GPUs in Kokkos. Our evaluation shows that the new Kokkos solution provides good scalability for many GPUs and a faster and simpler solution (from a programming productivity perspective) than MPI + Kokkos.

Valero Lara, Pedro↗

Correlations in alpha-alpha scattering and semi-classical optical models

We show the equivalence of semi-classical solutions to optical model coupled-channel equations derived from Watson's form of the nucleus-nucleus multiple-scattering series to the Glauber multiple-scattering series. A second-order solution to the semi-classical coupled-channel elastic amplitude is shown to be nearly equivalent to a second-order optical-phase-shift approximation to the Glauber amplitude if the densities of all nuclear excited states are approximated by the ground-state density. Using the Jastrow method to model the two-body density we find an average excited-state density to be of negligible importance in the double-scattering region of alpha-alpha scattering.

NASA Discipline Number 04-10↗

Evaluating Portable Parallelization Strategies for Heterogeneous Architectures in High Energy Physics

High-energy physics (HEP) experiments have developed millions of lines of code over decades that are optimized to run on traditional x86 CPU systems. However, we are seeing a rapidly increasing fraction of floating point computing power in leadership-class computing facilities and traditional data centers coming from new accelerator architectures, such as GPUs. HEP experiments are now faced with the untenable prospect of rewriting millions of lines of x86 CPU code, for the increasingly dominant architectures found in these computational accelerators. This task is made more challenging by the architecture-specific languages and APIs promoted by manufacturers such as NVIDIA, Intel and AMD. Producing multiple, architecture-specific implementations is not a viable scenario, given the available person power and code maintenance issues. The Portable Parallelization Strategies team of the HEP Center for Computational Excellence is investigating the use of Kokkos, SYCL, OpenMP, std::execution::parallel and alpaka as potential portability solutions that promise to execute on multiple architectures from the same source code, using representative use cases from major HEP experiments, including the DUNE experiment of the Long Baseline Neutrino Facility, and the ATLAS and CMS experiments of the Large Hadron Collider. This cross-cutting evaluation of portability solutions using real applications will help inform and guide the HEP community when choosing their software and hardware suites for the next generation of experimental frameworks. We present the outcomes of our studies, including performance metrics, porting challenges, API evaluations, and build system integration.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Numerical Simulations of STOVL Hot Gas Ingestion in Ground Proximity Using a Multigrid Solution Procedure

A multi grid solution procedure for the numerical simulation of turbulent flows in complex geometries has been developed. A Full Multigrid-Full Approximation Scheme (FMG-FAS) is incorporated into the continuity and momentum equations, while the scalars are decoupled from the multi grid V-cycle. A standard kappa-Epsilon turbulence model with wall functions has been used to close the governing equations. The numerical solution is accomplished by solving for the Cartesian velocity components either with a traditional grid staggering arrangement or with a multiple velocity grid staggering arrangement. The two solution methodologies are evaluated for relative computational efficiency. The solution procedure with traditional staggering arrangement is subsequently applied to calculate the flow and temperature fields around a model Short Take-off and Vertical Landing (STOVL) aircraft hovering in ground proximity.

Wang, Gang↗

Parallelized Quadrupole Simulations of Thermographic Responses of Composites

Thermography has been shown to be a viable technique for inspection of composites. Model inversion of the thermography data requires a fast method for performing the forward problem. Viable numerical methods for the thermal response forward problem are finite element, finite difference and the quadrupole method. Normally both the finite element and finite difference methods solve for the thermal response in the time domain which limits one’s ability to increase the speed of the simulation by parallelization. In contrast, the quadrupole method solves for the Laplace transform of the thermal response. One of the features of the Laplace transform methodology is the solution at any discrete time is independent of the solution at all other times. Therefore, it is easy to separate into a set of independent calculations with each of the times of interest being performed in parallel. Additionally, the numeric inversion of the Laplace transform typically involves numerically solving for the Laplace transform at multiple Laplace frequencies. Each of those solutions are also independent of solutions at other frequencies and can be calculated in parallel. By parallelization of this method, it is possible to perform the simulations of three-dimensional configurations in seconds. When the input stimulus for thermal response is a delta function heat flux (a reasonable approximation for flash heating), the thermal response is smooth. For this case, it is possible to accurately estimate the thermal response at any time within a given time interval from a set of simulations separated by exponentially increasing time steps. From these simulations, it is possible to accurately interpolate to find the response at intermediate times by a spline interpolation of the logarithm of time versus logarithm of temperature. The thermal response with exponential time stepping is shown to produce values for the thermal response which are within 1% of values within the time interval. The simulations are compared to finite element simulations of the same inspection configurations. The simulations are also compared to the thermographic measurements on composites where shape and depth of the delaminations are obtained from other inspection methods.

Thermography↗

Parallelized Quadrupole Simulations of Thermographic Responses of Composites

Thermography has been shown to be a viable technique for inspection of composites. Model inversion of the thermography data requires a fast method for performing the forward problem. Viable numerical methods for the thermal response forward problem are finite element, finite difference and the quadrupole method. Normally both the finite element and finite difference methods solve for the thermal response in the time domain which limits one’s ability to increase the speed of the simulation by parallelization. In contrast, the quadrupole method solves for the Laplace transform of the thermal response. One of the features of the Laplace transform methodology is the solution at any discrete time is independent of the solution at all other times. Therefore, it is easy to separate into a set of independent calculations with each of the times of interest being performed in parallel. Additionally, the numeric inversion of the Laplace transform typically involves numerically solving for the Laplace transform at multiple Laplace frequencies. Each of those solutions are also independent of solutions at other frequencies and can be calculated in parallel. By parallelization of this method, it is possible to perform the simulations of three-dimensional configurations in seconds. When the input stimulus for thermal response is a delta function heat flux (a reasonable approximation for flash heating), the thermal response is smooth. For this case, it is possible to accurately estimate the thermal response at any time within a given time interval from a set of simulations separated by exponentially increasing time steps. From these simulations, it is possible to accurately interpolate to find the response at intermediate times by a spline interpolation of the logarithm of time versus logarithm of temperature. The thermal response with exponential time stepping is shown to produce values for the thermal response which are within 1% of values within the time interval. The simulations are compared to finite element simulations of the same inspection configurations. The simulations are also compared to the thermographic measurements on composites where shape and depth of the delaminations are obtained from other inspection methods.

Thermography↗

High-throughput approach for investigating interdiffusion in medium- and high-entropy alloys

Interdiffusion experiments are usually time-consuming and tedious since diffusion couples must be annealed at several temperatures for a long time. The efforts required to study interdiffusion in multicomponent alloys increase dramatically as multiple diffusion couples are required to cover broad composition ranges and determine the diffusivities of individual elements in different chemical environments. To circumvent this challenge, we present a high-throughput approach applicable to single-phase and compositionally complex alloys, which are assumed to approximate ideal solid solutions. Here, a simple diffusion-multiple experiment combined with a physically based kinetic model is proposed to efficiently determine the diffusion coefficients of the constituent elements in quaternary CrFeCoNi alloys. Compared with tracer diffusivities reported in the literature, the results, thus, obtained do not differ by more than a factor of 2 and were obtained from a single interdiffusion experiment. In contrast, the diffusivities simulated with commercial mobility and thermodynamic databases are strongly overestimated by a factor ranging from 1 to 16. Therefore, our approach enables high-throughput determination of diffusivities and can help in the design of alloys for high-temperature applications where diffusion plays a key role.

36 MATERIALS SCIENCE↗

Chemical Robotics Enabled Exploration of Stability in Multicomponent Lead Halide Perovskites via Machine Learning

Metal halide perovskites have attracted immense interest as a promising material for a variety of optoelectronic and sensing applications. However, issues regarding long-term stability have emerged as the key bottleneck for commercialization. Here, we develop an automated experimental workflow based on combinatorial synthesis and rapid throughput characterization to explore long-term stability of these materials in ambient conditions. We apply it to four model perovskite systems: MA x FA y Cs 1–x–y PbBr 3 , MA x FA y Cs 1–x–y PbI 3 , Cs x FA y MA 1–x–y Pb(Br x+y I 1–x–y ) 3 , and Cs x MA y FA 1–x–y Pb(I x+y Br 1–x–y ) 3 . Non-negative matrix factorization and Gaussian process regression are used to interpolate the photoluminescent behavior of the phase diagram. This interpolative regression analysis helps to distinguish mixtures that form solid solutions from those that segregate into multiple materials, pointing out the most stable regions of the phase diagram. Here, we find stability dependence on composition to be nonuniform within the composition space, suggesting the presence of potential preferential compositional regions. This proposed workflow is universal and can be applied to other solution-processable materials.

36 MATERIALS SCIENCE↗

Crack Offset Measurement With the Projected Laser Target Device

The device and associated analysis methodology summarized in this report were developed for the purpose of estimating the size of discontinuities in the surface of the foam that covers the Space Shuttle External Tank. These surface offsets are thought to be due to subsurface cracks in the foam insulation. The mathematical analysis and procedure described here provide a method to quantity the dimensions of the crack offset in a direction perpendicular to the surface, making use of the projected laser target device (PLTD) tool and a laser line projector. The keys to the construction and use of the PLTD are the following geometrical design requirements: Laser dots are on a square grid: length on a side. Laser beams are perpendicular to projected surface. Beams are parallel out to the distance being projected. The PLTD can be used to (1) calibrate fixed cameras of unknown magnification and orientation (far-field solution); (2) provide equivalent calibration to multiple cameras, previously achieved only by the use of known target points, for example, in 3.D foreign-object debris tracking on a fixed launch platform; (3) compute scaling for conventional 2.D images, and depth of field for 3.D images (near-field solution); and (4) in conjunction with a laser line projector, achieve accurate measurements of surface discontinuity (cracks) in a direction perpendicular to the surface.

Source record↗

Hydrologic connectivity and dynamics of solute transport in a mountain stream: Insights from a long-term tracer test and multiscale transport modeling informed by machine learning

The movement of solutes in a watershed is a complex process with multiple interactions and feedbacks across spatial and temporal scales. Modeling the dynamics of solute transport along diverse hydrologic pathways within watersheds – from hillslopes to stream channels and in and out of the hyporheic zones – is challenging but critically important, as these processes integrate and contribute to the biogeochemical functioning of the river corridor up to the river network scale. Here we use results from a long-term network-scale tracer test at the H.J. Andrews experimental forest in western Cascade Mountains, Oregon, USA to inform a multiscale framework for transport in stream corridors. The framework uses a Lagrangian-based subgrid model to represent the effects of hyporheic exchange flow and advective transport at stream network scales. The spatially and temporally resolved stream discharge needed for the transport model is imputed across the river system by an entity-aware long short-term memory network. Modeled concentrations show good agreements with the observations and exhibit power scaling laws indicative of a very wide range of timescales over which hyporheic exchange flow occurs. Our results demonstrate a data-informed modeling framework that links dynamical processes occurring at small scales to a network context to help understand how changes at reach scale cascade into network-scale effects, providing a useful tool for sustainable river basin management.

54 ENVIRONMENTAL SCIENCES↗

A simple analytical function for bidirectional reflectance

A simple bidirectional reflectance model based on physical scattering laws is developed to calculate the bidirectional reflectance of a wide variety of surfaces. Chandrasekhar's radiative transfer solution, obtained for the anisotropically scattering semi-infinite medium, is used to compute the multiple-scattered radiances. However, since the radiative transfer solution does not take architectural effects into consideration, we adopt Hapke's (1986) approach and add an empirical term to explain the hot-spot phenomenon and use Cox and Munk's (1954) formulation to take into account the specular reflection. The physical parameters of the model are retrieved from bidirectional reflectance measurements. For each surface, only one set of model parameters is needed for application to all illumination and viewing geometries. The validity of the model is established by comparing the computed and measured reflectances for sets of viewing and illumination angles that were not included in the inversion algorithm. Good agreement is shown between the model-computed and the observed reflectances for dense prairie vegetation canopies, a sparse desert scrub community, a plowed agricultural field, and an alkali flat.

Ahmad, Suraiya P.↗

The Effect of Stabilization Heat Treatments on the Tensile and Creep Behavior of an Advanced Nickel-Based Disk Alloy

As part of NASA s Advanced Subsonic Technology Program, a study of stabilization heat treatment options for an advanced nickel-base disk alloy, ME 209, was performed. Using a simple, physically based approach, the effect of stabilization heat treatments on tensile and creep properties was analyzed in this paper. Solutions temperature, solution cooling rate, and stabilization temperature/time were found to have a significant impact on tensile and creep properties. These effects were readily quantified using the following methodology. First, the effect of solution cooling rate was assessed to determine its impact on a given property. The as-cooled property was then modified by using two multiplicative factors which assess the impact of solution temperature and stabilization parameters. Comparison of experimental data with predicted values showed this physically based analysis produced good results that rivaled the statistical analysis employed, which required numerous changes in the form of the regression equation depending on the property and temperature in question. As this physically based analysis uses the data for input, it should be noted that predictions which attempt to extrapolate beyond the bounds of the data must be viewed with skepticism. Future work aimed at expanding the range of the stabilization/aging parameters explored in this study would be highly desirable, especially at the higher solution cooling rates.

Gayda, John↗

Design Optimization of a Variable-Speed Power Turbine

NASA's Rotary Wing Project is investigating technologies that will enable the development of revolutionary civil tilt rotor aircraft. Previous studies have shown that for large tilt rotor aircraft to be viable, the rotor speeds need to be slowed significantly during the cruise portion of the flight. This requirement to slow the rotors during cruise presents an interesting challenge to the propulsion system designer as efficient engine performance must be achieved at two drastically different operating conditions. One potential solution to this challenge is to use a transmission with multiple gear ratios and shift to the appropriate ratio during flight. This solution will require a large transmission that is likely to be maintenance intensive and will require a complex shifting procedure to maintain power to the rotors at all times. An alternative solution is to use a fixed gear ratio transmission and require the power turbine to operate efficiently over the entire speed range. This concept is referred to as a variable-speed power-turbine (VSPT) and is the focus of the current study. This paper explores the design of a variable speed power turbine for civil tilt rotor applications using design optimization techniques applied to NASA's new meanline tool, the Object-Oriented Turbomachinery Analysis Code (OTAC).

design optimization↗

Effective Albedo of Vegetated Terrain at L-Band

This paper derives an explicit expression for an effective albedo of vegetated terrain from the zero- and multiple- order radiative transfer (RT) model comparison. The formulation establishes a direct physical link between the effective vegetation parameterization and the theoretical description of absorption and scattering within the canopy. The paper will present an evaluation of the derived albedo for corn canopies with data taken during an experiment at Alabama A&M Winfield A. Thomas Agricultural Research Station near Huntsville, Alabama in June, 1998. The test site consisted of two 50-m x 60-m plots - one with a bare surface and the other with grass cover - and four 30-m x 50-m plots of corn at different planting densities. One corn field was planted at a full density of 9.5 plants/sq m while the others were planted at 1/3, 1/2 and 2/3 of the full density. The fields were observed with a truck-mounted L-band radiometer at incident angle of 15 degree for the period of two weeks. Soil moisture (SM) changed daily due to irrigation and natural rainfall. Variations in gravimetric SM from 18 % to 34 % were seen during this period. Ground truth data, including careful characterization of the corn size and orientation statistics, and its dielectric, was also collected and used to simulate the effective albedo for the vegetation. The single-scattering albedo is defined as the fractional power scattered from individual vegetation constituents with respect to canopy extinction. It represents single-scattering properties of vegetation elements only, and is independent of ground properties. The values of the albedo get higher when there is dense vegetation (i.e. forest, mature corn, etc.) with scatterers, such as branches and trunks (or stalks in the case of corn), which are large with respect to the wavelength. This large albedo leads to a reduction in brightness temperature in the zero-order RT solution (known as tau-omega model). Higher-order multiple-scattering RT solutions are required for proper representation of scattering within vegetation. In this paper, an expression for an effective albedo for the whole canopy including the ground is derived for use in the zero-order RT model-based SM retrieval. This effective albedo takes into account of all the processes taking place within the canopy, including multiple-scattering. This new formulation will be presented and its importance for microwave SM retrieval will be evaluated for corn canopies in conjunction with the detailed ground truth data obtained during the experiment at Alabama in 1998. Emphasis will be placed on examining how the radiometer response to SM is modified by the corn canopy scattering under different field conditions. A semi-empirical parameterization of the effective albedo will be investigated through analysis of SM and vegetation water content effects on the effective albedo.

Kurum, Mehmet↗

Data Storage and Visualization Solutions for Real Time Simulations and Experiments

There are multiple challenges in performing real time simulations and real time experiments. Data storage and visualization are challenges that could be solved easily and when solved provide rich and valuable insights into the experiments. In this presentation, we will talk about the National Renewable Energy Laboratory solutions to tackle these issues while running real time experiments and simulations. The NREL team has made these open-source for the community to utilize with real time simulations and experiments.

data storage↗

A generalized distillation network synthesis model

In this report we propose a generalized superstructure-based distillation network synthesis model with improved modeling capabilities resulting in an extended solution space. The model can assign multiple mixtures to be separated (e.g., reactor network effluents) to different columns of the network while considering the interactions among different separation steps. In terms of outlets, products with general specifications, including pure components and multi-component mixtures, as well as streams without strict specifications (e.g., reactor recycle streams) can be readily handled. Furthermore, stream bypass is considered to avoid unnecessary separations, and thermal coupling is considered to yield novel solutions. The proposed model enables seamless integration with reactor network synthesis models thereby leading to solutions that are superior to the ones obtained by conventional sequential approaches.

42 ENGINEERING↗

Driving Force Analysis in an Infinite Anisotropic Plate with Multiple Crack Interactions

The methodology and a rigorous solution formulation are presented for stress intensity factors (SIF's, k) and total strain energy release rates (SERR, G(sub T)) of a multicracked plate, that has fully interacting cracks and is subjected to a far-field arbitrary stress state. The fundamental perturbation problem is derived, and the steps needed to formulate the system of singular integral equations whose solution gives rise to the evaluation of the SIF's are identified. Parametric studies are conducted for two, three, and four crack problems. The sensitivity and characteristics of the model is demonstrated.

Binienda, Wieslaw K.↗