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 1,297 records · Page 72

Enhancing In-Flight Structural Health Monitoring of Vertical Lift Vehicles Operating in an Urban Environment

In-situ airframe sensors have long been considered a potential solution for structural health monitoring that could change the design, certification, operation, and maintenance paradigms of flight vehicles. In this approach, large networks of sensors covering the entire, or most of, an airframe throughout its operational lifetime would support real-time decisions on airworthiness and obviate the need to overbuild components or perform multiple cycles of structural qualification testing and inspections. This concept would go beyond the current practice of placing select sensors in strategic locations or relying on such sensors only during airframe qualification flights and inspections. For vehicles in the emerging urban air mobility space, reducing weight associated with overbuilds and shortening down-time associated with inspections are critical for improving safety and affordability. In practice, the wide-scale use of in-situ sensors as primary assurance for structural health has not been demonstrated to be feasible or the best solution. A sensor integration testbed was developed as a platform to evaluate the potential of multiple sensor types to enable decision making on airworthiness. We report on the initial runs of this testbed with multiple sensors attached to a common test article. Metal foil strain gauges, eddy current, fiber optic, guided wave (acoustic emission and ultrasonic) and carbon nanotube roving sensors were affixed to the test article. Baseline as well as post damage initiation and fatiguing data are presented and discussed. Despite the relative simplicity of the test article, the interpretation of the as-captured test data was generally not conclusive or did not have wide enough coverage. This result emphasizes the challenges and current limitations both in testing and the practical broad application of embedded sensors as the determinative elements in critical decision making on wide-scale structural health.

Vertical Lift Vehicles↗

Two-epoch Orbit Estimation for Wide Binaries Resolved in Hipparcos and Gaia

The Hipparcos catalog and its Double and Multiple System Annex (DMSA) lists 4099 components with individual proper motions and coordinates on the epoch 1991.25. Many of these long-period binary stars are also present in Gaia Data Release 2 (DR2). Using the available relative positions and proper motions separated by 24.25 yr, the equations of relative orbital motion can be solved for the two epoch eccentric anomalies, orbital period, and eccentricity. This method employs elimination of the linear Thiele–Innes unknowns and nonlinear optimization of the remaining condition equations. The quality of these solutions is compromised by the insufficient condition and modest precision of the Hipparcos astrometric data, as revealed by Monte Carlo simulations with artificially perturbed data points. The presence of multiple systems and optical pairs can also perturb the results. Limited experiments with artificial data indicate that useful estimates can be obtained with a 25 yr epoch difference for wide binaries with orbital periods up to ~500 yr. The prospects of this method dramatically improve with the proposed next-generation space astrometry missions such as Gaia-NIR and Theia, especially when additional conditions are included from astrometric or spectroscopic measurements. An ancillary catalog of cross-identification and astrometric information for 1295 double-star pairs cross-matched in Gaia DR2 and Hipparcos is also published.

79 ASTRONOMY AND ASTROPHYSICS↗

Automated Unstructured Grid Adaptation on a Strut Fuel Injector at Hypervelocity Flow Conditions

Computational fluid dynamics (CFD) analysis is presented with the use of an automated unstructured grid adaptation tool on a strut fuel injector at hypervelocity flow conditions. The analysis was carried out with the VULCAN-CFD solver using Reynolds-averaged simulations (RAS). The hypervelocity flow conditions match the high Mach number flow of the experiments conducted as part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center (LaRC). The current work utilizes an automated grid adaptation tool recently implemented into VULCAN-CFD, and explores this tool’s ability to solve high-speed mixing problems. Simulation results obtained using the unstructured adaptive grid approach are compared to those on a user generated structured grid. These results are evaluated by analyzing how efficiently comparable fidelity results are obtained from both adapted and structured simulations. In addition, two adaptation strategies were used to explore the impact on the final solution. In the current work, the unstructured grid adaptation tool automatically generates unstructured grids and performs adaptation of the grid based on a Hessian error estimate of a specified flow field parameter. Multiple adaptations were executed using each run strategy with the one-dimensional values of the mixing efficiency used to determine grid convergence and for comparison with the structured grid simulation results. It was found that the unstructured adaptive grid simulations were able to produce results that matched closely with those on structured grids using far fewer grid cells, and thus, requiring far less computational time to reach the solution. It was also discovered that the adaptation run strategy influenced the total number of grid cells and the efficiency with which a final grid-adapted solution was reached. Overall, the investigation demonstrated that the automated unstructured grid adaptation tool implemented in VULCAN-CFD is capable of accurately and efficiently solving complex highspeed mixing problems using only a fraction of the grid cells required to obtain comparable results using a user-generated structured grid.

