SEARCH · Engineering Papers

Results for “cylindrical coordinates”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Electron Dynamics Within a MITL Containing a Load

In this paper, we derive the vacuum electric fields within specific cylindrically symmetric magnetically insulated transmission lines (MITLs) in the limit of an infinite speed of light for an arbitrary time-dependent current. We focus our attention on two types of MITLs: the radial MITL and a spherically curved MITL. We then simulate the motion of charged particles, such as electrons, present in these MITLs due to the vacuum fields. In general, the motion of charged particles due to the vacuum fields is highly nonlinear since the fields are nonlinear functions of spatial coordinates and depend on an arbitrary time-dependent current drive. Using guiding center theory, however, one can describe the gross particle kinetics using a combination of $\textbf {E} \times \textbf {B}$ and $\nabla B$ drifts. In addition, we compare our approximate inner MITL field models and particle kinetics with those from a fully electromagnetic simulation code. We find that the agreement between the approximate model and the electromagnetic simulations is excellent.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1992

Structural optimization of thin shells using finite element method

The objective of the present work was the structural optimization of thin shell structures that are subjected to stress and displacement constraints. In order to accomplish this, the structural optimization computer program DESAP1 was modified and improved. In the static analysis part of the DESAP1 computer program the torsional spring elements, which are used to analyze thin, shallow shell structures, were eliminated by modifying the membrane stiffness matrix of the triangular elements in the local coordinate system and adding a fictitious rotational stiffness matrix. This simplified the DESAP1 program input, improved the accuracy of the analysis, and saved computation time. In the optimization part of the DESAP1 program the stress ratio formula, which redesigns the thickness of each finite element of the structure, was solved by an analytical method. This scheme replaced the iterative solution that was previously used in the DESAP1 program, thus increasing the accuracy and speed of the design. The modified program was used to design a thin, cylindrical shell structure with optimum weight, and the results are reported in this paper.

Gotsis, Pascal K.↗

NASA NTRS1991

Large Deformation Behavior of Long Shallow Cylindrical Composite Panels

An exact solution is presented for the large deformation response of a simply supported orthotropic cylindrical panel subjected to a uniform line load along a cylinder generator. The cross section of the cylinder is circular and deformations up to the fully snapped through position are investigated. The orthotropic axes are parallel to the generator and circumferential directions. The governing equations are derived using laminated plate theory, nonlinear strain-displacement relations, and applying variational principles. The response is investigated for the case of a panel loaded exactly at midspan and for a panel with the load offset from midspan. The mathematical formulation is one dimensional in the circumferential coordinate. Solutions are obtained in closed-form. An experimental apparatus was designed to load the panels. Experimental results of displacement controlled tests performed on graphite-epoxy curved panels are compared with analytical predictions.

Carper, Douglas M.↗

NASA NTRS1985

Nonlinear theory for laminated and thick plates and shells including the effects of transverse shearing

Nonlinear strain displacement relations for three-dimensional elasticity are determined in orthogonal curvilinear coordinates. To develop a two-dimensional theory, the displacements are expressed by trigonometric series representation through-the-thickness. The nonlinear strain-displacement relations are expanded into series which contain all first and second degree terms. In the series for the displacements only the first few terms are retained. Insertion of the expansions into the three-dimensional virtual work expression leads to nonlinear equations of equilibrium for laminated and thick plates and shells that include the effects of transverse shearing. Equations of equilibrium and buckling equations are derived for flat plates and cylindrical shells. The shell equations reduce to conventional transverse shearing shell equations when the effects of the trigonometric terms are omitted and to classical shell equations when the trigonometric terms are omitted and the shell is assumed to be thin.

Stein, M.↗

NASA NTRS1993

A numerical study of confined turbulent jets

