Search NASA⌕ Search

SEARCH · Search NASA

Results for “Slip localization”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2005

Development of a Real-Time GPS/Seismic Displacement Meter: Applications to Civilian Infrastructure in Orange and Western Riverside Counties, California

We propose a three-year applications project that will develop an Integrated Real-Time GPS/Seismic System and deploy it in Orange and Western Riverside Counties, spanning three major strike-slip faults in southern California (San Andreas, San Jacinto, and Elsinore) and significant populations and civilian infrastructure. The system relying on existing GPS and seismic networks will collect and analyze GPS and seismic data for the purpose of estimating and disseminating real-time positions and total ground displacements (dynamic, as well as static) during all phases of the seismic cycle, from fractions of seconds to years. Besides its intrinsic scientific use as a real-time displacement meter (transducer), the GPS/Seismic System will be a powerful tool for local and state decision makers for risk mitigation, disaster management, and structural monitoring (dams, bridges, and buildings). Furthermore, the GPS/Seismic System will become an integral part of California's spatial referencing and positioning infrastructure, which is complicated by tectonic motion, seismic displacements, and land subsidence. Finally, the GPS/Seismic system will also be applicable to navigation in any environment (land, sea, or air) by combining precise real-time instantaneous GPS positioning with inertial navigation systems. This development will take place under the umbrella of the California Spatial Reference Center, in partnership with local (Counties, Riverside County Flood and Water Conservation District, Metropolitan Water District), state (Caltrans), and Federal agencies (NGS, NASA, USGS), the geophysics community (SCIGN/SCEC2), and the private sector (RBF Consulting). The project will leverage considerable funding, resources, and R&D from SCIGN, CSRC and two NSF-funded IT projects at UCSD and SDSU: RoadNet (Real-Time Observatories, Applications and Data Management Network) and the High Performance Wireless Research and Education Network (HPWREN). These two projects are funded to develop both the wireless networks and the integrated, seamless, and transparent information management system that will deliver seismic, geodetic, oceanographic, hydrological, ecological, and physical data to a variety of end users in real-time in the San Diego region. CSRC is interested in providing users access to real-time, accurate GPS data for a wide variety of applications including RTK surveying/GIS and positioning of moving platforms such as aircraft and emergency vehicles. SCIGN is interested in upgrading sites to high-frequency real-time operations for rapid earthquake response and GPS seismology. The successful outcome of the project will allow the implementation of similar systems elsewhere, particularly in plate boundary zones with significant populations and civilian infrastructure. CSRC would like to deploy the GPS/Seismic System in other parts of California, in particular San Diego, Los Angeles County and the San Francisco Bay Area.

Bock, Yehuda↗
NASA NTRS1994

Vortex distortion during vortex-surface interaction in a Mach 3 stream

An experimental study has been conducted in Mach 3 wind tunnel to investigate the behavior of supersonic vortices as they interact with a wedge surface placed in their passage. The experimental setup was arranged so that interactions resulted in a close encounter of the vortex core and the wedge leading edge. Spark shadow photographs of the flowfield along the pressure measurements on the wedge surface were used to study the interaction problem. In their most organized form, distortion of streamwise vortices upon interacting with the wedge was found to result in formation of symmetric detached shock fronts far upstream of the wedge leading edge followed by an apparent slip surface separating a subsonic region from a supersonic zone. Interaction experiments leading to substantial changes in the structure of vortices revealed that the supersonic vortex distortion has strong resemblances to the incompressible 'B-breakdown' reported in the literature. Experimental results also indicate that the interaction strongly depends on the vortex strength and vortex proximity to the wedge leading edge, and the generated flowfield was found to be highly unsteady. Interaction of concentrated streamwise vortices with the oblique shock formed over the wedge surface resulted in formation of a locally three-dimensional shock wave with a limited subsonic region downstream of the shock.

Kalkhoran, Iraj M.↗
NASA NTRS1978

Evaluation of a locally homogeneous model of spray evaporation