grid adaptation↗

Characterization of the Sulfur Saturated Versions of the High-Chromium Matrix Glasses

In this report, the Savannah River National Laboratory provides chemical analysis of a series of sulfur saturated melt versions of simulated nuclear waste glasses, and chemical analysis of the wash solutions resulting from the preparation of these glasses. The glasses were selected and fabricated by the Pacific Northwest National Laboratory as part of a broader study of the influence of glass composition on chemical durability, sulfur retention, and other properties. The resulting data will be used in the development of improved property/composition models to support operation of the Hanford Tank Waste Treatment and Immobilization Plant. Some degree of scatter among the Al₂O₃, B₂O₃, CaO, K₂O, Na₂O, and SiO₂ measurements was noted for the study glasses. The measured concentrations of B₂O₃, chlorine, and fluorine were below the targeted values for most of the study glasses, likely because of volatility during the multiple melting steps. The measured concentrations of Cr₂O₃ and K₂O were generally low relative to the targeted values. These components were identified in the wash solutions, indicating partitioning to the excess sulfate phase that was washed from the glasses prior to analyses. As expected, the measured concentrations of SO₃ in the glasses were higher than targeted due to the use of the sulfur saturation method in fabricating these glasses. Minor scatter among the triplicate measurements of some of the analytes of the wash solutions was noted. The wash solutions contained mainly sulfur and sodium, with moderate concentrations of chromium and potassium. Revision 1 of this report corrects the omission of Al₂O₃ values from Table A-4 in Appendix A.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Decarbonizing the grid: Utilizing demand-side flexibility for carbon emission reduction through locational marginal emissions in distribution networks

Decarbonization of the electric grid has become an important world-wide priority and is actively happening in many ways by introducing innovations and new technologies from the generation sectors to the demand sectors. In particular, one promising pathway toward such net-zero carbon emissions is to utilize the demand-side flexibility with the increasing number of flexible loads in distribution networks. In this paper, we explore a load shifting strategy with the emerging concept of location marginal emissions (LMEs) to reduce carbon emissions. LMEs measure the impact of carbon emissions including the locational aspect in more granular way and thus provide a novel mechanism for the system operator and load aggregators to design the LME-based load shifting strategy, which can efficiently guide consumers and thus adjust their consumption behaviors. Simulation case studies on the IEEE test networks are performed to validate the capability of the proposed load shifting method to reduce carbon emissions. We also compare this with other relevant strategies to discuss multiple scenarios and corresponding results. Finally, while each provides a different level of flexibility, all the explored strategies tested have led to solutions that have lower carbon emissions, indicating the great potential of demand-side flexibility in reducing carbon emissions for future distribution networks.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

An integrated numerical model for coupled poro-hydro-mechanics and fracture propagation using embedded meshes

Integrated models for fluid-driven fracture propagation and general multiphase flow in porous media are valuable to the study and engineering of several systems, including hydraulic fracturing, underground disposal of waste, and geohazard mitigation across such applications. This work extends the coupled model multiphase flow and poromechanical model of Ren et al. (2018) to admit fracture propagation (FP). The coupled XFEM-EDFM scheme utilizes a separate fracture mesh that is embedded on a static background mesh. The onset and dynamics of fracture propagation are governed by the equivalent stress intensity factor (SIF) criterion. A domain-integral method (J integral) is applied to compute this information. Additionally, an adaptive time-marching scheme is proposed to rapidly restrict and grow temporal resolution to match the underlying time-scales. The proposed model is verified with analytical solutions, and shows the capability to accurately and adaptively co-simulate fluid transport and deformation as well as the propagation of multiple fractures.

42 ENGINEERING↗

Can N95 Respirators Be Reused after Disinfection? How Many Times?