A numerical investigation is reported of turbulent incompressible jets confined in two ducts, one cylindrical and the other conical with a 5 degree divergence. In each case, three Craya-Curtet numbers are considered which correspond, respectively, to flow situations with no moderate and strong recirculation. Turbulence closure is achieved by using the k-epsilon model and a recently proposed realizable Reynolds stress algebraic equation model that relates the Reynolds stresses explicitly to the quadratic terms of the mean velocity gradients and ensures the positiveness of each component of the turbulent kinetic energy. Calculations are carried out with a finite-volume procedure using boundary-fitted curvilinear coordinates. A second-order accurate, bounded convection scheme and sufficiently fine grids are used to prevent the solutions from being contaminated by numerical diffusion. The calculated results are compared extensively with the available experimental data. It is shown that the numerical methods presented are capable of capturing the essential flow features observed in the experiments and that the realizable Reynolds stress algebraic equation model performs much better than the k-epsilon model for this class of flows of great practical importance.

Zhu, J.↗

NASA NTRS1994

Behavior of a Polydisperse Cluster of Drops Evaporating in a Vortex

A model has been developed to describe the dynamics and evaporation of polydisperse collections of liquid drops in an axisymmetric, infinite, cylindrical vortex. This formulation is valid both in the dense regime, where interactions between particles are important, and in the dilute regime, where interactions between particles are not important. In contrast to the standard way of discretizing polydisperse collections of particles, where the initial size distribution is partitioned into a finite number of bins and these bins are fixed during the calculation, here it is only the initial size distribution which is partitioned into a finite number of sizes. Each initial size class thus defined is followed dynamically and thermodynamically in its own system of coordinates which moves with the drops. Therefore, each initial size class develops a continuum of sizes as the drops are centrifuged towards the hotter ambient and evaporate. The gas phase is followed in its own system of coordinates.

liquid evaporation liquid dispersal particle flow↗

DOE OSTI2021

Diaboloidal mirrors: algebraic solution and surface shape approximations

A new type of optical element that can focus a cylindrical wave to a point focus (or vice versa) is analytically described. Such waves are, for example, produced in a beamline where light is collimated in one direction and then doubly focused by a single optic. A classical example in X-ray optics is the collimated two-crystal monochromator, with toroidal mirror refocusing. The element here replaces the toroid, and in such a system provides completely aberration free, point-to-point imaging of rays from the on-axis source point. We present an analytic solution for the mirror shape in its laboratory coordinate system with zero slope at the centre, and approximate solutions, based on bending an oblique circular cone and a bent right circular cylinder, that may facilitate fabrication and metrology.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

NASA NTRS1999

A Mixed Multi-Field Finite Element Formulation for Thermopiezoelectric Composite Shells

Analytical formulations are presented which account for the coupled mechanical, electrical, and thermal response of piezoelectric composite shell structures. A new mixed multi-field laminate theory is developed which combines "single layer" assumptions for the displacements along with layerwise fields for the electric potential and temperature. This laminate theory is formulated using curvilinear coordinates and is based on the principles of linear thermopiezoelectricity. The mechanics have the inherent capability to explicitly model both the active and sensory responses of piezoelectric composite shells in thermal environment. Finite element equations are derived and implemented for an eight-noded shell element. Numerical studies are conducted to investigate both the sensory and active responses of piezoelectric composite shell structures subjected to thermal loads. Results for a cantilevered plate with an attached piezoelectric layer are com- pared with corresponding results from a commercial finite element code and a previously developed program. Additional studies are conducted on a cylindrical shell with an attached piezoelectric layer to demonstrate capabilities to achieve thermal shape control on curved piezoelectric structures.

Lee, Ho-Jun↗

NASA NTRS1983

Electric field measurements at the magnetopause. I Observation of large convective velocities at rotational magnetopause discontinuities

