Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large displacement”

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 469 records · Page 26

Global, quasi-steady dynamics of the distant solar wind 2: Deformation of the heliospheric current sheet

A three-dimensional (3-D) magnetohydrodynamic (MHD) numerical model is used to trace global deformations of the heliospheric current sheet (HCS) caused by large-scale dynamical interactions associated with corotating solar wind flows. Configurations incorporating the tilted-dipole geometry are investigated out past 30 AU for a variety of dipole tilt angles, alpha. Inclusion of the full, 3-D interplanetary dynamics allows north-south displacements and the east-west warping of the HCS by advective corotational effects to be accurately assessed for the first time. It is found that large-scale spatial correlations between velocity and density imposed at the coronal source (i.e., the geometric arrangement whereby regions of high-speed low-density material lying adjacent to areas of slow, dense flow interact obliquely under the influence of solar rotation) result in a distinctive pattern of deformation of the HCS. For an alpha = 30 deg tilted-dipole example, it is shown that typical zonal variations in radial velocity lead to significant folding of the HCS within about 5 AU of the Sun. By 10 AU, additional sharp bends appear near the latitudinal extremes of the HCS surface, where it is overtaken by shock fronts driven by 3-D corotating interaction regions (CIRs). Moreover, the model suggests that inside about 20 AU, major plasma structures are systematically organized, about the HCS, such that the greatest concentrations of material and magnetic field (the centroids of the 3-D CIR structures) are coincident with the folded crests of the HCS (near heliographic latitudes lambda = +/- alpha). Thus, in these circumstances many of the more interesting dynamical features inconveniently lie well away from the heliographic equator. At larger heliocentric distances, where neighboring CIRs begin to interact strongly, the warping of the HCS abates dramatically and the association between folds in the HCS and major field and density concentrations is weakened and ultimately breaks down. On the basis of the model calculations, quasi-steady deformations and related phenomena should figure prominently in the interplanetary medium, except during those periods when the magnetic dipole and solar spin axes are nearly aligned (alpha is approximately less than 10 deg). Even then, they may be of significance if there are any substantial local warps or kinks in the streamer belt.

Pizzo, V. J.↗

Intermediate time sub-diffusion and stress relaxation in ring polymer melts

The slow dynamics of non-concatenated ring melts remains a frontier problem in polymer science with implications for many soft material environments including cellular biophysics. Here, in this work, we report large-scale simulations of model ring melts that analyze the monomer and center-of-mass (CM) mean square displacements (MSD) and stress relaxation function on intermediate time and length scales. The degree of dynamical slowing down is characterized by the maximally sub-diffusive fractional time scaling exponents. The data span an exceptionally wide range of ring degrees of polymerization and stiffnesses and are not successfully organized based on the classic measure linear chain entanglement, N/N e . Rather, we find that the crossover degree of polymerization, N D , based on ring macromolecular caging that successfully allows master curves to be constructed for the long-time CM self-diffusion constant also collapses these temporal dynamic scaling exponents. Different properties display different exponents and exhibit one or two regimes of linear variation with the logarithm of N D / N . A distinct crossover of the CM-MSD and stress relaxation exponents emerges at sufficiently large N or stiffness that is not found for the monomer MSD, indicating a novel form of dynamic decoupling. This crossover aligns with the predicted critical degree of polymerization for transitioning from a weak to strong caging regime, indicative of activated transport. The latter may reflect the emergence of an intermolecular collective contribution to stress in analogy with dense soft colloidal matter. Suggestions are made for future theoretical work to address the rich patterns of behavior discovered.

Anomalous diffusion↗

Optimizing trigger-level track reconstruction for sensitivity to exotic signatures

Many compelling beyond the Standard Model scenarios predict signals that result in unconventional charged particle trajectories. Signatures for which unusual tracks are the most conspicuous feature of the event pose significant challenges for experiments at the Large Hadron Collider (LHC), particularly for the trigger. This article presents a study of track-based triggers for a representative set of long-lived and unconventional signatures at the upcoming High Luminosity LHC, as well as resulting recommendations for the target parameters of a hardware-based tracking system. Scenarios studied include large multiplicities of low-p T tracks produced in a soft-unclustered-energy-pattern model, displaced leptons and anomalous prompt tracks predicted in a Supersymmetry model with long-lived staus, and displaced hadrons predicted in a Higgs portal scenario with long-lived scalars.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Interaction of SPI pellets with plasma on JET and associated disruptions