Measurements were conducted on an evaporating spray in a stagnant environment. The spray was formed using an air-atomizing injector to yield a Sauter mean diameter of the order of 30 microns. The region where evaporation occurred extended approximately 1 m from the injector for the test conditions. Profiles of mean velocity, temperature, composition, and drop size distribution, as well as velocity fluctuations and Reynolds stress, were measured. The results are compared with a locally homogeneous two-phase flow model which implies no velocity difference and thermodynamic equilibrium between the phases. The flow was represented by a k-epsilon-g turbulence model employing a clipped Gaussian probability density function for mixture fraction fluctuations. The model provides a good representation of earlier single-phase jet measurements, but generally overestimates the rate of development of the spray. Using the model predictions to represent conditions along the centerline of the spray, drop life-history calculations were conducted which indicate that these discrepancies are due to slip and loss of thermodynamic equilibrium between the phases.

Shearer, A. J.↗
NASA NTRS2020

A ROS-based Simulator for Testing the Enhanced Autonomous Navigation of the Mars 2020 Rover

In order to achieve the ambitious objectives of the Mars 2020 (M2020) mission, in particular the ability to autonomously traverse more challenging terrains more efficiently, new surface mobility software was developed for Enhanced Navigation (ENav). That decision was made early in the project, before most of the new surface flight software (FSW) existed, which created a need for a separate framework where the new navigation algorithms could be quickly prototyped and tested, before more realistic FSW-based testbeds became available. The JPL robotics team chose the Robot Operating System [1] (ROS) as the environment in which to test the new ENav algorithms. This made it possible to write the algorithms in the C language required by the FSW, so they could be directly ported over to the flight module later on, while leveraging all the C++ libraries and tools provided by ROS for simulation and testing. The ENav algorithms were developed as a separate C library, and stubs were used to replace any FSW-specific code, such as Event Reporting (EVRs) and data products (DPs). A ROS simulator was developed to generate a rich set of varied 3D terrains representative of the candidate Mars landing sites and simulate the physics of the rover motion, the point cloud perceived by the rover’s stereo vision system, and the new thinking-while-driving (TWD) navigation logic which directs the rover to drive autonomously to user-specified waypoints. To simulate the rover motion and perception, a ROS node was developed that uses a software library called HyperDrive Sim (HDSim), which is a wrapper for the Rover Sequencing and Visualization Program [2] (RSVP). That library provides roverterrain settling, realistic slip modelling, and camera rendering capability based on the rover’s NavCam machine vision models. To simulate the navigation logic, a ROS node was created that initializes and runs the ENav algorithms in a way that mimics the FSW execution, while also providing the capability to load and replay data products, including re-running the recorded inputs through the ENav algorithms for testing. An engineering Graphical User Interface (GUI) was also developed to visualize various elements, such as the rover pose during the drive, the simulated and perceived terrain, the selected local and global paths to the goal, the evaluated candidate paths and the reasons why they were rejected, the keep-in and keep-out zones (KIOZs), etc. Finally, an advanced Monte Carlo (MC) framework that can run many simulations in parallel on the Cloud and automatically generate reports that capture the key ENav performance metrics was developed to evaluate the system in a statisticallymeaningful way. This paper provides an overview of the ROSbased simulator used for testing the M2020 ENav algorithms.

Toupet, Olivier↗
NASA NTRS1993

Tectonic controls on rift basin morphology: Evolution of the northern Malawi (Nyasa) rift

