Search NASA⌕ Search

SEARCH · Search NASA

Results for “free boundary problem”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Methodology for simulation of unsteady flows to determine the time-dependent interference between stationary and moving boundaries

A new methodology is developed to simulate unsteady flows about prescribed and aerodynamically determined moving boundary problems. The method couples the fluid dynamics and rigid-body dynamics equations to capture the time-dependent interference between stationary and moving boundaries. The unsteady, compressible, inviscid (Euler) equations are solved on dynamic, unstructured grids by an explicit, finite-volume, upwind method. For efficiency, the grid adaptation is performed within a window around the moving object. The Eulerian equations of the rigid-body dynamics are solved by a Runge-Kutta method in a non-inertial frame of reference. The two-dimensional flow solver is validated by computing the flow past a sinusoidally-pitching airfoil and comparing these results with the experimental data. The overall methodology is used for two two-dimensional examples: the flow past an airfoil which is performing a three-degrees-of-freedom motion in a transonic freestream, and the free-fall of a store after separation from a wing-section. Then the unstructured mesh methodology is extended to three-dimensions to simulate unsteady flow past bodies in relative motion, where the trajectory is determined from the instantaneous aerodynamics. The flow solver and the adaptation scheme in three dimensions are validated by simulating the transonic, unsteady flow around a wing undergoing a forced, periodic, pitching motion, and comparing the results with the experimental data. To validate the trajectory code, the six-degrees-of-freedom motion of a store separating from a wing was computed using the experimentally determined force and moment fields, then comparing with an independently generated trajectory. Finally, the overall methodology was demonstrated by simulating the unsteady flowfield and the trajectory of a store dropped from a wing. The methodology, its computational cost notwithstanding, has proven to be accurate, automated, easy for dynamic gridding, and relatively efficient for the required man-hours.

Singh, Kamakhya Prasad↗
NASA NTRS1988

Free vibration of rectangular plates with a small initial curvature

The method of matched asymptotic expansions is used to solve the transverse free vibration of a slightly curved, thin rectangular plate. Analytical results for natural frequencies and mode shapes are presented in the limit when the dimensionless bending rigidity, epsilon, is small compared with in-plane forces. Results for different boundary conditions are obtained when the initial deflection is: (1) a polynomial in both directions, and (2) the product of a polynomial and a trigonometric function, and arbitrary. For the arbitrary initial deflection case, the Fourier series technique is used to define the initial deflection. The results obtained show that the natural frequencies of vibration of slightly curved plates are coincident with those of perfectly flat, prestressed rectangular plates. However, the eigenmodes are very different from those of initially flat prestressed rectangular plates. The total deflection is found to be the sum of the initial deflection, the deflection resulting from the solution of the flat plate problem, and the deflection resulting from the static problem.

Adeniji-Fashola, A. A.↗
NASA NTRS2000

Magnetic Field Suppression of Flow in Semiconductor Melt

One of the most promising approaches for the reduction of convection during the crystal growth of conductive melts (semiconductor crystals) is the application of magnetic fields. Current technology allows the experimentation with very intense static fields (up to 80 KGauss) for which nearly convection free results are expected from simple scaling analysis in stabilized systems (vertical Bridgman method with axial magnetic field). However, controversial experimental results were obtained. The computational methods are, therefore, a fundamental tool in the understanding of the phenomena accounting during the solidification of semiconductor materials. Moreover, effects like the bending of the isomagnetic lines, different aspect ratios and misalignments between the direction of the gravity and magnetic field vectors can not be analyzed with analytical methods. The earliest numerical results showed controversial conclusions and are not able to explain the experimental results. Although the generated flows are extremely low, the computational task is a complicated because of the thin boundary layers. That is one of the reasons for the discrepancy in the results that numerical studies reported. Modeling of these magnetically damped crystal growth experiments requires advanced numerical methods. We used, for comparison, three different approaches to obtain the solution of the problem of thermal convection flows: (1) Spectral method in spectral superelement implementation, (2) Finite element method with regularization for boundary layers, (3) Multiquadric method, a novel method with global radial basis functions, that is proven to have exponential convergence. The results obtained by these three methods are presented for a wide region of Rayleigh and Hartman numbers. Comparison and discussion of accuracy, efficiency, reliability and agreement with experimental results will be presented as well.

