Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed parallelization”

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,693 records · Page 94
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Parallelization of Rocket Engine Simulator Software (PRESS)

We have outlined our work in the last half of the funding period. We have shown how a demo package for RESSAP using MPI can be done. However, we also mentioned the difficulties with the UNIX platform. We have reiterated some of the suggestions made during the presentation of the progress of the at Fourth Annual HBCU Conference. Although we have discussed, in some detail, how TURBDES/PUMPDES software can be run in parallel using MPI, at present, we are unable to experiment any further with either MPI or PVM. Due to X windows not being implemented, we are also not able to experiment further with XPVM, which it will be recalled, has a nice GUI interface. There are also some concerns, on our part, about MPI being an appropriate tool. The best thing about MPr is that it is public domain. Although and plenty of documentation exists for the intricacies of using MPI, little information is available on its actual implementations. Other than very typical, somewhat contrived examples, such as Jacobi algorithm for solving Laplace's equation, there are few examples which can readily be applied to real situations, such as in our case. In effect, the review of literature on both MPI and PVM, and there is a lot, indicate something similar to the enormous effort which was spent on LISP and LISP-like languages as tools for artificial intelligence research. During the development of a book on programming languages [12], when we searched the literature for very simple examples like taking averages, reading and writing records, multiplying matrices, etc., we could hardly find a any! Yet, so much was said and done on that topic in academic circles. It appears that we faced the same problem with MPI, where despite significant documentation, we could not find even a simple example which supports course-grain parallelism involving only a few processes. From the foregoing, it appears that a new direction may be required for more productive research during the extension period (10/19/98 - 10/18/99). At the least, the research would need to be done on Windows 95/Windows NT based platforms. Moreover, with the acquisition of Lahey Fortran package for PC platform, and the existing Borland C + + 5. 0, we can do work on C + + wrapper issues. We have carefully studied the blueprint for Space Transportation Propulsion Integrated Design Environment for the next 25 years [13] and found the inclusion of HBCUs in that effort encouraging. Especially in the long period for which a map is provided, there is no doubt that HBCUs will grow and become better equipped to do meaningful research. In the shorter period, as was suggested in our presentation at the HBCU conference, some key decisions regarding the aging Fortran based software for rocket propellants will need to be made. One important issue is whether or not object oriented languages such as C + + or Java should be used for distributed computing. Whether or not "distributed computing" is necessary for the existing software is yet another, larger, question to be tackled with.

Cezzar, Ruknet↗
NASA NTRS1997

Root-Raised Cosine Filter Implementation That Uses Canonical Signed Digits for High-Speed Digital Filter Applications

NASA Lewis Research Center's Space Communications Division has been investigating high-speed digital filters that can operate at a higher speed than those in current use for a digital modulator and demodulator (modem). Using the Canonical Signed Digits (CSD) number representation for filter coefficients is a very effective way to increase the filter's speed while reducing complexity in the digital filter hardware design. This approach is a good alternative to using an expensive parallel-processing design technique or custom, application-specific integrated circuits. Such integrated circuits may not be suitable for applications that require filter speeds faster than what application-specific integrated circuits digital signal processors can offer for a dedicated channel. When a communication channel is a dedicated, multiplication process--a costly, time-consuming process--it can be greatly simplified by a replacement of the filter coefficients with CSD numbers. A computer code written with the MATLAB software package runs the program and generates CSD-represented filter coefficients that are based on minimizing minimum mean square errors. Also, the Alta Group of Cadence's Signal Processing Workstation is used to simulate and analyze the CSD filter responses. The impulse response of the root-raised cosine filter that is used as a base model is defined. From this filter, a set of coefficients is sampled and stored in a file. For the all coefficients, the optimal CSD number for each coefficient is searched on the basis of the minimum-mean-square-errors criterion. Because the distribution of CSD numbers is not uniform, quantization errors tend to be bigger for coefficients greater than 1/2. To offset errors that occur in a region of coefficients between 1/2 to 1 and to better represent fractions with CSD numbers, an extra nonzero digit is allowed for any coefficients exceeding 1/2. This will greatly improve frequency response as well as intersymbol interference at the receiver. The frequency response of a set of collected CSD-represented filter coefficients was compared with the same filter that was conventionally implemented. Analyses show CSD-implemented filters perform as well as conventional filters. Comparison of eye diagrams and bit-error-rate curves between CSD filters and traditionally implemented filters are almost indistinguishable. However, filter complexity was reduced from almost 3.5 to 1 for CSD filters. Complete computer simulation results are available. In the near future, work will focus on building actual working digital filter hardware in a field programmable gate array (FPGA).

