Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “consistent scheme”

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 649 records · Page 36

Conical Euler solution for a highly-swept delta wing undergoing wing-rock motion

Modifications to an unsteady conical Euler code for the free-to-roll analysis of highly-swept delta wings are described. The modifications involve the addition of the rolling rigid-body equation of motion for its simultaneous time-integration with the governing flow equations. The flow solver utilized in the Euler code includes a multistage Runge-Kutta time-stepping scheme which uses a finite-volume spatial discretization on an unstructured mesh made up of triangles. Steady and unsteady results are presented for a 75 deg swept delta wing at a freestream Mach number of 1.2 and an angle of attack of 30 deg. The unsteady results consist of forced harmonic and free-to-roll calculations. The free-to-roll case exhibits a wing rock response produced by unsteady aerodynamics consistent with the aerodynamics of the forced harmonic results. Similarities are shown with a wing-rock time history from a low-speed wind tunnel test.

Lee, Elizabeth M.↗

Structure of confined laminar spray diffusion flames: Numerical investigation

The structure of confined laminar spray diffusion flames is investigated numerically by solving the gas-phase conservation equations for mass species, continuity, momentum, and energy and the liquid-phase equations for droplet position, velocity, size, and temperature. A one-step global reaction scheme along with six equilibrium reactions are employed to model the flame chemistry. Monodisperse as well as polydisperse sprays are considered. The numerical results demonstrate that liquid spray flames substantially differ from gaseous flames in their structure, i.e., temperature, concentration, and velocity fields, shape, and dimensions under the same conditions. Spray flames are predicted to be taller and narrower than their counterpart gaseous ones and their shapes are almost cylindrical. This is in agreement with experimental observations. The numerical computations also show that the use of the equilibrium reactions with the one-step reaction scheme decreases the flame temperature compared to the one-step reaction scheme without the equilibrium reactions and more importantly increases the surface area of the flame zone due to a phenomenon termed 'equilibrium broadening.' The spray flames also possess a finite thickness with minimal overlap of the fuel and oxygen species. A case for which a fuel-mixture consisting of 20 to 80 percent gas-liquid by mass is introduced into the combustor is also investigated and compared with predictions using only gaseous or liquid fuel.

Mawid, M. A.↗

RTE: A computer code for Rocket Thermal Evaluation

The numerical model for a rocket thermal analysis code (RTE) is discussed. RTE is a comprehensive thermal analysis code for thermal analysis of regeneratively cooled rocket engines. The input to the code consists of the composition of fuel/oxidant mixture and flow rates, chamber pressure, coolant temperature and pressure. dimensions of the engine, materials and the number of nodes in different parts of the engine. The code allows for temperature variation in axial, radial and circumferential directions. By implementing an iterative scheme, it provides nodal temperature distribution, rates of heat transfer, hot gas and coolant thermal and transport properties. The fuel/oxidant mixture ratio can be varied along the thrust chamber. This feature allows the user to incorporate a non-equilibrium model or an energy release model for the hot-gas-side. The user has the option of bypassing the hot-gas-side calculations and directly inputting the gas-side fluxes. This feature is used to link RTE to a boundary layer module for the hot-gas-side heat flux calculations.

Naraghi, Mohammad H. N.↗

Experimental Simulation of Active Control With On-line System Identification on Sound Transmission Through an Elastic Plate

An adaptive control algorithm with on-line system identification capability has been developed. One of the great advantages of this scheme is that an additional system identification mechanism such as an additional uncorrelated random signal generator as the source of system identification is not required. A time-varying plate-cavity system is used to demonstrate the control performance of this algorithm. The time-varying system consists of a stainless-steel plate which is bolted down on a rigid cavity opening where the cavity depth was changed with respect to time. For a given externally located harmonic sound excitation, the system identification and the control are simultaneously executed to minimize the transmitted sound in the cavity. The control performance of the algorithm is examined for two cases. First, all the water was drained, the external disturbance frequency is swept with 1 Hz/sec. The result shows an excellent frequency tracking capability with cavity internal sound suppression of 40 dB. For the second case, the water level is initially empty and then raised to 3/20 full in 60 seconds while the external sound excitation is fixed with a frequency. Hence, the cavity resonant frequency decreases and passes the external sound excitation frequency. The algorithm shows 40 dB transmitted noise suppression without compromising the system identification tracking capability.