Abstract The presented data refer to the Shattered Pellet Injector (SPI) experiments carried out at JET in 2019–2020. This paper is a full journal version of the data originally presented as posters at TMPDM_2020 and EPS_2021. This paper presents various aspects of the interaction of pellets with plasma and associated disruptions. The experiment was performed with I p = (1.1–3.1) MA plasmas and mainly with Ne + D 2 pellet composition, but also with Ar pellets. The Current Quench (CQ) time, τ 80−20 , is the key characteristic of mitigation effectiveness. A pellet with a high content of Ne or Ar can reduce the CQ duration below the upper required JET threshold. Plasmas with high (thermal + internal poloidal magnetic) pre-disruptive plasma energy require a high content of Ne pellets to obtain a short CQ duration. Pellets with a small amount of Ne (and accordingly large amount of D), instead of causing a mitigated CQ, create the conditions for a ‘cold’ Vertical Displacement Events (VDE). The SPI was applied to plasma with different status: mainly to normal (‘healthy’) plasma, i.e. not prone to disruption, post-disruptive and VDE plasma. This study shows that SPI effectiveness in terms of CQ duration and, accordingly, EM loads does not depend on the state of the plasma, whether it is ‘healthy’ or post-disruptive plasma. SPI has been shown to reduce the axisymmetric vertical vessel reaction forces by about (30–40) % compared to unmitigated disruptions. On JET, the VDE, whether ‘hot’ or ‘cold’, always creates the conditions for a toroidal asymmetry in the plasma, so the VDE on the JET is referred to as Asymmetric VDE (AVDE). The interrupting of VDE and prevention of AVDE with SPI has been demonstrated. Thus, the effectiveness of disruption mitigation using SPI has been confirmed.

Gerasimov, S. N. (ORCID:0009000237937211)↗

Dynamic Deformation in Nuclear Graphite and Underlying Mechanisms

Time-dependent deformation in nuclear graphite is influenced by the creation and migration of radiation-induced defects in the reactor environment. This study investigates the role of pre-existing defects such as point defect clusters and Mrozowski cracks in nuclear graphite IG-110. Separate specimens were irradiated with a 2.8 MeV Au 2+ beam with a fluence of 4.38 × 10 14 cm −2 and an 8 MeV C 2+ beam with a fluence of 1.24 × 10 16 cm −2 . Microscopic specimens were either mechanically loaded inside a transmission electron microscope (TEM) or subjected to ex situ indentation-based creep loading. In situ TEM tests showed significant plasticity in regions highly localized around the Mrozowski cracks, resembling slip or ripplocation bands. Slip bands were also seen near regions without pre-existing defects but at very high stresses. Ex situ self-ion irradiation embrittled the specimens and decreased the creep displacement and rate, while heavy ion irradiation resulted in the opposite behavior. We hypothesize that the large-sized gold ions (compared to the carbon atoms) induced interplanar swelling as well as cross-plane channels for increased defect mobility. These findings illustrate the role of pre-existing defects in the dynamic relaxation of stresses during irradiation and the need for more studies into the radiation environment’s impact on the mechanical response of nuclear graphite.

creep deformation↗

Dynamic analysis of large structures by modal synthesis techniques.

Several criteria that may be used to evaluate the merits of some of the existing techniques for the dynamic analysis of large structures which involve division into substructures or components are examined. These techniques make use of component displacement modes to synthetize global systems of generalized coordinates and, for that reason, they have come to be known as modal synthesis or component mode methods. Two techniques have been found to be particularly useful - i.e., the modal synthesis method with fixed attachment modes, and the modal synthesis method with free attachment modes. These two methods are treated in detail, and general flow charts are presented for guidance in computer programming.

Hurty, W. C.↗

Computational fracture mechanics

