Search NASA⌕ Search

SEARCH · Search NASA

Results for “AXISYMMETRIC DEFORMATION”

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 NTRS1993

Time-Accurate Unsteady Flow Simulations Supporting the SRM T+68-Second Pressure Spike Anomaly Investigation (STS-54B)

Time-accurate unsteady flow simulations are being performed supporting the SRM T+68sec pressure 'spike' anomaly investigation. The anomaly occurred in the RH SRM during the STS-54 flight (STS-54B) but not in the LH SRM (STS-54A) causing a momentary thrust mismatch approaching the allowable limit at that time into the flight. Full-motor internal flow simulations using the USA-2D axisymmetric code are in progress for the nominal propellant burn-back geometry and flow conditions at T+68-sec--Pc = 630 psi, gamma = 1.1381, T(sub c) = 6200 R, perfect gas without aluminum particulate. In a cooperative effort with other investigation team members, CFD-derived pressure loading on the NBR and castable inhibitors was used iteratively to obtain nominal deformed geometry of each inhibitor, and the deformed (bent back) inhibitor geometry was entered into this model. Deformed geometry was computed using structural finite-element models. A solution for the unsteady flow has been obtained for the nominal flow conditions (existing prior to the occurrence of the anomaly) showing sustained standing pressure oscillations at nominally 14.5 Hz in the motor IL acoustic mode that flight and static test data confirm to be normally present at this time. Average mass flow discharged from the nozzle was confirmed to be the nominal expected (9550 lbm/sec). The local inlet boundary condition is being perturbed at the location of the presumed reconstructed anomaly as identified by interior ballistics performance specialist team members. A time variation in local mass flow is used to simulate sudden increase in burning area due to localized propellant grain cracks. The solution will proceed to develop a pressure rise (proportional to total mass flow rate change squared). The volume-filling time constant (equivalent to 0.5 Hz) comes into play in shaping the rise rate of the developing pressure 'spike' as it propagates at the speed of sound in both directions to the motor head end and nozzle. The objectives of the present analysis are to: (1) capture the dynamic responses of the motor combustion gas flow to correlate with available low-frequency (less than 12.5 sample/sec) data and (2) observe the high-frequency (up to 50 Hz) characteristics of the response to determine any potentials for dynamic coupling.

Dougherty, N. S.↗
NASA NTRS1976

The axisymmetric jet counterflow problem

A numerical method for computing the flow field of an underexpanded axisymmetric supersonic jet flowing counter to a uniform supersonic stream is described. The strong shock waves and shear layers which occur in the flow are treated as normal and tangential discontinuities bounding regions of inviscid flow. Within the regions, the unsteady gasdynamic equations are solved in integral form on a moving computational mesh. As the numerical flow field evolves, the mesh is deformed according to the jump relations applied at the discontinuities, thus determining their positions as part of the solution. Also presented are computed results obtained for a wide range of jet total pressure. The standoff distances of the bow shock and jet shock are in good agreement with those measured experimentally, thus demonstrating the accuracy of the numerical method.

Schiff, L. B.↗
NASA NTRS1981

Application of the line-spring model to a cylindrical shell containing a circumferential or axial part-through crack

An approximate solution was obtained for a cylindrical shell containing a part-through surface crack. It was assumed that the shell contains a circumferential or axial semi-elliptic internal or external surface crack and was subjected to a uniform membrane loading or a uniform bending moment away from the crack region. A Reissner type theory was used to account for the effects of the transverse shear deformations. The stress intensity factor at the deepest penetration point of the crack was tabulated for bending and membrane loading by varying three dimensionless length parameters of the problem formed from the shell radius, the shell thickness, the crack length, and the crack depth. The upper bounds of the stress intensity factors are provided by the results of the elasticity solution obtained from the axisymmetric crack problem for the circumferential crack, and that found from the plane strain problem for a circular ring having a radial crack for the axial crack. The line-spring model gives the expected results in comparison with the elasticity solutions. Results also compare well with the existing finite element solution of the pressurized cylinder containing an internal semi-elliptic surface crack.