Kim, Heechul↗
NASA NTRS1992

Physical properties of neutral gas in M31 and the Galaxy

The present study analyzes, in parallel with published data for the Galaxy, neutral hydrogen (H I) absorption and deduced emission detected along seven lines of sight through the disk of M31. It is shown that the brightness temperature of H I emission is coupled to the opacity of the gas. The Galactic relationship shows asymptotic trends at both large and small opacities. A simple yet effective physical model which accounts for this behavior consists of only two independent components: a high-opacity cool component of fixed mean temperature, and a low-opacity warm component of fixed mean brightness. A lower mean gas pressure by a factor of about 2 is argued to be the most plausible mechanism for accounting for a higher cool-component H I temperature in M31. Deduced volume filling factors of the Galactic H I are about 1 and 15 percent, respectively, for the cool and warm components, while for M31 they are 8 and 30 percent. The large ratio of surface to volume filling factors for both cool and warm H I suggests that these components are distributed predominantly as large sheet- or shell-like structures.

Braun, Robert↗
NASA NTRS2014

Structural Characterization of Lateral-grown 6H-SiC am-plane Seed Crystals by Hot Wall CVD Epitaxy

The performance of commercially available silicon carbide (SiC) power devices is limited due to inherently high density of screw dislocations (SD), which are necessary for maintaining polytype during boule growth and commercially viable growth rates. The NASA Glenn Research Center (GRC) has recently proposed a new bulk growth process based on axial fiber growth (parallel to the c-axis) followed by lateral expansion (perpendicular to the c-axis) for producing multi-faceted m-plane SiC boules that can potentially produce wafers with as few as one SD per wafer. In order to implement this novel growth technique, the lateral homoepitaxial growth expansion of a SiC fiber without introducing a significant number of additional defects is critical. Lateral expansion is being investigated by hot wall chemical vapor deposition (HWCVD) growth of 6H-SiC am-plane seed crystals (0.8mm x 0.5mm x 15mm) designed to replicate axially grown SiC single crystal fibers. The post-growth crystals exhibit hexagonal morphology with approximately 1500 m (1.5 mm) of total lateral expansion. Preliminary analysis by synchrotron white beam x-ray topography (SWBXT) confirms that the growth was homoepitaxial, matching the polytype of the respective underlying region of the seed crystal. Axial and transverse sections from the as grown crystal samples were characterized in detail by a combination of SWBXT, transmission electron microscopy (TEM) and Raman spectroscopy to map defect types and distribution. X-ray diffraction analysis indicates the seed crystal contained stacking disorders and this appears to have been reproduced in the lateral growth sections. Analysis of the relative intensity for folded transverse acoustic (FTA) and optical (FTO) modes on the Raman spectra indicate the existence of stacking faults. Further, the density of stacking faults is higher in the seed than in the grown crystal. Bundles of dislocations are observed propagating from the seed in m-axis lateral directions. Contrast extinction analysis of these dislocation lines reveals they are edge type basal plane dislocations that track the growth direction. Polytype phase transition and stacking faults were observed by high-resolution TEM (HRTEM), in agreement with SWBXT and Raman scattering.

Single Crystal↗
NASA NTRS1994

Shock propagation and the generation of magnetohydrodynamic wave fields in inhomogeneous molecular clouds

