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 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Inverse aspect-ratio expanded tokamak equilibria

Following Greene et al. [Phys. Fluids 14, 671 (1971)] and Connor et al. [Phys. Plasmas 31, 577 (1988); Plasma Phys. Control. Fusion 34, 161 (1992); and Nucl. Fusion 33, 1533 (1993)], the Grad-Shafranov equation for an axisymmetric tokamak plasma equilibrium is solved via an expansion in the, supposedly small, inverse aspect-ratio of the plasma, ϵ. The displacements of equilibrium magnetic flux-surfaces due to plasma shaping are assumed to be $\mathcal{O}$(ϵ) smaller than the minor radii of the surfaces, but no other restriction is placed on the nature of the shaping. The solution of the Grad-Shafranov equation is matched to a vacuum solution that extends to infinity, and consists of an expansion in toroidal functions. The external poloidal magnetic field generated by a finite set of discrete external poloidal magnetic field-coils is calculated, and incorporated into the toroidal function expansion. In this manner, the shape of a large aspect-ratio tokamak plasma is directly related to the currents flowing in the external poloidal field-coils. Finally, a pedestal in the plasma pressure, and the associated spike in the bootstrap current, are incorporated into the model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

An analytic implementation of the IR-resummation for the BAO peak

We develop an analytic method for implementing the IR-resummation of [1], which allows one to correctly and consistently describe the imprint of baryon acoustic oscillations (BAO) on statistical observables in large-scale structure. We show that the final IR-resummed correlation function can be computed analytically without relying on numerical integration, thus allowing for an efficient and accurate use of these predictions on real data in cosmological parameter fitting. In this work we focus on the one-loop correlation function and the BAO peak. We show that, compared with the standard numerical integration method of IR-resummation, the new method is accurate to better than 0.2%, and is quite easily improvable. We also give an approximate resummation scheme which is based on using the linear displacements of a fixed fiducial cosmology, which when combined with the method described above, is about six times faster than the standard numerical integration. Finally, we show that this analytic method is generalizable to higher loop computations.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1986

Large Deformation Dynamic Bending of Composite Beams

Studies were conducted on the large deformation response of composite beams subjected to a dynamic axial load. The beams were loaded with a moderate eccentricity to promote bending. The study was primarily experimental but some finite element results were obtained. Both the deformation and the failure of the beams were of interest. The static response of the beams was also studied to determine potential differences between the static and dynamic failure. Twelve different laminate types were tested. The beams were loaded dynamically with a gravity driven impactor traveling at 19.6 ft/sec and quasi-static tests were conducted on identical beams in a displacement controlled manner. For laminates of practical interest, the failure modes under static and dynamic loadings were identical. Failure in most of the laminate types occurred in a single event involving 40% to 50% of the plies. However, failure in laminates with 30 deg or 15 deg off-axis plies occured in several events. All laminates exhibited bimodular elastic properties. Using empirically determined flexural properties, a finite element analysis was reasonably accurate in predicting the static and dynamic deformation response.

Derian, E. J.↗
NASA NTRS1994

Electromagnetically levitated vibration isolation system for the manufacturing process of silicon monocrystals

This paper introduces a study on an Electromagnetically Levitated Vibration Isolation System (ELVIS) for isolation control of large-scale vibration. This system features no mechanical contact between the isolation table and the installation floor, using a total of four electromagnetic actuators which generate magnetic levitation force in the vertical and horizontal directions. The configuration of the magnet for the vertical direction is designed to prevent any generation of restoring vibratory force in the horizontal direction. The isolation system is set so that vibration control effects due to small earthquakes can be regulated to below 5(gal) versus horizontal vibration levels of the installation floor of up t 25(gal), and those in the horizontal relative displacement of up to 30 (mm) between the floor and levitated isolation table. In particular, studies on the relative displacement between the installation floor and the levitated isolation table have been made for vibration control in the horizontal direction. In case of small-scale earthquakes (Taft wave scaled: max. 25 gal), the present system has been confirmed to achieve a vibration isolation to a level below 5 gal. The vibration transmission ratio of below 1/10 has been achieved versus continuous micro-vibration (approx. one gal) in the horizontal direction on the installation floor.

Kanemitsu, Yoichi↗
NASA NTRS2016

