Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional properties”

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 541 records · Page 30

Parametric Studies of Square Solar Sails Using Finite Element Analysis

Parametric studies are performed on two generic square solar sail designs to identify parameters of interest. The studies are performed on systems-level models of full-scale solar sails, and include geometric nonlinearity and inertia relief, and use a Newton-Raphson scheme to apply sail pre-tensioning and solar pressure. Computational strategies and difficulties encountered during the analyses are also addressed. The purpose of this paper is not to compare the benefits of one sail design over the other. Instead, the results of the parametric studies may be used to identify general response trends, and areas of potential nonlinear structural interactions for future studies. The effects of sail size, sail membrane pre-stress, sail membrane thickness, and boom stiffness on the sail membrane and boom deformations, boom loads, and vibration frequencies are studied. Over the range of parameters studied, the maximum sail deflection and boom deformations are a nonlinear function of the sail properties. In general, the vibration frequencies and modes are closely spaced. For some vibration mode shapes, local deformation patterns that dominate the response are identified. These localized patterns are attributed to the presence of negative stresses in the sail membrane that are artifacts of the assumption of ignoring the effects of wrinkling in the modeling process, and are not believed to be physically meaningful. Over the range of parameters studied, several regions of potential nonlinear modal interaction are identified.

Sleight, David W.↗

Dielectric Properties of Piezoelectric Polyimides

Molecular modeling and dielectric measurements are being used to identify mechanisms governing piezoelectric behavior in polyimides such as dipole orientation during poling, as well as degree of piezoelectricity achievable. Molecular modeling on polyimides containing pendant, polar nitrile (CN) groups has been completed to determine their remanent polarization. Experimental investigation of their dielectric properties evaluated as a function of temperature and frequency has substantiated numerical predictions. With this information in hand, we are then able to suggest changes in the molecular structures, which will then improve upon the piezoelectric response.

Ounaies, Z.↗

Damage Initiation and Ultimate Tensile Strength of Scaled [0 deg n/90 deg n/0 deg n]sub T Graphite-Epoxy Coupons

Previous research on scaling effects in composite materials has demonstrated that the stress levels at first ply failure and ultimate failure of composite laminates are dependent on the size of the laminate. In particular, the thickness dimension has been shown to be the most influential parameter in strength scaling of composite coupons loaded in tension. Geometrically and constitutively scaled laminates exhibit decreasing strength with increasing specimen size, and the magnitude of the strength-size effect is a function of both material properties and laminate stacking sequence. Some of the commonly used failure criteria for composite materials such as maximum stress, maximum strain, and tensor polynomial (e.g., Tsai-Wu) cannot account for the strength-size effect. In this paper, three concepts are developed and evaluated for incorporating size dependency into failure criteria for composite materials. An experimental program of limited scope was performed to determine the first ply failure stress in scaled cross-ply laminates loaded in tension. Test specimens were fabricated of AS-4/3502 graphite-epoxy composite material with laminate stacking sequences of [0 deg n/90 deg n/o deg n]subT where n=1-6. Two experimental techniques were used to determine first ply failure, defined as a transverse matrix crack in the 90 deg ply: (1) step loading with dye penetrant x-ray of the specimen at each load interval, and (2) acoustic emission. The best correlation between first ply failure analysis and experimental data was obtained using a modified Weibull approach which incorporated the residual thermal stress and the outer ply constraint, as well as the ply thickness effect. Finally, a second set of experiments was performed to determine the tensile response and ultimate failure of the scaled cross-ply laminates. The results of these experiments indicated no influence of specimen size on tensile response or ultimate strength.

Jackson, Karen E.↗

On the Motions of an Oscillating System Under the Influence of Flip-Flop Controls

So-called flip-flop controls (also called "on-off-course controls") are frequently preferred to continuous controls because of their simple construction. Thus they are used also for the steering control of airplanes. Such a body possesses-even if one thinks, for instance, only of the symmetric longitudinal motion - three degrees of freedom so that a study of its motions under the influence of an intermittent control is at least lengthy. Thus, it is suggested that an investigation of the basic effect of such a control first be made on a system with one degree of freedom. Furthermore, we limit ourselves in the resent report to the investigation of an "ideal" control where the control surface immediately obeys the command given by the "steering control function". Thus the oscillation properties of the control surface and the defects in linkage, sensing element, and mixing device are, at first, neglected. As long as the deviations from the "ideal" control may be neglected in practice, also the motion of the control surface takes place at the heat of the motion of the principal system. The aim of our investigation is to obtain a survey of the influence of the system and control coefficients on the damping behavior which is to be attained.