Fedoseyev, A. I.↗
NASA NTRS2016

Optical Extinction Measurements of Dust Density in the GMRO Regolith Test Bin

A regolith simulant test bin was constructed and completed in the Granular Mechanics and Regolith Operations (GMRO) Lab in 2013. This Planetary Regolith Test Bed (PRTB) is a 64 sq m x 1 m deep test bin, is housed in a climate-controlled facility, and contains 120 MT of lunar-regolith simulant, called Black Point-1 or BP-1, from Black Point, AZ. One of the current uses of the test bin is to study the effects of difficult lighting and dust conditions on Telerobotic Perception Systems to better assess and refine regolith operations for asteroid, Mars and polar lunar missions. Low illumination and low angle of incidence lighting pose significant problems to computer vision and human perception. Levitated dust on Asteroids interferes with imaging and degrades depth perception. Dust Storms on Mars pose a significant problem. Due to these factors, the likely performance of telerobotics is poorly understood for future missions. Current space telerobotic systems are only operated in bright lighting and dust-free conditions. This technology development testing will identify: (1) the impact of degraded lighting and environmental dust on computer vision and operator perception, (2) potential methods and procedures for mitigating these impacts, (3) requirements for telerobotic perception systems for asteroid capture, Mars dust storms and lunar regolith ISRU missions. In order to solve some of the Telerobotic Perception system problems, a plume erosion sensor (PES) was developed in the Lunar Regolith Simulant Bin (LRSB), containing 2 MT of JSC-1a lunar simulant. PES is simply a laser and digital camera with a white target. Two modes of operation have been investigated: (1) single laser spot - the brightness of the spot is dependent on the optical extinction due to dust and is thus an indirect measure of particle number density, and (2) side-scatter - the camera images the laser from the side, showing beam entrance into the dust cloud and the boundary between dust and void. Both methods must assume a mean particle size in order to extract a number density. The optical extinction measurement yields the product of the 2nd moment of the particle size distribution and the extinction efficiency Qe. For particle sizes in the range of interest (greater than 1 micrometer), Qe approximately equal to 2. Scaling up of the PES single laser and camera system is underway in the PRTB, where an array of lasers penetrate a con-trolled dust cloud, illuminating multiple targets. Using high speed HD GoPro video cameras, the evolution of the dust cloud and particle size density can be studied in detail.

lunar simulant↗
NASA NTRS1994

Adaptive-wall wind-tunnel research at Ames Research Center: A retrospective

This paper reviews adaptive-wall wind-tunnel research conducted at Ames Research Center between 1978 and 1988. This research focused on developing ways to apply the concept of adaptive walls in transonic test sections with ventilated walls. In the approach pursued at Ames, local mass flow through slotted test-section walls was controlled by adjusting the pressure in compartments of a segmented plenum; the flow measurements required to test compatibility of the wind-tunnel flow with free-air boundary conditions were made using a laser velocimeter; empirical influence coefficients were used to predict how pressure changes in the plenum would effect the flow. Both two- and three-dimensional proof-of-concept experiments were conducted in a small indraft wind tunnel. Subsequently, a two-dimensional test section was demonstrated in the Ames 2x2 ft. Transonic Wind Tunnel-a wind tunnel representative of a production facility. The 2x2 ft.tests showed large reductions of wall interference; however, the time required to make the necessary flow measurements and plenum pressure adjustments far exceeded what would be acceptable for production testing. This was primarily because the imaginary surface where free-air compatibility was tested had to be separated from the test section walls to avoid the complex, viscous flow adjacent to the walls. It is unlikely that the Ames approach to adaptive walls will be applied in production wind tunnels. The two biggest unresolved problems are how to quickly and accurately make the necessary flow measurements and how to predict the effects of wall adjustments. In contrast, flexible-wall technology for two-dimensional testing is ready for application in production wind tunnels. However, neither the ventilated- nor flexible-wall approach has been shown to be a technically viable and cost-effective solution to the three-dimensional adaptive wall problem.