A Continuum Damage Mechanics Model to Predict Kink-Band Propagation Using Deformation Gradient Tensor Decomposition

A new model is proposed that represents the kinematics of kink-band formation and propagation within the framework of a mesoscale continuum damage mechanics (CDM) model. The model uses the recently proposed deformation gradient decomposition approach to represent a kink band as a displacement jump via a cohesive interface that is embedded in an elastic bulk material. The model is capable of representing the combination of matrix failure in the frame of a misaligned fiber and instability due to shear nonlinearity. In contrast to conventional linear or bilinear strain softening laws used in most mesoscale CDM models for longitudinal compression, the constitutive response of the proposed model includes features predicted by detailed micromechanical models. These features include: 1) the rotational kinematics of the kink band, 2) an instability when the peak load is reached, and 3) a nonzero plateau stress under large strains.

Bergan, Andrew C.↗
NASA NTRS1993

Oceanic transform earthquakes with unusual mechanisms or locations - Relation to fault geometry and state of stress in the adjacent lithosphere

Results are presented of a search for transform earthquakes departing from the pattern whereby they occur on the principal transform displacement zone (PTDZ) and have strike-slip mechanisms consistent with transform-parallel motion. The search was conducted on the basis of source mechanisms and locations taken from the Harvard centroid moment tensor catalog and the bulletin of the International Seismological Center. The source mechanisms and centroid depths of 10 such earthquakes on the St. Paul's, Marathon, Owen, Heezen, Tharp, Menard, and Rivera transforms are determined from inversions of long-period body waveforms. Much of the anomalous earthquake activity on oceanic transforms is associated with complexities in the geometry of the PTDZ or the presence of large structural features that may influence slip on the fault.

Wolfe, Cecily J.↗
NASA NTRS2012

Deployment Analysis of a Simple Tape-Spring Hinge Using Probabilistic Methods

Acceptance of new deployable structures architectures and concepts requires validated design methods to minimize the expense involved with technology validation flight testing. Deployable concepts for large lightweight spacecraft include booms, antennae, and masts. This paper explores the implementation of probabilistic methods in the design process for the deployment of a strain-energy mechanism, specifically a simple tape-spring hinge. Strain-energy mechanisms are attractive for deployment in very lightweight systems because they do not require the added mass and complexity associated with motors and controllers. However, designers are hesitant to include free deployment, strain-energy mechanisms because of the potential for uncontrolled behavior. In the example presented here, the tapespring cross-sectional dimensions have been varied and a target displacement during deployment has been selected as the design metric. Specifically, the tape-spring should reach the final position in the shortest time with the minimal amount of overshoot and oscillations. Surrogate models have been used to reduce computational expense. Parameter values to achieve the target response have been computed and used to demonstrate the approach. Based on these results, the application of probabilistic methods for design of a tape-spring hinge has shown promise as a means of designing strain-energy components for more complex space concepts.

Lyle, Karen H.↗
DOE OSTI2026

High-throughput computation of electric polarization in solids via Berry flux diagonalization

Electric polarization in the absence of an externally applied electric field is a key property of polar materials, but the standard interpolation-based ab initio approach to compute polarization differences within the modern theory of polarization presents challenges for automated high-throughput calculations. Berry flux diagonalization [J. Bonini et al., Phys. Rev. B 102, 045141 (2020)] has been proposed as an efficient and reliable alternative, though it has yet to be widely deployed. Here, we assess Berry flux diagonalization using ab initio calculations of a large set of materials, introducing and validating heuristics that ensure branch alignment with a minimal number of intermediate interpolated structures. Our automated implementation of Berry flux diagonalization succeeds in cases where prior interpolation-based workflows fail due to band-gap closures or branch ambiguities. Benchmarking with ab initio calculations of 176 candidate ferroelectrics, we demonstrate the efficacy of the approach on a broad range of insulating materials and obtain accurate effective polarization values with fewer interpolated structures than prior automated interpolation-based workflows. Our real-space heuristics that can predict gauge stability a priori from ionic displacements enable a general automated framework for reliable polarization calculations and efficient high-throughput screening of chemically and structurally diverse polar insulators. These results establish Berry flux diagonalization as a robust and efficient method to compute the effective polarization of solids and to accelerate the data-driven discovery of functional polar materials.