Source record↗

Use of Combined A-Train Observations to Validate GEOS Model Simulated Dust Distributions During NAMMA

During August 2006, the NASA African Multidisciplinary Analyses Mission (NAMMA) field experiment was conducted to characterize the structure of African Easterly Waves and their evolution into tropical storms. Mineral dust aerosols affect tropical storm development, although their exact role remains to be understood. To better understand the role of dust on tropical cyclogenesis, we have implemented a dust source, transport, and optical model in the NASA Goddard Earth Observing System (GEOS) atmospheric general circulation model and data assimilation system. Our dust source scheme is more physically based scheme than previous incarnations of the model, and we introduce improved dust optical and microphysical processes through inclusion of a detailed microphysical scheme. Here we use A-Train observations from MODIS, OMI, and CALIPSO with NAMMA DC-8 flight data to evaluate the simulated dust distributions and microphysical properties. Our goal is to synthesize the multi-spectral observations from the A-Train sensors to arrive at a consistent set of optical properties for the dust aerosols suitable for direct forcing calculations.

Nowottnick, E.↗

Assembly and Integration Process of the First High Density Detector Array for the Atacama Cosmology Telescope

The Advanced ACTPol (AdvACT) upgrade on the Atacama Cosmology Telescope (ACT) consists of multichroicTransition Edge Sensor (TES) detector arrays to measure the Cosmic Microwave Background (CMB) polarization anisotropies in multiple frequency bands. The first AdvACT detector array, sensitive to both 150 and 230 GHz, is fabricated on a 150 mm diameter wafer and read out with a completely different scheme compared to ACTPol. Approximately 2000 TES bolometers are packed into the wafer leading to both a much denser detector density and readout circuitry. The demonstration of the assembly and integration of the AdvACT arrays is important for the next generation CMB experiments, which will continue to increase the pixel number and density. We present the detailed assembly process of the first AdvACT detector array.

Li, Yaqiong↗

Interplanetary shock normals

A scheme has been developed to improve the estimate of an interplanetary shock normal by using both magnetic field and plasma data from a single spacecraft. Calculation of the basic shock model employs a subset of the eight magnetohydrodynamic conservation relations for a shock in an isotropic medium. This subset consists of six equations that are devoid of pressure and temperature terms. A sigma-weighted least squares loss function technique is used to best fit the overdetermination equations with respect to the eleven parameters of the system. This procedure usually yields a normal much more accurate than the one obtained using magnetic field values alone. An example is given of the use of the technique on the data for the Pioneer 7 shock of 29 August 1966.

Lepping, R. P.↗

A study and experiment plan for digital mobile communication via satellite

The viability of mobile communications is examined within the context of a frequency division multiple access, single channel per carrier satellite system emphasizing digital techniques to serve a large population of users. The intent is to provide the mobile users with a grade of service consistant with the requirements for remote, rural (perhaps emergency) voice communications, but which approaches toll quality speech. A traffic model is derived on which to base the determination of the required maximum number of satellite channels to provide the anticipated level of service. Various voice digitalization and digital modulation schemes are reviewed along with a general link analysis of the mobile system. Demand assignment multiple access considerations and analysis tradeoffs are presented. Finally, a completed configuration is described.

Jones, J. J.↗

Elasticity Solution of an Adhesively Bonded Cover Plate of Various Geometries

The plane strain of adhesively bonded structures consisting of two different isotropic adherends is considered. By expressing the x-y components of the displacements in terms of Fourier integrals and using the corresponding boundary and continuity conditions, the integral equations for the general problem are obtained and solved numerically by applying Gauss-Chebyshev integration scheme. The shear and the normal stresses in the adhesive are calculated for various geometries and material properties for a stiffened plate under uniaxial tension. Numerical results involving the stress intensity factors and the strain energy release rate are presented. The closed-form expressions for the Fredholm kernels are provided to obtain the solution for an arbitrary geometry and material properties. For the general geometry, the contribution of the normal stress is quite significant, while for symmetric geometries, the shear stress is dominant, the normal stress vanishes if the adherends are of the same material and the same thickness.