Schairer, Edward T.↗
NASA NTRS1972

An analysis of the relaxation of laminar boundary layer on a flat plate after passage of an interface with application to expansion-tube flows

The relaxation of the accelerating-gas boundary layer to the test-gas boundary layer over a flat plate in an expansion tube is analyzed. Several combinations of test gas and acceleration gas are considered. The problem is treated in two conically similar limits: (1) when the time lag between the arrival of the shock and the interface at the leading edge of the plate is very large, and (2) when this lag is negligible. The time-dependent laminar-boundary-layer equations of a binary mixture of perfect gases are taken as the flow-governing equations. This coupled set of differential equations, written in terms of the Lam-Crocco variables, has been solved by a line-relaxation finite-difference techniques. The results presented include the Stanton number and the local skin-friction coefficient as functions of shock Mach number and the nondimensional distance-time variable. The results indicate that more than 95 percent of the test-gas boundary layer exists over a length, measured from the leading edge of the plate, equal to about three-tenths of the distance traversed by the interface in the free stream.

Gupta, R. N.↗
NASA NTRS2001

Thermal Effects on the Bearing Behavior of Composite Joints

Thermal effects on the pin-bearing behavior of an IM7/PET15 composite laminate are studied comprehensively. A hypothesis presents factors influencing a change in pin-bearing strength with a change in temperature for a given joint design. The factors include the change in the state of residual cure stress, the material properties, and the fastener fit with a change in temperature. Experiments are conducted to determine necessary lamina and laminate material property data for the IM7/PET15 being utilized in this study. Lamina material properties are determined between the assumed stress free temperature of 460 F down to -200 F. Laminate strength properties are determined for several lay-ups at the operating temperatures of 350 F, 70 F, and -200 F. A three-dimensional finite element analysis model of a composite laminate subject to compressive loading is developed. Both the resin rich layer located between lamina and the thermal residual stresses present in the laminate due to curing are determined to influence the state of stress significantly. Pin-bearing tests of several lay-ups were conducted to develop an understanding on the effect of temperature changes on the pin-bearing behavior of the material. A computational study investigating the factors influencing pin-bearing strength was performed. A finite element model was developed and used to determine the residual thermal cure stresses in the laminate containing a hole. Very high interlaminar stress concentrations were observed two elements away from the hole boundary at all three operating temperatures. The pin-bearing problem was modeled assuming a rigid frictionless pin and restraining only radial displacements at the hole boundary. A uniform negative pressure load was then applied to the straight end of the model. A solution, where thermal residual stresses were combined with the state of stress due to pin-bearing loads was evaluated. The presence of thermal residual stresses intensified the interlaminar stresses predicted at the hole boundary in the pin-bearing problem. This dissertation shows that changes in material properties drives pin-bearing strength degradation with increasing temperature.

Walker, Sandra Polesky↗
NASA NTRS1979

A model for the noise radiated by turbulent boundary layers and its interaction with laminar layers in supersonic flow

Motivated by the noise-boundary layer transition problem in supersonic flow, the interaction between incoming aerodynamic noise and the laminar boundary layer has been studied. The noise field has been modeled as a Mach wave radiation field consisting of discrete waves emanating from coherent turbulent entities moving downstream within the supersonic turbulent boundary layer. The individual disturbances are likened to miniature sonic booms and the laminar boundary layer is strafed by the waves as the sources move downstream. The mean, autocorrelation, and power spectral density of the field are expressed in terms of the wave shapes and their average arrival rates. The emphasis in the interaction analysis is on the behavior of the shocklets in the noise field. The shocklets are shown to be focused by the laminar boundary layer in its outer region. Borrowing wave-propagation terminology, this region is termed the caustic region. Using scaling laws from sonic boom work, focus factors at the caustic are estimated to vary from 2 to 6 for incoming shocklet strengths of 1 to .01 percent of the free stream pressure level.

Schopper, M. R.↗
NASA NTRS2001