Radiometric (K-Ar and Ar-40/Ar-39) age determinations of volcanic and volcaniclastic rocks, combined with structural, gravity, and seismic reflection data, are used to constrain the age of sedimentary strata contained within the seismically and volcanically active northern Malawi (Nyasa) rift and to characterize changes in basin and flank morphologies with time. Faulting and volcanism within the Tukuyu-Karonga basin began at approximately 8.6 Ma, when sediments were deposited in abroad, initially asymmetric lake basin bounded on its northeastern side by a border fault system with minor topographic relief. Extensions, primarily by a slip along the border fault, and subsequent regional isostatic compensation led to the development of a 5-km-deep basin bounded by broad uplifted flanks. Along the low-relief basin margin opposite border fault, younger stratigraphic sequences commonly onlap older wedge-shaped sequences, although their internal geometry is often progradational. Intrabasinal faulting, flankuplift, and basaltic and felsic volcanism from centers at the northern end of the basin became more important at about 2.5 Ma when cross-rift transfer faults developed to link the Tukuyu-Karonga basin to the Rukwa basin. Local uplift and volcanic construction at the northern end of the basin led to a southeastward shift in the basin's depocenter. Sequence boundaries are commonly erosional along this low-relief (hanging wall) margin and conformable in the deep lake basin. The geometry of stratigraphic sequences and the distribution of the erosion indicate that horizontal and vertical crustal movements both across and along the length of the rift basin led to changes in levels of the lake, irrespective of paleoclimatic fluctuations.

Ebinger, C. J.↗
NASA NTRS1993

Drag reduction at a plane wall

The objective is to determine by analytical means how drag on a plane wall may be modified favorably using a minimal amount of flow information - preferably only information at the wall. What quantities should be measured? How should that information be assimilated in order to arrive at effective control? As a prototypical problem, incompressible, viscous flow, governed by the Navier-Stokes equations, past a plane wall at which the no-slip condition was modified was considered. The streamwise and spanwise velocity components are required to be zero, but the normal component is to be specified according to some control law. The challenge is to choose the wall-normal velocity component based on flow conditions at the wall so that the mean drag is as small as possible. There can be no net mass flux through the wall, and the total available control energy is constrained. A turbulent flow is highly unsteady and has detailed spatial structure. The mean drag on the wall is the integral over the wall of the local shear forces exerted by the fluid, which is then averaged in time; it is a 'macroscopic' property of the flow. It is not obvious how unsteady boundary control is to be applied in order to modify the mean flow most effectively, especially in view of the non- self-adjoint nature of the governing equations. An approximate analytical solution to the suboptimal scheme is pursued.

Hill, D. C.↗
NASA NTRS1989

Strain accumulation and surface deformation along the San Andreas, California

The goal of this project remains to be the achievement of a better understanding of the regional and local deformation and crustal straining processes in western North America, particularly the effect of the San Andreas and nearby faults on the spatial and temporal crustal deformation behavior. Construction of theoretical models based on the mechanics of coupled elastic plate/viscoelastic foundation and large scale crack mechanics provide a rational basis for the interpretation of seismic and aseismic anomalies and expedite efforts in forecasting the stability of plate boundary deformations. In the present period, special focus is placed on the 3-D effect of irregular fault locked patches on the ground measured deformation fields. Specifically, use is made of a newly developed 3-D boundary element program to analyze the fault slip and vertical ground motion in the Parkfield area on the San Andreas.

Li, Victor C.↗
NASA NTRS1994

Modeling of periodic great earthquakes on the San Andreas fault: Effects of nonlinear crustal rheology

We analyze the cycle of great earthquakes along the San Andreas fault with a finite element numerical model of deformation in a crust with a nonlinear viscoelastic rheology. The viscous component of deformation has an effective viscosity that depends exponentially on the inverse absolute temperature and nonlinearity on the shear stress; the elastic deformation is linear. Crustal thickness and temperature are constrained by seismic and heat flow data for California. The models are for anti plane strain in a 25-km-thick crustal layer having a very long, vertical strike-slip fault; the crustal block extends 250 km to either side of the fault. During the earthquake cycle that lasts 160 years, a constant plate velocity v(sub p)/2 = 17.5 mm yr is applied to the base of the crust and to the vertical end of the crustal block 250 km away from the fault. The upper half of the fault is locked during the interseismic period, while its lower half slips at the constant plate velocity. The locked part of the fault is moved abruptly 2.8 m every 160 years to simulate great earthquakes. The results are sensitive to crustal rheology. Models with quartzite-like rheology display profound transient stages in the velocity, displacement, and stress fields. The predicted transient zone extends about 3-4 times the crustal thickness on each side of the fault, significantly wider than the zone of deformation in elastic models. Models with diabase-like rheology behave similarly to elastic models and exhibit no transient stages. The model predictions are compared with geodetic observations of fault-parallel velocities in northern and central California and local rates of shear strain along the San Andreas fault. The observations are best fit by models which are 10-100 times less viscous than a quartzite-like rheology. Since the lower crust in California is composed of intermediate to mafic rocks, the present result suggests that the in situ viscosity of the crustal rock is orders of magnitude less the rock viscosity determined in the laboratory.