Numerical procedures for accurate determination of elastic stress intensity factors for the general two-dimensional crack problem are reviewed. The elastic perfectly plastic state of crack tip deformation is studied by a finite element procedure. Elastic-plastic fields in the immediate vicinity of a crack tip are determined numerically by finite element procedures based on asymptotic studies of crack tip singularities in plastic materials. The small-scale yielding problem is modeled, and expressions for crack tip opening displacement, shear singularity amplitude, and plastic zone extent are derived. A finite element solution to the large-scale yielding of a circumferentially cracked round tension bar is obtained. The three-dimensional aspects of flawed structures and numerical methods of treating them are studied. Ductile fracture mechanisms, in particular crack tip fracture on the microscale, are discussed.

Rice, J. R.↗

Criterion for propagation of magnetohydrodynamic waves An equivalence

It is shown that the requirement that a magnetohydrodynamic wave propagate over a distance equal to a large number of wavelengths without substantial attenuation is equivalent to the requirement that the lateral velocity of displacement of the field lines far exceed a characteristic velocity. This latter velocity is the velocity of slippage of the field lines with respect to the medium.

De, B. R.↗

Multiple-mode nonlinear free and forced vibrations of beams using finite element method

Effects of large deflection geometric nonlinearity and multiple modes for free and forced vibrations of beams are investigated. Inplane displacement and inplane inertia are included in the formulation. The finite element method is employed. A harmonic force matrix is derived for forced vibration analysis. The relations of amplitude, frequency ratio and forcing intensity for beams of different boundary conditions and various load conditions are presented.

Mei, Chuh↗

A computational procedure for the dynamics of flexible beams within multibody systems

The dynamic analysis of three dimensional elastic beams which experience large rotational and large deformational motions are examined. The beam motion is modeled using an inertial reference for the translational displacements and a body-fixed reference for the rotational quantities. Finite strain rod theories are then defined in conjunction with the beam kinematic description which accounts for the effects of stretching, bending, torsion, and transverse shear deformations. A convected coordinate representation of the Cauchy stress tensor and a conjugate strain definition is introduced to model the beam deformation. To treat the beam dynamics, a two-stage modification of the central difference algorithm is presented to integrate the translational coordinates and the angular velocity vector. The angular orientation is then obtained from the application of an implicit integration algorithm to the Euler parameter/angular velocity kinematical relation. The combined developments of the objective internal force computation with the dynamic solution procedures result in the computational preservation of total energy for undamped systems. The present methodology is also extended to model the dynamics of deployment/retrieval of the flexible members. A moving spatial grid corresponding to the configuration of a deployed rigid beam is employed as a reference for the dynamic variables. A transient integration scheme which accurately accounts for the deforming spatial grid is derived from a space-time finite element discretization of a Hamiltonian variational statement. The computational results of this general deforming finite element beam formulation are compared to reported results for a planar inverse-spaghetti problem.

Downer, Janice Diane↗

On the structure of the Southern Oscillation

A realistic oceanic general circulation model is forced with winds observed over the tropical Pacific between 1967 and 1979. The structure of the simulated Southern Oscillation is strikingly different in the western and eastern sides of the basin, because the principal interannual zonal-wind fluctuations are confined to the west and are in the form of an equatorial jet. This causes thermocline displacements to have maxima off the equator in the west (where the curl of the wind is large) but on the equator in the east. Zonal phase propagation, both on and off the equator, is at different speeds in the west and east. The phase pattern is complex, and there is, on interannual time scale, no explicit evidence of individual equatorial waves. These results lead to a modification of the 'delayed oscillator' mechanism originally proposed by Schopf and Suarez to explain a continual Southern Oscillation. The results also permit an evaluation of the various coupled ocean-atmosphere models that simulate the Southern Oscillation and indicate which measurements are necessary to determine which models are most relevant to reality.

Chao, YI↗

Crustal Displacements Due to Continental Water Loading