Field Effects of Buoyancy on Lean Premixed Turbulent Flames

Buoyancy affects the entire flowfield of steady turbulent flames and this aspect of flame buoyancy coupling is largely unexplored by experiments or by theory. Open flames and flames within large confinements are free to expand and interact with the surrounding environment. In addition to fluid and combustion conditions, their aerodynamic flowfields are determined by the flame brush orientation and geometry, wake of the stabilizer, enclosure size, and of course, the gravitational field. Because the flowfield consists mainly of cold reactants (mostly in the nearfield) and hot products (mostly in the farfield), buoyancy effects are manifested in the farfield region. In upward pointing flames, an obvious effect is a favorable axial pressure gradient that accelerates the products thereby increasing the axial aerodynamic stretch rate. Intrinsic to turbulent flows, changes in mean aerodynamic stretch also couple to the fluctuating pressure field. Consequently, buoyancy can influence the turbulence intensities upstream and downstream of the flame. Flame wrinkling process, and heat release rate are also directly affected. This backward coupling mechanism is the so-called elliptic problem. To resolve the field effects of buoyancy would require the solution of three-dimensional non-linear Navier Stokes equations with full specification of the upstream, wall and downstream boundary conditions.

Cheng, R. K.↗
NASA NTRSDate not supplied

Coupled Material Response Simulations of Dragonfly Capsule and DrEAM Reconstruction

Icarus is a three-dimensional, unstructured, finite-volume material response solver developed at NASA Ames Research Center and has been recently used to analyze the material response of the Dragonfly capsule for a variety of problems. Since the Dragonfly capsule will be instrumented in a similar manner to the Mars 2020 and MSL capsules, it is important to assess how our current simulation tools can aid in understanding Dragonfly Entry Aerosciences Measurements (DrEAM). In this presentation, the Ares multi-physics tool that couples Icarus to the flow and radiation physics of the aerothermal environment will be used to better understand how modeling sensitivities might affect environment reconstructions and our understanding of measured data. Ares couples Icarus to the US3D flow solver and NERO, a reduced-order, finite-volume radiation transport solver and uses a customized time-scale management to enable coupled simulations for a large portion of the Dragonfly trajectory. As a result of the coupling, a more accurate and consistent assessment of modelling sensitivities to aerothermal reconstruction can be assessed. For example, the radiometer is sensitive to the quantity of CN in the flow, which is a byproduct of the free-stream methane but also the pyrolysis gas products injected into the boundary layer due to ablation of the heatshield. Ares will be used to conduct a simulation of the full Dragonfly capsule for 50 seconds of the trajectory. The different aerothermal and material response models will be discussed and the key results will be presented in terms of the simulated thermocouple and radiometer measurements on the backshell.

Ablation↗
NASA NTRS1991

Turbulization of shear flows in astrophysics

The problem of instability in shear flows is treated in a number of works which present a global analysis of irrotational (potential) modes in a 2-D compressible shear flow, implying that, in the presence of reflecting boundaries, those modes can increase and that the characteristic time of the increase is several orders as high as the dynamic time in the medium. Due to the slowness of the increase of irrotational modes and to the problematic character of the theory's applicability for thin accretion disks, an attempt was made to develop an alternative theory of turbulization shear flows. It is quite contrary to the above theory. Instead of the global analysis, a local one is made far from the flow boundaries. Instead of irrotational perturbations, vortical ones are examined. Lastly, the medium is considered to be incompressible. The flow in the accretion disk can be modeled by a plane Couette flow, without regard for boundary conditions: it can be modeled by a plane shear flow in infinite incompressible fluid. According to this scenario, there is a critical perturbation level in the free shear flow. If the level of initial perturbations exceeds the critical one, the flow is turbulized.

Chagelishvili, G. D.↗
NASA NTRS1975

Advances in Application of Fast Semidirect Computational Methods in Transonic Flow