Reches, Ze'ev↗
NASA NTRS1995

Hypersonic Aerodynamic Characteristics of a Proposed Single-Stage-To-Orbit Vehicle

The hypersonic aerodynamic characteristics of a winged body concept representing a candidate single- stage-to-orbit vehicle which features wing tip fin controllers and elevon/body flap control surfa'Fs are predicted at points along a nominal trajectory for Mach numbers from 5 to 27 and angles of attack from 19 to 32 degrees. Predictions are derived from surface properties based on flow solvers for inviscid and viscous, laminar flows acting as a perfect gas, as a gas in chemical equilibrium and as a gas in chemical non- equilibrium. At a Mach number of 22, the lateral aerodynamic characteristics of the vehicle are determined based on an inviscid analysis at side slip angles of 2 and 4 degrees and 32 degrees angle of attack; a viscous analysis was carried out to determine the effect of gas chemistry model on surface pressure and to determine the incremental aerodynamics for control surface deflections. The results show that the longitudinal pitch characteristics of the baseline configuration, i.e., zero control surface deflections, are significantly altered by real gas chemistry at angles of attack greater than 30 degrees and Mach numbers greater than 9; and, that aerodynamics derived from inviscid solutions are of sufficient accuracy for preliminary analysis. Also, it is shown that a Mach number of 22, the choice of gas chemistry model has a large impact on surface pressure levels at highly localized regions on the vehicle and that the vehicle can be trimmed at control surface deflections less than 11 degrees.

Weilmuenster, K. James↗
NASA NTRS2005

Gears Based on Carbon Nanotubes

Gears based on carbon nanotubes (see figure) have been proposed as components of an emerging generation of molecular- scale machines and sensors. In comparison with previously proposed nanogears based on diamondoid and fullerene molecules, the nanotube-based gears would have simpler structures and are more likely to be realizable by practical fabrication processes. The impetus for the practical development of carbon-nanotube- based gears arises, in part, from rapid recent progress in the fabrication of carbon nanotubes with prescribed diameters, lengths, chiralities, and numbers of concentric shells. The shafts of the proposed gears would be made from multiwalled carbon nanotubes. The gear teeth would be rigid molecules (typically, benzyne molecules), bonded to the nanotube shafts at atomically precise positions. For fabrication, it may be possible to position the molecular teeth by use of scanning tunneling microscopy (STM) or other related techniques. The capability to position individual organic molecules at room temperature by use of an STM tip has already been demonstrated. Routes to the chemical synthesis of carbon-nanotube-based gears are also under investigation. Chemical and physical aspects of the synthesis of molecular scale gears based on carbon nanotubes and related molecules, and dynamical properties of nanotube- based gears, have been investigated by computational simulations using established methods of quantum chemistry and molecular dynamics. Several particularly interesting and useful conclusions have been drawn from the dynamical simulations performed thus far: The forces acting on the gears would be more sensitive to local molecular motions than to gross mechanical motions of the overall gears. Although no breakage of teeth or of chemical bonds is expected at temperatures up to at least 3,000 K, the gears would not work well at temperatures above a critical range from about 600 to about 1,000 K. Gear temperature could probably be controlled by use of coolant gases. For a given application, the gears would work well at temperatures below the critical range, provided that the rotational energy was less than the energy required to tilt the teeth through an angle of 20 . The predominant mechanism of gear failure would be slippage caused by tilting of teeth. Gears would resume functioning if the slipping gears were decelerated sufficiently.

