Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuum modeling”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Motion of Drops on Surfaces with Wettability Gradients

A liquid drop present on a solid surface can move because of a gradient in wettability along the surface, as manifested by a gradient in the contact angle. The contact angle at a given point on the contact line between a solid and a liquid in a gaseous medium is the angle between the tangent planes to the liquid and the solid surfaces at that point and is measured within the liquid side, by convention. The motion of the drop occurs in the direction of increasing wettability. The cause of the motion is the net force exerted on the drop by the solid surface because of the variation of the contact angle around the periphery. This force causes acceleration of an initially stationary drop, and leads to its motion in the direction of decreasing contact angle. The nature of the motion is determined by the balance between the motivating force and the resisting hydrodynamic force from the solid surface and the surrounding gaseous medium. A wettability gradient can be chemically induced as shown by Chaudhury and Whitesides who provided unambiguous experimental evidence that drops can move in such gradients. The phenomenon can be important in heat transfer applications in low gravity, such as when condensation occurs on a surface. Daniel et al have demonstrated that the velocity of a drop on a surface due to a wettability gradient in the presence of condensation can be more than two orders of magnitude larger than that observed in the absence of condensation. In the present research program, we have begun to study the motion of a drop in a wettability gradient systematically using a model system. Our initial efforts will be restricted to a system in which no condensation occurs. The experiments are performed as follows. First, a rectangular strip of approximate dimensions 10 x 20 mm is cut out of a silicon wafer. The strip is cleaned thoroughly and its surface is exposed to the vapor from an alkylchlorosilane for a period lasting between one and two minutes inside a desiccator. This is done using an approximate line source of the vapor in the form of a string soaked in the alkylchlorosilane. Ordinarily, many fluids, including water, wet the surface of silicon quite well. This means that the contact angle is small. But the silanized surface resists wetting, with contact angles that are as large as 100 degs. Therefore, a gradient of wettability is formed on the silicon surface. The region near the string is highly hydrophobic, and the contact angle decreases gradually toward a small value at the hydrophilic end away from this region. The change in wettability occurs over a distance of several mm. The strip is placed on a platform within a Plexiglas cell. Drops of a suitable liquid are introduced on top of the strip near the hydrophobic end. An optical system attached to a video camera is trained on the drop so that images of the moving drop can be captured on videotape for subsequent analysis. We have performed preliminary experiments with water as well as ethylene glycol drops. Results from these experiments will be presented in the poster. Future plans include the refinement of the experimental system so as to permit images to be recorded from the side as well as the top, and the conduct of a systematic study in which the drop size is varied over a good range. Experiments will be conducted with different fluids so as to obtain the largest possible range of suitably defined Reynolds and Capillary numbers. Also, an effort will be initiated on theoretical modeling of this motion. The challenges in the development of the theoretical description lie in the proper analysis of the region in the vicinity of the contact line, as well as in the free boundary nature of the problem. It is known that continuum models assuming the no slip condition all the way to the contact line fail by predicting that the stress on the solid surface becomes singular as the contact line is approached. One approach for dealing with this issue has been to relax the no-slip boundary condition using the Navier model. Molecular dynamics simulations of the contact line region show that for a non-polar liquid on a solid surface, the no-slip boundary condition is in fact incorrect near the contact line. Furthermore, the same simulations also show that the usual relationship between stress and the rate of deformation breaks down in the vicinity of the contact line. In developing continuum theoretical models of the system, we shall accommodate this knowledge to the extent possible.

Subramanian, R. Shankar↗
NASA NTRS1996

Robust Stabilizing Compensators for Flexible Structures with Collocated Controls

For flexible structures with collocated rate and attitude sensors/actuators, we characterize compensator transfer functions which guarantee modal stability even when stiffness/inertia parameters are uncertain. While the compensators are finite-dimensional, the structure models are allowed to be infinite-dimensional (continuum models), with attendant complexity of the notion of stability; thus exponential stability is not possible and the best we can obtain is strong stability. Robustness is interpreted essentially as maintaining stability in the worst case. The conditions require that the compensator transfer functions be positive real and use is made of the Kalman-Yakubovic lemma to characterize them further. The concept of positive realness is shown to be equivalent to dissipativity in infinite dimensions. In particular we show that for a subclass of compensators it is possible to make the system strongly stable as well as dissipative in an appropriate energy norm.

Balakrishman, A. V.↗
NASA NTRS1984