The effects of long-wavelength (> 100 km), seasonal variability in continental water storage on vertical crustal motions are assessed. The modeled vertical displacements (delta-r(sub M)) have root-mean-square (RMS) values for 1994-1998 as large as 8 mm with ranges up to 30 mm, and are predominantly annual in character. Regional strains are on the order of 20 nanostrain for tilt and 5 nanostrain for horizontal deformation. We compare delta-r(sub M) with observed Global Positioning System (GPS) heights (delta-r(sub O)) (which include adjustments to remove estimated effects of atmospheric pressure and annual tidal and non-tidal ocean loading) for 147 globally distributed sites. When the delta-r(sub O) time series are adjusted by delta-r(sub M), their variances are reduced, on average, by an amount equal to the variance of the delta-r(sub M). Of the delta-r(sub O) time series exhibiting a strong annual signal, more than half are found to have an annual harmonic that is in phase and of comparable amplitude with the annual harmonic in the delta-r(sub M). The delta-r(sub M) time series exhibit long-period variations that could be mistaken for secular tectonic trends or post-glacial rebound when observed over a time span of a few years.

vanDam, T.↗

Testing and Analyses of Advanced Composite Tow-Steered Shells with Cutouts

The structural response of two composite tow-steered shells with small cutouts in end compression is assessed using analyses and experimental measurements. The cylindrical shells were manufactured without cutouts using an automated fiber placement system, where the shells’ fiber orientation angles vary continuously around the shell circumference from ±10 degrees on the axially stiff crown and keel, to ±45 degrees on the shear-stiff sides. The first shell with overlaps has laminate thickness variations on the crown and keel that result from application of all 24 tows during each pass of the fiber placement system. The second shell without overlaps uses the fiber placement system’s tow drop/add capability to achieve a more uniform shell wall thickness. These shells without cutouts were tested in axial compression and buckled elastically. An unreinforced cutout representative of a passenger door on a commercial aircraft fuselage is machined into one side of each shell. Results from nonlinear finite element analysis are compared with their corresponding load, displacement, and strain measurements in prebuckling through stable postbuckling. Results from tests and analyses of shells with large cutouts are also discussed. The very good agreement observed in this study increases confidence in applying tow-steered composites to operational vehicles.

Tow-Steered↗

Thermal and chemical convection in planetary mantles

Melting of the upper mantle and extraction of melt result in the formation of a less dense depleted mantle. This paper describes series of two-dimensional models that investigate the effects of chemical buoyancy induced by these density variations. A tracer particles method has been set up to follow as closely as possible the chemical state of the mantle and to model the chemical buoyant force at each grid point. Each series of models provides the evolution with time of magma production, crustal thickness, surface heat flux, and thermal and chemical state of the mantle. First, models that do not take into account the displacement of plates at the surface of Earth demonstrate that chemical buoyancy has an important effect on the geometry of convection. Then models include horizontal motion of plates 5000 km wide. Recycling of crust is taken into account. For a sufficiently high plate velocity which depends on the thermal Rayleigh number, the cell's size is strongly coupled with the plate's size. Plate motion forces chemically buoyant material to sink into the mantle. Then the positive chemical buoyancy yields upwelling as depleted mantle reaches the interface between the upper and the lower mantle. This process is very efficient in mixing the depleted and undepleted mantle at the scale of the grid spacing since these zones of upwelling disrupt the large convective flow. At low spreading rates, zones of upwelling develop quickly, melting occurs, and the model predicts intraplate volcanism by melting of subducted crust. At fast spreading rates, depleted mantle also favors the formation of these zones of upwelling, but they are not strong enough to yield partial melting. Their rapid displacement toward the ridge contributes to faster large-scale homogenization.

Dupeyrat, L.↗

Modeling and Simulation of Variable Mass, Flexible Structures