Jaffe, Richard↗
NASA NTRS1994

Neotectonics of Asia: Thin-shell finite-element models with faults

As India pushed into and beneath the south margin of Asia in Cenozoic time, it added a great volume of crust, which may have been (1) emplaced locally beneath Tibet, (2) distributed as regional crustal thickening of Asia, (3) converted to mantle eclogite by high-pressure metamorphism, or (4) extruded eastward to increase the area of Asia. The amount of eastward extrusion is especially controversial: plane-stress computer models of finite strain in a continuum lithosphere show minimal escape, while laboratory and theoretical plane-strain models of finite strain in a faulted lithosphere show escape as the dominant mode. We suggest computing the present (or neo)tectonics by use of the known fault network and available data on fault activity, geodesy, and stress to select the best model. We apply a new thin-shell method which can represent a faulted lithosphere of realistic rheology on a sphere, and provided predictions of present velocities, fault slip rates, and stresses for various trial rheologies and boundary conditions. To minimize artificial boundaries, the models include all of Asia east of 40 deg E and span 100 deg on the globe. The primary unknowns are the friction coefficient of faults within Asia and the amounts of shear traction applied to Asia in the Himalayan and oceanic subduction zones at its margins. Data on Quaternary fault activity prove to be most useful in rating the models. Best results are obtained with a very low fault friction of 0.085. This major heterogeneity shows that unfaulted continum models cannot be expected to give accurate simulations of the orogeny. But, even with such weak faults, only a fraction of the internal deformation is expressed as fault slip; this means that rigid microplate models cannot represent the kinematics either. A universal feature of the better models is that eastern China and southeast Asia flow rapidly eastward with respect to Siberia. The rate of escape is very sensitive to the level of shear traction in the Pacific subduction zones, which is below 6 MPa. Because this flow occurs across a wide range of latitudes, the net eastward escape is greater than the rate of crustal addition in the Himalaya. The crustal budget is balanced by extension and thinning, primarily within the Tibetan plateau and the Baikal rift. The low level of deviation stresses in the best models suggests that topographic stress plays a major role in the orogeny; thus, we have to expect that different topography in the past may have been linked with fundamentally different modes of continental collision.

Kong, Xianghong↗
NASA NTRS2011

Fracture Mechanics Analyses of Subsurface Defects in Reinforced Carbon-Carbon Joggles Subjected to Thermo-Mechanical Loads

Coating spallation events have been observed along the slip-side joggle region of the Space Shuttle Orbiter wing-leading-edge panels. One potential contributor to the spallation event is a pressure build up within subsurface voids or defects due to volatiles or water vapor entrapped during fabrication, refurbishment, or normal operational use. The influence of entrapped pressure on the thermo-mechanical fracture-mechanics response of reinforced carbon-carbon with subsurface defects is studied. Plane-strain simulations with embedded subsurface defects are performed to characterize the fracture mechanics response for a given defect length when subjected to combined elevated-temperature and subsurface-defect pressure loadings to simulate the unvented defect condition. Various subsurface defect locations of a fixed-length substrate defect are examined for elevated temperature conditions. Fracture mechanics results suggest that entrapped pressure combined with local elevated temperatures have the potential to cause subsurface defect growth and possibly contribute to further material separation or even spallation. For this anomaly to occur, several unusual circumstances would be required making such an outcome unlikely but plausible.

Knight, Norman F., Jr.↗
NASA NTRS1994

Topography and tectonics of the central New Madrid seismic zone: Results of numerical experiements using a three-dimensional boundary element program