Aksel, G. N.↗

Enhanced Schapery Theory Software Development for Modeling Failure of Fiber-Reinforced Laminates

Progressive damage and failure analysis (PDFA) tools are needed to predict the nonlinear response of advanced fiber-reinforced composite structures. Predictive tools should incorporate the underlying physics of the damage and failure mechanisms observed in the composite, and should utilize as few input parameters as possible. The purpose of the Enhanced Schapery Theory (EST) was to create a PDFA tool that operates in conjunction with a commercially available finite element (FE) code (Abaqus). The tool captures the physics of the damage and failure mechanisms that result in the nonlinear behavior of the material, and the failure methodology employed yields numerical results that are relatively insensitive to changes in the FE mesh. The EST code is written in Fortran and compiled into a static library that is linked to Abaqus. A Fortran Abaqus UMAT material subroutine is used to facilitate the communication between Abaqus and EST. A clear distinction between damage and failure is imposed. Damage mechanisms result in pre-peak nonlinearity in the stress strain curve. Four internal state variables (ISVs) are utilized to control the damage and failure degradation. All damage is said to result from matrix microdamage, and a single ISV marks the micro-damage evolution as it is used to degrade the transverse and shear moduli of the lamina using a set of experimentally obtainable matrix microdamage functions. Three separate failure ISVs are used to incorporate failure due to fiber breakage, mode I matrix cracking, and mode II matrix cracking. Failure initiation is determined using a failure criterion, and the evolution of these ISVs is controlled by a set of traction-separation laws. The traction separation laws are postulated such that the area under the curves is equal to the fracture toughness of the material associated with the corresponding failure mechanism. A characteristic finite element length is used to transform the traction-separation laws into stress-strain laws. The ISV evolution equations are derived in a thermodynamically consistent manner by invoking the stationary principle on the total work of the system with respect to each ISV. A novel feature is the inclusion of both pre-peak damage and appropriately scaled, post-peak strain softening failure. Also, the characteristic elements used in the failure degradation scheme are calculated using the element nodal coordinates, rather than simply the square root of the area of the element.

Pineda, Evan J.↗

Evaluation of CALIOP 532-nm Aerosol Optical Depth Over Opaque Water Clouds

With its height-resolved measurements and near global coverage, the CALIOP lidar onboard the CALIPSO satellite offers a new capability for aerosol retrievals in cloudy skies. Validation of these retrievals is difficult, however, as independent, collocated and co-temporal data sets are generally not available. In this paper, we evaluate CALIOP aerosol products above opaque water clouds by applying multiple retrieval techniques to CALIOP Level 1 profile data and comparing the results. This approach allows us to both characterize the accuracy of the CALIOP above-cloud aerosol optical depth (AOD) and develop an error budget that quantifies the relative contributions of different error sources. We focus on two spatial domains: the African dust transport pathway over the tropical North Atlantic and the African smoke transport pathway over the southeastern Atlantic. Six years of CALIOP observations (2007-2012) from the northern hemisphere summer and early fall are analyzed. The analysis is limited to cases where aerosol layers are located above opaque water clouds so that a constrained retrieval technique can be used to directly retrieve 532 nm aerosol optical depth and lidar ratio. For the moderately dense Sahara dust layers detected in the CALIOP data used in this study, the mean/median values of the lidar ratios derived from a constrained opaque water cloud (OWC) technique are 45.1/44.4 +/- 8.8 sr, which are somewhat larger than the value of 40 +/- 20 sr used in the CALIOP Level 2 (L2) data products. Comparisons of CALIOP L2 AOD with the OWC-retrieved AOD reveal that for nighttime conditions the L2 AOD in the dust region is underestimated on average by approx. 26% (0.183 vs. 0.247). Examination of the error sources indicates that errors in the L2 dust AOD are primarily due to using a lidar ratio that is somewhat too small. The mean/median lidar ratio retrieved for smoke is 70.8/70.4 +/- 16.2 sr, which is consistent with the modeled value of 70 +/- 28 sr used in the CALIOP L2 retrieval. Smoke AOD is found to be underestimated, on average, by approx. 39% (0.191 vs. 0.311). The primary cause of AOD differences in the smoke transport region is the tendency of the CALIOP layer detection scheme to prematurely assign layer base altitudes and thus underestimate the geometric thickness of smoke layers.

