Search NASA⌕ Search

SEARCH · Search NASA

Results for “stochastic geometry”

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.

83 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

TURBOMAT: A Probabilistic Turbomachinery Aeroelastic Analysis Tool

An integration of aeroelastic analysis procedures with probabilistic analysis methods enables us to design safe reliable engines with quantified reliability. Towards this goal, a graphical user interface (GUI) based tool that integrates the codes Aeroelastic analysis of propfans (ASTROP2) and Numerical Evaluation of Stochastic Structures Under Stress (NESSUS) is developed. The tool entitled TURBOMachinery Aeroelastic Analysis Tool (TURBOMAT), is developed utilizing the MATrix Laboratory (Matlab) Guide (Graphical User Interface Development) tool box. TURBOMAT provides a user friendly computational environment for rapid assessment of Turbomachinery blades flutter characteristics, subjected to uncertain loading conditions with variability in material and aerodynamic properties. The tool is seen as an education tool for new students and young engineers starting their careers in structural Aeroelasticity who want to learn and understand aeroelastic aspects of turbomachinery components, fans, compressors and turbines, including uncertainties in loading and material properties.A typical fan blade configuration geometry was chosen to demonstrate the tool. The results are presented in the form of probabilistic density function (PDF), the cumulative distribution function (CDF) and sensitivity factors. Both first order fast probability integration (FPI) and the Monto Carlo (MC) techniques are used in the analysis and compared. The tool enabled us to quantify blade flutter reliability as well as the ranking of uncertain variables and their importance to blade flutter response.

Aeroelastic Analysis↗
NASA NTRS2004

Scatter in Carbon/Silicon Carbide (C/SiC) Composites Quantified

Carbon-fiber-reinforced silicon carbide matrix (C/SiC) composites processed by chemical vapor infiltration are candidate materials for aerospace thermal structures. Carbon fibers can retain properties at very high temperatures, but they are known to have poor oxidation resistance in adverse, high-temperature environments. Nevertheless, the combination of CVI-SiC matrix with higher stiffness and oxidation resistance, the interfacial coating, and additional surface-seal coating provides the necessary protection to the carbon fibers, and makes the material viable for high-temperature space applications operating under harsh environments. Furthermore, C/SiC composites, like other ceramic matrix composites (CMCs), exhibit graceful non-catastrophic failure because of various inherent energy dissipating mechanisms. The material exhibits nonlinearity in deformation even at very low stress levels. This is the result of the severe matrix microcracking present in the as processed composite because of large differences between the coefficients of thermal expansion of the fiber and the matrix. Utilization of these advanced composites in next generation space vehicles will require innovative structural configurations, updated materials, and refined analyses. Structural safety issues for these vehicles are in direct competition with performance and cost. One would have to quantify the uncertainties associated with the design using formal probabilistic methods. Specifically four fundamental aspects on which analyses are based-- (1) loading conditions, (2) material behavior, (3) geometrical configurations, and (4) structural connections between the composite components and baseline structure--are stochastic in nature. A direct way to formally account for uncertainties is to develop probabilistic structural analysis methods where all participating variables are described by appropriate probability density functions. The present work, however, focuses on analyzing the stochastic material behavior of these advanced composites using formal probabilistic analysis methods. Often, some of the desirable property characteristics that allow composites to offer advantages over conventional structural materials (like tailoring of composite properties) and the complexity are in fact responsible for their greater statistical variability and the requirements for more characterization tests. Composite properties are anisotropic as well, having different properties in different directions. This means that characterization of a property such as stiffness--which will vary greatly depending on the orientation of the fiber relative to the direction of the testing--must be repeated for several different directions and loading conditions. The fabrication process for composites also introduces statistical variations in properties and geometry. A composite part is produced in a number of steps, each of which introduces statistical variability. The matrix is usually produced from a combination of raw materials; and the fiber, which has its own set of properties, is often coated or surface treated, introducing yet another source of variability.

Murthy, Pappu L. N.↗
NASA NTRS2006

The Effect of Cumulus Cloud Field Anisotropy on Domain-Averaged Solar Fluxes and Atmospheric Heating Rates