Large convective electric fields of the order of 10 mV/m (sometimes as high as 22 mV/m) are observed at rotational magnetopause discontinuities. These observations were made with the long cylindrical (179-m base line) probes carried on the ISEE 1 satellite. These electric field observations yield convective velocity magnitudes (equal to the cross product of the vector E and the vector B, the latter divided by the square of the magnitude of B) of the order of 150 km/s. In this format for the convective velocity magnitudes, some of these observations are similar to the high speed plasma velocity observations that were made at the magnetopause with the plasma experiment carried on the ISEE 1 satellite. It is shown that, for many of these magnetopause crossings, there exists a special moving coordinate system where the observed electric fields vanish. Such a unique reference system is often used in theoretical studies of magnetic discontinuities. This special coordinate system does not move at the local plasma velocity but moves instead at a velocity intermediate between the convective velocity and the local Alfven velocity. It is used here as a diagnostic tool for the experimental investigation of rotational discontinuities at the magnetopause.

Aggson, T. L.↗

DOE OSTI2020

Analytical expressions of the surface shape of ‘diaboloid’ mirrors

Modern deterministic polishing processes allow fabrication of x-ray optics with almost any arbitrary aspherical surface shape. Among these optics, the so called "diaboloid"mirror is of special interest. The diaboloid mirror that converts a cylindrical wave to a spherical wave would improve focusing in x-ray beamlines implementing a diffraction element between a parabolic cylinder and a toroidal mirror. The replacement of the toroidal mirror in existing beamlines by the diaboloid mirror would mitigate aberrations. The shape of the diaboloid mirror is usually calculated numerically based on a truncated polynomial solution of the optical path problem. Here, we present an exact analytical solution for the shape of a diaboloid mirror as a function of the conjugate parameters of the mirror placed in a beamline. The derived analytical expressions for the diaboloid mirror in both the canonical and mirror-based coordinate systems are implemented in ray-tracing simulations to verify the beamline performances.

Yashchuk, Valeriy V.↗

NASA NTRS2000

Flow in the Proximity of the Pin-Tool in Friction Stir Welding and Its Relation to Weld Homogeneity

In the Friction Stir Welding (FSW) process a rotating pin inserted into a seam literally stirs the metal from each side of the seam together. It is proposed that the flow in the vicinity of the pin-tool comprises a primary rapid shear over a cylindrical envelope covering the pin-tool and a relatively slow secondary flow taking the form of a ring vortex about the tool circumference. This model is consistent with a plastic characterization of metal flow, where discontinuities in shear flow are allowed but not viscous effects. It is consistent with experiments employing several different kinds of tracer: atomic markers, shot, and wire. If a rotating disc with angular velocity w is superposed on a translating continuum with linear velocity omega, the trajectories of tracer points become circular arcs centered upon a point displaced laterally a distance v/omega from the center of rotation of the disc in the direction of the advancing side of the disc. In the present model a stream of metal approaching the tool (taken as the coordinate system of observation) is sheared at the slip surface, rapidly rotated around the tool, sheared again on the opposite side of the tool, and deposited in the wake of the tool. Local shearing rates are high, comparable to metal cutting in this model. The flow patterns in the vicinity of the pin-tool determine the level of homogenization and dispersal of contaminants that occurs in the FSW process. The approaching metal streams enfold one another as they are rotated around the tool. Neglecting mixing they return to the same lateral position in the wake of the tool preserving lateral tracer positions as if the metal had flowed past the tool like an extrusion instead of being rotated around it. (The seam is, however, obliterated.) The metal stream of thickness approximately that of the tool diameter D is wiped past the tool at elevated temperatures drawn out to a thickness of v/2(omega) in the wiping zone. Mixing distances in the wiping zone are multiplied in the unfolded metal. Inhomogeneities on a smaller scale than the mixing length are obliterated, but structure on a larger scale may be transmitted to the wake of a FSW weld.

Nunes, Arthur C., Jr.↗

NASA NTRS1994

Interlaminar stress analysis of dropped-ply laminated plates and shells by a mixed method