We develop a simple one-dimensional model for the interaction of a steady, thin, planar shock wave with a nonrigid cloud which may be in motion relative to the surrounding medium, and we apply the model to shocks impinging on, and propagating through, molecular clouds. Both 'adiabatic' (gamma = 5/3) and radiative (gamma = 1) shocks are considered and we allow for the presence of a uniform magnetic field directed either parallel or perpendicular to the shock normal. The former field orientation is equivalent to the hydrodynamic case, and the latter involves only fast MHD shocks. We focus on the manner in which such shocks can generate internal kinetic motions in the cloud on a range of size and density scales through the direct acceleration of cores and clumps by shocks transmitted into them and through the generation of an MHD wavefield via the reflection of the incident shock at clump boundaries. We find that stronger incident Mach numbers and smaller density contrasts lead to more efficient cloud acceleration, as do isothermal intercloud shocks and small intercloud magnetic field strengths. The acceleration efficiency is insensitive to the adiabatic index and the magnetic field strength in the cloud itself. For typical parameter choices, the direct acceleration of clouds and clumps by strong shocks is found to be substantial and could at least in part account for their observed velocity dispersions. If the shocks are moderately weak, the final velocity of the cloud is linearly related to its initial velocity, with higher acceleration giving shallower slopes (i.e., final velocity distributions which are less sensitive to the initial distribution). Compared to the kinetic energy of the postshock cloud, the energy given to the wavefield at each encounter is small, and the heating of the interclump medium by the dissipation of this wavefield is found to be insufficient to balance the cooling rate in the cloud as a whole (although it may be important in particular regions), even if this medium is warm, unless it is also extremely tenuous (n approximately less than 0.1/cu cm). Nevertheless, the correction for the velocity imparted to the cloud leads to a substantial increase in the critical incident Mach numner for wave emission over that reported by Spitzer for the rigid case. The implications of our model for shock-induced star formation are discussed briefly.

Miesch, Mark S.↗
NASA NTRS2011

Catalytic Microtube Rocket Igniter

Devices that generate both high energy and high temperature are required to ignite reliably the propellant mixtures in combustion chambers like those present in rockets and other combustion systems. This catalytic microtube rocket igniter generates these conditions with a small, catalysis-based torch. While traditional spark plug systems can require anywhere from 50 W to multiple kW of power in different applications, this system has demonstrated ignition at less than 25 W. Reactants are fed to the igniter from the same tanks that feed the reactants to the rest of the rocket or combustion system. While this specific igniter was originally designed for liquid methane and liquid oxygen rockets, it can be easily operated with gaseous propellants or modified for hydrogen use in commercial combustion devices. For the present cryogenic propellant rocket case, the main propellant tanks liquid oxygen and liquid methane, respectively are regulated and split into different systems for the individual stages of the rocket and igniter. As the catalyst requires a gas phase for reaction, either the stored boil-off of the tanks can be used directly or one stream each of fuel and oxidizer can go through a heat exchanger/vaporizer that turns the liquid propellants into a gaseous form. For commercial applications, where the reactants are stored as gases, the system is simplified. The resulting gas-phase streams of fuel and oxidizer are then further divided for the individual components of the igniter. One stream each of the fuel and oxidizer is introduced to a mixing bottle/apparatus where they are mixed to a fuel-rich composition with an O/F mass-based mixture ratio of under 1.0. This premixed flow then feeds into the catalytic microtube device. The total flow is on the order of 0.01 g/s. The microtube device is composed of a pair of sub-millimeter diameter platinum tubes connected only at the outlet so that the two outlet flows are parallel to each other. The tubes are each approximately 10 cm long and are heated via direct electric resistive heating. This heating brings the gasses to their minimum required ignition temperature, which is lower than the auto-thermal ignition temperature, and causes the onset of both surface and gas phase ignition producing hot temperatures and a highly reacting flame. The combustion products from the catalytic tubes, which are below the melting point of platinum, are injected into the center of another combustion stage, called the primary augmenter. The reactants for this combustion stage come from the same source but the flows of non-premixed methane and oxygen gas are split off to a secondary mixing apparatus and can be mixed in a near-stoichiometric to highly lean mixture ratio. The primary augmenter is a component that has channels venting this mixed gas to impinge on each other in the center of the augmenter, perpendicular to the flow from the catalyst. The total crosssectional area of these channels is on a similar order as that of the catalyst. The augmenter has internal channels that act as a manifold to distribute equally the gas to the inward-venting channels. This stage creates a stable flame kernel as its flows, which are on the order of 0.01 g/s, are ignited by the combustion products of the catalyst. This stage is designed to produce combustion products in the flame kernel that exceed the autothermal ignition temperature of oxygen and methane.