Delale, F.↗
NASA NTRS1984

Theoretical model of impact damage in structural ceramics

This paper presents a mechanistically consistent model of impact damage based on elastic failures due to tensile and shear overloading. An elastic axisymmetric finite element model is used to determine the dynamic stresses generated by a single particle impact. Local failures in a finite element are assumed to occur when the primary/secondary principal stresses or the maximum shear stress reach critical tensile or shear stresses, respectively. The succession of failed elements thus models macrocrack growth. Sliding motions of cracks, which closed during unloading, are resisted by friction and the unrecovered deformation represents the 'plastic deformation' reported in the literature. The predicted ring cracks on the contact surface, as well as the cone cracks, median cracks, radial cracks, lateral cracks, and damage-induced porous zones in the interior of hot-pressed silicon nitride plates, matched those observed experimentally. The finite element model also predicted the uplifting of the free surface surrounding the impact site.

Liaw, B. M.↗
NASA NTRS1984

Stress-derivative control of keystoning deformation in finite element codes

The stress derivative technique for control of keystoning deformation in under-integrated finte elements is based on expansion of the stress in a Taylor series about the element center and retention of additional terms beyond the constant stress term. It has the advantage over other control techniques that keystoning resistance is provided by actual rather than artificial material properties. Application of this technique to the quadrilateral ring elements used for modelling solids of revolution subjected to axisymmetric loads is described. In a cylindrical coordinate system additional terms appear in the formulation which must be dealt with in arriving at a workable keystoning control scheme.

Schulz, J. C.↗
NASA NTRS2015

Modeling the Effects of Spaceflight on the Posterior Eye in VIIP

Purpose: Visual Impairment and Intracranial Pressure (VIIP) syndrome is a new and significant health concern for long-duration space missions. Its etiology is unknown, but is thought to involve elevated intracranial pressure (ICP)that induces connective tissue changes and remodeling in the posterior eye (Alexander et al. 2012). Here we study the acute biomechanical response of the lamina cribrosa (LC) and optic nerve to elevations in ICP utilizing finite element (FE) modeling. Methods: Using the geometry of the posterior eye from previous axisymmetric FE models (Sigal et al. 2004), we added an elongated optic nerve and optic nerve sheath, including the pia and dura. Tissues were modeled as linear elastic solids. Intraocular pressure and central retinal vessel pressures were set at 15 mmHg and 55 mmHg, respectively. ICP varied from 0 mmHg (suitable for standing on earth) to 30 mmHg (representing severe intracranial hypertension, thought to occur in space flight). We focused on strains and deformations in the LC and optic nerve (within 1 mm of the LC) since we hypothesize that they may contribute to vision loss in VIIP. Results: Elevating ICP from 0 to 30 mmHg significantly altered the strain distributions in both the LC and optic nerve (Figure), notably leading to more extreme strain values in both tension and compression. Specifically, the extreme (95th percentile) tensile strains in the LC and optic nerve increased by 2.7- and 3.8-fold, respectively. Similarly, elevation of ICP led to a 2.5- and 3.3-fold increase in extreme (5th percentile) compressive strains in the LC and optic nerve, respectively. Conclusions: The elevated ICP thought to occur during spaceflight leads to large acute changes in the biomechanical environment of the LC and optic nerve, and we hypothesize that such changes can activate mechanosensitive cells and invoke tissue remodeling. These simulations provide a foundation for more comprehensive studies of microgravity effects on human vision, e.g. to guide biological studies in which cells and tissues are mechanically loaded in a ranger elevant for microgravity conditions.

gravitational physiology↗
NASA NTRS1983

Dynamics of Tisserand's frame for an elastic spacecraft with stored angular momentum