Shear deformation plate continua of large double-layered space structures

The energy equivalence to construct equivalent continua for the actual lattice structure is described. An energy equivalent continuum is defined as that which has the same amount of strain and kinetic energies stored in it as the original lattice structure when both are subjected to the same loading conditions. The equivalent continuum is characterized by its strain and kinetic energies from which the constitutive relations and the equations of motion can be derived. A simple method to model large rigid jointed lattices as continuous media with couple stresses is presented. The transition from the discrete system to the continuous medium is achieved by expanding the displacements and the rotations of the nodal points in a Taylor series about a suitable chosen origin. Basic kinematic assumptions are introduced to insure that the assumptions used in deriving the governing equations of the modeled continuum are satisfied. The number of terms retained in the Taylor series expansion will depend upon the properties to be evaluated. This implies that the kind of continuum needed to model from the discrete lattice, before the actual properties are derived was to predetermined.

Hefzy, M. S.↗
NASA NTRS1983

State and model error estimation for elliptic systems: Applications to large antenna static shape determination

This paper outlines the application of various estimation approaches to the problem of static shape determination for large antenna systems. The problem consists of estimating the shape of an antenna surface from measurements of its static deflection. The estimation schemes are based on any one of the modeling options of a single PDE for early insight and understanding, coarse-resolution multiple-PDE models for parametric studies and fine-resolution piecewise-continuum models for detailed design. For any one of these three models, estimator design can be developed using an infinite-dimensional approach, where the necessary finite-element truncation and approximation is conducted after the analytical design has taken place, or it can be based on a finite-dimensional approach, where the model is truncated before the estimation problem is formulated. One of the main objectives of the paper is to develop both approaches while simultaneously investigating their differences and similarities. Simulation results of an application of the finite-dimensional approach to a large parabolic reflector are presented.

Rodriguez, G.↗
NASA NTRS1992

Highlights of the theoretical work reported at the workshop

A brief overview is given of the theoretical work done on gamma-ray bursts (GRBs) done in this workshop. Papers on energy source models, optical transients, X-ray and gamma-ray continuum models, and cyclotron line models are reviewed.

Lamb, D. Q.↗
NASA NTRS1997

Coupling of Sph and Finite Element Codes for Multi-Layer Orbital Debris Shield Design

Particle-based hydrodynamics models offer distinct advantages over Eulerian and Lagrangian hydrocodes in particular shock physics applications. Particle models are designed to avoid the mesh distortion and state variable diffusion problems which can hinder the effective use of Lagrangian and Eulerian codes respectively. However conventional particle-in-cell and smooth particle hydrodynamics methods employ particles which are actually moving interpolation points. A new particle-based modeling methodology, termed Hamiltonian particle hydrodynamics, was developed by Fahrenthold and Koo (1997) to provide an alternative, fully Lagrangian, energy-based approach to shock physics simulations. This alternative formulation avoids the tensile and boundary instabilities associated with standard smooth particle hydrodynamics formulations and the diffusive grid- to-particle mapping schemes characteristic of particle-in-cell methods. In the work described herein, the method of Fahrenthold and Koo has been extended, by coupling the aforementioned hydrodynamic particle model to a hexahedral finite element based description of the continuum dynamics. The resulting continuum model retains all of the features (including general contact-impact effects) of Hamiltonian particle hydrodynamics, while in addition accounting for tensile strength, plasticity, and damage effects important in the simulation of hypervelocity impact on orbital debris shielding. A three dimensional, vectorized, and autotasked implementation of the extended particle method described here has been coded for application to orbital debris shielding design. Source code for the pre-processor (PREP), analysis code (EXOS), post-processor (POST), and rezoner (ZONE), have been delivered separately, along with a User's Guide describing installation and application of the software.

Fahrenthold, Eric P.↗
NASA NTRSDate not supplied

Modeling Deformable Linear Objects for Autonomous Robotic Outfitting of Lunar Surface Systems