Schneider, Steven J.↗
NASA NTRS2000

Distribution of Slip at the Northern Sumatran Fault System

We model spatial variations in horizontal displacements of 117 geodetic sites measured during annual surveys in 1989-1996 with the Global Positioning System (GPS) as elastic strain across a locked strike-slip fault to infer the contemporary slip rate, locking depth, and location of the Sumatran fault (SF) in northern Sumatra (1 S-3 N). GPS-derived slip rate estimates increase slightly northward from 23 plus or minus 3 mm/yr at 0.8 deg S to 26 plus or minus 2mm/yr at 2.7 N. They agree with geologic estimates north of the Equator, but at 0.5 S they are about 10 mm/yr higher. Strain appears to be distributed asymmetrically about the fault. South of 2 N, about 5 mm/yr of shear is required within the offshore forearc, west of the fault, to achieve a closer agreement of fault locations inferred from GPS velocities with geologically identified traces of the SF. Locking depth estimates are on the order of 10-20 km. The western branch of the major fault bifurcation near 1 N slips at a rate five times higher than the eastern branch. The two main strands of the fault at the northwestern tip of Sumatra (5.5 N) appear to be nearly free of horizontal strain; significant slip must occur away from the two strands, probably further east at two other geologically active branches. The Banda Aceh embayment is extruded to the northwest at a rate of 5 plus or minus 2 mm/yr. Within the estimated velocity uncertainties of several mm/yr, fault-normal deformation along the SF is insignificant. Almost strain free, the northern part of the back-arc basin is part of a rigid Sunda shelf, while the northern forearc is subjected to 8 plus or minus 5 x 10 (exp -8)/yr of extension nearly parallel to the arc.

Genrich, J. F.↗
NASA NTRS2005

ICESat Observations of Topographic Change in the Northern Segment of the 2004 Sumatra-Andaman Islands Earthquake Rupture Zone

The Andaman Islands are located 120 km east of the Sunda trench in the northern quarter of the 1300 km long rupture zone of the 2004 Sumatra-Andaman Islands earthquake inferred from the distribution of aftershocks. Initial field reports indicate that several meters of uplift and up to a meter of submergence occurred on the western and eastern shorelines of the Andaman Islands, respectively, associated with the earthquake (Bilham, 2005). Satellite images also document uplift of western shoreline coral reef platforms above sea level. Body-wave (Ji, 2005; Yamamaka, 2005) and tide-gauge (Ortiz, 2005) slip inversions only resolve coseismic slip in the southern one-third to one-half of the rupture zone. The amount of coseismic slip in the Andaman Islands region is poorly constrained by these inversions. The Ice, Cloud, and land Elevation Satellite (ICESat), a part of the NASA Earth Observing System, is being used to document the spatial pattern of Andaman Islands vertical displacements in order to constrain models of slip distribution in the northern part of the rupture zone. ICESat carries the Geoscience Laser Altimeter System (GLAS) that obtains elevation measurements from 80 m diameter footprints spaced 175 m apart along profiles. For surfaces of low slope, single-footprint absolute elevation and horizontal accuracies of 10 cm and 6 m (1 sigma), respectively, referenced to the ITRF 2002 TOPEX/Poseidon ellipsoid are being obtained. Laser pulse backscatter waveforms enable separation of ground topography and overlying vegetation cover. During each 33-day observing period ICESat acquires three profiles crossing the Andaman Islands. A NNE-SSW oriented track consists of 1600 laser footprints along the western side of North, Middle, and South Andaman Islands and 240 laser footprints across the center of Great Andaman Island. Two NNW-SSE tracks consist of 440 footprints across Middle Andaman Island and 25 footprints across the west side of Sentinel Island. Cloud-free profiles were acquired in the fall of 2003 and 2004. During February-March, 2005 ICESat's precise pointing capability will be used to exactly repeat these three profiles, with a cross-track accuracy of better than 100 m, providing trench- parallel and -perpendicular observations of topographic change of the Andaman Islands that will compliment geodetic field surveys. The observed elevation changes will be compared to models of coseismic deformation associated with the mainshock and large aftershocks in the Andaman Islands region.