A mixed method of approximation based on Reissner's variational principle is developed for the linear analysis of interlaminar stresses in laminated composites, with special interest in laminates that contain terminated internal plies (dropped-ply laminates). Two models are derived, one for problems of generalized plane deformation and the other for the axisymmetric response of shells of revolution. A layerwise approach is taken in which the stress field is assumed with an explicit dependence on the thickness coordinate in each layer. The dependence of the stress field on the thickness coordinate is determined such that the three-dimensional equilibrium equations are satisfied by the approximation. The solution domain is reduced to one dimension by integration through the thickness. Continuity of tractions and displacements between layers is imposed. The governing two-point boundary value problem is composed of a system of both differential and algebraic equations (DAE's) and their associated boundary conditions. Careful evaluation of the system of DAE's was required to arrive at a form that allowed application of a one-step finite difference approximation. A two-stage Gauss implicit Runge-Kutta finite difference scheme was used for the solution because of its relatively high degree of accuracy. Patch tests of the two models revealed problems with solution accuracy for the axisymmetric model of a cylindrical shell loaded by internal pressure. Parametric studies of dropped-ply laminate characteristics and their influence on the interlaminar stresses were performed using the generalized plane deformation model. Eccentricity of the middle surface of the laminate through the ply drop-off was found to have a minimal effect on the interlaminar stresses under longitudinal compression, transverse tension, and in-plane shear. A second study found the stiffness change across the ply termination to have a much greater influence on the interlaminar stresses.

Harrison, Peter N.↗

NASA NTRS2001

Ceiling Fires Studied to Simulate Low-Gravity Fires

A unique new way to study low-gravity flames in normal gravity has been developed. To study flame structure and extinction characteristics in low-stretch environments, a normal gravity low-stretch diffusion flame was generated using a cylindrical PMMA sample of varying large radii, as shown in the photograph. These experiments have demonstrated that low-gravity flame characteristics can be generated in normal gravity through the proper use of scaling. On the basis of this work, it is feasible to apply this concept toward the development of an Earth-bound method of evaluating material flammability in various gravitational environments from normal gravity to microgravity, including the effects of partial gravity low-stretch rates such as those found on the Moon (1/6g) or Mars (1/3g). During these experiments, the surface regression rates for PMMA were measured for the first time over the full range of flammability in air, from blowoff at high stretch, to quenching at low stretch, as plotted in the graph. The solid line drawn through the central portion of the data (3<a<100 sec-1) has a slope of unity, which indicates regression is proportional to stretch. The figure coordinates assume that the values of stretch are equivalent, whether derived from forced stretch or from buoyant stretch. The excellent correlation of the regression-rate data over the two-order-of-magnitude variation of stretch shows the reasonableness of this assumption.

Olson, Sandra L.↗

NASA NTRS1993

System identification of the Large-Angle Magnetic Suspension Test Facility (LAMSTF)

The Large-Angle Magnetic Suspension Test Facility (LAMSTF), a laboratory-scale research project to demonstrate the magnetic suspension of objects over wide ranges of attitudes, has been developed. This system represents a scaled model of a planned Large-Gap Magnetic Suspension System (LGMSS). The LAMSTF system consists of a planar array of five copper electromagnets which actively suspend a small cylindrical permanent magnet. The cylinder is a rigid body and can be controlled to move in five independent degrees of freedom. Five position variables are sensed indirectly by using infra-red light-emitting diodes and light-receiving phototransistors. The motion of the suspended cylinder is in general nonlinear and hence only the linear, time-invariant perturbed motion about an equilibrium state is considered. One of the main challenges in this project is the control of the suspended element over a wide range of orientations. An accurate dynamic model plans an essential role in controller design. The analytical model of the LAMSTF system includes highly unstable real poles (about 10 Hz) and low-frequency flexible modes (about 0.16 Hz). Projection filters are proposed to identify the state space model from closed-loop test data in time domain. A canonical transformation matrix is also derived to transform the identified state space model into the physical coordinate. The LAMSTF system is stabilized by using a linear quadratic regulator (LQR) feedback controller. The rate information is obtained by calculating the back difference of the sensed position signals. The reference inputs contain five uncorrelated random signals. This control input and the system reponse are recorded as input/output data to identify the system directly from the projection filters. The sampling time is 4 ms and the model is fairly accurate in predicting the step responses for different controllers while the analytical model has a deficiency in the pitch axis.