Tisserand's frame represents a frame in which the deformations in elastic members and the motions in a central rigid body add zero translational and angular momenta to the overall momentum of the vehicle. The considered relations have been known for some time in the case of nongyroscopic spacecraft. The present investigation is concerned with an extension of the relations to gyroscopic systems, taking into account an entirely elastic spacecraft with significant stored angular momentum. A two-dimensional free, rectangular, elastic structure with axisymmetric rotors is considered. Equations of motion are obtained with the aid of the Newton-Eulerian technique and Hamilton's principle. An analysis is conducted of the dynamics of an elastic gyroscopic free vehicle, and the dynamics of the Tisserand's frame in response to external stimuli is examined.

Hablani, H. B.↗
NASA NTRS1972

Hydrodynamic water impact

The hydrodynamic impact of a falling body upon a viscous incompressible fluid was investigated by numerically solving the equations of motion. Initially the mathematical model simulated the axisymmetric impact of a rigid right circular cylinder upon the initially quiescent free surface of a fluid. A compressible air layer exists between the falling cylinder and the liquid free surface. The mathematical model was developed by applying the Navier-Stokes equations to the incompressible air layer and the incompressible fluid. Assuming the flow to be one dimensional within the air layer, the average velocity, pressure and density distributions were calculated. The liquid free surface was allowed to deform as the air pressure acting on it increases. For the liquid the normalized equations were expressed in two-dimensional cylindrical coordinates. The governing equations for the air layer and the liquid were expressed in finite difference form and solved numerically. For the liquid a modified version of the Marker-and-Cell method was used. The mathematical model has been reexamined and a new approach has recently been initiated. Essentially this consists of examining the impact of an inclined plate onto a quiesent water surface with the equations now formulated in cartesian coordinates.

Kettleborough, C. F.↗
NASA NTRS1993

Application of CFD Analysis to Design Support and Problem Resolution for ASRM and RSRM

The use of Navier-Stokes CFD codes to predict the internal flow field environment in a solid rocket motor is a very important analysis element during the design phase of a motor development program. These computational flow field solutions uncover a variety of potential problems associated with motor performance as well as suggesting solutions to these problems. CFD codes have also proven to be of great benefit in explaining problems associated with operational motors such as in the case of the pressure spike problem with the STS-54B flight motor. This paper presents results from analyses involving both motor design support and problem resolution. The issues discussed include the fluid dynamic/mechanical stress coupling at field joints relative to significant propellant deformations, the prediction of axial and radial pressure gradients in the motor associated with motor performance and propellant mechanical loading, the prediction of transition of the internal flow in the motor associated with erosive burning, the accumulation of slag at the field joints and in the submerged nozzle region, impingement of flow on the nozzle nose, and pressure gradients in the nozzle region of the motor. The analyses presented in this paper have been performed using a two-dimensional axisymmetric model. Fluent/BFC, a three dimensional Navier-Stokes flow field code, has been used to make the numerical calculations. This code utilizes a staggered grid formulation along with the SIMPLER numerical pressure-velocity coupling algorithm. Wall functions are used to represent the character of the viscous sub-layer flow, and an adjusted k-epsilon turbulence model especially configured for mass injection internal flows, is used to model the growth of turbulence in the motor port. Conclusions discussed in this paper consider flow field effects on the forward, center, and aft propellant grains except for the head end star grain region of the forward propellant segment. The field joints and the submerged nozzle are discussed as well. Conclusions relative to both the design evaluation of the ASRM and the RSRM scenarios explaining the pressure spikes were based on the flow field solutions presented in this paper.

Dill, Richard A.↗
NASA NTRS2005

Parametric Investigation of Liquid Jets in Low Gravity

An axisymmetric phase field model is developed and used to model surface tension forces on liquid jets in microgravity. The previous work in this area is reviewed and a baseline drop tower experiment selected for model comparison. This paper uses the model to parametrically investigate the influence of key parameters on the geysers formed by jets in microgravity. Investigation of the contact angle showed the expected trend of increasing contact angle increasing geyser height. Investigation of the tank radius showed some interesting effects and demonstrated the zone of free surface deformation is quite large. Variation of the surface tension with a laminar jet showed clearly the evolution of free surface shape with Weber number. It predicted a breakthrough Weber number of 1.