This paper presents structural models of deformable linear objects (DLOs). DLOs are a subclass of deformable objects that encompasses common outfitting elements such as cables and ropes. Models are validated through hardware experiments, and integration in a robotic autonomy architecture for space environments is discussed. A persistent human presence on the lunar surface is one of the next major milestones in space exploration. This requires the development of robust extraplanetary construction technologies including structures and materials modeling and robotic systems. Previous robotic construction technology development has primarily focused on structural assembly, with significantly less focus on robotically performed outfitting tasks to instantiate subsystems providing power, data, life support, etc. These tasks involve manipulation of highly flexible elements, which are difficult to model, such as cable harnesses, ropes, and hoses. Robotic manipulation of DLOs, especially cable harnesses, is an active area of research as cable harnesses are essential for providing power and data to space assets. DLO models that can be used for robot manipulator trajectory generation are necessary for autonomous operation of lunar infrastructure. There are many proposed methods for modeling DLOs, and they primarily fall into three types: 1) discrete model-based, 2) continuum model-based, and 3) Neural Network-based. These types each have pros and cons, and the tradeoff between model accuracy and computational speed informs which type should be used. An understanding of this trade-off is imperative for real-time control of autonomous systems. High computational requirements reduce the speed of the model, making real-time control difficult, while accuracy is critical to preventing collisions. Discrete models, such as a mass-spring multibody representation, require relatively few calculations, and accuracy is directly tied to the step size of the discretization. Continuum models, such as a B-spline representation or a Cosserat rod model (a mix of continuous and discrete), are more informed of the structural properties of the cable and are much more accurate than a rigid body mass-spring model, but at significant computational cost. A Neural Network approach can provide an online solution with very few computational steps, but properly generating training data can be difficult and validation for an in-space application is not trivial. This paper explores the trade-off between different modeling approaches and compares accuracy and computational speed/complexity of the three types mentioned above. Model accuracy is evaluated using a cable in a static configuration. True cable shape is obtained using a depth camera for RGB images and point-cloud segmentation. The purpose of this experiment is to evaluate the trade-offs of different approaches to the DLO modeling problem. Understanding the tradeoffs between different cable modeling techniques paves the way for developing robotic control and planning architectures necessary for real-time manipulation of DLOs for lunar infrastructure outfitting. Real-time control is required for robotic systems to be able to actively manipulate a cable in a harsh environment where model and sensor errors compound, and environmental conditions can cause significant disturbances. Cable routing must be performed in areas with high density of objects/obstacles: through truss structures, near solar panels or mirror arrays, next to bundles of electrical equipment. Understanding the best way to plan and manipulate a cable without disrupting the environment or damaging the cable is imperative to robotic outfitting operations on the lunar surface.

Amy M Quartaro↗
NASA NTRS1983

Simplified Modeling of Tetrahedral Trusses

Tetrahedral-lattice structure exhibits great complexity in mathematical analysis, and complexity quickly increases with number of elements. Solutions for static and dynamic behavior of structure obtained with model consisting of analogous structure with fewer repeating elements. Advantage of continuum models is analytic solutions can be obtained, either in closed form or in relatively-tractable Taylor series.

Leondis, A.↗
NASA NTRS2004

Cross-Linked Nanotube Materials with Variable Stiffness Tethers

The constitutive properties of a cross-linked single-walled carbon nanotube material are predicted with a multi-scale model. The material is modeled as a transversely isotropic solid using concepts from equivalent-continuum modeling. The elastic constants are determined using molecular dynamics simulation. Some parameters of the molecular force field are determined specifically for the cross-linker from ab initio calculations. A demonstration of how the cross-linked nanotubes may affect the properties of a nanotube/polyimide composite is included using a micromechanical analysis.

Frankland, Sarah-Jane V.↗
NASA NTRS2009

Suzaku observation of IGR J16318-4848

We report on the first Suzaku observation of IGR J16318-4848, the most extreme example of a new group of highly absorbed X-ray binaries that have recently been discovered by the International Gamma-Ray Astrophysics Laboratory (INTEGRAL). The Suzaku observation was carried out between 2006 August 14 and 17, with a net exposure time of 97 ks. The average X-ray spectrum of the source can be well described (chi-square (sub red)= 0.99) with a continuum model typical for neutron stars i.e., a strongly absorbed power law continuum with a photon index of 0.676(42) and an exponential cutoff at 20.5(6) keY. The absorbing column is N(sub H) = 1.95(3) X 10(exp 24)/square cm. Consistent with earlier work, strong fluorescent emission lines of Fe K-alpha, Fe K-beta, and Ni K-alpha are observed. Despite the large N(sub H), no Compton shoulder is seen in the lines, arguing for a non-spherical and inhomogeneous absorber. Seen at an average 5-60 keV absorbed flux of 3.4 x 10(exp -10) erg/square cm/second, the source exhibits significant variability on timescales of hours.