Harding, David↗
NASA NTRS2014

The Role of Higher-Order Modes on the Electromagnetic Whistler-Cyclotron Wave Fluctuations of Thermal and Non-Thermal Plasmas

Two fundamental challenging problems of laboratory and astrophysical plasmas are the understanding of the relaxation of a collisionless plasmas with nearly isotropic velocity distribution functions and the resultant state of nearly equipartition energy density with electromagnetic plasma turbulence. Here, we present the results of a study which shows the role that higher-order-modes play in limiting the electromagnetic whistler-like fluctuations in a thermal and non-thermal plasma. Our main results show that for a thermal plasma the magnetic fluctuations are confined by regions that are bounded by the least-damped higher order modes. We further show that the zone where the whistler-cyclotron normal modes merges the electromagnetic fluctuations shifts to longer wavelengths as the beta(sub e) increases. This merging zone has been interpreted as the beginning of the region where the whistler-cyclotron waves losses their identity and become heavily damped while merging with the fluctuations. Our results further indicate that in the case of nonthermal plasmas, the higher-order modes do not confine the fluctuations due to the effective higher-temperature effects and the excess of suprathermal plasma particles. The analysis presented here considers the second-order theory of fluctuations and the dispersion relation of weakly transverse fluctuations, with wave vectors parallel to the uniform background magnetic field, in a finite temperature isotropic bi-Maxwellian and Tsallis-kappa-like magnetized electron-proton plasma. Our results indicate that the spontaneously emitted electromagnetic fluctuations are in fact enhanced over these quasi modes suggesting that such modes play an important role in the emission and absorption of electromagnetic fluctuations in thermal or quasi-thermal plasmas.

Vinas, Adolfo F.↗
NASA NTRS2021

An Efficient and Accurate Algorithm for Computing Grid-Averaged Solar Fluxes for Horizontally Inhomogeneous Clouds

A computationally efficient method is presented to account for the horizontal cloud inhomogeneity by using a radiatively equivalent plane parallel homogeneous (PPH) cloud. The algorithm can accurately match the calculations of the reference (rPPH) independent column approximation (ICA) results, but use only the same computational time required for a single plane parallel computation. The effective optical depth of this synthetic sPPH cloud is derived by exactly matching the direct transmission to that of the inhomogeneous ICA cloud. The ffective9 scattering asymmetry factor is found from a pre-calculated albedo inverse look-up-table that is allowed to vary over the range from -1.0 to 1.0. In the special cases of conservative scattering and total absorption, the synthetic method is exactly equivalent to the ICA, with only a small bias (about 0.2% in flux) relative to ICA due to imperfect interpolation in using the look-up tables. In principle, the ICA albedo can be approximated accurately regardless of cloud inhomogeneity. For a more complete comparison, the broadband shortwave albedo and transmission calculated from the synthetic sPPH cloud and averaged over all incident directions, have the RMS biases of 0.26% and 0.76%, respectively, for inhomogeneous clouds over a wide variation of particle size. The advantages of the synthetic PPH method are that (1) it is not required that all the cloud subcolumns have uniform microphysical characteristic, (2) it is applicable to any 1D radiative transfer scheme, and (3) it can handle arbitrary cloud optical depth distributions and an arbitrary number of cloud subcolumns with uniform computational efficiency.

cloud inhomogeneity↗
NASA NTRS2006

"One-Stop Shopping" for Ocean Remote-Sensing and Model Data