Chato, David J.↗
NASA NTRS1999

Dynamic Nucleation of Supercooled Melts and Measurement of the Surface Tension and Viscosity

We investigate the phenomenon of acoustic pressure-induced nucleation by using a novel approach involving the large amplitude resonant radial oscillations and collapse of a single bubble intentionally injected into a supercooled liquid. Using a combination of previously developed and proven techniques, the bubble is suspended in a fluid host by an ultrasonic field which supplies both the levitation capability as well as the forcing of the radial oscillations. We observe the effects of an increase in pressure (due to bubble collapse) in a region no larger than 100 microns within the supercooled melt to rigorously probe the hypothesis of pressure-induced nucleation of the solid phase. The use of single bubbles operating in narrow temporal and spatial scales will allow the direct and unambiguous correlation between the origin and location of the generation of the disturbance and the location and timing of the nucleation event. In a companion research effort, we are developing novel techniques for the non-contact measurements of the surface tension and viscosity of highly viscous supercooled liquids. Currently used non-invasive methods of surface tension measurement for the case of undercooled liquids generally rely of the quantitative determination of the resonance frequencies of drop shape oscillations, of the dynamics of surface capillary waves, or of the velocity of streaming flows. These methods become quickly ineffective when the liquid viscosity rises to a significant value. An alternate and accurate method which would be applicable to liquids of significant viscosity is therefore needed. We plan to develop such a capability by measuring the equilibrium shape of levitated undercooled melt droplets as they undergo solid-body rotation. The experimental measurement of the characteristic point of transition (bifurcation point) between axisymmetric and two-lobed shapes will be used to calculate the surface tension of the liquid. Such an approach has already been validated through the experimental verification of numerical modeling results. The experimental approach involves levitation, melting, and solidification of undercooled droplets using a hybrid ultrasonic-electrostatic technique in both a gaseous as well as a vacuum environment. A shape relaxation method will be investigated in order to derive a reliable method to measure the viscosity of undercooled melts. The analysis of the monotonic relaxation to equilibrium shape of a drastically deformed and super-critically damped free drop has been used to derive interfacial tension of immiscible liquid combinations where one of the component has high viscosity. A standard approach uses the initial elongation of a droplet through shear flows, but an equivalent method could involve the initial deformation of a drop levitated in a gas by ultrasonic radiation pressure, electric stresses, or even solid body rotation. The dynamic behavior of the free drop relaxing back to equilibrium shape will be modeled, and its characteristic time dependence should provide a quantitative means to evaluate the liquid viscosity.

Trinh, E. H.↗
NASA NTRS1993

A finite-element-based perturbation model for the rotordynamic analysis of shrouded pump impellers: Part 2: User's guide

This report describes the computational steps involved in executing a finite-element-based perturbation model for computing the rotor dynamic coefficients of a shrouded pump impeller or a simple seal. These arise from the fluid/rotor interaction in the clearance gap. In addition to the sample cases, the computational procedure also applies to a separate category of problems referred to as the 'seal-like' category. The problem, in this case, concerns a shrouded impeller, with the exception that the secondary, or leakage, passage is totally isolated from the primary-flow passage. The difference between this and the pump problem is that the former is analytically of the simple 'seal-like' configuration, with two (inlet and exit) flow-permeable stations, while the latter constitutes a double-entry / double-discharge flow problem. In all cases, the problem is that of a rotor clearance gap. The problem here is that of a rotor excitation in the form of a cylindrical whirl around the housing centerline for a smooth annular seal. In its centered operation mode, the rotor is assumed to give rise to an axisymmetric flow field in the clearance gap. As a result, problems involving longitudinal or helical grooves, in the rotor or housing surfaces, go beyond the code capabilities. Discarding, for the moment, the pre- and post-processing phases, the bulk of the computational procedure consists of two main steps. The first is aimed at producing the axisymmetric 'zeroth-order' flow solution in the given flow domain. Detailed description of this problem, including the flow-governing equations, turbulence closure, boundary conditions, and the finite-element formulation, was covered by Baskharone and Hensel. The second main step is where the perturbation model is implemented, with the input being the centered-rotor 'zeroth-order' flow solution and a prescribed whirl frequency ratio (whirl frequency divided by the impeller speed). The computational domain, in the latter case, is treated as three dimensional, with the number of computational planes in the circumferential direction being specified a priori. The reader is reminded that the deformations in the finite elements are all infinitesimally small because the rotor eccentricity itself is a virtual displacement. This explains why we have generically termed the perturbation model the 'virtually' deformable finite-element category. The primary outcome of implementing the perturbation model is the tangential and radial components, F(sub theta)(sup *) and F(sub r)(sup *) of the fluid-exerted force on the rotor surface due to the whirling motion. Repetitive execution of the perturbation model subprogram over a sufficient range of whirl frequency ratios, and subsequent interpolation of these fluid forces, using the least-square method, finally enable the user to compute the impeller rotor dynamic coefficients of the fluid/rotor interaction. These are the direct and cross-coupled stiffness, damping, and inertia effects of the fluid/rotor interaction.