Poteshman, Abigail N. [University of Chicago, IL (↗
DOE OSTI2023

An adjoint-based optimization method for jointly inverting heterogeneous material properties and fault slip from earthquake surface deformation data

SUMMARY Analysis of tectonic and earthquake-cycle associated deformation of the crust can provide valuable insights into the underlying deformation processes including fault slip. How those processes are expressed at the surface depends on the lateral and depth variations of rock properties. The effect of such variations is often tested by forward models based on a priori geological or geophysical information. Here, we first develop a novel technique based on an open-source finite-element computational framework to invert geodetic constraints directly for heterogeneous media properties. We focus on the elastic, coseismic problem and seek to constrain variations in shear modulus and Poisson’s ratio, proxies for the effects of lithology and/or temperature and porous flow, respectively. The corresponding nonlinear inversion is implemented using adjoint-based optimization that efficiently reduces the cost function that includes the misfit between the calculated and observed displacements and a penalty term. We then extend our theoretical and numerical framework to simultaneously infer both heterogeneous Earth’s structure and fault slip from surface deformation. Based on a range of 2-D synthetic cases, we find that both model parameters can be satisfactorily estimated for the megathrust setting-inspired test problems considered. Within limits, this is the case even in the presence of noise and if the fault geometry is not perfectly known. Our method lays the foundation for a future reassessment of the information contained in increasingly data-rich settings, for example, geodetic GNSS constraints for large earthquakes such as the 2011 Tohoku-oki M9 event, or distributed deformation along plate boundaries as constrained from InSAR.

Geochemistry & Geophysics↗
NASA NTRS1989

Minimizing distortion and internal forces in truss structures by simulated annealing

Inaccuracies in the length of members and the diameters of joints of large truss reflector backup structures may produce unacceptable levels of surface distortion and member forces. However, if the member lengths and joint diameters can be measured accurately it is possible to configure the members and joints so that root-mean-square (rms) surface error and/or rms member forces is minimized. Following Greene and Haftka (1989) it is assumed that the force vector f is linearly proportional to the member length errors e(sub M) of dimension NMEMB (the number of members) and joint errors e(sub J) of dimension NJOINT (the number of joints), and that the best-fit displacement vector d is a linear function of f. Let NNODES denote the number of positions on the surface of the truss where error influences are measured. The solution of the problem is discussed. To classify, this problem was compared to a similar combinatorial optimization problem. In particular, when only the member length errors are considered, minimizing d(sup 2)(sub rms) is equivalent to the quadratic assignment problem. The quadratic assignment problem is a well known NP-complete problem in operations research literature. Hence minimizing d(sup 2)(sub rms) is is also an NP-complete problem. The focus of the research is the development of a simulated annealing algorithm to reduce d(sup 2)(sub rms). The plausibility of this technique is its recent success on a variety of NP-complete combinatorial optimization problems including the quadratic assignment problem. A physical analogy for simulated annealing is the way liquids freeze and crystallize. All computational experiments were done on a MicroVAX. The two interchange heuristic is very fast but produces widely varying results. The two and three interchange heuristic provides less variability in the final objective function values but runs much more slowly. Simulated annealing produced the best objective function values for every starting configuration and was faster than the two and three interchange heuristic.

Kincaid, Rex K.↗
NASA NTRS1990

Validation of parameterization scheme for eddy diffusion from satellite data

The eddy diffusion coefficient K(yy) has been calculated usign LIMS for the months of December 1978 and January and February 1979. Two methods have been used. The first implements the suggestion made by Tung (1987) to parameterize the eddy transport as a diffusive process along isentropes. The second method integrates the equation relating the parcel displacements to the eddy velocity fields. The latter method uses a filtering on both space and time domains to isolate transients and is referred to as the 'spectral method'. Results from the first method are shown to be reliable only for quiescent periods, breaking down when the meridional gradient of potential vorticity is negligible. Results from the two methods are in agreement only for very disturbed conditions, when transience is readily isolated. It is concluded that the parameterizations suggested for eddy transport and calculated in this paper may be meaningful for quiet periods, but are not reliable for unsteady and very large amplitude disturbances.

Sassi, F.↗
DOE OSTI2020

ELM control optimization for various ITER scenarios based on linear and quasi-linear figures of merit

For the purpose of better understanding type-I edge localized mode (ELM) control in ITER with resonant magnetic perturbation (RMP) fields, the plasma response to RMP is computed by a resistive full magneto-hydrodynamic model in toroidal geometry. Five scenarios designed for ITER are considered, ranging from the pre-nuclear to nuclear phases. The plasma response to RMP is quantified by the plasma surface displacement near the X-point of the divertor plasma and at the outboard mid-plane. The optimal coil configurations between two high- Q deuterium-tritium (DT) scenarios (at the same plasma current of 15 MA and the same magnetic field of 5.3 T but different fusion gains, Q = 5 and 10) are predicted to be similar. For the other ITER scenarios with similar edge safety factor q 95 ~ 3 to that of the baseline scenario, the optimal coil phasing is also similar. The optimization results are different for a half-current full-field (7.5 MA/5.3 T) scenario, largely due to the difference in q 95 . The RMP coil currents are also optimized to tailor the core vs edge toroidal torques exerted by the 3D RMP fields on the plasma column. Torque optimization, with various objective functions proposed in the study, is useful for minimizing the side effects of RMP on the plasma core flow in ITER, while still maintaining the ELM control capability. Full utilization of three rows of ELM control coils in ITER is found to be essential to ensure both flexibility and robustness of ELM control, in terms of both linear and quasilinear plasma responses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1987

Simultaneous VLA observations of a flare at 6 and 20 centimeter wavelengths

Using the Very Large Array at 6 and 20 cm wavelengths, a May 15, 1980 solar flare was mapped. The 1B flare, as observed in H alpha at Mees Solar Observatory, Maui, Hawaii, appeared as two sequential flares occurring at different neutral lines. The peaks of the time profile at 20 cm were delayed with respect to the 6 cm counterparts, but they were related to each other and to the H alpha activity. At 20 cm, precursor activity occurred, and an oppositely polarized source an arcminute away from the main burst may have triggered the flare. The main 20 cm emission appeared to be displaced limbward from the 6 cm burst. If both the 6 and 20 cm emission originated in the same system of loops, it is inferred that the apparent lateral displacement was caused by a height difference of 33,000 km in the sources of emission.

Velusamy, T.↗
NASA NTRS2004

Finite Element Analysis in Concurrent Processing: Computational Issues

The purpose of this research is to investigate the potential application of new methods for solving large-scale static structural problems on concurrent computers. It is well known that traditional single-processor computational speed will be limited by inherent physical limits. The only path to achieve higher computational speeds lies through concurrent processing. Traditional factorization solution methods for sparse matrices are ill suited for concurrent processing because the null entries get filled, leading to high communication and memory requirements. The research reported herein investigates alternatives to factorization that promise a greater potential to achieve high concurrent computing efficiency. Two methods, and their variants, based on direct energy minimization are studied: a) minimization of the strain energy using the displacement method formulation; b) constrained minimization of the complementary strain energy using the force method formulation. Initial results indicated that in the context of the direct energy minimization the displacement formulation experienced convergence and accuracy difficulties while the force formulation showed promising potential.