This paper is intended as a review and summary of the advances made in a recently developed approach for rapid numerical solution of the equations of inviscid transonic aerodynamics. The investigation has been limited to two-dimensional, steady, inviscid flow over airfoils in a subsonic free stream, with emphasis on development of a rapid computational technique, rather than on generality of application. The approach uses finite-difference algorithms called "fast direct elliptic solvers" within an iteration scheme. "Direct" means that the entire computation field is solved at once, rather than in successive traverses over the field as in a point- or line-relaxation method. Such an iterative method is referred to as "semidirect." The iterative convergence can be faster than in other relaxation methods because changes are felt simultaneously at all points in each succeeding iteration. Direct elliptic solvers and semidirect methods have restrictions, but these are gradually being removed. Direct solvers were first developed for solving Poisson's equation on a rectangle without interior boundaries. A method to treat first-order systems, a direct Cauchy-Riemann solver has also been developed. Numerical treatment of part of a system of nonlinear equations by a Poisson solver has been reported. Also Poisson solvers in semidirect methods were used for nonseparable elliptic equations. The semidirect method was extended to the solution of a problem of mixed type, where the improved Murman-Cole transonic small-disturbance difference equations were solved. A slightly supercritical flow over a biconvex airfoil was treated successfully, but the iterations did not converge for more strongly supercritical conditions In another work the addition of terms ot both sides of the difference equations stabilized the iteration for supercritical conditions with large supersonic zones. For this, the Cauchy-Riemann solver was revised to incl,ude the needed terms. Most recently, the evaluation of parameters for rapid convergence and comparisons, with Murman's line-relaxation method was described. The method was extended to full second order accuracy in a fully conservative formulation in another work.

Martin, E. Dale↗
NASA NTRS1992

High speed transition prediction

The main objective of this work period was to develop, acquire and apply state-of-the-art tools for the prediction of transition at high speeds at NASA Ames. Although various stability codes as well as basic state codes were acquired, the development of a new Parabolized Stability Equation (PSE) code was minimal. The time that was initially allocated for development was used on other tasks, in particular for the Leading Edge Suction problem, in acquiring proficiency in various graphics tools, and in applying these tools to evaluate various Navier-Stokes and Euler solutions. The second objective of this work period was to attend the Transition and Turbulence Workshop at NASA Langley in July and August, 1991. A report on the Workshop follows. From July 8, 1991 to August 2, 1991, the author participated in the Transition and Turbulence Workshop at NASA Langley. For purposes of interest here, analysis can be said to consist of solving simplified governing equations by various analytical methods, such as asymptotic methods, or by use of very meager computer resources. From the composition of the various groups at the Workshop, it can be seen that analytical methods are generally more popular in Great Britain than they are in the U.S., possibly due to historical factors and the lack of computer resources. Experimenters at the Workshop were mostly concerned with subsonic flows, and a number of demonstrations were provided, among which were a hot-wire experiment to probe the boundary layer on a rotating disc, a hot-wire rake to map a free shear layer behind a cylinder, and the use of heating strips on a flat plate to control instability waves and consequent transition. A highpoint of the demonstrations was the opportunity to observe the rather noisy 'quiet' supersonic pilot tunnel in operation.

Gasperas, Gediminis↗
NASA NTRS1984

Progress in wind tunnel wall interference assessment/correction procedures at the NAE

Wall corrections as a function of wall porosity in the transonic wall interference problem was assessed. Effective porosities primarily for the two dimensional case were established as follows: (1) comparison of experimental data for two geometrically similar models of different chord/height ratio, an overall value of wall porosity could be deduced; (2) theoretical development which allows for unequal porosity for the floor and ceiling and wall boundary pressure measurements, porosities for floor and ceiling could be deduced; (3) a scheme was developed which allowed unequal porosity of floor and ceiling and streamwise varying porosity. The boundary layer development along the perforated floor and ceiling under the influence of the model pressure field, variations in boundary layer thickness underlining the difficulties in deducing meaningful values of wall porosity were determined. Wall boundary pressure measurement, in combination with singularity modelling of the airfoil, was sufficient to yield required information on the wall interference flow without having to establish some value for wall porosity. The singularity modelling of the airfoil initially covered only lift and volume but was extended to include drag and pitching moment, and second order volume term. It is shown by asymptotic transonic small disturbance analysis, that the derived corrections to angle of attack and free stream Mach number are correct to the first order.