The advent of the new Ares I launch vehicle has highlighted the need for advanced dynamic analysis tools for variable mass, flexible structures. This system is composed of interconnected flexible stages or components undergoing rapid mass depletion through the consumption of solid or liquid propellant. In addition to large rigid body configuration changes, the system simultaneously experiences elastic deformations. In most applications, the elastic deformations are compatible with linear strain-displacement relationships and are typically modeled using the assumed modes technique. The deformation of the system is approximated through the linear combination of the products of spatial shape functions and generalized time coordinates. Spatial shape functions are traditionally composed of normal mode shapes of the system or even constraint modes and static deformations derived from finite element models of the system. Equations of motion for systems undergoing coupled large rigid body motion and elastic deformation have previously been derived through a number of techniques [1]. However, in these derivations, the mode shapes or spatial shape functions of the system components were considered constant. But with the Ares I vehicle, the structural characteristics of the system are changing with the mass of the system. Previous approaches to solving this problem involve periodic updates to the spatial shape functions or interpolation between shape functions based on system mass or elapsed mission time. These solutions often introduce misleading or even unstable numerical transients into the system. Plus, interpolation on a shape function is not intuitive. This paper presents an approach in which the shape functions are held constant and operate on the changing mass and stiffness matrices of the vehicle components. Each vehicle stage or component finite element model is broken into dry structure and propellant models. A library of propellant models is used to describe the distribution of mass in the fuel tank or Solid Rocket Booster (SRB) case for various propellant levels. Based on the mass consumed by the liquid engine or SRB, the appropriate propellant model is coupled with the dry structure model for the stage. Then using vehicle configuration data, the integrated vehicle model is assembled and operated on by the constant system shape functions. The system mode shapes and frequencies can then be computed from the resulting generalized mass and stiffness matrices for that mass configuration. The rigid body mass properties of the vehicle are derived from the integrated vehicle model. The coupling terms between the vehicle rigid body motion and elastic deformation are also updated from the constant system shape functions and the integrated vehicle model. This approach was first used to analyze variable mass spinning beams and then prototyped into a generic dynamics simulation engine. The resulting code was tested against Crew Launch Vehicle (CLV-)class problems worked in the TREETOPS simulation package and by Wilson [2]. The Ares I System Integration Laboratory (SIL) is currently being developed at the Marshall Space Flight Center (MSFC) to test vehicle avionics hardware and software in a hardware-in-the-loop (HWIL) environment and certify that the integrated system is prepared for flight. The Ares I SIL utilizes the Ares Real-Time Environment for Modeling, Integration, and Simulation (ARTEMIS) tool to simulate the launch vehicle and stimulate avionics hardware. Due to the presence of vehicle control system filters and the thrust oscillation suppression system, which are tuned to the structural characteristics of the vehicle, ARTEMIS must incorporate accurate structural models of the Ares I launch vehicle. The ARTEMIS core dynamics simulation models the highly coupled nature of the vehicle flexible body dynamics, propellant slosh, and vehicle nozzle inertia effects combined with mass and flexible body properties that vary significant with time during the flight. All forces that act on the vehicle during flight must be simulated, including deflected engine thrust force, spatially distributed aerodynamic forces, gravity, and reaction control jet thrust forces. These forces are used to excite an integrated flexible vehicle, slosh, and nozzle dynamics model for the vehicle stack that simulates large rigid body translations and rotations along with small elastic deformations. Highly effective matrix math operations on a distributed, threaded high-performance simulation node allow ARTEMIS to retain up to 30 modes of flex for real-time simulation. Stage elements that separate from the stack during flight are propagated as independent rigid six degrees of freedom (6DOF) bodies. This paper will present the formulation of the resulting equations of motion, solutions to example problems, and describe the resulting dynamics simulation engine within ARTEMIS.

Tobbe, Patrick A.↗

Efficiency Improvements to the Displacement Based Multilevel Structural Optimization Algorithm