Baskharone, Erian A.↗
NASA NTRS1996

Forced Oscillations of Supported Drops

Oscillations of supported liquid drops are the subject of wide scientific interest, with applications in areas as diverse as liquid-liquid extraction, synthesis of ceramic powders, growing of pure crystals in low gravity, and measurement of dynamic surface tension. In this research, axisymmetric forced oscillations of arbitrary amplitude of viscous liquid drops of fixed volume which are pendant from or sessile on a rod with a fixed or moving contact line and surrounded by an inviscid ambient gas are induced by moving the rod in the vertical direction sinusiodally in time. In this paper, a preliminary report is made on the computational analysis of the oscillations of supported drops that have 'clean' interfaces and whose contact lines remain fixed throughout their motions. The relative importance of forcing to damping can be increased by either increasing the amplitude of rod motion A or Reynolds number Re. It is shown that as the ratio of forcing to damping rises, for drops starting from an initial rest state a sharp increase in deformation can occur when they are forced to oscillate in the vicinity of their resonance frequencies, indicating the incipience of hysteresis. However, it is also shown that the existence of a second stable limit cycle and the occurrence of hysteresis can be observed if the drop is subjected to a so-called frequency sweep, where the forcing frequency is first increased and then decreased over a suitable range. Because the change in drop deformation response is abrupt in the vicinity of the forcing frequencies where hysteresis occurs, it should be possible to exploit the phenomenon to accurately measure the viscosity and surface tension of the drop liquid.

Wilkes, Edward D.↗
NASA NTRS1988

Modeling of combustion processes of stick propellants via combined Eulerian-Lagrangian approach

This research is motivated by the improved ballistic performance of large-caliber guns using stick propellant charges. A comprehensive theoretical model for predicting the flame spreading, combustion, and grain deformation phenomena of long, unslotted stick propellants is presented. The formulation is based upon a combined Eulerian-Lagrangian approach to simulate special characteristics of the two phase combustion process in a cartridge loaded with a bundle of sticks. The model considers five separate regions consisting of the internal perforation, the solid phase, the external interstitial gas phase, and two lumped parameter regions at either end of the stick bundle. For the external gas phase region, a set of transient one-dimensional fluid-dynamic equations using the Eulerian approach is obtained; governing equations for the stick propellants are formulated using the Lagrangian approach. The motion of a representative stick is derived by considering the forces acting on the entire propellant stick. The instantaneous temperature and stress fields in the stick propellant are modeled by considering the transient axisymmetric heat conduction equation and dynamic structural analysis.

Kuo, K. K.↗
NASA NTRS1989

The structural response of unsymmetrically laminated composite cylinders