Fluegge-Lotz, I.↗

Finite-Difference Time Domain Solution of Light Scattering and Absorption by Particles in an Absorbing Medium

The three-dimensional (3-D) finite-difference time-domain (FDTD) technique has been extended to simulate light scattering and absorption by nonspherical particles embedded in an absorbing dielectric medium. A uniaxial perfectly matched layer (UPML) absorbing boundary condition is used to truncate the computational domain. When computing the single-scattering properties of a particle in an absorbing dielectric medium, we derive the single-scattering properties including scattering phase functions, extinction, and absorption efficiencies using a volume integration of the internal field. A Mie solution for light scattering and absorption by spherical particles in an absorbing medium is used to examine the accuracy of the 3-D UPML FDTD code. It is found that the errors in the extinction and absorption efficiencies from the 3-D UPML FDTD are less than similar to 2%. The errors in the scattering phase functions are typically less than similar to 5%. The errors in the asymmetry factors are less than similar to 0.l%. For light scattering by particles in free space, the accuracy of the 3-D UPML FDTD scheme is similar to a previous model.

Sun, W.↗

Subscale Test Methods for Combustion Devices

Stated goals for long-life LRE s have been between 100 and 500 cycles: 1) Inherent technical difficulty of accurately defining the transient and steady state thermochemical environments and structural response (strain); 2) Limited statistical basis on failure mechanisms and effects of design and operational variability; and 3) Very high test costs and budget-driven need to protect test hardware (aversion to test-to-failure). Ambitious goals will require development of new databases: a) Advanced materials, e.g., tailored composites with virtually unlimited property variations; b) Innovative functional designs to exploit full capabilities of advanced materials; and c) Different cycles/operations. Subscale testing is one way to address technical and budget challenges: 1) Prototype subscale combustors exposed to controlled simulated conditions; 2) Complementary to conventional laboratory specimen database development; 3) Instrumented with sensors to measure thermostructural response; and 4) Coupled with analysis

Anderson, W. E.↗

An Axisymmetric Hydrodynamical Model for the Torus Wind in AGN: X-ray Excited Funnel Flow - 2

We have calculated a series of models of outflows from the obscuring torus in active galactic nuclei (AGN). Our modeling assumes that the inner face of a rotationally supported torus is illuminated and heated by the intense X-rays from the inner accretion disk and black hole. As a result of such heating a strong biconical outflow is observed in our simulations. We calculate 3-dimensional hydrodynamical models, assuming axial symmetry, and including the effects of X-ray heating, ionization, and radiation pressure. We discuss the behavior of a large family of these models, their velocity fields, mass fluxes and temperature, as functions of the torus properties and X-ray flux. Synthetic warm absorber spectra are calculated, assuming pure absorption, for sample models at various inclination angles and observing times. We show that these models have mass fluxes and flow speeds which are comparable to those which have been inferred from observations of Seyfert 1 warm absorbers, and that they can produce rich absorption line spectra.

Dorodnitsyn, A.↗

Optimization of Neutral Atom Imagers

The interactions between plasma structures and neutral atom populations in interplanetary space can be effectively studied with energetic neutral atom imagers. For neutral atoms with energies less than 1 keV, the most efficient detection method that preserves direction and energy information is conversion to negative ions on surfaces. We have examined a variety of surface materials and conversion geometries in order to identify the factors that determine conversion efficiency. For chemically and physically stable surfaces smoothness is of primary importance while properties such as work function have no obvious correlation to conversion efficiency. For the noble metals, tungsten, silicon, and graphite with comparable smoothness, conversion efficiency varies by a factor of two to three. We have also examined the way in which surface conversion efficiency varies with the angle of incidence of the neutral atom and have found that the highest efficiencies are obtained at angles of incidence greater then 80deg. The conversion efficiency of silicon, tungsten and graphite were examined most closely and the energy dependent variation of conversion efficiency measured over a range of incident angles. We have also developed methods for micromachining silicon in order to reduce the volume to surface area over that of a single flat surface and have been able to reduce volume to surface area ratios by up to a factor of 60. With smooth micro-machined surfaces of the optimum geometry, conversion efficiencies can be increased by an order of magnitude over instruments like LENA on the IMAGE spacecraft without increase the instruments mass or volume.