We present results of a series of numerical experiments designed to test hypothetical mechanisms that derive deformation in the New Madrid seismic zone. Experiments are constrained by subtle topography and the distribution of seismicity in the region. We use a new boundary element algorithm that permits calcuation of the three-dimensional deformation field. Surface displacement fields are calculated for the New Madrid zone under both far-field (plate tectonics scale) and locally derived driving strains. Results demonstrate that surface displacement fields cannot distinguish between either a far-field simple or pure shear strain field or one that involves a deep shear zone beneath the upper crustal faults. Thus, neither geomorphic nor geodetic studies alone are expected to reveal the ultimate driving mechanism behind the present-day deformation. We have also tested hypotheses about strain accommodation within the New Madrid contractional step-over by including linking faults, two southwest dipping and one vertical, recently inferred from microearthquake data. Only those models with step-over faults are able to predict the observed topography. Surface displacement fields for long-term, relaxed deformation predict the distribution of uplift and subsidence in the contractional step-over remarkably well. Generation of these displacement fields appear to require slip on both the two northeast trending vertical faults and the two dipping faults in the step-over region, with very minor displacements occurring during the interseismic period when the northeast trending vertical faults are locked. These models suggest that the gently dippling central step-over fault is a reverse fault and that the steeper fault, extending to the southeast of the step-over, acts as a normal fault over the long term.

Gomberg, Joan↗
NASA NTRS1946

Altitude-Wind-Tunnel Investigation of Performance of Several Propellers on YP-47M Airplane at High Blade Loading: Curtiss 838-1C2-18R1 Four-Blade Propeller - 2

An investigation was conducted in the Cleveland altitude wind tunnel to determine the performance of a Curtiss propeller with four 838-lC2-lSRl blades on a YP-47M airplane at high blade loadings and engine powers. The study was made for a range of power coefficients between 0.30 and 1.00 at free-stream Mach numbers of 0.40 and 0.50. The results of the force measurements indicate primarily the trend of propeller efficiency for changes in power coefficient or advance-diameter ratio, inasmuch as corrections for the effects of tunnel-wall constriction on the installation have not been applied. Slip-stream pressure surveys across the propeller disk are presented to illustrate blade thrust load distribution for several operating conditions. At a free-stream Mach number of 0.40, nearly constant peak efficiencies were obtained at power coefficients from 0.30 to 0.70. A change in power coefficient from 0.70 to 0.90 reduced the peak efficiency about 5 percent. Blade stall at the tip sections became evident for a power coefficient of 0.91 when the advance-diameter ratio was reduced to 1.87. At a free-stream Mach number of 0.50, the highest propeller efficiencies were obtained for power coefficients from 0.80 to 1.00 at advance-diameter ratios above 2.90. At advance-diameter ratios below 2.90, the highest efficiencies were obtained for power coefficients of 0.60 and 0.70. The envelope of the efficiency curves decreased about 12 percent between advance-diameter ratios of 2.60 and 4.20. Local compressibility effects became evident for a power coefficient of 0.40 when the advance-diameter ratio was decreased to 1.75.

Wallner, Lewis E.↗
NASA NTRS1989

Eastern Ishtar Terra: Tectonic evolution derived from recognized features