The responses of an unsymmetrically laminated fiber-reinforced composite cylinder to an axial compressive load, a torsional load, and the temperature change associated with cooling from the processing temperature to the service temperature are investigated. These problems are considered axisymmetric and the response is studied in the context of linear elastic material behavior and geometrically linear kinematics. Four different laminates are studied: a general unsymmetric laminate; two unsymmetric but more conventional laminates; and a conventional quasi-isotropic symmetric laminate. The responses based on closed-form solutions for different boundary conditions are computed and studied in detail. Particular emphasis is directed at understanding the influence of elastic couplings in the laminates. The influence of coupling decreased from a large effect in the general unsymmetric laminate, to practically no effect in the quasi-isotropic laminate. For example, the torsional loading of the general unsymmetric laminate resulted in a radial displacement. The temperature change also caused a significant radial displacement to occur near the ends of the cylinder. On the other hand, the more conventional unsymmetric laminate and the quasi-isotropic cylinder did not deform radially when subjected to a torsional load. From the results obtained, it is clear the degree of elastic coupling can be controlled and indeed designed into a cylinder, the degree and character of the coupling being dictated by the application.

Butler, T. A.↗
NASA NTRS1999

Waves in Radial Gravity Using Magnetic Fluid

Terrestrial laboratory experiments studying various fluid dynamical processes are constrained, by being in an Earth laboratory, to have a gravitational body force which is uniform and unidirectional. Therefore fluid free-surfaces are horizontal and flat. Such free surfaces must have a vertical solid boundary to keep the fluid from spreading horizontally along a gravitational potential surface. In atmospheric, oceanic, or stellar fluid flows that have a horizontal scale of about one-tenth the body radius or larger, sphericity is important in the dynamics. Further, fluids in spherical geometry can cover an entire domain without any sidewall effects, i.e. have truly periodic boundary conditions. We describe spherical body-force laboratory experiments using ferrofluid. Ferrofluids are dilute suspensions of magnetic dipoles, for example magnetite particles of order 10 nm diameter, suspended in a carrier fluid. Ferrofluids are subject to an additional body force in the presence of an applied magnetic field gradient. We use this body force to conduct laboratory experiments in spherical geometry. The present study is a laboratory technique improvement. The apparatus is cylindrically axisymmetric. A cylindrical ceramic magnet is embedded in a smooth, solid, spherical PVC ball. The geopotential field and its gradient, the body force, were made nearly spherical by careful choice of magnet height-to-diameter ratio and magnet size relative to the PVC ball size. Terrestrial gravity is eliminated from the dynamics by immersing the "planet" and its ferrofluid "ocean" in an immiscible silicone oil/freon mixture of the same density. Thus the earth gravity is removed from the dynamics of the ferrofluid/oil interface and the only dynamically active force there is the radial magnetic gravity. The entire apparatus can rotate, and waves are forced on the ferrofluid surface by exterior magnets. The biggest improvement in technique is in the wave visualization. Fluorescing dye is added to the oil/freon mixture and an argon ion laser generates a horizontal light that can be scanned vertically. Viewed from above, the experiment is a black circle with wave deformations surrounded by a light background. A contour of the image intensity at any light sheet position gives the surface of the ferrofluid "ocean" at that "latitude". Radial displacements of the waves as a function of longitude are obtained by subtracting the contour line positions from a no-motion contour at that laser sheet latitude. The experiments are run by traversing the forcing magnet with the laser sheet height fixed and images are frame grabbed to obtain a time-series at one latitude. The experiment is then re-run with another laser-sheet height to generate a full picture of the three-dimensional wave structure in the upper hemisphere of the ball as a function of time. We concentrate here on results of laboratory studies of waves that are important in Earth's atmosphere and especially the ocean. To get oceanic scaling in the laboratory, the experiment must rotate rapidly (4-second rotation period) so that the wave speed is slow compared to the planetary rotation speed as in the ocean. In the Pacific Ocean, eastward propagating Kelvin waves eventually run into the South American coast. Theory predicts that some of the wave energy should scatter into coastal-trapped Kelvin waves that propagate north and south along the coast. Some of this coastal wave energy might then scatter into mid-latitude Rossby waves that propagate back westward. Satellite observations of the Pacific Ocean sea-surface temperature and height seem to show signatures of westward propagating mid-latitude Rossby waves, 5 to 10 years after the 1982-83 El Nino. The observational data is difficult to interpret unambiguously owing to the large range of motions that fill the ocean at shorter timescales. This series of reflections giving eastward, north- ward, and then westward traveling waves is observed cleanly in the laboratory experiments, confirming the theoretical expectations.