Shappirio, M.↗

SDO/AIA Light Curves and Implications for Coronal Heating: Model Predictions

It seems largely agreed that many coronal loops---those observed at a temperature of about 1 MK---are bundles of unresolved strands that are heated by storms of impulsive nanoflares. The nature of coronal heating in hotter loops and in the very important but largely ignored diffuse component of active regions is much less clear. Is it also impulsive or is it quasi steady? The spectacular new data from the Atmospheric Imaging Assembly (AIA) telescopes on the Solar Dynamics Observatory (SDO) offer an excellent opportunity to address this question. We analyze the light curves of coronal loops and the diffuse corona in 6 different AIA channels and compare them with the predicted light curves from theoretical models. Light curves in the different AIA channels reach their peak intensities with predictable orderings as a function the nanoflare storm properties. We show that while some sets of light curves exhibit clear evidence of cooling after nanoflare storms, other cases are less straightforward to interpret. Complications arise because of line-of-sight integration through many different structures, the broadband nature of the AIA channels, and because physical properties can change substantially depending on the magnitude of the energy release. Nevertheless, the light curves exhibit predictable and understandable patterns. This presentation emphasizes the modeling aspects of our study. A companion presentation emphasizes the observations.

Klimchuk, James A.↗

Energetic Electrons in Dipolarization Events: Spatial Properties and Anisotropy

Using the electromagnetic fields of an MHD simulation of magnetotail reconnection, flow bursts, and dipolarization, we further investigate the acceleration of electrons to suprathermal energies. Particular emphasis is on spatial properties and anisotropies as functions of energy and time. The simulation results are compared with Time History of Events and Macroscale Interactions during Substorms observations. The test particle approach successfully reproduces several observed injection features and puts them into a context of spatial maps of the injection region(s): a dominance of perpendicular anisotropies farther down the tail and closer to the equatorial plane, an increasing importance of parallel anisotropy closer to Earth and at higher latitudes, a drop in energy fluxes at energies below approximately 10 keV, coinciding with the plasma density drop, together with increases at higher energy, a triple peak structure of flux increases near 0 deg, 90 deg, and 180 deg, and a tendency of flux increases to extend to higher energy closer to Earth and at lower latitudes. We identified the plasma sheet boundary layers and adjacent lobes as a main source region for both increased and decreased energetic electron fluxes, related to the different effects of adiabatic acceleration at high and low energies. The simulated anisotropies tend to exceed the observed ones, particularly for perpendicular fluxes at high energies. The most plausible reason is that the MHD simulation lacks the effects of anisotropy-driven microinstabilities and waves, which would reduce anisotropies.

anisotropy↗

Verification of Unstructured Grid Adaptation Components

Adaptive unstructured grid techniques have made limited impact on production analysis workflows where the control of discretization error is critical to obtaining reliable simulation results. Recent progress has matured a number of independent implementations of flow solvers, error estimation methods, and anisotropic grid adaptation mechanics. Known differences and previously unknown differences in grid adaptation components and their integrated processes are identified here for study. Unstructured grid adaptation tools are verified using analytic functions and the Code Comparison Principle. Three analytic functions with different smoothness properties are adapted to show the impact of smoothness on implementation differences. A scalar advection-diffusion problem with an analytic solution that models a boundary layer is adapted to test individual grid adaptation components. Laminar flow over a delta wing and turbulent flow over an ONERA M6 wing are verified with multiple, independent grid adaptation procedures to show consistent convergence to fine-grid forces and a moment. The scalar problems illustrate known differences in a grid adaptation component implementation and a previously unknown interaction between components. The wing adaptation cases in the current study document a clear improvement to existing grid adaptation procedures. The stage is set for the infusion of verified grid adaptation into production fluid flow simulations.

Park, Michael A.↗

Verification of Unstructured Grid Adaptation Components