OurOcean Portal 2.0 (http:// ourocean.jpl.nasa.gov) is a software system designed to enable users to easily gain access to ocean observation data, both remote-sensing and in-situ, configure and run an Ocean Model with observation data assimilated on a remote computer, and visualize both the observation data and the model outputs. At present, the observation data and models focus on the California coastal regions and Prince William Sound in Alaska. This system can be used to perform both real-time and retrospective analyses of remote-sensing data and model outputs. OurOcean Portal 2.0 incorporates state-of-the-art information technologies (IT) such as MySQL database, Java Web Server (Apache/Tomcat), Live Access Server (LAS), interactive graphics with Java Applet at the Client site and MatLab/GMT at the server site, and distributed computing. OurOcean currently serves over 20 real-time or historical ocean data products. The data are served in pre-generated plots or their native data format. For some of the datasets, users can choose different plotting parameters and produce customized graphics. OurOcean also serves 3D Ocean Model outputs generated by ROMS (Regional Ocean Model System) using LAS. The Live Access Server (LAS) software, developed by the Pacific Marine Environmental Laboratory (PMEL) of the National Oceanic and Atmospheric Administration (NOAA), is a configurable Web-server program designed to provide flexible access to geo-referenced scientific data. The model output can be views as plots in horizontal slices, depth profiles or time sequences, or can be downloaded as raw data in different data formats, such as NetCDF, ASCII, Binary, etc. The interactive visualization is provided by graphic software, Ferret, also developed by PMEL. In addition, OurOcean allows users with minimal computing resources to configure and run an Ocean Model with data assimilation on a remote computer. Users may select the forcing input, the data to be assimilated, the simulation period, and the output variables and submit the model to run on a backend parallel computer. When the run is complete, the output will be added to the LAS server for

Li, P. Peggy↗
NASA NTRS1999

Study of Magnetic Damping Effect on Convection and Solidification Under G-Jitter Conditions

As shown by NASA resources dedicated to measuring residual gravity (SAMS and OARE systems), g-jitter is a critical issue affecting space experiments on solidification processing of materials. This study aims to provide, through extensive numerical simulations and ground based experiments, an assessment of the use of magnetic fields in combination with microgravity to reduce the g-jitter induced convective flows in space processing systems. We have so far completed asymptotic analyses based on the analytical solutions for g-jitter driven flow and magnetic field damping effects for a simple one-dimensional parallel plate configuration, and developed both 2-D and 3-D numerical models for g-jitter driven flows in simple solidification systems with and without presence of an applied magnetic field. Numerical models have been checked with the analytical solutions and have been applied to simulate the convective flows and mass transfer using both synthetic g-jitter functions and the g-jitter data taken from space flight. Some useful findings have been obtained from the analyses and the modeling results. Some key points may be summarized as follows: (1) the amplitude of the oscillating velocity decreases at a rate inversely proportional to the g-jitter frequency and with an increase in the applied magnetic field; (2) the induced flow approximately oscillates at the same frequency as the affecting g-jitter, but out of a phase angle; (3) the phase angle is a complicated function of geometry, applied magnetic field, temperature gradient and frequency; (4) g-jitter driven flows exhibit a complex fluid flow pattern evolving in time; (5) the damping effect is more effective for low frequency flows; and (6) the applied magnetic field helps to reduce the variation of solutal distribution along the solid-liquid interface. Work in progress includes numerical simulations and ground-based measurements. Both 2-D and 3-D numerical simulations are being continued to obtain further information on g-jitter driven flows and magnetic field effects. A physical model for ground-based measurements is completed and some measurements of the oscillating convection are being taken on the physical model. The comparison of the measurements with numerical simulations is in progress. Additional work planned in the project will also involve extending the 2-D numerical model to include the solidification phenomena with the presence of both g-jitter and magnetic fields.

Li, Ben Q.↗
NASA NTRS2002

Mechanisms of Gravity-Evoked Neuronal Plasticity

The grant focuses on a gene we identified called, serum and glucocorticoid- induced kinase (SGK), during a previously funded NASA project. The abundance of SGK messenger RNA (mRNA) and protein is increased in CNS tissues from animals reared in microgravity in comparison with 1G reared animals. In the funded proposal we had three aims: 1) characterize the distribution of SGK mRNA in the developing and adult rat CNS, 2) determine if expression of enzymatically active or inactive forms of SGK in cells influenced cell morphology (neurite growth), and 2) determine if SGK is a CREB kinase - that is, a protein kinase that adds phosphate groups to the transcription factor CREB. Over the past year we have made strong progress in the two most difficult parts of the project, namely specific aims 2 and 3. In specific aim #2 we planned to express a dominant negative or a constitutively active form of SGK in PC12 cells and assay the effects on neurite growth. Several methods are available for examining the effects of a transgene on PC12 neurite growth. Relevant variables include the performance of the assay +/- serum, +/- NGF, substratum for growth, timing between transfection and assay. Over the past 8 months we have customized the assay to enable us to most readily determine the effects of transgene expression on neurite growth. We have also compared the relative utility of transfecting DNA as opposed to protein itself. We are now well positioned to study the effects of SGK on neurite growth. We have also made progress in parallel studies in primary neurons. We have made constructs which will lead to transgene expression in cultures of spinal cord neurons. Co-transfection of a reporter and the SGK constructs can now be performed.