Barragan, Laura↗
NASA NTRS2018

Evidence for an Optical Low-Frequency Quasi-Periodic Oscillation in the Kepler Light Curve of an Active Galaxy

We report evidence for a quasi-periodic oscillation (QPO) in the optical light curve of KIC 9650712, a narrow-line Seyfert 1 galaxy in the original Kepler field. After the development and application of a pipeline for Kepler data specific to active galactic nuclei (AGNs), one of our sample of 21 AGNs selected by infrared photometry and X-ray flux demonstrates a peak in the power spectrum at log ν = −6.58 Hz, corresponding to a temporal period of t = 44 days. We note that although the power spectrum is well fit by a model consisting of a Lorentzian and a single power law, alternative continuum models cannot be ruled out. From optical spectroscopy, we measure the black hole mass of this AGN as log (M(sub BH)/solar mass) = 8.17. We find that this frequency lies along a correlation between low-frequency QPOs and black hole mass from stellar and intermediate mass black holes to AGNs, similar to the known correlation in high-frequency QPOs.

quasi-periodic oscillation (QPO)↗
NASA NTRS1988

X-ray, optical, and IR luminosity correlations in quasars

For a sample of 88 radio-quiet quasars that includes 12 X-ray upper limits, it is found that the X-ray and IR luminosities follow a significantly nonlinear relation that is consistent with that found for the X-ray and optical luminosities. The redshift dependence is weak and may be nonexistent. Since the optical through IR continuum of quasars is much the same shape as a function of luminosity, this implies that the nonlinear relation between the X-ray luminosity and the optical and IR luminosities is due to a systematic change in the X-ray continuum of quasars with increasing luminosity. Either the X-ray continuum must flatten with increasing luminosity, or the overall normalization relative to the optical and IR must decrease. Simple continuum models which link the IR and X-ray portions of the continuum by a single power law arising from the same synchrotron mechanism can be ruled out.

Kriss, Gerard A.↗
NASA NTRS1990

Recent literature on structural modeling, identification, and analysis

The literature on the mathematical modeling of large space structures is first reviewed, with attention given to continuum models, model order reduction, substructuring, and computational techniques. System identification and mode verification are then discussed with reference to the verification of mathematical models of large space structures. In connection with analysis, the paper surveys recent research on eigensolvers and dynamic response solvers for large-order finite-element-based models.

Craig, Roy R., Jr.↗
NASA NTRS2006

Precision Fe K-Alpha and Fe K-Beta Line Spectroscopy of the Seyfert 1.9 Galaxy NGC 2992 with Suzaku

We present detailed time-averaged X-ray spectroscopy in the 0.5-10 keV band of the Seyfert 1.9 galaxy NGC 2992 with the Suzaku X-ray Imaging Spectrometers (XIS). The source had a factor approximately 3 higher 2-10 keV flux (approximately 1.2 x l0(exp -11) erg per square cm per s) than the historical minimum and a factor approximately 7 less than the historical maximum. The XIS spectrum of NGC 2992 can be described by several components. There is a primary continuum, modeled as a power-law with a photon index of Gamma = 1.57(sup +0.06) (sup -0.03) that is obscured by a Compton-thin absorber with a column density of 8.01(sup +0.6) (sup -0.5)x l0 (exp 21) per square cm. . There is another, weaker, unabsorbed power-law component (modeled with the same slope as the primary), that is likely to be due to the primary continuum being electron-scattered into our line-of-sight by a region extended on a scale of hundreds of parsecs. We measure the Thomson depth of the scattering zone to be Tau = 0.072 +/- 0.021. An optically-thin thermal continuum emission component, which probably originates in the same extended region, is included in the model and yields a temperature and luminosity of KT = 0.656(sup +0.088) (sup -0.0.61) keV and approximately 1.2 +/- 0.4 x l0 (exp 40) erg per s respectively. We detect an Fe K emission complex which we model with broad and narrow lines and we show that the intensities of the two components are decoupled at a confidence level > 3 sigma. The broad Fe K alpha line has an equivalent width of 118(sup +32) (sup -61) eV and could originate in an accretion disk (with inclination angle greater than approximately 30 deg) around the putative central black hole. The narrow Fe K alpha line has an equivalent width of 1632(sup +47) (sup -26) eV and is unresolved (FWHM < 4630 km per s) and likely originates in distant matter. The absolute flux in the narrow line implies that the column density out of the line-of-sight could be much higher than measured in the line-of-sight, and that the mean (historically-averaged) continuum luminosity responsible for forming the line could be a factor of several higher than that measured from the data. We also detect the Fe K Beta line (corresponding to the narrow Fe K alpha line) with a high signal-to-noise ratio and describe a new robust method to constrain the ionization state of Fe responsible for the Fe K alpha and Fe K Beta lines that does not require any knowledge of possible gravitational and Doppler energy shifts affecting the line energies. For the distant line-emitting matter (e. g. the putative obscuring torus) we deduce that the predominant ionization state is lower than Fe VIII (at 99% confidence), conservatively taking into account residual calibration uncertainties in the XIS energy scale and theoretical and experimental uncertainties in the Fe K fluorescent line energies. From the limits on a possible Compton-reflection continuum it is likely that the narrow Fe K alpha and Fe K Beta lines originate in a Compton-thin structure.