Adaptive unstructured grid techniques have made limited impact on production analysis workflows where the control of discretization error is critical to obtaining reliable simulation results. Recent progress has matured a number of independent implementations of flow solvers, error estimation methods, and anisotropic grid adaptation mechanics. Known differences and previously unknown differences in grid adaptation components and their integrated processes are identified here for study. Unstructured grid adaptation tools are verified using analytic functions and the Code Comparison Principle. Three analytic functions with different smoothness properties are adapted to show the impact of smoothness on implementation differences. A scalar advection-diffusion problem with an analytic solution that models a boundary layer is adapted to test individual grid adaptation components. The scalar problems illustrate known differences in a grid adaptation component implementation and a previously unknown interaction between components. Laminar flow over a delta wing is verified with multiple, independent grid adaptation procedures to show consistent convergence to fine-grid forces and pitching moment.

Park, Michael A.↗

Residual Stresses Induced due to Curing of the Bulk Matrix in a Simplified Three-Dimensional (3D) Woven Repeating Unit Cell

It has been observed that 3D woven materials are highly susceptible to processes induced defects including tow misalignment, intratow and intertow cracking and/or voids. Adequate process modeling of 3D woven Polymer Matrix Composites (PMCs) is necessary to predict and estimate the effects of the manufacturing process on these defects. A preliminary step towards this is achieved with the commercial Finite Element Analysis (FEA) tool Abaqus complemented with written user subroutines to account for the effect of shrinkage and thermoelastic properties evolution as a function of the degree of cure of the bulk matrix phase in a 3D woven composite repeating unit cell (RUC). As a pathfinder, the process modeling framework is demonstrated on a simplified 3D woven RUC of an AS4/RIMR 135 system. Curing simulations are performed on the 3D woven RUC to observe the processing induced residual stresses. The results of the analysis show that high stress concentrations, present within the 3D intertow matrix bordering the binder tow, coincide with X-Ray computed tomography (X-Ray-CT) scan data depicting cracking in the same location

3D Woven↗

Implementation of Charged Particle Behavior in Discrete Element Method (DEM) Simulations

Lunar landers will agitate the surface of the Moon with an exhaust plume during descent which will, due to the particulate nature of the lunar regolith, loosen and eject grains from the surface. This ejection is not only coupled with the charged plume gas, but also results in significant particle-particle interactions. Settling of these grains after plume effects have subsided takes much longer than expected in a ballistic trajectory. The prevailing hypothesis attributes this behavior to the accumulated charge on the particles. We are thus developing a discrete element method (DEM) approach to explore these charged particle interactions on the lunar surface. The Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) Improved for General Granular and Granular Heat Transfer Simulations (LIGGGHTS) software package provides a DEM modeling framework for granular interactions. It includes many complexities such as non-spherical particle shapes, cohesion and frictional forces, and heat transfer, but has no provision for inter-particle electrostatic forces and charge transfer that are important to examine in the lunar environment. In this work, a standard Coulomb potential and a Yukawa potential are integrated into the LIGGGHTS framework to provide a basis for particle-particle electrostatic interactions, as well as a gravitational potential to enable inter-grain gravitational interactions. A preliminary approach to charge transfer between grains incorporating properties such as work function and electrical conductivity to the library of available material characteristics will be presented. Several scenarios have been simulated that include charged particle interactions within a diffuse granular gas, settling of charged grains into a regolith bed, sliding of granular material along an incline, and vibration of settled grains to produce a behavior similar to fluidization. There are numerous challenges to incorporate realistic interactions between complex lunar particles. Currently, grains are modeled to behave as if the entirety of the charge acts at the center of mass, such as conductors with spherical symmetry and insulators with homogeneously distributed charge. We are developing improvements that will include the use of non-spherical particle geometries, as well as reasonable approximations of insulating/dielectric materials that have non-uniform charge distributions. The cases simulated thus far will be examined in a relevant environment within a vacuum chamber to validate the simulations. These simulations will be bounded by experiments utilizing high-speed camera observations of the motion for validation. The grains in the experiment will exchange charge during their motion and this can be quantified by collection within a charge measurement device such as a Faraday cup. Such a device may be modeled within the software by defining an integration region and computing the contained charge as a function of simulation time, allowing for side-by-side comparison of simulated and measured bulk charging results. Any differences will be reconciled by updating the mathematical mechanisms described within the simulation suite. Successfully combining results from experiments within a relevant environment into the LIGGGHTS framework will improve modeling of the charged grain dynamics experienced on the Moon to provide insights into dust behavior for future lunar exploration missions.