Kalb, Robert↗
NASA NTRS2017

Simulation Propulsion System and Trajectory Optimization

A number of new aircraft concepts have recently been proposed which tightly couple the propulsion system design and operation with the overall vehicle design and performance characteristics. These concepts include propulsion technology such as boundary layer ingestion, hybrid electric propulsion systems, distributed propulsion systems and variable cycle engines. Initial studies examining these concepts have typically used a traditional decoupled approach to aircraft design where the aerodynamics and propulsion designs are done a-priori and tabular data is used to provide inexpensive look ups to the trajectory ana-ysis. However the cost of generating the tabular data begins to grow exponentially when newer aircraft concepts require consideration of additional operational parameters such as multiple throttle settings, angle-of-attack effects on the propulsion system, or propulsion throttle setting effects on aerodynamics. This paper proposes a new modeling approach that eliminated the need to generate tabular data, instead allowing an expensive propulsion or aerodynamic analysis to be directly integrated into the trajectory analysis model and the entire design problem optimized in a fully coupled manner. The new method is demonstrated by implementing a canonical optimal control problem, the F-4 minimum time-to-climb trajectory optimization using three relatively new analysis tools: Open M-DAO, PyCycle and Pointer. Pycycle and Pointer both provide analytic derivatives and Open MDAO enables the two tools to be combined into a coupled model that can be run in an efficient parallel manner that helps to cost the increased cost of the more expensive propulsion analysis. Results generated with this model serve as a validation of the tightly coupled design method and guide future studies to examine aircraft concepts with more complex operational dependencies for the aerodynamic and propulsion models.

Optimization↗
NASA NTRS2011

Fire Return Interval Within the Northern Boundary of the Larch Forest

Larch (Larix spp.) dominant forests compose a large proportion of the forests of Russia (i.e., about 40% of forested areas). These forests range from the Yenisei ridge on the west to the Pacific Ocean on the east, and from Lake Baikal on the south to the 73rd parallel in the north. Larch stands comprise the world s northern most forest at Ary-Mas (72 deg 28' N, 102 deg 15' E). Larch dominated forests occupy about 70% of the permafrost areas in Siberia. Larch forms high closure stands as well as open forests, and is found mainly over permafrost, where other tree species barely survive. Wildfires are typical for this territory with the majority occurring as ground fires due to low crown closure. Due to the thin active layer in permafrost soils and a dense lichen-moss cover, ground fires may cause stand mortality. The vast areas of larch-dominant forests is generally considered as a "carbon sink"; however, positive long-term temperature trends at higher latitudes are expected to result in an increase of fire frequency, and thus may convert this area to a source for greenhouse gases. There are recent observations regarding the increase of fire frequency within non-protected territories. Surprisingly, there are few publications on fire chronoseqences for the huge forested territory between the Ural Mountains and the Pacific Ocean. Also there is a general understanding that bimodal (late spring -- early summer and late summer-beginning of fall) fire seasonal distribution in the south becomes uni-modal (late spring -- early summer) in the north. The purpose of this study is to investigate the wildfire history at the northern edge of the zone of larch dominance.

Kharuk, V. I.↗
NASA NTRS1999

Investigation of the Influence of Microgravity on Transport Mechanisms in a Virtual Spaceflight Chamber: A Ground Based Program