Cumulus clouds can become tilted or elongated in the presence of wind shear. Nevertheless, most studies of the interaction of cumulus clouds and radiation have assumed these clouds to be isotropic. This paper describes an investigation of the effect of fair-weather cumulus cloud field anisotropy on domain-averaged solar fluxes and atmospheric heating rate profiles. A stochastic field generation algorithm was used to produce twenty three-dimensional liquid water content fields based on the statistical properties of cloud scenes from a large eddy simulation. Progressively greater degrees of x-z plane tilting and horizontal stretching were imposed on each of these scenes, so that an ensemble of scenes was produced for each level of distortion. The resulting scenes were used as input to a three-dimensional Monte Carlo radiative transfer model. Domain-average transmission, reflection, and absorption of broadband solar radiation were computed for each scene along with the average heating rate profile. Both tilt and horizontal stretching were found to significantly affect calculated fluxes, with the amount and sign of flux differences depending strongly on sun position relative to cloud distortion geometry. The mechanisms by which anisotropy interacts with solar fluxes were investigated by comparisons to independent pixel approximation and tilted independent pixel approximation computations for the same scenes. Cumulus anisotropy was found to most strongly impact solar radiative transfer by changing the effective cloud fraction, i.e., the cloud fraction when the field is projected on a surface perpendicular to the direction of the incident solar beam.

Hinkelman, Laura M.↗
NASA NTRS2008

Light Echos in Kerr Geometry: A Source of High Frequency QPOs from Random X-ray Bursts

We propose that high frequency quasi-periodic oscillations (HFQPOs) can be produced from randomly-formed X-ray bursts (flashes) by plasma interior to the ergosphere of a rapidly-rotating black hole. We show by direct computation of their orbits that the photons comprising the observed X-ray light curves, if due to a multitude of such flashes, are affected significantly by the black hole's dragging of inertial frames; the photons of each such burst arrive to an observer at infinity in multiple (double or triple), distinct 'bunches' separated by a roughly constant time lag of t/M approximately equal to 14, regardless of the bursts' azimuthal position. We argue that every other such 'bunch' represents photons that follow trajectories with an additional orbit around the black hole at the photon circular orbit radius (a photon 'echo'). The presence of this constant lag in the response function of the system leads to a QPO feature in its power density spectra, even though the corresponding light curve consists of a totally stochastic signal. This effect is by and large due to the black hole spin and is shown to gradually diminish as the spin parameter a decreases or the radial position of the burst moves outside the static limit surface (ergosphere). Our calculations indicate that for a black hole with Kerr parameter of a/M=0.99 and mass of M=10*Msun the QPO is expected at a frequency of approximately 1.3-1.4 kHz. We discuss the plausibility and observational implications of our model/results as well as its limitations.

Fukumura, K.↗
NASA NTRS2008

Light Echoes in Kerr Geometry: A Source of High Frequency QPOs from Random X-ray Bursts

We propose that high frequency quasi-periodic oscillations (HFQPOs) can be produced from randomly-formed X-ray bursts (flashes) by plasma interior to the ergosphere of a rapidly-rotating black hole. We show by direct computation of their orbits that the photons comprising the observed X-ray light curves, if due to a multitude of such flashes, are affected significantly by the black hole's dragging of inertial frames; the photons of each such burst arrive to an observer at infinity in multiple (double or triple), distinct "bunches" separated by a roughly constant time lag of Deltat(t(sub lag))/M approx. 14, regardless of the bursts' azimuthal position. We argue that every other such "bunch" represents photons that follow trajectories with an additional orbit around the black hole at the photon circular orbit radius (a photon "echo"). The presence of this constant lag in the response function of the system leads to a QPO feature in its power density spectra, even though the corresponding light curve consists of a totally stochastic signal. This effect is by and large due to the black hole spin and is shown to gradually diminish as the spin parameter a decreases or the radial position of the burst moves outside the static limit surface (ergosphere). Our calculations indicate that for a black hole with Kerr parameter of a/M = 0.99 and mass of M = 10Stellar Mass the QPO is expected at a frequency of v(sub QPO) approx. 1.3 - 1.4 kHz. We discuss the plausibility and observational implications of our model/results as well as its limitations. Subject headings: accretion, accretion disks - black hole physics - X-rays: galaxies - stars: oscillations

Fukumura, Keigo↗
NASA NTRSDate not supplied

A Parallelized Oxidation-Driven Surface Recession Framework in DSMC Code, SPARTA

Spacecrafts rely on ablative thermal protection systems (TPS) made of composites consisting of a carbon-based reinforcement and a polymeric matrix. These materials are designed to withstand high-temperature oxidation and surface recession during re-entry into the Earth's atmosphere. However, ablation occurs due to a complex interplay of thermal, mechanical, and chemical factors, making it challenging to determine the individual impact of each on the TPS's overall degradation. In this study, we have developed an ablation model that can leverage a finite rate carbon oxidation model to predict material recession and surface states more accurately. Stochastic PArallel Rarified-gas Time-accurate Analyzer (SPARTA), a direct-simulation Monte Carlo (DSMC) code, is modified to allow oxidation-driven ablation of implicitly defined carbon surfaces. In SPARTA, implicit surfaces are generated from the grid corner point values via a marching cubes algorithm, therefore creating a new set of surface elements every time ablation is performed. The finite-rate oxidation model developed by Gopalan et. al can perform both gas-surface and pure-surface reactions and is now adapted to tally surface data on a per grid cell basis. The ablation functionality was also adjusted so once the reactions have occurred, the number of reactions leading to CO formation can be converted to corner point reduction values; therefore, carbon removal is directly proportional to surface recession. We also briefly discuss some unique challenges associated with parallelizing this dynamic surface state and geometry. Finally, we analyze the performance of this parallelized implicit chemistry model with simple 2D and 3D benchmark cases by producing surface state statistics, area changes over time, and visualization across a range of surface temperatures and processors with and without load-balancing.