Electrostatics↗

Comparison of Multisine Peak Factor Minimization Algorithms for Aircraft System Identification

Two phase-optimized multisine peak factor minimization algorithms are presented and evaluated. The first algorithm minimizes peak factor by iteratively clipping the peaks of generated multisine signals. The second algorithm optimizes peak factor indirectly through minimization of an approximation of the infinity norm of the multisine. Algorithm performance was evaluated as a function of different signal properties, including the number of harmonics, harmonic spacing, and number of snow harmonics (extra harmonics included for further reduction of the peak factor). The two algorithms are compared against results obtained by minimizing peak factor directly using a simplex algorithm, which has been a common approach when designing phase-optimized multisines for system identification flight tests. Sample results show that the clipping and infinity norm algorithms produced multisine signals with comparable peak factors that were lower than that of the simplex algorithm. However, the clipping algorithm runs an order of magnitude faster than the other two algorithms, which also makes it practical to repeat the algorithm multiple times to achieve even lower peak factors.

system identification↗

Comparison of Multisine Peak Factor Minimization Algorithms for Aircraft System Identification(Presentation)

Two phase-optimized multisine peak factor minimization algorithms are presented and evaluated. The first algorithm minimizes peak factor by iteratively clipping the peaks of generated multisine signals. The second algorithm optimizes peak factor indirectly through minimization of an approximation of the infinity norm of the multisine. Algorithm performance was evaluated as a function of different signal properties, including the number of harmonics, harmonic spacing, and number of snow harmonics (extra harmonics included for further reduction of the peak factor). The two algorithms are compared against results obtained by minimizing peak factor directly using a simplex algorithm, which has been a common approach when designing phase-optimized multisines for system identification flight tests. Sample results show that the clipping and infinity norm algorithms produced multisine signals with comparable peak factors that were lower than that of the simplex algorithm. However, the clipping algorithm runs an order of magnitude faster than the other two algorithms, which also makes it practical to repeat the algorithm multiple times to achieve even lower peak factors.

flight test↗

Wicking of liquids in screens

An investigation was conducted to determine the magnitude of the wicking rates of liquids in various screens. Evaluation of the parameters characterizing the wicking process resulted in the development of an expression which defined the wicking velocity in terms of screen and system geometry, liquid properties, and gravitational effects. Experiment data obtained both in normal gravity and in weightlessness demonstrated that the model successfully predicted the functional relation of the liquid properties and the distance from the liquid source to the wicking velocity. Because the pore geometry in the screens was complex, several screen geometric parameters were lumped into a single constant which was determined experimentally for each screen.

Symons, E. P.↗

Adjoint Methods for Adjusting Three-Dimensional Atmosphere and Surface Properties to Fit Multi-Angle/Multi-Pixel Polarimetric Measurements

This paper derives an efficient procedure for using the three-dimensional (3D) vector radiative transfer equation (VRTE) to adjust atmosphere and surface properties and improve their fit with multi-angle/multi-pixel radiometric and polarimetric measurements of scattered sunlight. The proposed adjoint method uses the 3D VRTE to compute the measurement misfit function and the adjoint 3D VRTE to compute its gradient with respect to all unknown parameters. In the remote sensing problems of interest, the scalar-valued misfit function quantifies agreement with data as a function of atmosphere and surface properties, and its gradient guides the search through this parameter space. Remote sensing of the atmosphere and surface in a three-dimensional region may require thousands of unknown parameters and millions of data points. Many approaches would require calls to the 3D VRTE solver in proportion to the number of unknown parameters or measurements. To avoid this issue of scale, we focus on computing the gradient of the misfit function as an alternative to the Jacobian of the measurement operator. The resulting adjoint method provides a way to adjust 3D atmosphere and surface properties with only two calls to the 3D VRTE solver for each spectral channel, regardless of the number of retrieval parameters, measurement view angles or pixels. This gives a procedure for adjusting atmosphere and surface parameters that will scale to the large problems of 3D remote sensing. For certain types of multi-angle/multi-pixel polarimetric measurements, this encourages the development of a new class of three-dimensional retrieval algorithms with more flexible parametrizations of spatial heterogeneity, less reliance on data screening procedures, and improved coverage in terms of the resolved physical processes in the Earth's atmosphere.

linearization↗