In January 1992, the IML-1 FES experiment produced a set of classic experimental data and a 40 hour holographic "movie" of an ensemble of spheres in a fluid in microgravity. Because the data are in the form of holograms, we can study the three-dimensional distribution of particles with unprecedented detail by a variety of methods and for a wide variety of interests. The possession of the holographic movie is tantamount to having a complex experiment in space while working in an easily accessible laboratory on earth. The movie contains a vast amount of useful data, including residual g, g-jitter, convection and transport data, and particle fluid interaction data. The information content in the movie is so great that we have scarcely begun to tap into the data that is actually available in the more than 1000 holograms, each containing as much as 1000 megabytes of information. This ground-based project is exploiting this data and the concept of holographic storage of spaceflight data to provide an understanding of the effects of microgravity in materials processing. This paper provides the foundation, objectives, and status of the ground based project. The primary objective of this project is to advance the understanding of microgravity effects on crystal growth, convection in materials processing in the space environment, and complex transport phenomena at low Reynolds numbers. This objective is being achieved both experimentally and theoretically. Experiments are making use of existing holographic data recorded during the IML- I spaceflight. A parallel theoretical effort is providing the models for understanding the particle fields and their physics in the microgravity environment.

Trolinger, James D.↗
NASA NTRS1991

Dynamic modelling and analysis of a magnetically suspended flexible rotor

A 12-state lumped-element model is presented for a flexible rotor supported by two attractive force electromagnetic journal bearings. The rotor is modeled as a rigid disk with radial mass unbalance mounted on a flexible, massless shaft with internal damping (Jeffcott rotor). The disk is offset axially from the midspan of the shaft. Bearing dynamics in each radial direction are modeled as a parallel combination of a negative (unstable) spring and a linear current-to-force actuator. The model includes translation and rotation of the rigid mass and the first and second bending models of the flexible shaft, and it simultaneously includes internal shaft damping, gyroscopic effects, and the unstable nature of the attractive force magnetic bearings. The model is used to analyze the dependence of the system transmission zeros and open-loop poles on system parameters. The dominant open-loop poles occur in stable/unstable pairs with bandwidth dependent on the ratios of bearing (unstable) stiffnesses to rotor mass and damping dependent on the shaft spin rate. The zeros occur in complex conjugate pairs with bandwidth dependent on the ratios of shaft stiffness to rotor mass and damping dependent on the shaft spin rate. Some of the transmission zeros are non-minimum phase when the spin rate exceeds the shaft critical speed. The transmission zeros and open-loop poles impact the design of magnetic bearing control systems. The minimum loop cross-over frequency of the closed-loop system is the speed of the unstable open-loop poles. For the supercritical shaft spin rates, the presence of non-minimum phase zeros limits the distribution rejection achievable at frequencies near or above the shaft critical speed. Since non-minimum phase transmission zeros can only be changed by changing the system inputs and/or outputs, closed-loop performance is limited for supercritical spin rates unless additional force or torque actuators are added.

Mccallum, Duncan C.↗
NASA NTRS2019

Dynamic Smagorinsky Modeled Large-Eddy Simulations of Turbulence Using Tetrahedral Meshes

Eddy-resolving numerical computations of turbulent flows are emerging as viable alternatives to Reynolds Averaged Navier-Stokes (RANS) calculations for flows with an intrinsically steady mean state due to the advances in large-scale parallel computing. In these computations, medium to large turbulent eddies are resolved by the numerics while the smaller or subgrid scales are either modeled or taken care of by the inherent numerical dissipation. To advance the state of the art of unstructured-mesh turbulence simulation capabilities, large eddy simulations (LES) using the dynamic Smagorinsky model (DSM) on tetrahedral meshes are carried out with the space-time conservation element, solution element (CESE) method. In contrast to what has been reported in the literature, the present implementation of dynamic models allows for active backscattering without any ad-hoc limiting of the eddy viscosity calculated from the subgrid-scale model. For the benchmark problems involving compressible isotropic turbulence decay as well as the shock/turbulent boundary layer interaction benchmark problems, no numerical instability associated with kinetic energy growth is observed and the volume percentage of the backscattering portion accounts for about 38-40% of the simulation domain. A slip-wall model in conjunction with the implemented DSM is used to simulate a relatively high Reynolds number Mach 2.85 turbulent boundary layer over a 30° ramp with several tetrahedral meshes and a wall-normal spacing of either Δγ& = 10 or Δγ& = 20. The computed mean wall pressure distribution, separation region size, mean velocity profiles, and Reynolds stress agree reasonably well with experimental data.

CESE method↗