DSMC↗
NASA NTRS1998

Estimation of 3-D Cloud Effects on TOMS Satellite Retrieval of Surface UV Irradiance

To improve surface UV irradiance retrieval from the Total Ozone Mapping Spectrometer (TOMS) we simulate errors of the TOMS cloud correction algorithm for summertime broken cloud conditions. Cloud scenes (50 km by 50 km) are modeled by a normal random (Gaussian) field with a fixed lower boundary and conservative scattering. The model relates stochastic field characteristics with the cloud amount, mean cloud diameter and aspect ratio. Clouds are embedded into Rayleigh atmosphere with standard ozone profile. Radiative transfer calculations of the radiance at the top of the atmosphere and irradiance at the surface were performed using 3-D Monte Carlo (MC) code. The results are averaged over the satellite field of view on the surface (50 km by 50 km) and compared with TOMS predicted surface irradiance for the same scene reflectance. The TOMS algorithm assumes horizontally homogeneous Cl-type cloud between 3 km and 5.5 km. The effective optical depth is determined by fitting observed (MC) radiance at 380 nm. Having the same radiance at the satellite the homogeneous and broken cloud models predict different average irradiances at the surface. This is due to the differences in Bidirectional Reflection Distribution Function (BRDF) for homogeneous and broken cloud scenes with the same hemispherical albedo. For typical TOMS observational geometry at mid-latitudes the simulated single pixels errors may be as large as +/- 20%. Qualitatively these errors are due to the dominance of the non-horizontal cloud surfaces, which are not accounted for in the homogeneous cloud model. However, due to high variability of the real cloud shapes and types it is unclear how these single pixel errors would affect TOMS time-integrated UV exposure over extended periods (weeks to months) for different regions.

Krotkov, Nickolay A.↗
NASA NTRS2022

Coupling Carbon Oxidation and Surface Recession in Direct-Simulation Monte Carlo Code, SPARTA

Ablative thermal protection system (TPS) materials for spacecraft are composites that are often made out of carbon-based reinforcement and a polymeric matrix. They endure high-temperature oxidation and surface recession when re-entering Earth’s atmosphere. Ablation is the result of many coupled and competing thermal, mechanical, and chemical phenomena, and it is difficult to isolate the role of each on the overall degradation of the TPS. Here we develop an ablation model for material recession coupled explicitly to finite rate carbon oxidation in complex microstructures. In this work, Stochastic PArallel Rarified-gas Time-accurate Analyzer (SPARTA), a direct-simulation Monte Carlo (DSMC) code, is modified to allow oxidation-driven ablation of implicitly defined carbon surfaces. In SPARTA, implicit surfaces are generated from the grid corner point values via a marching cubes algorithm, therefore creating a new set of surface elements every time ablation is performed. The finite-rate oxidation model developed by Gopalan et. al, was adapted to tally surface reactions and other surface data on a per-grid cell basis. The ablation functionality was also adjusted so once the reactions have occurred, the number of reactions leading to CO formation can be converted to corner point reduction values; therefore, carbon removal is directly proportional to surface recession. We also develop robust algorithms which handle the evolution of the flow cells and solid material regions, including split cells (flow cell divided in two by a solid surface). Finally, we demonstrate our implicit chemistry model for 2D and 3D geometries by producing reaction statistics and detailed visualization of oxidation-induced material recession at the microscale.

V Arias↗
NASA NTRS2017

Reconfigurable Cellular Composite Structures for Lighter than Air Vehicles with Scalable Size and Endurance