Ohlsen, D. R.↗
NASA NTRS1996

LES on unstructured deforming meshes: Towards reciprocating IC engines

A variable explicit/implicit characteristics-based advection scheme that is second-order accurate in space and time has been developed recently for unstructured deforming meshes (O'Rourke & Sahota 1996a). To explore the suitability of this methodology for Large-Eddy Simulation (LES), three subgrid-scale turbulence models have been implemented in the CHAD CFD code (O'Rourke & Sahota 1996b): a constant-coefficient Smagorinsky model, a dynamic Smagorinsky model for flows having one or more directions of statistical homogeneity, and a Lagrangian dynamic Smagorinsky model for flows having no spatial or temporal homogeneity (Meneveau et al. 1996). Computations have been made for three canonical flows, progressing towards the intended application of in-cylinder flow in a reciprocating engine. Grid sizes were selected to be comparable to the coarsest meshes used in earlier spectral LES studies. Quantitative results are reported for decaying homogeneous isotropic turbulence, and for a planar channel flow. Computations are compared to experimental measurements, to Direct-Numerical Simulation (DNS) data, and to Rapid-Distortion Theory (RDT) where appropriate. Generally satisfactory evolution of first and second moments is found on these coarse meshes; deviations are attributed to insufficient mesh resolution. Issues include mesh resolution and computational requirements for a specified level of accuracy, analytic characterization of the filtering implied by the numerical method, wall treatment, and inflow boundary conditions. To resolve these issues, finer-mesh simulations and computations of a simplified axisymmetric reciprocating piston-cylinder assembly are in progress.

Haworth, D. C.↗
NASA NTRS1993

Patterns of deformation and volcanic flows associated with lithospheric loading by large volcanoes on Venus

Magellan radar imaging and topography data are now available for a number of volcanoes on Venus greater than 100 km in radius. These data can be examined to reveal evidence of the flexural response of the lithosphere to the volcanic load. On Venus, erosion and sediment deposition are negligible, so tectonic evidence of deformation around large volcanoes should be evident except where buried by very young flows. Radar images reveal that most tectonic features and flow units on the flanks of these volcanoes have predominantly radial orientations. However, both Tepev Mons in Bell Regio and Sapas Mons in Atla Regio exhibit circumferential graben on their flanks. In addition, images reveal several flow units with an annular character around the north and west flanks of Tepev Mons. This pattern most likely results from ponding of flows in an annular flexural moat. Maat Mons in Atla Regio and Sif Mons in Eistla Regio are examples of volcanoes that lack circumferential graben and annular flows; discernible flow units and fractures on these constructs appear to be predominantly radial. Altimetry data can also provide evidence of flexural response. Tepev Mons is partially encircled by depressions that may be sections of a flexural moat that has not been completely filled. The locations of these depressions generally coincide with the annular flows described above. There is weaker evidence for such depressions around Maat Mons as well. The lack of circumferential tectonic features around most volcanoes on Venus might be explained by gradual moat filling and coverage by radial flows. The depressions around Tepev (and possible Maat) may indicate that this process is currently continuing. We use analytic models of plate flexure in an axisymmetric geometry to constrain the elastic plate thickness supporting Tepev Mons. If we consider the outer radius of the ponded flows to be the edge of a moat, we find that models with elastic plate thickness of 10-20 km fit best. Finite element models of a volcanic load detached from the underlying lithosphere predict overthrusting and radial normal faulting at the volcano's edge. Such a mechanism for the formation of radial rift zones on Venus volcanoes would make such features analogous to structures on the flanks of volcanoes on Earth.

Mcgovern, Patrick J.↗