Multilevel Structural Optimization (MSO) continues to be an area of research interest in engineering optimization. In the present project, the weight optimization of beams and trusses using Displacement based Multilevel Structural Optimization (DMSO), a member of the MSO set of methodologies, is investigated. In the DMSO approach, the optimization task is subdivided into a single system and multiple subsystems level optimizations. The system level optimization minimizes the load unbalance resulting from the use of displacement functions to approximate the structural displacements. The function coefficients are then the design variables. Alternately, the system level optimization can be solved using the displacements themselves as design variables, as was shown in previous research. Both approaches ensure that the calculated loads match the applied loads. In the subsystems level, the weight of the structure is minimized using the element dimensions as design variables. The approach is expected to be very efficient for large structures, since parallel computing can be utilized in the different levels of the problem. In this paper, the method is applied to a one-dimensional beam and a large three-dimensional truss. The beam was tested to study possible simplifications to the system level optimization. In previous research, polynomials were used to approximate the global nodal displacements. The number of coefficients of the polynomials equally matched the number of degrees of freedom of the problem. Here it was desired to see if it is possible to only match a subset of the degrees of freedom in the system level. This would lead to a simplification of the system level, with a resulting increase in overall efficiency. However, the methods tested for this type of system level simplification did not yield positive results. The large truss was utilized to test further improvements in the efficiency of DMSO. In previous work, parallel processing was applied to the subsystems level, where the derivative verification feature of the optimizer NPSOL had been utilized in the optimizations. This resulted in large runtimes. In this paper, the optimizations were repeated without using the derivative verification, and the results are compared to those from the previous work. Also, the optimizations were run on both, a network of SUN workstations using the MPICH implementation of the Message Passing Interface (MPI) and on the faster Beowulf cluster at ICASE, NASA Langley Research Center, using the LAM implementation of UP]. The results on both systems were consistent and showed that it is not necessary to verify the derivatives and that this gives a large increase in efficiency of the DMSO algorithm.

Plunkett, C. L.↗

Mesh Deformation Based on Fully Stressed Design: The Method and Two-Dimensional Examples

Mesh deformation in response to redefined boundary geometry is a frequently encountered task in shape optimization and analysis of fluid-structure interaction. We propose a simple and concise method for deforming meshes defined with three-node triangular or four-node tetrahedral elements. The mesh deformation method is suitable for large boundary movement. The approach requires two consecutive linear elastic finite-element analyses of an isotropic continuum using a prescribed displacement at the mesh boundaries. The first analysis is performed with homogeneous elastic property and the second with inhomogeneous elastic property. The fully stressed design is employed with a vanishing Poisson s ratio and a proposed form of equivalent strain (modified Tresca equivalent strain) to calculate, from the strain result of the first analysis, the element-specific Young s modulus for the second analysis. The theoretical aspect of the proposed method, its convenient numerical implementation using a typical linear elastic finite-element code in conjunction with very minor extra coding for data processing, and results for examples of large deformation of two-dimensional meshes are presented in this paper. KEY WORDS: Mesh deformation, shape optimization, fluid-structure interaction, fully stressed design, finite-element analysis, linear elasticity, strain failure, equivalent strain, Tresca failure criterion

Hsu, Su-Yuen↗

Hydrokinetic energy conversion using flow induced oscillations of single-cylinder with large passive turbulence control

Various types of flow-induced oscillations (FIOs) have been implemented in development of marine hydrokinetic (MHK) energy converters. With passive turbulence control (PTC), energy harvesting starts at a flow speed of about 0.5 m/s. However, there is worldwide MHK energy available in even slower currents. In the present study, the effect of damping on FIO and power extraction is investigated for a converter with large turbulence stimulation (PTC) consisting of straight strips with a height of 15% of the cylinder diameter and placed symmetrically on the cylinder surface. The oscillating amplitude decreases, as the damping ratio increases, with unchanged sinusoidal pattern of the displacement time-history. The frequency ratio is also affected by damping especially in the VIV initial branch and transition region between VIV and galloping. An important flow characteristic of the large-PTC cylinder is that a recirculation region is formed behind the PTC, causing appreciable disturbance to the flow past the cylinder. Power can be harvested in the whole FIO range and the harnessed power maximum appears at the largest inflow velocity tested. However, the optimum of harnessing efficiency is located at the beginning of the VIV upper branch. The gap between VIV and galloping is bridged when large PTC is used, eliminating the drop in power and efficiency even at higher damping, which would be a weakness of regular-PTC cylinder for energy harvesting. Finally, the mechanism behind the variation of harnessing efficiency with inflow velocity and damping ratio is revealed, and the optimality criterion for the converter design is discussed.

16 TIDAL AND WAVE POWER↗