Sobieszczanski-Sobieski, Jaroslaw↗
DOE OSTI2026

Modeling the effect of MHD activity on runaway electron generation during SPARC disruptions

Magnetohydrodynamic (MHD) instabilities and runaway electrons (REs) interact in several ways, making it important to self-consistently model these interactions for accurate predictions of RE generation and the design of mitigation strategies, such as massive gas injection (MGI). Using M3D-C1 – an extended MHD code with a RE fluid model – we investigate the effects of 3-D nonlinear MHD activity, material injection, and 2-D axisymmetric vertical displacement events (VDEs) on RE evolution during disruptions on SPARC – a high-field, high-current tokamak designed to achieve a fusion gain Q > 1. Several cases, comprising different combinations of neon (Ne) and deuterium (D 2 ) injection, are considered. Here, our results demonstrate key effects that arise from the self-consistent RE + MHD coupling, such as an initial increase in RE generation due to MHD instability growth, decreased saturation energies of the m/n = 1/1 mode driving sawteeth-like activity, RE losses in stochastic magnetic fields, and subsequent RE confinement and plateau formation due to re-healing of flux surfaces. Large RE plateaus (>5 MA) are obtained with Ne-only injection (2-5 × 10 21 atoms), while combined D 2 + Ne injection (2 × 10 21 Ne atoms; 1.8 × 10 22 D 2 molecules) produces a lower RE current (<2 MA). With D 2 + Ne injection, a post thermal quench "cold" VDE terminates the RE beam, preventing a steady plateau. These simulations couple REs, 3-D MHD instabilities, MGI, and axisymmetric VDEs for the first time in SPARC disruption simulations and represent a crucial step in understanding RE generation and mitigation in high-current devices like SPARC.