Engineered non-stochastic cellular materials show promising characteristics on the laboratory scale,with nearly ideal specific stiffness and strength scaling at ultralight mass density. These propertiessuggest performance benefits in any application with combined stiffness and mass constraints, suchas air vehicles. We investigate here the application of re-configurable cellular composite materialsand structures to lighter than air vehicles. We describe the properties and applicability of these materials,provide an example analysis of governing loading conditions associated with airships, showan example optimization method for navigating the design space, and describe how recent advancesin cellular material manufacturing and reconfiguration enable system performance benefits includingnew concepts of operation. Lastly, we propose lighter than air vehicles that are assembled andmaintained in-flight, eliminating structural compromises associated with transitional flight modesand ground handling.Engineered non-stochastic cellular material properties suggest performance benefits in lighter than air vehicles due tostiffness and mass constraints that are intrinsic to the airship design problem. Recent advances in cellular materialmanufacturing and reconfiguration enable system performance benefits including new concepts of operation, such aslighter than air vehicles that are assembled and maintained in-flight, eliminating structural compromises associatedwith transitional flight modes and ground handling. Existing engineered cellular materials display properties allowinglarge large scale airships design as monocoque cellular solids. Inevitable improvements in cellular material propertiesand manufacturing will improve feasibility even further. Given the suggestion that the two most significant technologygaps exist across all current airship projects are manufacturing and assembly processes and ground handling [7],a strategy that encompasses construction and maintenance in flight could provide critical rephrasing of the systemdesign problem through these new concepts of operation. Refactoring of traditional manufacturing, operation, andservice process constraints could extend to other domains in aerospace systems and manufacturing in general.In future work, the complexity of the design task would benefit from a form of optimization in order to find themost suitable geometry for a chosen application. For example, the Sequential Least SQuares Programming (SLSQP)function from within the SciPy Minimize library is a multiobjective constrained optimization method that has beenapplied to fixed wing aircraft design. [17] In this situation it would allow for several objective functions such as drag,bending stiffness, buoyancy and cost of transport to be incorporated into a composite objective function.

Copplestone, Grace↗
NASA NTRS2003

A Eulerian-Lagrangian Model to Simulate Two-Phase/Particulate Flows

Figure 1 shows a snapshot of liquid fuel spray coming out of an injector nozzle in a realistic gas-turbine combustor. Here the spray atomization was simulated using a stochastic secondary breakup model (Apte et al. 2003a) with point-particle approximation for the droplets. Very close to the injector, it is observed that the spray density is large and the droplets cannot be treated as point-particles. The volume displaced by the liquid in this region is significant and can alter the gas-phase ow and spray evolution. In order to address this issue, one can compute the dense spray regime by an Eulerian-Lagrangian technique using advanced interface tracking/level-set methods (Sussman et al. 1994; Tryggvason et al. 2001; Herrmann 2003). This, however, is computationally intensive and may not be viable in realistic complex configurations. We therefore plan to develop a methodology based on Eulerian-Lagrangian technique which will allow us to capture the essential features of primary atomization using models to capture interactions between the fluid and droplets and which can be directly applied to the standard atomization models used in practice. The numerical scheme for unstructured grids developed by Mahesh et al. (2003) for incompressible flows is modified to take into account the droplet volume fraction. The numerical framework is directly applicable to realistic combustor geometries. Our main objectives in this work are: Develop a numerical formulation based on Eulerian-Lagrangian techniques with models for interaction terms between the fluid and particles to capture the Kelvin- Helmholtz type instabilities observed during primary atomization. Validate this technique for various two-phase and particulate flows. Assess its applicability to capture primary atomization of liquid jets in conjunction with secondary atomization models.

Apte, S. V.↗
NASA NTRS2019

Testbed Results for Scalar and Vector Radiative Transfer Computations of Light in Atmosphere-Ocean Systems

We generate and tabulate reflectance values of the Stokes parameters I, Q, and U of upwelling radiance just above a rough ocean surface and at the top of the atmosphere (TOA) for 100 scattering geometries, four atmosphere-ocean systems, and four wavelengths. The atmosphere-ocean systems increase in complexity from (a) a molecular atmosphere above a rough ocean surface (AOS-I model); to (b) a pure water body below a rough ocean surface (AOS-II model); to (c) a fully-coupled simple atmosphere-ocean system (AOS-III model) containing a molecular atmosphere, rough ocean surface, and pure water; to (d) a fully-coupled complex atmosphere-ocean system (AOS-IV model) that includes scattering by molecules, rough ocean surface, pure water, and hydrosols. Our wavelengths (350, 450, 550, and 650 nm) capture the ultraviolet-visible range. Our tables provide radiative transfer (RT) testbed results for atmosphere-ocean systems with an accuracy that surpasses the measurement accuracy of state-of-the-art polarimeters. To validate the accuracy of these tables we performed computations using three independent RT codes that provide deterministic numerical solutions for the RT equation. The agreement is 10(exp –5) for AOS-IV model, and 10(exp –6) for the other models. The degree of linear polarization computed by these RT codes differs by ≤0.2% for 15 isolated cases of tabulated reflectance values, and by ≤0.1% for all remaining cases. We also provide comparisons with results obtained by a stochastic RT code for AOS-I model. The agreement between the deterministic and stochastic results for this model is 10(exp –5) at TOA, and 10(exp –6) above the ocean surface.

Polarimeters↗