Ohman, L. H.↗
NASA NTRS2010

Trajectory Optimization: OTIS 4

The latest release of the Optimal Trajectories by Implicit Simulation (OTIS4) allows users to simulate and optimize aerospace vehicle trajectories. With OTIS4, one can seamlessly generate optimal trajectories and parametric vehicle designs simultaneously. New features also allow OTIS4 to solve non-aerospace continuous time optimal control problems. The inputs and outputs of OTIS4 have been updated extensively from previous versions. Inputs now make use of objectoriented constructs, including one called a metastring. Metastrings use a greatly improved calculator and common nomenclature to reduce the user s workload. They allow for more flexibility in specifying vehicle physical models, boundary conditions, and path constraints. The OTIS4 calculator supports common mathematical functions, Boolean operations, and conditional statements. This allows users to define their own variables for use as outputs, constraints, or objective functions. The user-defined outputs can directly interface with other programs, such as spreadsheets, plotting packages, and visualization programs. Internally, OTIS4 has more explicit and implicit integration procedures, including high-order collocation methods, the pseudo-spectral method, and several variations of multiple shooting. Users may switch easily between the various methods. Several unique numerical techniques such as automated variable scaling and implicit integration grid refinement, support the integration methods. OTIS4 is also significantly more user friendly than previous versions. The installation process is nearly identical on various platforms, including Microsoft Windows, Apple OS X, and Linux operating systems. Cross-platform scripts also help make the execution of OTIS and post-processing of data easier. OTIS4 is supplied free by NASA and is subject to ITAR (International Traffic in Arms Regulations) restrictions. Users must have a Fortran compiler, and a Python interpreter is highly recommended.

Riehl, John P.↗
NASA NTRS1993

A photostationary state analysis of the NO2-NO system based on airborne observations from the subtropical/tropical North and South Atlantic

The Chemical Instrumentation Test and Evaluation 3 (CITE 3) NO-NO2 database has provided a unique opportunity to examine important aspects of tropospheric photochemistry as related to the rapid cycling between NO and NO2. Our results suggest that when quantitative testing of this photochemical system is based on airborne field data, extra precautions may need to be taken in the analysis. This was particularly true in the CITE 3 data analysis where different regional environments produced quite different results when evaluating the photochemical test ratio (NO2)(sub expt)/(NO2)(sub calc), designated here as R(sub E)/R(sub C). The quantity (NO2)(sub Calc) was evaluated using the following photostationary state expression: (NO2)(sub Calc) = k(sub 1)(O3) + k(sub 4)(HO2) + k(sub 5)(CH3O2) + k(sub 6)(RO2))(NO)(sub Expt)/J(sub 2). The four most prominent regional environmental data sets identified in this analysis were those labeled here as free-tropospheric northern hemisphere (FTNH), free-tropospheric tropical northern hemisphere (FTTNH), free-tropospheric southern hemisphere (FTSH), and tropical-marine boundary layer (plume) (TMBL(P)). The respective R(sub E)/R(sub C) mean and median values for these four data subsets were 1.74, 1.69; 3.00, 2.79; 1.01, 0.97; and 0.99, 0.94. Of the four data subsets listed, the two that were statistically the most robust were FTNH and FTSH; for these the respective R(sub E)/R(sub C) mean and standard deviation of the mean values were 1.74 +/- 0.07 and 1.01 +/- 0.04. The FTSH observations were in good agreement with theory, whereas those from the FTNH data set were in significant disagreement. An examination of the critical photochemical parameters O3, UV(zenith), NO, NO2, and non-methane hydrocarbons (NMHCs) for these two databases indicated that the most likely source of the R(sub E)/R(sub C) bias in the FTNH results was the presence of a systematic error in the observational data rather than a shortening in our understanding of fundamental photochemical processes. Although neither a chemical nor meteorological analyses of these data identified this error with complete certainty, they did point to the three most likely possibilities: (1) an NO2 interference from a yet unidentified NO(y) species: (2) the presence of unmeasured hydrocarbons, the integrated reactivity of which would be equivalent to approximately 2.7 parts per billion by volume (ppbv) of toluene; or (3) some combination of points (1) and (2). Details concerning hypotheses (1) and (2) as well as possible ways to minimize these problems in future airborne missions are discussed.