Liu, Z.↗

The Darkening of the Greenland Ice Sheet: Trends, Drivers and Projections (1981-2100)

The surface energy balance and meltwater production of the Greenland ice sheet (GrIS) are modulated by snow and ice albedo through the amount of absorbed solar radiation. Here we show, using space-borne multispectral data collected during the 3 decades from 1981 to 2012, that summertime surface albedo over the GrIS decreased at a statistically significant (99 %) rate of 0.02 decade(sup -1) between 1996 and 2012. Over the same period, albedo modelled by the Modele Atmospherique Regionale (MAR) also shows a decrease, though at a lower rate (approximately -0.01 decade(sup -1)) than that obtained from space-borne data. We suggest that the discrepancy between modelled and measured albedo trends can be explained by the absence in the model of processes associated with the presence of light-absorbing impurities. The negative trend in observed albedo is confined to the regions of the GrIS that undergo melting in summer, with the dry snow zone showing no trend. The period 1981-1996 also showed no statistically significant trend over the whole GrIS. Analysis of MAR outputs indicates that the observed albedo decrease is attributable to the combined effects of increased near-surface air temperatures, which enhanced melt and promoted growth in snow grain size and the expansion of bare ice areas, and to trends in light-absorbing impurities (LAI) on the snow and ice surfaces. Neither aerosol models nor in situ and remote sensing observations indicate increasing trends in LAI in the atmosphere over Greenland. Similarly, an analysis of the number of fires and BC emissions from fires points to the absence of trends for such quantities. This suggests that the apparent increase of LAI in snow and ice might be related to the exposure of a "dark band" of dirty ice and to increased consolidation of LAI at the surface with melt, not to increased aerosol deposition. Albedo projections through to the end of the century under different warming scenarios consistently point to continued darkening, with albedo anomalies averaged over the whole ice sheet lower by 0.08 in 2100 than in 2000, driven solely by a warming climate. Future darkening is likely underestimated because of known underestimates in modelled melting (as seen in hindcasts) and because the model albedo scheme does not currently include the effects of LAI, which have a positive feedback on albedo decline through increased melting, grain growth, and darkening.

surface energy↗

Retrieval of H2O Abundance in Titan's Stratosphere: A (Re)Analysis of CIRS/Cassini and PACS/Herschel Observations

Since its first measurement 20 years ago by the Infrared Space Observatory (ISO), the water (H2O) mole fraction in Titan’s stratosphere remains uncertain due to large differences between the determinations from available measurements. More particularly, the recent measurements made from the Herschel observatory (PACS (Photodetector Array Camera and Spectrometer) and HIFI (Heterodyne Instrument for the Far Infrared)) estimated the H2O mole fraction to be 0.023 parts per billion at 12.1 millibars. A mixing ratio of 0.14 parts per billion at 10.7 millibars was, however, retrieved from nadir spatially-resolved observations of Cassini/CIRS (Composite Infrared Spectrometer). At the same pressure level (10.7 millibars), this makes a difference of a factor of 5.5 between PACS and CIRS measurements, and this has notably prevented current models from fully constraining the oxygen flux flowing into Titan’s atmosphere. In this work, we try to understand the differences between the H2O mole fractions estimated from Herschel/PACS and Cassini/CIRS observations. The strategy for this is to 1) analyse recent disc-averaged observations of CIRS to investigate if the observation geometry could explain the previous observed differences, and 2) (re)analyse the three types of observation with the same retrieval scheme to assess if previous differences in retrieval codes/methodology could be responsible for the previous discrepancies. With this analysis, we show that using the same retrieval method better reconcile the previous measurements of these instruments. However, the addition of the disc-averaged CIRS observations, instead of confirming the consistency between the different datasets, reveals discrepancies between one of the CIRS disc-averaged set of observations and PACS measurements. This raises new questions regarding the possibility of latitudinal variations of H2O, which could be triggered by seasonal changes of the meridional circulation. As it has already been shown for nitriles and hydrocarbons, this circulation could potentially impact the latitudinal distribution of H2O through the subsidence or up- welling of air rich in H2O. The possible influence of spatial/time variations of the OH/H2O input flux in Titan’s atmosphere is also discussed. The analysis of more observations will be needed in future work to address the questions arising from this work and to improve the understanding of the sources of H2O in Titan’s atmosphere.