The coronavirus disease 2019 (COVID-19) pandemic has led to a major shortage of N95 respirators, which are essential for protecting healthcare professionals and the general public who may come into contact with the virus. Thus, it is essential to determine how we can reuse respirators and other personal protective equipment in these urgent times. We investigated multiple commonly used disinfection schemes on media with particle filtration efficiency of 95%. Heating was recently found to inactivate the virus in solution within 5 min at 70 °C and is among the most scalable, user-friendly methods for viral disinfection. We found that heat (≤85 °C) under various humidities (≤100% relative humidity, RH) was the most promising, nondestructive method for the preservation of filtration properties in meltblown fabrics as well as N95-grade respirators. At 85 °C, 30% RH, we were able to perform 50 cycles of heat treatment without significant changes in the filtration efficiency. At low humidity or dry conditions, temperatures up to 100 °C were not found to alter the filtration efficiency significantly within 20 cycles of treatment. Ultraviolet (UV) irradiation was a secondary choice, which was able to withstand 10 cycles of treatment and showed small degradation by 20 cycles. However, UV can potentially impact the material strength and subsequent sealing of respirators. Finally, treatments involving liquids and vapors require caution, as steam, alcohol, and household bleach all may lead to degradation of the filtration efficiency, leaving the user vulnerable to the viral aerosols.

60 APPLIED LIFE SCIENCES↗

High-Throughput Computing: Case Study of Medical Image Processing Applications

HPC is designed for large-scale simulations using monolithic codes of tightly coupled processes highly optimized to deliver decreased time to solution. Medical image processing is not a traditional field of HPC. Similar to AI applications, medical image processing parses large datasets, typically multiple times, to support a variety of studies for classification, diagnosis or monitoring purposes. The convergence of AI, HPC and Big Data encouraged more fields using image processing to transition to HPC. However, not all applications benefit from the same optimizations. In this paper we focus on high throughput medical image processing applications that analyze a huge dataset of small MRI images and that require HPC systems to decrease the time of parsing the entire dataset and not individual MRIs. We show in this research the performance of running SLANT, an image processing application for a whole brain segmentation, on large-scale systems and highlight performance limitations. We present optimizations prioritizing throughput that exhibit a 3.5x speed-up on the Summit Supercomputer that can be used as a baseline for building a high-throughput execution framework for other HPC systems.

Predescu, Maria↗

All-Polyester Multilayer Plastics (‘All-Polyester MLPs’): A Redesign for Inherently Recyclable Plastics

Multilayer plastics (MLPs) are widely used in modern packaging because they combine multiple functions such as oxygen and moisture barriers, mechanical strength, puncture resistance, and heat sealability into lightweight and cost-effective packaging solutions. These attributes are essential for food, beverage, pharmaceutical, and consumer goods packaging. However, conventional commercial MLPs typically consist of five to twelve layers made from chemically incompatible materials, including polyolefins, polyethylene terephthalate, nylon, ethylene-vinyl alcohol (EVOH), adhesives, and tie layers. This complexity makes MLPs extremely difficult to recycle. As a result, the vast majority of multilayer plastics (MLPs) are disposed of through landfilling or incineration. This disposal pathway perpetuates demand for virgin material production, thereby driving high industrial energy consumption, increasing greenhouse gas (GHG) emissions, and contributing to the long-term accumulation of plastic waste. The scope of work included: (1) Design and fabrication of all-polyester multilayer structures using commercially relevant processing methods; (2) Experimental validation of barrier, mechanical, and sealing performance; (3) Demonstration of both mechanical and chemical recycling pathways; and (4) Comprehensive techno-economic analysis (TEA) and life-cycle assessment (LCA) to quantify cost, energy, and environmental impacts.

36 MATERIALS SCIENCE↗

Metrics for Intercomparison of Remapping Algorithms (MIRA) protocol applied to Earth system models

Abstract. Strongly coupled nonlinear phenomena such as those described by Earth system models (ESMs) are composed of multiple component models with independent mesh topologies and scalable numerical solvers. A common operation in ESMs is to remap or interpolate component solution fields defined on their computational mesh to another mesh with a different combinatorial structure and decomposition, e.g., from the atmosphere to the ocean, during the temporal integration of the coupled system. Several remapping schemes are currently in use or available for ESMs. However, a unified approach to compare the properties of these different schemes has not been attempted previously. We present a rigorous methodology for the evaluation and intercomparison of remapping methods through an independently implemented suite of metrics that measure the ability of a method to adhere to constraints such as grid independence, monotonicity, global conservation, and local extrema or feature preservation. A comprehensive set of numerical evaluations is conducted based on a progression of scalar fields from idealized and smooth to more general climate data with strong discontinuities and strict bounds. We examine four remapping algorithms with distinct design approaches, namely ESMF Regrid (Hill et al., 2004), TempestRemap (Ullrich and Taylor, 2015), generalized moving least squares (GMLS) (Trask and Kuberry, 2020) with post-processing filters, and WLS-ENOR (Li et al., 2020). By repeated iterative application of the high-order remapping methods to the test fields, we verify the accuracy of each scheme in terms of their observed convergence order for smooth data and determine the bounded error propagation using challenging, realistic field data on both uniform and regionally refined mesh cases. In addition to retaining high-order accuracy under idealized conditions, the methods also demonstrate robust remapping performance when dealing with non-smooth data. There is a failure to maintain monotonicity in the traditional L2-minimization approaches used in ESMF and TempestRemap, in contrast to stable recovery through nonlinear filters used in both meshless GMLS and hybrid mesh-based WLS-ENOR schemes. Local feature preservation analysis indicates that high-order methods perform better than low-order dissipative schemes for all test cases. The behavior of these remappers remains consistent when applied on regionally refined meshes, indicating mesh-invariant implementations. The MIRA intercomparison protocol proposed in this paper and the detailed comparison of the four algorithms demonstrate that the new schemes, namely GMLS and WLS-ENOR, are competitive compared to standard conservative minimization methods requiring computation of mesh intersections. The work presented in this paper provides a foundation that can be extended to include complex field definitions, realistic mesh topologies, and spectral element discretizations, thereby allowing for a more complete analysis of production-ready remapping packages.