Yaqoob, Tahir↗
NASA NTRS1992

A piecewise continuous Timoshenko beam model for the dynamic analysis of tapered beam-like structures

Distributed parameter modeling offers a viable alternative to the finite element approach for modeling large flexible space structures. The introduction of the transfer matrix method into the continuum modeling process provides a very useful tool to facilitate the distributed parameter model applied to some more complex configurations. A uniform Timoshenko beam model for the estimation of the dynamic properties of beam-like structures has given comparable results. But many aeronautical and aerospace structures are comprised of non-uniform sections or sectional properties, such as aircraft wings and satellite antennas. This paper proposes a piecewise continuous Timoshenko beam model which is used for the dynamic analysis of tapered beam-like structures. A tapered beam is divided into several segments of uniform beam elements. Instead of arbitrarily assumed shape functions used in finite element analysis, the closed-form solution of the Timoshenko beam equation is used. Application of the transfer matrix method relates all the elements as a whole. By corresponding boundary conditions and compatible conditions a characteristic equation for the global tapered beam has been developed, from which natural frequencies can be derived. A computer simulation is shown in this paper, and compared with the results obtained from the finite element analysis. While piecewise continuous Timoshenko beam model decreases the number of elements significantly; comparable results to the finite element method are obtained.

Shen, Ji Yao↗
NASA NTRS2024

The Giant Out Burst of EXO 2030+375 II: Broadband Spectroscopy and Evolution

In 2021, the High-mass X-ray Binary EXO 2030+375 underwent a giant X-ray outburst, the first since 2006, which reached a peak flux of ∼600 mCrab (3–50 keV). The goal of this work is to study the spectral evolution over the course of the outburst, search for possible cyclotron resonance scattering features (CRSFs), and to associate spectral components with the emission pattern of the accretion column. We use broadband spectra taken with NuSTAR, NICER, and Chandra near the peak and during the decline phase of the outburst. We describe the data with established empirical continuum models and perform pulse-phase-resolved spectroscopy. We compare the spectral evolution with pulse phase with a proposed geometrical emission model. We find a significant spectral hardening toward lower luminosity, a behavior that is expected for super-critical sources. The continuum shape and evolution cannot be described by a simple power law model with exponential cutoff but requires additional absorption or emission components. We can confirm the presence of a narrow absorption feature at ∼10 keV in both NuSTAR observations. The absence of harmonics question the interpretation of this feature as a CRSF. The empirical spectral components cannot be directly associated with identified emission components from the accretion column.

R. Ballhausen↗
NASA NTRS1992

An ionized accretion disk in Cygnus X-1

Results of an analysis of archival Exosat and HEAO 1-A2 observations of Cyg X-1 are presented. It is found, in agreement with Ebisawa (1991), that the X-ray spectrum of Cyg X-1 exhibits a 'high-energy excess' above 8 keV, similar to that seen in Ginga spectra of AGN. It is shown, on the basis of a likelihood ratio test, that the data are significantly better fit by a Compton reflection model than by a partial covering model. Even with the energy resolution of the Exosat GSPC, it is found that the derived iron line energy and shape are very dependent on the continuum model used. The best-fit temperature for a Suniaev-Titarchuk Comptonization model is found to be significantly larger than that found in previous analyses based primarily on data in the 30-140 keV band. This suggests that the 5-50 keV continuum is not well fit by the same spectrum form applicable at high energies and that a multiple temperature thermal Compton model is necessary to fit the broadband data.

Done, C.↗