Bauduin, S.↗

HARM3D+NUC: A New Method for Simulating the Post-merger Phase of Binary Neutron Star Mergers with GRMHD, Tabulated EOS, and Neutrino Leakage

The first binary neutron star merger has already been detected in gravitational waves. The signal was accompanied by an electromagnetic counterpart including a kilonova component powered by the decay of radioactive nuclei, as well as a short γ-ray burst. In order to understand the radioactively powered signal, it is necessary to simulate the outflows and their nucleosynthesis from the post-merger disk. Simulating the disk and predicting the composition of the outflows requires general relativistic magnetohydrodynamical (GRMHD) simulations that include a realistic, finite-temperature equation of state (EOS) and self-consistently calculating the impact of neutrinos. In this work, we detail the implementation of a finite-temperature EOS and the treatment of neutrinos in the GRMHD code HARM3D+NUC, based on HARM3D. We include formal tests of both the finite-temperature EOS and the neutrino-leakage scheme. We further test the code by showing that, given conditions similar to those of published remnant disks following neutron star mergers, it reproduces both recombination of free nucleons to a neutron-rich composition and excitation of a thermal wind.

Ariadna Murguia-Berthier↗

The composite sequential clustering technique for analysis of multispectral scanner data

The clustering technique consists of two parts: (1) a sequential statistical clustering which is essentially a sequential variance analysis, and (2) a generalized K-means clustering. In this composite clustering technique, the output of (1) is a set of initial clusters which are input to (2) for further improvement by an iterative scheme. This unsupervised composite technique was employed for automatic classification of two sets of remote multispectral earth resource observations. The classification accuracy by the unsupervised technique is found to be comparable to that by traditional supervised maximum likelihood classification techniques. The mathematical algorithms for the composite sequential clustering program and a detailed computer program description with job setup are given.

Su, M. Y.↗

Unsupervised classification of earth resources data.

A new clustering technique is presented. It consists of two parts: (a) a sequential statistical clustering which is essentially a sequential variance analysis and (b) a generalized K-means clustering. In this composite clustering technique, the output of (a) is a set of initial clusters which are input to (b) for further improvement by an iterative scheme. This unsupervised composite technique was employed for automatic classification of two sets of remote multispectral earth resource observations. The classification accuracy by the unsupervised technique is found to be comparable to that by existing supervised maximum liklihood classification technique.

Su, M. Y.↗

Predictions for partial and monolayer coverages of O2 on graphite

Monolayer properties of O2 on graphite are calculated using a pattern recognition, optimization scheme. Equilibrium monolayers are predicted at two different densities with properties in agreement with recent X-ray diffraction, specific heat, and neutron scattering data. Properties of the extremely low density regime are calculated using a model based upon a distribution of two-dimensional O2 clusters. The results are consistent with experimental evidence.

Pan, R. P.↗

Nonlinear control of an experimental beam by IMSC

Results are reported from an experiment designed to control the vibratory motion of a beam at NASA, Langley Research Center. The experimental setup consists of a free-free uniform beam acted upon by four electromagnetic force actuators, with the motion being measured by nine displacement sensors. The entire assembly is linked to a CDC Cyber 175 computer which permits on-line real-time computation of the control forces. The control scheme is based on Independent Modal-Space Control (IMSC), whereby the modal forces are computed using a nonlinear, on-off control law. The actuator forces are then synthesized using a linear transformation, resulting in quantized forces. The sensors data is processed by modal filters. It is observed that the controls designed on the basis of the IMSC method are very effective in suppressing the vibratory motion of the beam, even though there is about a 50 percent differene between the actual and the analytically computed frequencies.

Meirovitch, L.↗