Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex loading”

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 451 records · Page 25

Virtual diffraction simulations using the quasi-coarse-grained dynamics method to understand and interpret plasticity contributions during in situ shock experiments

The dynamic deformation response of metallic materials has contributions from dislocations, deformation twinning, and plastic deformation. The current state-of-art techniques can detail the complex mechanistic history of deformation modes under shock loading in real-time using in situ X-ray diffraction (XRD). However, the capability of these experiments to unravel plasticity contributions is challenging due to limitations in interpreting results and the lack of validation from atomistic simulations. Molecular dynamics (MD) simulations can successfully capture various deformation modes in metals and complement experiments using simulated diffractograms at various stages of evolution. However, the difference in length and time scales of MD simulations and experiments is a substantial obstacle in completing and interpreting in situ diffractograms. Therefore, the existing modeling methods require various approximations to model defect evolution and interaction at the mesoscales. However, while using approximations to correlate the peak broadening behavior to the density of dislocations or the shifts/splitting due to the presence of twins, the interpretations of the plasticity contributions from diffractograms are non-trivial, especially when multiple modes of deformation may be operating. This viewpoint discusses combining a mesoscale modeling method called quasi-coarse-grained dynamics and virtual XRD to characterize the plasticity contributions in BCC metals from slip, twinning, and phase transformation behavior. Furthermore, the combined approach shows promise in bridging the mesoscale gap between the capabilities of atomic-scale simulations and in situ experiments to characterize the dynamic deformation of materials.

36 MATERIALS SCIENCE↗

In-situ observations of cyclic deformation in an extruded Mg-2Nd-1Y-0.1Zr-0.1Ca alloy

In this study, the evolution of deformation mechanisms during cyclic loading in an extruded, solution-treated Mg–2Nd–1Y–0.1Zr–0.1Ca alloy was investigated using a combination of in-situ loading, scanning electron microscopy (SEM), electron backscatter diffraction (EBSD), and focused ion beam (FIB) nanofabrication. The initial microstructure exhibited a random crystallographic texture with no preferred grain orientation. Flat, rectangular dog-bone specimens were subjected to load-controlled, fully reversed fatigue for 50 cycles, during which the same region was sequentially mapped to track microstructural changes. After 10 cycles of loading deformation twins were observed. During tensile reloading detwinning or narrowing of those twinned regions occurred. After 20 cycles, detwinning ceased and residual twins remained in the material. SEM imaging revealed numerous surface slip traces after cyclic loading. EBSD-assisted slip trace analysis identified the activation of prismatic and pyramidal < c+a> slip systems during low-cycle fatigue. Site-specific scanning transmission electron microscopy (STEM) further revealed that deformation was also accommodated by basal < a> slip and the dissociation of < c+a> dislocations. Center-of-symmetry (COS) analysis confirmed that the dissociation of < c+a> dislocations resulted in the formation of I₁ intrinsic stacking faults after cyclic loading. These findings provide new insights into the complex interplay of dislocation mechanisms governing fatigue deformation in rare-earth-containing Mg alloys.

Cyclic deformation↗

Confinement effects and acid strength in zeolites

Chemical reactivity and sorption in zeolites are coupled to confinement and—to a lesser extent—to the acid strength of Brønsted acid sites (BAS). In presence of water the zeolite Brønsted acid sites eventually convert into hydronium ions. The gradual transition from zeolite Brønsted acid sites to hydronium ions in zeolites of varying pore size is examined by ab initio molecular dynamics combined with enhanced sampling based on Well-Tempered Metadynamics and a recently developed set of collective variables. While at low water content (1–2 water/BAS) the acidic protons prefer to be shared between zeolites and water, higher water contents (n > 2) invariably lead to solvation of the protons within a localized water cluster adjacent to the BAS. At low water loadings the standard free energy of the formed complexes is dominated by enthalpy and is associated with the acid strength of the BAS and the space around the site. Conversely, the entropy increases linearly with the concentration of waters in the pores, favors proton solvation and is independent of the pore size/shape.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Interleaved dual-species arrays of single atoms using a passive optical element and one trapping laser