Previous analyses have recognized several styles and orientations of compressional deformation, crustal convergence, and crustal thickening in Eastern Ishtar Terra. An east to west sense of crustal convergence through small scale folding, thrusting, and buckling is reflected in the high topography and ridge-and-valley morphology of Maxwell Montes and the adjacent portion of Fortuna Tessera. This east to west convergence was accompanied by up to 1000 km of lateral motion and large scale strike-slip faulting within two converging shear zones which has resulted in the present morphology of Maxwell Montes. A more northeast to southwest sense of convergence through large scale buckling and imbrication is reflected in large, northwest-trending scarps along the entire northern boundary of Ishtar Terra, with up to 2 km of relief present at many of the scarps. It was previously suggested that both styles of compression have occurred at the expense of pre-existing tessera regions which have then been overprinted by the latest convergence event. The difference in style is attributed mostly to differences in the properties of the crust converging with the tessera blocks. If one, presumably thick, tessera block converges with another tessera region, then the widespread, distributed style of deformation occurs, as observed in western Fortuna Tessera. However, if relatively thin crust (such as suggested for the North Polar Plains converges with thicker tessera regions, then localized deformation occurs, as reflected in the scarps along Northern Ishtar Terra. The purpose is to identify the types of features observed in Eastern Ishtar Terra. Their potential temporal and spatial relationships, is described, possible origins for them is suggested, and how the interpretation of some of these features has led to the multiple-style tectonic evolution model described is shown.

Vorderbruegge, R. W.↗
NASA NTRS1993

Photogeological analysis of Europan tectonic features

Preliminary photogeological analyses of the Pelorus Linea and Sidon Flexus regions of Europa were conducted to explore the proposal by Schenk that lateral motion of crustal blocks has occurred in a 'rift zone' including possible strike-slip, tension fracturing, and geometric plate rotation about an Euler pole. These analyses revealed features interpreted as tensional structures and block rotation in a strike-slip regime consistent with the Schenk hypotheses and implied the presence of at least two stages of crustal deformation consistent with a chronology developed by Lucchitta. Confirmation of regional scale Euler pole rotation was ambiguous, however. Up to 80 kilometers of possible extension was identified in the rift zone; to accommodate this, 'cryosubduction' is speculatively proposed as a mechanism for recycling Europan 'ice lithosphere'. The cumulative width of wedge-shaped bands included in the rift zone was measured and plotted versus distance from the inferred rotation pole. Three sharp decreases in the total width were noted. These occur roughly where certain triple bands cross the rift zone suggesting that the bands are structural features that predate and influence the zone. While the curve hints at one or more sinusoidal relationships consistent with rotation geometry, given the low photographic resolution and the preliminary nature of this examination the question of whether the observations represent coherent regional rotation modified by crosscutting structures or instead imply independent local rotations separated by these structures is unanswered by this analysis.

Tufts, B. R.↗
NASA NTRS2015

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL Multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. This work builds upon previous efforts to verify the use of the acoustic velocity potential equation (AVPE) laid out by Campos. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, del squared psi - (lambda/c) squared psi - M x [M x del((del)(psi))] - 2((lambda)(M)/c + M x del(M) x (del)(psi) - 2(lambda)(psi)[M x del(1/c)] = 0. with M as the Mach vector, c as the speed of sound, and λ as the complex eigenvalue. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, axial velocity, exit velocity). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The one way coupled analysis is perform four times utilizing geometries determined through traditional SRM modeling procedures. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, S. R.↗
NASA NTRSDate not supplied

Uncertainty Reduction With Multi-Model Monte Carlo for Crystal Plasticity Simulations of Additively Manufactured Metals

In this work, multi-model Monte Carlo estimators are developed to reduce uncertainty in quantities of interest (QoIs) extracted from crystal plasticity simulations of additively manufactured (AM) metals. A significant concern in AM parts is uncertainty in mechanical properties caused in part by complex microstructures that arise from the AM process. Quantifying uncertainty in microstructure-sensitive behavior using experiments alone is costly, especially when mechanical allowables must be established. Quantitative relationships among microstructure, micromechanical metrics like slip accumulation, crack initiation, and failure are also difficult to capture with limited experiments. Crystal plasticity material models instead enable computational prediction of micromechanical stress and strain fields given a discretized microstructure. However, high-fidelity finely discretized crystal plasticity simulations are computationally expensive, while lower-fidelity models are less accurate and generally biased, making uncertainty quantification and reduction computationally difficult as well. Multi-model Monte Carlo methods leverage correlations between high- and low-fidelity models to produce unbiased estimators for QoIs with reduced uncertainty relative to standard Monte Carlo. Crystal plasticity QoIs considered in this work include yield strength and the mean and extreme values of micromechanical fields that are relevant to crack initiation. Multi-model Monte Carlo estimators are developed for each individual QoI and several groups of QoIs. The results of this work establish relationships among model correlations, sample allocation, and uncertainty reduction for different combinations of QoIs and demonstrate a trend of less uncertainty reduction as QoIs become more sensitive to local microstructure. Limitations from using pilot samples to estimate model covariances and train low-fidelity models are also addressed. The uncertainty reduction achieved by multi-model Monte Carlo is an important step toward using computational mechanics models to predict microstructure-sensitive crack initiation and failure in AM parts.

uncertainty quantification↗