Huang, Jen-Kuang↗

NASA NTRS2003

Strength Differential Measured in Inconel 718: Effects of Hydrostatic Pressure Studied

Aeropropulsion components, such as disks, blades, and shafts, are commonly subjected to multiaxial stress states at elevated temperatures. Experimental results from loadings as complex as those experienced in service are needed to help guide the development of accurate viscoplastic, multiaxial deformation models that can be used to improve the design of these components. During a recent study on multiaxial deformation (ref. 1) on a common aerospace material, Inconel 718, it was shown that the material in the aged state exhibits a strength differential effect (SDE), whereby the uniaxial compressive yield and subsequent flow behavior are significantly higher than those in uniaxial tension. Thus, this material cannot be described by a standard von Mises yield formulation. There have been other formulations postulated (ref. 2) that involve other combinations of the stress invariants, including the effect of hydrostatic stress. The question remained as to which invariants are necessary in the flow model. To capture the physical mechanisms occurring during deformation and reflect them in the plasticity formulation, researchers examined the flow of Inconel 718 under various amounts of hydrostatic stress to determine whether or not hydrostatic stress is needed in the formulation. Under NASA Grant NCC3-464, monitored by the NASA Glenn Research Center, a series of tensile tests were conducted at Case Western Reserve University on aged (precipitation hardened) Inconel 718 at 650 C and with superimposed hydrostatic pressure. Dogbone shaped tensile specimens (3-mm-diameter gauge by 16-mm gauge length) and cylindrical compression specimens (3-mm-diameter gauge by 6-mm gauge length) were strain gauged and loaded in a high-pressure testing apparatus. Hydrostatic pressures were obtained with argon and ranged from 210 to 630 MPa. The aged Inconel 718 showed a pronounced difference in the tension and compression yield strength (i.e., an SDE), as previously observed. Also, there were no significant effects of hydrostatic pressure on either the tensile and compressive yield strength (see the graph) or on the magnitude of the SDE. This behavior is not consistent with the pressure-dependent theory of the SDE, which postulates that the SDE is associated with pressure-dependent and/or internal friction dependent deformation associated with non-Schmid effects at the crystal level (refs. 3 and 4). Flow in Inconel 718 appears to be independent of hydrostatic pressure, suggesting that this invariant may be removed from the phenomenological constitutive model. As part of an ongoing effort to develop advanced constitutive models, Glenn s Life Prediction Branch coordinated this work with that of research on the multiaxial deformation behavior of Inconel 718 being conducted at Pennsylvania State University under NASA Grant NCC597.

Lewandowski, John J.↗

DOE OSTI2019

Chapter Six - Derivation, cross-validation, and comparison of analytic formulas for electrostatic deflector aberrations

Here, we derived first and second order analytic aberration formulas in the horizontal transverse plane for an electrostatic deflector specified by the reference orbit radius, the central angle spanning the deflector, and its inhomogeneity coefficients. The derivation was performed using an iterative order-by-order perturbation method in a Frenet–Serret beamline coordinate system. We produced a C program edabrt for calculation of the first and second order aberrations using the formulas that we derived. The electrostatic deflector aberration formulas from Wollnik (1965) disagree with COSY INFINITY and the aberration formulas that we derived here, and they also deviate from the first and second order symplecticity conditions. The reason for this discrepancy is that the aberration formulas from Wollnik (1965) consider only motion in the main field of the deflector and do not properly account for fringe field effects, in particular the necessarily occurring change of kinetic energy due to the change in potential. The code GIOS uses electrostatic aberration formulas from Wollnik (1965) and exhibits the same discrepancy. We performed a comparison of first and second order electrostatic deflector aberrations for (1) the analytic formulas derived in this work, (2) differential-algebraic (DA) numerical integration of the equations of motion using the code COSY INFINITY, (3) the aberration formulas from Wollnik (1965), (4) aberrations computed using the code GIOS, and (5) the aberration formulas from Wollnik (1965) adjusted to account for the occurring change in potential. An electrostatic spherical deflector and an electrostatic cylindrical deflector were used as test cases for this comparison. There is excellent agreement between methods (1), (2), and (5), and these three methods satisfy the first and second order symplecticity conditions.