We demonstrate trapping of individual rubidium (Rb) and cesium (Cs) atoms in an interleaved array of bright tweezers and dark bottle-beam traps, using a microfabricated optical element illuminated by a single-laser beam and a 4f system with spatial filtering. Our approach exploits the opposite-sign dynamic polarizabilities of Rb and Cs, ensuring that each species is exclusively trapped in either bright or dark sites. The passive optical mask creates optimal trap depths for both species using three transmittance levels while minimizing the optical phase difference, implemented using a variable-thickness absorbing layer of amorphous germanium. This trapping architecture achieves atom loading rates close to 50% while reducing system complexity compared to conventional methods using active optoelectronic components and/or multiple-laser wavelengths.

Fang, Chengyu [Univ. of Wisconsin, Madison, WI (Un↗

Machine Learning Solutions for a Stable Grid Recovery

Grid operating security studies are typically employed to establish operating boundaries, ensuring secure and stable operation for a range of operation under NERC guidelines. However, if these boundaries are severely violated, existing system security margins will be largely unknown, as would be a secure incremental dispatch path to higher security margins while continuing to serve load. As an alternative to the use of complex optimizations over dynamic conditions, this work employs the use of machine learning to identify a sequence of secure state transitions which place the grid in a higher degree of operating security with greater static and dynamic stability margins. Several reinforcement learning solution methods were developed using deep learning neural networks, including Deep Q-learning, Mu-Zero, and the continuous algorithms Proximal Reinforcement Learning, and Advantage Actor Critic Learning. The work is demonstrated on a power grid with three control dimensions but can be scaled in size and dimensionality, which is the subject of ongoing research.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Structure parameters in rotating Couette-Poiseuille channel flow

It is well-known that a number of steady state problems in fluid mechanics involving systems of nonlinear partial differential equations can be reduced to the problem of solving a single operator equation of the form: v + lambda Av + lambda B(v) = 0, v is the summation of H, lambda is the summation of one-dimensional Euclid space, where H is an appropriate (real or complex) Hilbert space. Here lambda is a typical load parameter, e.g., the Reynolds number, A is a linear operator, and B is a quadratic operator generated by a bilinear form. In this setting many bifurcation and stability results for problems were obtained. A rotating Couette-Poiseuille channel flow was studied, and it showed that, in general, the superposition of a Poiseuille flow on a rotating Couette channel flow is destabilizing.

Knightly, George H.↗

Response of discretely stiffened structures

This paper presents a theoretical study on dynamic response of complex discretely stiffened and interconnected structures to random loads. Modal and transfer matrix techniques are utilized to analyze these multidimensional structures. Application of this formulation for response and noise generation studies in space station and fuselage-wing design is discussed. The response spectral densities are calculated for several typical discretely stiffened and interconnected structural systems. It is found that this approach allows for easy parametric evaluation and that by proper selection of structural parameters and loading conditions, the vibration response can be reduced.

Vaicaitis, R.↗

Fuzzy logic based robotic controller

Existing Proportional-Integral-Derivative (PID) robotic controllers rely on an inverse kinematic model to convert user-specified cartesian trajectory coordinates to joint variables. These joints experience friction, stiction, and gear backlash effects. Due to lack of proper linearization of these effects, modern control theory based on state space methods cannot provide adequate control for robotic systems. In the presence of loads, the dynamic behavior of robotic systems is complex and nonlinear, especially where mathematical modeling is evaluated for real-time operators. Fuzzy Logic Control is a fast emerging alternative to conventional control systems in situations where it may not be feasible to formulate an analytical model of the complex system. Fuzzy logic techniques track a user-defined trajectory without having the host computer to explicitly solve the nonlinear inverse kinematic equations. The goal is to provide a rule-based approach, which is closer to human reasoning. The approach used expresses end-point error, location of manipulator joints, and proximity to obstacles as fuzzy variables. The resulting decisions are based upon linguistic and non-numerical information. This paper presents a solution to the conventional robot controller which is independent of computationally intensive kinematic equations. Computer simulation results of this approach as obtained from software implementation are also discussed.

Attia, F.↗

Meshless Local Petrov-Galerkin (MLPG) Method with Orthogonal Polynomials for Euler-Bernoulli Beam Problems

In this paper, the feasibility of orthogonal polynomials in the meshless local Petrov Galerkin method (MLPG) method is studied. The orthogonal polynomials, Chebyshev and Legendre polynomials, are used in this MLPG method as trial functions. The test functions used were power functions with smooth derivatives at their ends. The performance of these methods is studied by applying these methods to Euler-Bernoulli beam problems. The MLPG-Galerkin and Legendre methods passed all the patch tests for simple beam problems. Next the formulations are tested on complex beam problems such as beams with partial loadings and continuous beam problems. Problems with load discontinuities and additional supports require special attention. Near discontinuities, judicious choice of number of nodes and nodal placements are needed to obtain accurate deflections, slopes, moments and shear forces. As polynomial functions are used, the large number of nodes can create a transformation matrix that is ill-conditioned, resulting in problems with the inversion of the matrix. The conditioning worsens as the number of nodes are increased beyond 20. Quadruple precision was needed for models to obtain accurate solutions. Even with quadruple precision the accuracy of the method suffers as the number of nodes is increased beyond 20. This appears to be a drawback of the MLPG-Chebyshev and MLPG-Legendre methods.

Raju, Ivatury S.↗

Toughening of Boron Carbide Composites by Hierarchical Microstructuring

Due to a unique combination of properties including high hardness, low density, chemical and thermal stability, semi-conductivity, and high neutron absorption, boron carbide (B4C) is a potential candidate for various applications involving extreme environment. However, B4C’s current application is limited because of its low fracture toughness. In this study, a hierarchical microstructure design with features including TiB2 grains and graphite platelets was used to toughen B4C by simultaneously utilizing multiple toughening mechanisms including crack deflection, bridging, and micro-crack toughening. Using field-assisted sintering technology (FAST), B4C composites with dense and hierarchical microstructure were fabricated. Previously, the fracture toughness of fabricated B4C composites was measured at micro-scale using micro-indentation to have up to 56% improvement. In this work, the B4C composites’ fracture toughness was characterized at macro-scale using four-point bending methods and compared with previous results obtained at micro-scale. Micromechanics modeling of fracture behaviors for B4C-TiB2 composites was also performed to evaluate the contributions from experimentally observed toughening mechanisms. From four-point bending tests, B4C composites reinforced with both TiB2 grains (~15 vol%) and graphite platelets (~8.7 vol%) exhibited the highest fracture toughness enhancement from 2.38 to 3.65 MPa∙m^(1/2). The measured values were lower than those obtained using micro-indentation but maintained the general trends. The discrepancy between the indentation and four-point bending test results originated from the complex deformation behaviors triggered by the high contact load during indentation tests. Through micromechanics modeling, introduced thermal residual stress due to thermal expansion mismatch between B4C and TiB2, and weak interphases at B4C-TiB2 boundaries were identified as the main causes for experimentally observed toughness enhancement. These results proved the effectiveness of hierarchical microstructure designs for B4C toughening and can provide reference for the future design of B4C composites with optimized microstructures for further fracture toughness enhancement.

composites↗

Design of Hypervelocity Flow Generator (HFG) and Its Diagnostics

Ground facilities for hypersonic research are a key element for successful development of aerodynamically proven hypervelocity vehicles. Design concepts and diagnostics of a hypervelocity flow generator (HFG) were made as a test platform for hypersonic/hypervelocity spacecraft models at the NASA Langley Research Center. The HFG is a hypersonic flow field generator using optically heated gas which is blown into an 80 cu m vacuum chamber. The vacuum chamber is kept at a stable vacuum pressure with a combination of three large vacuum pumps, while the HFG is in the test mode. The HFG provides a relatively small test section with approximately a 20 cm window. This facility was designed to generate 2.45 km/sec of flow speed, and potentially generate a continuous flow with the nozzle and vacuum system. The window in test section provides a direct view of the shock wave around a model in order to measure temperature, pressure, and density profile within the shock layer. One of the key test goals under this project is to understand why the emission spectra from the standing shockwave plasma predicted by Lora-Loran codes are significantly different from the measured emission spectra from the Flight Investigation Reentry Environment (FIRE II) Flight. The correct estimate of the thermal loading on the leading edge of hypersonic vehicles greatly affects the aerodynamic design, the material selection for the vehicle, and the cooling requirement and can be obtained by the precise modeling of emission spectra from shockwave plasma. However, the estimation of thermal loading is not an easy task due to complex non-equilibrium radiative process within high temperature shock layers that still falls into a category of cold plasma. Direct flight experiments are the most desirable, but not a cost-effective approach. Analysis by computational fluid dynamics (CFD) offers many test flexibilities. However, the CFD codes must be fully tested and validated with experimental data before the codes are effectively used for practical design. Large discrepancies between experiments and codes appear in hypersonic/hypervelocity flow regimes at high altitudes of 60 km ~ 90 km. This HFG facility offers some important parameters for CFD code validation, such as collision cross-sections, relaxation times, reaction rate coefficients and transportation coefficients. The HFG test facility is based on the ejection flow of high temperature gas heated over to 3500 K through a nozzle. The tungsten gas chamber of the HFG is heated up to a desired temperature by a 60 kW optical power beam source. This system consists of an optical power source, a thermal chamber, an expansion nozzle, a test section, and an 80 cu m vacuum tank. 60 kW optical input power is obtained from the 150 kW Vortek arc lamp system (by Vortek Industries, Vancouver, Canada). This optical beam is focused to heat the gas chamber within which a flow media is heated. The maximum achievable temperature of the flow medium reached approximately 3500 K or even higher but is limited by the melting point of the chamber material used. The exhaust velocity through the nozzle was determined by the stagnation temperature and the molecular weight of the working medium at the test section. To provide design parameters, a NASA Chemical Equilibrium with Application (CEA) computer program is used for the simulation of aerothermal data. This CEA program can calculate chemical equilibrium and properties of complex mixtures using shock tube parameters. For nitrogen gas at 2666 K stagnation temperature, the maximum achievable velocity at the test section is approximately 2.45 km/sec which is within the range of the thermal velocity of 8000 K shock layer. Based on the calculation through the CEA program, the design parameters of a HFG were determined and implemented for the test section that includes an expansion nozzle. The installed test section of the experimental facility can sustain a condition of a re-entry vehicle from the Space at an altitude of 60-90 km. After installation of the HFG, the system was fully tested and its operational parameters were measured. An 80 cu m vacuum chamber of HFG was set at 1 torr level to keep a stable downstream condition. The pumping time to reach the minimum vacuum pressure (~ 1 torr) at the test chamber from atmospheric pressure was approximately an hour using Kinney (MBV-14000/MB -1600/KT-300) and Stokes (Model 1772 and 412) pumps. Such a setup condition allowed a continued stable operation of the HFG experiment with flowrates through 1-, 3-, or 5-mm diameters of nozzle throat. The flow characteristics of the HFG for various operating conditions were performed using a focused Schlieren method. At 0.8 torr chamber pressure, a barrel shock was observed at the test section. With a sphere obstacle of 12 mm diameter in the flow, a bow shock (~ 2-mm thickness) was observed by a focused Schlieren visualization method.

Hypervelocity flow↗

Topological Machine Learning Methods for Power System Responses to Contingencies: Preprint

While deep learning tools, coupled with the emerging machinery of topological data analysis, are proven to deliver various performance gains in a broad range of applications, from image classification to biosurveillance to blockchain fraud detection, their utility in areas of high societal importance such as power system modeling and, particularly, resilience quantification in the energy sector yet remain untapped. To provide fast acting synthetic regulation and contingency reserve services to the grid while having minimal disruptions on customer quality of service, we propose a new topologybased system that depends on neural network architecture for impact metrics classification and prediction in power systems. This novel topology-based system allows one to evaluate the impact of three power system contingency types, namely, in conjunction with transmission lines, transformers, and transmission lines combined with transformers. We show that the proposed new neural network architecture equipped with local topological measures facilitates both more accurate classification of unserved load as well as the amount of unserved load. In addition, we are able to learn complex relationships between electrical properties and local topological measurements on the simulated response to contingencies for NREL-SIIP power system.

graph theory↗

A Deformed Geometry Synthesis Technique for Determining Stacking and Cryogenically Induced Preloads for the Space Launch System

The Space Launch System (SLS) stacking and Core Stage (CS) fueling induce significant preloads that contribute to the liftoff pad separation “twang”. To accurately capture this, an approach is required that can replicate the physics of all SLS physical stacking steps, CS cryogenic shrinkage, associated geometric nonlinearities, and the transient behavior and decay of the preloads with changing boundary conditions as the vehicle separates from the pad. The Deformed Geometry Synthesis (DGS) approach presented here satisfies the above requirements. DGS determines induced preloads by modeling components in their deformed geometry states and then enforcing compatibility by closing the resulting “deadbands”. DGS seamlessly integrates into the multibody modal synthesis framework and does not require the use of artificial external loads to enforce preloads or post-processing steps to remove their influence. Since DGS iterates to solve for the deformed state inclusive of geometric nonlinearities, running linearized parametrics to exercise different potential orientations of ball jointed struts that connect the CS to Boosters for cryo-shrinkage analyses is entirely avoided. Relative to the transient behavior and decay of stacking and cryo-induced preloads with SLS liftoff pad separation, this is an area of considerable interest to the SLS program. To capture this in the most accurate way possible, DGS algorithms are designed to work with Henkel-Mar nonlinear pad separation algorithms which operate on the separating longitudinal and lateral degrees of freedom (DoFs) between the vehicle and the pad. As the separating DoFs release, in whatever manner as dictated by the interface geometries, interface loads and interface flexibilities as well as the external loading on the vehicle, the subject preloads generate a complex twang/decay time-trace as dictated by the physics of the problem. This paper presents DGS numerical verification against the closed-form solution for Timoshenko’s 3 ball-jointed strut preload problem. This problem is then extended by the authors to the geometric nonlinear case where DGS is compared to the Newton-Raphson solution of the nonlinear equations. Next, DGS is utilized to solve the SLS stacking and cryogenic shrinkage coupled loads analyses. Finally, Henkel-Mar pad separation simulations are executed that isolate the impact of the induced preloads’ twang and decay characteristics.

Space Launch System↗

Shuttle Debris Impact Tool Assessment Using the Modern Design of Experiments

Computational tools have been developed to estimate thermal and mechanical reentry loads experienced by the Space Shuttle Orbiter as the result of cavities in the Thermal Protection System (TPS). Such cavities can be caused by impact from ice or insulating foam debris shed from the External Tank (ET) on liftoff. The reentry loads depend on cavity geometry and certain Shuttle state variables, among other factors. Certain simplifying assumptions have been made in the tool development about the cavity geometry variables. For example, the cavities are all modeled as shoeboxes , with rectangular cross-sections and planar walls. So an actual cavity is typically approximated with an idealized cavity described in terms of its length, width, and depth, as well as its entry angle, exit angle, and side angles (assumed to be the same for both sides). As part of a comprehensive assessment of the uncertainty in reentry loads estimated by the debris impact assessment tools, an effort has been initiated to quantify the component of the uncertainty that is due to imperfect geometry specifications for the debris impact cavities. The approach is to compute predicted loads for a set of geometry factor combinations sufficient to develop polynomial approximations to the complex, nonparametric underlying computational models. Such polynomial models are continuous and feature estimable, continuous derivatives, conditions that facilitate the propagation of independent variable errors. As an additional benefit, once the polynomial models have been developed, they require fewer computational resources to execute than the underlying finite element and computational fluid dynamics codes, and can generate reentry loads estimates in significantly less time. This provides a practical screening capability, in which a large number of debris impact cavities can be quickly classified either as harmless, or subject to additional analysis with the more comprehensive underlying computational tools. The polynomial models also provide useful insights into the sensitivity of reentry loads to various cavity geometry variables, and reveal complex interactions among those variables that indicate how the sensitivity of one variable depends on the level of one or more other variables. For example, the effect of cavity length on certain reentry loads depends on the depth of the cavity. Such interactions are clearly displayed in the polynomial response models.

DeLoach, R.↗

Load-Following Power Timeline Analyses for the International Space Station

Spacecraft are typically complex assemblies of interconnected systems and components that have highly time-varying thermal communications, and power requirements. It is essential that systems designers be able to assess the capability of the spacecraft to meet these requirements which should represent a realistic projection of demand for these resources once the vehicle is on-orbit. To accomplish the assessment from the power standpoint, a computer code called ECAPS has been developed at NASA Lewis Research Center that performs a load-driven analysis of a spacecraft power system given time-varying distributed loading and other mission data. This program is uniquely capable of synthesizing all of the changing spacecraft conditions into a single, seamless analysis for a complete mission. This paper presents example power load timelines with which numerous data are integrated to provide a realistic assessment of the load-following capabilities of the power system. Results of analyses show how well the power system can meet the time-varying power resource demand.

Fincannon, James↗

Comparison of RTA Vibratory Loading in the 40- by 80-Foot Wind Tunnel and 80- by 120-Foot Wind Tunnel at the National Full-Scale Aerodynamics Complex

This paper presents data taken in the 80- by 120-Foot Wind Tunnel in 1992 to data taken in the 40- by 80-Foot Wind Tunnel in 2016 at the NFAC. In both sets of data, the RTA was equipped with an identical set of S-76 helicopter blades. These data sets were compared to determine the repeatability of the RTA as a testing device, and to quantify differences in 40- by 80-Foot Wind Tunnel vs. 80- by 120-Foot Wind Tunnel results. The contents include an overview of the RTA and test facility, as well as a discussion of the methodology used. Vibratory load data is presented for a variety of wind and RPM conditions in various test section configurations. Data presented in this report quantifies the repeatability of the RTA.

NFAC↗

Tri-Level Linear Programming Model for Automatic Load Shedding Using Spectral Clustering

Traditional load shedding schemes can be inadequate in grids with high renewable penetration, leading to unstable events and unnecessary grid islanding. Although for both manual and automatic operating modes load shedding areas have been predefined by grid operators, they have remained fixed, and may be sub-optimal due to dynamic operating conditions. In this work, a distributed tri-level linear programming model for automatic load shedding to avoid system islanding is presented. Preventing islanding is preferred because it reduces the need for additional load shedding besides the disconnection of transmission lines between islands. This is crucial as maintaining the local generation-demand balance is necessary to preserve frequency stability. Furthermore, uneven distribution of generation resources among islands can lead to increased load shedding, causing economic and reliability challenges. This issue is further compounded in modern power systems heavily dependent on non-dispatchable resources like wind and solar. The upper-level model uses complex power flow measurements to determine the system areas to shed load depending on actual operating conditions using a spectral clustering approach. The mid-level model estimates the area system state, while the lower-level model determines the locations and load values to be shed. The solution is practical and promising for real-world applications.

Baquedano-Aguilar, Mario D.↗

Multidisciplinary Probabilistic Heat Transfer/Structural Analysis Code Developed: NESTEM

High-Speed Civil Transport (HSCT) engine combustor liners are subjected to complex thermal environments and have to endure these for thousands of hours with assured reliability. In the past, several deterministic analyses have been performed, including detailed heat transfer analyses to obtain thermal profiles and deterministic stress analyses to identify critical locations of high stresses. Actual rig tests also have been performed for segments by simulating these loading situations as closely as possible. However, it is well known that many uncertainties exist in loading (primarily thermal loads due to heat transfer), boundary conditions (end fixity unknowns), and material properties (moduli, thermal-expansion coefficients, and conductivities). The present in-house effort at the NASA Lewis Research Center is directed toward accounting for these in a formal way to assess the performance of liner components under complex and uncertain loading conditions as well as subject to other geometry- and material-related uncertainties.

Murthy, Pappu L. N.↗