58 GEOSCIENCES↗

Numerical studies of three-dimensional breakdown in trailing vortex wakes

The development of dangerous vortex wakes trailing between aircraft for appreciable distance downstream, and posing hazards to following aircraft, is calculated using finite-element three-dimensional relaxation methods. Near-wake, middle-wake, and far-wake conditions are delineated. Vortex system finite lifetimes can be calculated with greater facility owing to the inclusion of a self-induction term in the solution dependent upon local curvature and vortex core radius. Looping and convergence phenomena are studied for single vortex pairs, twin pairs, and multiple pairs. Proportions of flap span and wing span are shown to affect the time to convergence of vortex pairs. The self-induction term renders the convergence time of vortices a sensitive function of wavelength. Wakes trailing between wide-body ('Jumbo') aircraft are also examined.

Hackett, J. E.↗

Role of the doubly stochastic Neyman type-A and Thomas counting distributions in photon detection

The sums of Meyman type-A and Thomas random variables are shown to retain their form under the constant multiplication parameter constraint. The conditions under which the two random variables converge in distribution to the fixed multiplicative Poisson, and to the Gaussian, are presented. It is shown that the latter result is important, in that it provides a solution to likelihood-ratio detection, estimation and discrimination problems in the presence of many kinds of signal and noise. Among the explicit applications presented are: (1) the photon-counter scintillation detection of nuclear particles when particle flux is low, (2) the photon-counting detection of weak optical signals in the presence of ionizing radiation, and (3) the design of a star-scanner spacecraft guidance system adequate for hostile space environments.

Teich, M. C.↗

Implementation of a block Lanczos algorithm for Eigenproblem solution of gyroscopic systems

The details of implementation of a general numerical procedure developed for the accurate and economical computation of natural frequencies and associated modes of any elastic structure rotating along an arbitrary axis are described. A block version of the Lanczos algorithm is derived for the solution that fully exploits associated matrix sparsity and employs only real numbers in all relevant computations. It is also capable of determining multiple roots and proves to be most efficient when compared to other, similar, exisiting techniques.

Gupta, Kajal K.↗

Implementation of a block Lanczos algorithm for eigenproblem solution of gyroscopic systems

This paper describes the details of implementation of a general numerical procedure developed for the accurate and economical computation of natural frequencies and associated modes of any elastic structure rotating along an arbitrary axis. A block version of the Lanczos algorithm is derived for the solution that fully exploits associated matrix sparsity and employs only real numbers in all relevant computations. It is also capable of determining multiple roots and proves to be most efficient when compared to other, similar, existing techniques.

Gupta, K. K.↗

Simulation of cooperating robot manipulators on a mobile platform

The dynamic equations of motion are presented for two or more cooperating manipulators on a freely moving mobile platform. The system of cooperating robot manipulators forms a closed kinematic chain where the force of interaction must be included in the formulation of robot and platform dynamics. The formulation includes the full dynamic interactions from arms to platform and arm tip to arm tip, and the possible translation and rotation of the platform. The equations of motion are shown to be identical in structure to the fixed-platform cooperative manipulator dynamics. The number of DOFs of the system is sufficiently large to make recursive dynamic calculation methods potentially more efficient than closed-form solutions. A complete simulation with two 6-DOF manipulators of a free-floating platform is presented along a with a multiple-arm controller to position the common load.

Murphy, Stephen H.↗

Coupled structural/thermal/electromagnetic analysis/tailoring of graded composite structures