43 PARTICLE ACCELERATORS↗

NASA NTRS1994

System identification of the Large-Angle Magnetic Suspension Test Fixture (LAMSTF)

The Large-Angle Magnetic Suspension Test Fixture (LAMSTF), a laboratory-scale research project to demonstrate the magnetic suspension of objects over wide ranges of attitudes, has been developed. This system represents a scaled model of a planned Large-Gap Magnetic Suspension System (LGMSS). The LAMSTF consists of a small cylindrical permanent magnet suspended element which is levitated above a planar array of five electromagnets mounted in a circular configuration. The cylinder is a rigid body and can be controlled to move in five independent degrees of freedom. Six position variables are sensed indirectly by using infrared light-emitting diodes and light-receiving phototransistors. The motion of the suspended cylinder is in general nonlinear and hence only the linear, time-invariant perturbed motion about an equilibrium state is considered. One of the main challenges in this project is the control of the suspended element over a wide range of orientations. An accurate dynamic model plays an essential role in controller design. The analytical model is first derived and open-loop characteristics discussed. The system is shown to be highly unstable and requires feedback control for system identification. Projection filters are first proposed to identify the state space model from closed-loop input/output test data in the time domain. This method is then extended to identify linear systems from the frequency test data. A canonical transformation matrix is also derived to transform the identified state space model into the physical coordinate. The LAMSTF system is stabilized by using a linear quadratic regulator (LQR) feedback controller for closed-loop identification. The rate information is obtained by calculating the back difference of the sensed position signals. Only the closed-loop random input/output data are recorded. Preliminary results from numerical simulations demonstrate that the identified system model is fairly accurate from either time domain or frequency-domain data. Experiments will be performed to validate the proposed closed-loop identification algorithms.

Huang, Jen-Kuang↗

DOE OSTI2022

Current Progress of the Final Design of a Subcritical Assembly at the Oak Ridge National Laboratory

During the past decade, the US Department of Energy (DOE)/National Nuclear Security Administration (NNSA) Nuclear Criticality Safety Program (NCSP) has been coordinating nuclear criticality safety (NCS) courses to provide training and qualification to new staff, supervisors, managers, regulators, and other professionals involved with NCS. To address the growing need for more NCS training possibilities, a subcritical assembly was designed at Oak Ridge National Laboratory (ORNL). The Oak Ridge Subcritical Assembly (ORSA), operated at ORNL, could be used in a new course for fissile material handlers/operators or as a backup training facility for currently existing NCSP hands-on NCS training. This paper describes some of the design steps necessary for obtaining a final ORSA design, beginning with a successful feasibility study performed in 2020. Numerous parameters within each main element of the assembly were considered, such as the reflector used, the type of fuel coating or plating used, and the number of fuel plates used. The design choices related to each parametric study described in this paper are detailed and justified at the end of each subsection. An overview of the current final design characteristics of ORSA and the experiments’ possibilities is also provided. ORSA will be built with fuel plates from legacy Aerojet General Nucleonics (AGN)-201M research reactors currently stored at the Y-12 National Security Complex (Y-12). The cylindrical fuel plates comprise UO 2 particles embedded in polyethylene and have a nominal enrichment of 19.5 ± 0.5 wt % 235 U. The basis for the final design was the feasibility study from 2020. In the feasibility study, ORSA achieved k eff = 0.95, corresponding to a multiplication factor of 20, with approximately 620 g of 235 U (four 4 cm plates, three 2 cm plates, and two 1 cm plates) and an 8.5 cm thick graphite reflector radially and on top of the assembly. The bottom of the assembly has a graphite thickness of 5 cm, and the chosen graphite density was 2.3 g/cm 3 . All design calculations were performed using SCALE 6.2.4 KENO-VI and with the ENDF/B VII.1 Continuous Energy cross-section library.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