Davis, D. D.↗
NASA NTRS1987

Transient response for interaction of two dynamic bodies

During the launch sequence of any space vehicle complicated boundary interactions occur between the vehicle and the launch stand. At the start of the sequence large forces exist between the two; contact is then broken in a short but finite time which depends on the release mechanism. The resulting vehicle response produces loads which are very high and often form the design case. It is known that the treatment of the launch pad as a second dynamic body is significant for an accurate prediction of launch response. A technique was developed for obtaining loads generated by the launch transient with the effect of pad dynamics included. The method solves uncoupled vehicle and pad equations of motion. The use of uncoupled models allows the simulation of vehicle launch in a single computer run. Modal formulation allows a closed-form solution to be written, eliminating any need for a numerical integration algorithm. When the vehicle is on the pad the uncoupled pad and vehicle equations have to be modified to account for the constraints they impose on each other. This necessitates the use of an iterative procedure to converge to a solution, using Lagrange multipliers to apply the required constraints. As the vehicle lifts off the pad the coupling between the vehicle and the pad is eliminated point by point until the vehicle flies free. Results obtained by this method were shown to be in good agreement with observed loads and other analysis methods. The resulting computer program is general, and was used without modification to solve a variety of contact problems.

Prabhakar, A.↗
NASA NTRS2010

Comments on "Atlantic Tropical Cyclogenetic Processes during SOP-3 NAMMA in the GEOS-5 Global Data Assimilation and Forecast System"

Considerable attention has been given to the potential negative impacts of the Saharan air layer (SAL) in recent years. Researchers recently raised questions about the negative impacts of Dunion and Velden and other studies in terms of storms that reached at least tropical storm strength and suggested that the SAL was an intrinsic part of the tropical cyclone environment for both storms that weaken after formation and those that intensify. Braun also suggested that several incorrect assumptions underlie many of the studies on the negative impacts of the SAL, including assumptions that most low-to-midlevel dry tropical air is SAL air, that the SAL is dry throughout its depth, and that the proximity of the SAL to storms struggling to intensify implies some role in that struggle. The recent paper by Reale et al.(RL1) is an example of the problems inherent in some of these assumptions. In their paper, RL1 analyze a simulation from the Global Earth Observing System (GEOS-5) global model and describe an extensive tongue of warm, dry air that stretches southward from at least 30 deg N (the northern limit of their plots) and wraps into a low pressure system during the period 26-29 August 2006, suppressing convection and possibly development of the African easterly wave associated with that low pressure system. They attributed the warm, dry tongue to the SAL (i.e., heating of the air mass during passage over the Sahara and radiative warming of the dust layer). Whether it was their intention, the implication is that this entire feature is due solely to the SAL and not to other possible sources of dry air or warmth. In addition, they suggested that a cool tongue of air in the boundary layer located directly beneath the elevated warm, dry tongue (forming a thermal dipole) was possibly the result of reduced solar radiation caused by an overlying dust layer. They stated that "the cool anomaly in the lower levels does not have any plausible explanation relying only on transport." In this comment, evidence from satellite and global meteorological analyses is presented that casts considerable doubt upon RL1 s interpretation of the GEOS-5 forecasts and their conclusion that the results implied a negative role of the SAL. We show that the major portion of the warm, dry air aloft was located in a nearly dust-free slot between two Saharan dust outbreaks, had a significant source from the midlatitudes (>30 N), and was likely driven by strong subsidence warming and drying. In addition, when wind fields are examined in a reference frame moving with the wave, National Centers for Environmental Prediction (NCEP) global meteorological analyses suggest that the cool tongue in the boundary layer can be readily explained by transport of cooler air from higher latitudes. At the very least, it offers a plausible alternative explanation for the cool tongue that does not rely on radiative impacts of the dust.

Braun, Scott A.↗