Accomplishments are described for the fourth years effort of a 5-year program to develop a methodology for coupled structural/thermal/electromagnetic analysis/tailoring of graded component structures. These accomplishments include: (1) demonstration of coupled solution capability; (2) alternate CSTEM electromagnetic technology; (3) CSTEM acoustic capability; (4) CSTEM tailoring; (5) CSTEM composite micromechanics using ICAN; and (6) multiple layer elements in CSTEM.

Mcknight, R. L.↗

Adjoint Techniques for Topology Optimization of Structures Under Damage Conditions

The objective of this cooperative agreement was to seek computationally efficient ways to optimize aerospace structures subject to damage tolerance criteria. Optimization was to involve sizing as well as topology optimization. The work was done in collaboration with Steve Scotti, Chauncey Wu and Joanne Walsh at the NASA Langley Research Center. Computation of constraint sensitivity is normally the most time-consuming step of an optimization procedure. The cooperative work first focused on this issue and implemented the adjoint method of sensitivity computation (Haftka and Gurdal, 1992) in an optimization code (runstream) written in Engineering Analysis Language (EAL). The method was implemented both for bar and plate elements including buckling sensitivity for the latter. Lumping of constraints was investigated as a means to reduce the computational cost. Adjoint sensitivity computation was developed and implemented for lumped stress and buckling constraints. Cost of the direct method and the adjoint method was compared for various structures with and without lumping. The results were reported in two papers (Akgun et al., 1998a and 1999). It is desirable to optimize topology of an aerospace structure subject to a large number of damage scenarios so that a damage tolerant structure is obtained. Including damage scenarios in the design procedure is critical in order to avoid large mass penalties at later stages (Haftka et al., 1983). A common method for topology optimization is that of compliance minimization (Bendsoe, 1995) which has not been used for damage tolerant design. In the present work, topology optimization is treated as a conventional problem aiming to minimize the weight subject to stress constraints. Multiple damage configurations (scenarios) are considered. Each configuration has its own structural stiffness matrix and, normally, requires factoring of the matrix and solution of the system of equations. Damage that is expected to be tolerated is local and represents a small change in the stiffness matrix compared to the baseline (undamaged) structure. The exact solution to a slightly modified set of equations can be obtained from the baseline solution economically without actually solving the modified system.. Shennan-Morrison-Woodbury (SMW) formulas are matrix update formulas that allow this (Akgun et al., 1998b). SMW formulas were therefore used here to compute adjoint displacements for sensitivity computation and structural displacements in damaged configurations.

Akgun, Mehmet A.↗

Optimization of Aerospace Structure Subject to Damage Tolerance Criteria

The objective of this cooperative agreement was to seek computationally efficient ways to optimize aerospace structures subject to damage tolerance criteria. Optimization was to involve sizing as well as topology optimization. The work was done in collaboration with Steve Scotti, Chauncey Wu and Joanne Walsh at the NASA Langley Research Center. Computation of constraint sensitivity is normally the most time-consuming step of an optimization procedure. The cooperative work first focused on this issue and implemented the adjoint method of sensitivity computation in an optimization code (runstream) written in Engineering Analysis Language (EAL). The method was implemented both for bar and plate elements including buckling sensitivity for the latter. Lumping of constraints was investigated as a means to reduce the computational cost. Adjoint sensitivity computation was developed and implemented for lumped stress and buckling constraints. Cost of the direct method and the adjoint method was compared for various structures with and without lumping. The results were reported in two papers. It is desirable to optimize topology of an aerospace structure subject to a large number of damage scenarios so that a damage tolerant structure is obtained. Including damage scenarios in the design procedure is critical in order to avoid large mass penalties at later stages. A common method for topology optimization is that of compliance minimization which has not been used for damage tolerant design. In the present work, topology optimization is treated as a conventional problem aiming to minimize the weight subject to stress constraints. Multiple damage configurations (scenarios) are considered. Each configuration has its own structural stiffness matrix and, normally, requires factoring of the matrix and solution of the system of equations. Damage that is expected to be tolerated is local and represents a small change in the stiffness matrix compared to the baseline (undamaged) structure. The exact solution to a slightly modified set of equations can be obtained from the baseline solution economically without actually solving the modified system. Sherrnan-Morrison-Woodbury (SMW) formulas are matrix update formulas that allow this. SMW formulas were therefore used here to compute adjoint displacements for sensitivity computation and structural displacements in damaged configurations.

Akgun, Mehmet A.↗