Datta, Rishabh [Massachusetts Inst. of Technology ↗
DOE OSTI2022

Price Formation and Grid Operation Impacts from Variable Renewable Energy Resources

Increasing amounts of Variable Renewable Energy Resources (VREs) impact electricity markets and their operation. VREs are intermittent, zero marginal cost resources that tend to displace emissions-intensive generators in electricity dispatch, reducing emissions, but impacting price formation, revenue sufficiency, reliability, and market power mitigation processes of electricity markets. But VREs are not the only factor that affects operational and financial challenges in electricity markets. Declining natural gas prices, changing resource mixes, as well as different electricity market designs and regulatory policies all factor into the challenges both electricity market participants and operators face in today’s electricity markets. With this in-depth examination of electricity markets and related literature review, we aim to inform on key challenges of market design and operation for successful integration of large amounts of zero marginal cost resources. We’ve identified several areas, including VREs impact on price formation, revenue sufficiency, reliability, market power monitoring and mitigation, as well as how state-level incentives and market design impact VREs and these challenges. With each key challenge, we survey the literature to answer the question: To what extent is the problem, and how has it evolved over time? We first conduct a thorough review of ongoing challenges in electricity markets to understand the problem and review the empirical literature to capture important findings on how VREs, specifically, impact the problem. From this review, we highlight metrics that are important to understanding VRE integration and how market designs and outcomes are evolving with increasing levels of VREs. We propose several empirical models for future research to determine the impact of VREs on these identified challenges.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1985

Advanced ion thruster research

A series of experiments conducted on a ring cusp magnetic field ion thruster; in which the anode, cathode and discharge chamber backplate were moved relative to the magnetic cusp; are described. Optimum locations for the anode, cathode and backplate which yield the lowest energy cost per plasma ion and highest extracted ion fraction are identified. The results are discussed in terms of simple physical models. The results of preliminary experiments into the operation of hollow cathodes on nitrogen and xenon over a large pressure range (0.1 to 100 Torr) are presented. They show that the cathode discharge transfers from the cathode insert to the exterior edge of the orifice plate as the interelectrode pressure is increased. Experimental evidence showing that a new ion extractor grid concept can be used to stabilize the plasma sheath at the screen grid is presented. This concept, identified by the term constrained sheath optics, is shown to hold ion beamlet divergence and impingement characteristics to stable values as the beamlet current and the net and total accelerating voltages are changed. The current status of a study of beamlet vectoring induced by displacing the accelerator and/or decelerator grids of a three grid ion extraction system relative to the screen grid is discussed.

Wilbur, P. J.↗
NASA NTRS1994

A survey of the core-congruential formulation for geometrically nonlinear TL finite elements

This article presents a survey of the core-congruential formulation (CCF) for geometrically nonlinear mechanical finite elements based on the total Lagrangian (TL) kinematic description. Although the key ideas behind the CCF can be traced back to Rajasekaran and Murray in 1973, it has not subsequently received serious attention. The CCF is distinguished by a two-phase development of the finite element stiffness equations. The initial phase developed equations for individual particles. These equations are expressed in terms of displacement gradients as degrees of freedom. The second phase involves congruential-type transformations that eventually binds the element particles of an individual element in terms of its node-displacement degrees of freedom. Two versions of the CCF, labeled direct and generalized, are distinguished. The direct CCF (DCCF) is first described in general form and then applied to the derivation of geometrically nonlinear bar, and plane stress elements using the Green-Lagrange strain measure. The more complex generalized CCF (GCCF) is described and applied to the derivation of 2D and 3D Timoshenko beam elements. Several advantages of the CCF, notably the physically clean separation of material and geometric stiffnesses, and its independence with respect to the ultimate choice of shape functions and element degrees of freedom, are noted. Application examples involving very large motions solved with the 3D beam element display the range of applicability of this formulation, which transcends the kinematic limitations commonly attributed to the TL description.

Felippa, Carlos A.↗