Search NASA⌕ Search

SEARCH · Search NASA

Results for “model simplification”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Comparison of Eulerian Bin and Lagrangian Particle-Based Schemes in Simulations of Pi Chamber Dynamics and Microphysics

This paper presents a comparison of simulations applying either a traditional Eulerian bin microphysics or a novel particle-based Lagrangian approach to represent CCN activation and cloud droplet growth. The Eulerian microphysics solve the evolution equation for the spectral density function, whereas the Lagrangian approach follows computational particles referred to as superdroplets. Each superdroplet represents a multiplicity of natural droplets that makes the Lagrangian approach computationally feasible. The two schemes apply identical representation of CCN activation and use the same droplet growth equation; these make direct comparison between the two schemes practical. The comparison, the first of its kind, applies an idealized simulation setup motivated by laboratory experiments with the Pi Chamber and previous model simulations of the Pi Chamber dynamics and microphysics. The Pi Chamber laboratory apparatus considers interactions between turbulence, CCN activation, and cloud droplet growth in moist Rayleigh–Bénard convection. Simulated steady-state droplet spectra averaged over the entire chamber are similar, with the mean droplet concentration, mean radius, and spectral width close in Eulerian and Lagrangian simulations. Small differences that do exist are explained by the inherent differences between the two schemes and their numerical implementation. The local droplet spectra differ substantially, again in agreement with the inherent limitations of the theoretical foundation behind each approach. There is a general agreement between simulations and Pi Chamber observations, with simplifications of the CCN activation and droplet growth equation used in the simulations likely explaining specific differences.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1985

Optimization studies for the development of adaptive wall wind tunnel

The development of adaptive wall wind tunnels makes it possible to establish interference-free boundaries in the region of the model. It is, therefore, feasible to simulate free air conditions. However, the amount of total testing time involved in establishing the required conditions and the complexity of the needed instrumentation prevent currently the use of such methods for routine production testing. Optimization studies have been proposed with the objective to find approaches which will make a more extended application of adaptive wall techniques possible. A new 25cm x 11cm adaptive wall test section was designed and installed to carry out some of these studies. Investigations envisaged are related to the development of a one step convergence method to reduce total testing time, and to measurement procedures which will permit a simplification of the instrumentation. Attention is given to the proposed concept, the test section, instrumentation, and the influence coefficients.

Bodapati, S.↗
NASA NTRS1991

Latitudinal dynamics of steady solar wind flows

In many previous studies, it has been assumed that the streamlines of the solar wind are in a purely radial direction beyond a few solar radii. Therefore, the nonradial modulations by the magnetic force and gas pressure is neglected at large heliocentric distances. The work reported in this paper includes a two-dimensional magnetohydrodynamic (MHD) study of the latitudinal structure of solar wind flows. The dynamical effects of the interplanetary magnetic field (IMF) are investigated by MHD solutions and by comparing them with both one-dimensional and two-dimensional hydrodynamic (HD) solutions. The corotational effect of the sun and the influence of transient events in the solar atmosphere are neglected. With such a simplification, the latitudinal structure of the steady solar wind is controlled by the dynamic balance of the magnetic force and the gas pressure. Several different patterns of the evolution of the solar wind latitudinal structure are exhibited in the calculation. The results indicate the existence of a proton number density maximum at the magnetic neutral line, whether or not there is a density maximum or minimum at the inner boundary. The drift motion of magnetic field lines toward the magnetic neutral line enhances the magnetic field strength around the neutral sheet, which may provide a possible explanation of the discrepancy between the measured IMF at 1 AU and that extrapolated from the photospheric magnetic field by the current source-surface modeling.

Yang, W.-H.↗
NASA NTRS1995

Development Status of the NSTAR Ion Propulsion System Power Processor

A 0.5-2.3 kW xenon ion propulsion system is presently being developed under the NASA Solar Electric Propulsion Technology Application Readiness (NSTAR) program. This propulsion system includes a 30 cm diameter xenon ion thruster, a Digital Control Interface Unit, a xenon feed system, and a power processing unit (PPU). The PPU consists of the power supply assemblies which operate the thruster neutralizer, main discharge chamber, and ion optics. Also included are recycle logic and a digital microcontroller. The neutralizer and discharge power supplies employ a dual use configuration which combines the functions of two power supplies into one, significantly simplifying the PPU. Further simplification was realized by implementing a single thruster control loop which regulates the beam current via the discharge current. Continuous throttling is possible over a 0.5-2.3 kW output power range. All three power supplies have been fabricated and tested with resistive loads, and have been combined into a single breadboard unit with the recycle logic and microcontroller. All line and load regulation test results show the power supplies to be within the NSTAR flight PPU specified power output of 1.98 kW. The overall efficiency of the PPU, calculated as the combined efficiencies of the power supplies and controller, at 2.3 kW delivered to resistive loads was 0.90. The component was 6.16 kg. Integration testing of the neutralizer and discharge power supplies with a functional model thruster revealed no issues with discharge ignition or steady state operation.

Hamley, John A.↗
NASA NTRS1989

Numerical Simulation of Flow Through an Artificial Heart

A solution procedure was developed that solves the unsteady, incompressible Navier-Stokes equations, and was used to numerically simulate viscous incompressible flow through a model of the Pennsylvania State artificial heart. The solution algorithm is based on the artificial compressibility method, and uses flux-difference splitting to upwind the convective terms; a line-relaxation scheme is used to solve the equations. The time-accuracy of the method is obtained by iteratively solving the equations at each physical time step. The artificial heart geometry involves a piston-type action with a moving solid wall. A single H-grid is fit inside the heart chamber. The grid is continuously compressed and expanded with a constant number of grid points to accommodate the moving piston. The computational domain ends at the valve openings where nonreflective boundary conditions based on the method of characteristics are applied. Although a number of simplifing assumptions were made regarding the geometry, the computational results agreed reasonably well with an experimental picture. The computer time requirements for this flow simulation, however, are quite extensive. Computational study of this type of geometry would benefit greatly from improvements in computer hardware speed and algorithm efficiency enhancements.

Rogers, Stuart E.↗
DOE OSTI2022

Metaplectic geometrical optics for ray-based modeling of caustics: Theory and algorithms

The optimization of radio frequency-wave (RF) systems for fusion experiments is often performed using ray-tracing codes, which rely on the geometrical-optics (GO) approximation. However, GO fails at caustics such as cutoffs and focal points, erroneously predicting the wave intensity to be infinite. This is a critical shortcoming of GO, since the caustic wave intensity is often the quantity of interest, e.g., RF heating. Full-wave modeling can be used instead, but the computational cost limits the speed at which such optimizations can be performed. Here, we have developed a less expensive alternative called metaplectic geometrical optics (MGO). Instead of evolving waves in the usual x (coordinate) or k (spectral) representation, MGO uses a mixed X$\equiv$Ax+Bk representation. By continuously adjusting the matrix coefficients A and B along the rays, one can ensure that GO remains valid in the X coordinates without caustic singularities. The caustic-free result is then mapped back onto the original x space using metaplectic transforms. Here, we overview the MGO theory and review algorithms that will aid the development of an MGO-based ray-tracing code. We show how using orthosymplectic transformations leads to considerable simplifications compared to previously published MGO formulas. We also prove explicitly that MGO exactly reproduces standard GO when evaluated far from caustics (an important property that until now has only been inferred from numerical simulations), and we relate MGO to other semiclassical caustic-removal schemes published in the literature. Finally this discussion is then augmented by an explicit comparison of the computed spectrum for a wave bounded between two cutoffs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1982

The CMC:3DPNS computer program for prediction of three-dimensional, subsonic, turbulent aerodynamic juncture region flow. Volume 1: Theoretical

An order-of-magnitude analysis of the subsonic three dimensional steady time averaged Navier-Stokes equations, for semibounded aerodynamic juncture geometries, yields the parabolic Navier-Stokes simplification. The numerical solution of the resultant pressure Poisson equation is cast into complementary and particular parts, yielding an iterative interaction algorithm with an exterior three dimensional potential flow solution. A parabolic transverse momentum equation set is constructed, wherein robust enforcement of first order continuity effects is accomplished using a penalty differential constraint concept within a finite element solution algorithm. A Reynolds stress constitutive equation, with low turbulence Reynolds number wall functions, is employed for closure, using parabolic forms of the two-equation turbulent kinetic energy-dissipation equation system. Numerical results document accuracy, convergence, and utility of the developed finite element algorithm, and the CMC:3DPNS computer code applied to an idealized wing-body juncture region. Additional results document accuracy aspects of the algorithm turbulence closure model.

Baker, A. J.↗
NASA NTRS2002

Partially Decentralized Control Architectures for Satellite Formations

In a partially decentralized control architecture, more than one but less than all nodes have supervisory capability. This paper describes an approach to choosing the number of supervisors in such au architecture, based on a reliability vs. cost trade. It also considers the implications of these results for the design of navigation systems for satellite formations that could be controlled with a partially decentralized architecture. Using an assumed cost model, analytic and simulation-based results indicate that it may be cheaper to achieve a given overall system reliability with a partially decentralized architecture containing only a few supervisors, than with either fully decentralized or purely centralized architectures. Nominally, the subset of supervisors may act as centralized estimation and control nodes for corresponding subsets of the remaining subordinate nodes, and act as decentralized estimation and control peers with respect to each other. However, in the context of partially decentralized satellite formation control, the absolute positions and velocities of each spacecraft are unique, so that correlations which make estimates using only local information suboptimal only occur through common biases and process noise. Covariance and monte-carlo analysis of a simplified system show that this lack of correlation may allow simplification of the local estimators while preserving the global optimality of the maneuvers commanded by the supervisors.

Carpenter, J. Russell↗
DOE OSTI2022

Measurement and modeling of methane diffusion in hydrocarbon mixtures

Methane (CH 4 ) dissolution and diffusive mass transfer in liquid hydrocarbon mixtures is of key interest in the context of enhanced oil recovery from tight shales. In this paper, we have studied CH 4 dissolution and diffusion in normal alkane mixtures at a temperature of 50 ⁰C and at pressures of ~8 MPa. For the measurement of CH 4 dissolution/diffusion in bulk liquid hydrocarbon mixtures, we have utilized a high pressure and temperature, constant-volume diffusion (CV-D) setup. Here, we have studied CH 4 diffusion in mixtures with three long-chain normal alkanes: decane (C 10 ), dodecane (C 12 ) and hexadecane (C 16 ). We have measured CH 4 solubility and diffusion in its binary mixtures with each of these three normal alkanes, as well as in ternary and quaternary mixtures. During the experiments, the swelling of the liquid mixture due to CH 4 dissolution was measured in situ via a cathetometer and was subsequently integrated into the data analysis. A key conclusion from this study is that the solubility and transport properties of the multicomponent mixtures can be predicted accurately from binary mixture measurements using an appropriate Equation of State (EOS) and Wilke’s simplification of the classical Maxwell-Stefan (MS) diffusion theory. This observation can facilitate accurate prediction of diffusive mass transfer in more complex liquid hydrocarbon mixtures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Nucleon-induced inelastic scattering with statistical strength functions and the ECIS direct reaction code

Modern theoretical descriptions of inelastic scattering make use of multi-step direct reaction approaches together with transition potentials obtained from sophisticated nuclear structure models. Here we demonstrate how the complexity of such calculations can be reduced to permit simpler ones, also using the ECIS code, but providing an almost equally precise alternative to a much more detailed calculation. We have studied the transition form factors within the random phase approximation (RPA), where these are obtained as linear combinations of particle–hole states. At moderate to high excitation energies, where interference effects tend to disappear, we have proposed an independent particle–hole formalism in which particle–hole states are spread in energy with an appropriate strength function obtained from the RPA. The effects of more complex modes, such as 2p–2h ones, are simulated with widths calculated in a semi-classical context. In this work, we verify the validity of our approximations for pre-equilibrium proton-induced reactions on 90 Zr target. Our calculations provide a good description of the reaction data and point toward a simplification of the description of nucleon-induced reactions based on averages of microscopic details of the projectile–target interaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1993

Receptivity in parallel flows: An adjoint approach

Linear receptivity studies in parallel flows are aimed at understanding how external forcing couples to the natural unstable motions which a flow can support. The vibrating ribbon problem models the original Schubauer and Skramstad boundary layer experiment and represents the classic boundary layer receptivity problem. The process by which disturbances are initiated in convectively-unstable jets and shear layers has also received attention. Gaster was the first to handle the boundary layer analysis with the recognition that spatial modes, rather than temporal modes, were relevant when studying convectively-unstable flows that are driven by a time-harmonic source. The amplitude of the least stable spatial mode, far downstream of the source, is related to the source strength by a coupling coefficient. The determination of this coefficient is at the heart of this type of linear receptivity study. The first objective of the present study was to determine whether the various wave number derivative factors, appearing in the coupling coefficients for linear receptivity problems, could be reexpressed in a simpler form involving adjoint eigensolutions. Secondly, it was hoped that the general nature of this simplification could be shown; indeed, a rather elegant characterization of the receptivity properties of spatial instabilities does emerge. The analysis is quite distinct from the usual Fourier-inversion procedures, although a detailed knowledge of the spectrum of the Orr-Sommerfeld equation is still required. Since the cylinder wake analysis proved very useful in addressing control considerations, the final objective was to provide a foundation upon which boundary layer control theory may be developed.

Hill, D. Christopher↗
DOE OSTI2022

Influence of cementation on the yield surface of rocks numerically determined from digital microstructures

Digital Rock Physics has reached a level of maturity on the characterisation of primary properties that depend on the microstructure – such as porosity, permeability or elastic moduli – by numerically solving field equations on μCT scan images of rock. After small deformations or at depth though, most rocks eventually reach their limit of elasticity and the complementary plastic properties are needed to describe the full mechanical behaviour. Currently, determination of a rock’s yield surface from its microstructure is often restricted to semi-analytical criteria derived by limit analysis or numerical simulations performed on idealised geometries. Such simplification lacks representativeness, particularly for processes that affect directly the pore-grain interface such as the cementation phenomenon, happening during diagenesis. Eventually, only direct numerical simulation of elasto-plasticity performed on digitalised microstructures can be used to assess the strength of different cemented materials and its evolution with the alteration of the microstructure. In this study, we provide a comprehensive parametric study on the impact of cementation on rock strength for real microstructures of cemented granular materials. Compared to most previous studies, the whole yield surface is determined numerically (using Finite Element Method) in order to assess the influence of cementation for different stress-paths. The previously known tendency of rock to strengthen with increasing cementation volume is verified. New results on the influence of cement property namely Young’s modulus, friction and cohesion on the rock’s yield surface are explored. The envelopes obtained are compared to the ones obtained by experimental data and existing models. Finally, the framework presented in this study showcases the wider possibility of determining any rock’s or porous material’s yield surface from its microstructure.

36 MATERIALS SCIENCE↗
NASA NTRS2017

Simulation of the Ozone Monitoring Instrument Aerosol Index Using the NASA Goddard Earth Observing System Aerosol Reanalysis Products

We provide an analysis of the commonly used Ozone Monitoring Instrument (OMI) aerosol index (AI) product for qualitative detection of the presence and loading of absorbing aerosols. In our analysis, simulated top-of-atmosphere (TOA) radiances are produced at the OMI footprints from a model atmosphere and aerosol profile provided by the NASA Goddard Earth Observing System (GEOS-5) Modern-Era Retrospective Analysis for Research and Applications aerosol reanalysis (MERRAero). Having established the credibility of the MERRAero simulation of the OMI AI in a previous paper we describe updates in the approach and aerosol optical property assumptions. The OMI TOA radiances are computed in cloud-free conditions from the MERRAero atmospheric state, and the AI is calculated. The simulated TOA radiances are fed to the OMI aerosol retrieval algorithms, and its retrieved AI (OMAERUV AI) is compared to the MERRAero calculated AI. Two main sources of discrepancy are discussed: one pertaining the OMI algorithm assumptions of the surface pressure, which are generally different from what the actual surface pressure of an observation is, and the other related to simplifying assumptions in the molecular atmosphere radiative transfer used in the OMI algorithms. Surface pressure assumptions lead to systematic biases in the OMAERUV AI, particularly over the oceans. Simplifications in the molecular radiative transfer lead to biases particularly in regions of topography intermediate to surface pressures of 600hPa and 1013.25hPa. Generally, the errors in the OMI AI due to these considerations are less than 0.2 in magnitude, though larger errors are possible, particularly over land. We recommend that future versions of the OMI algorithms use surface pressures from readily available atmospheric analyses combined with high-spatial resolution topographic maps and include more surface pressure nodal points in their radiative transfer lookup tables.

remote sensing↗
DOE OSTI2022

The role of pre-existing heterogeneities in materials under shock and spall

There has been a challenge for many decades to understand how heterogeneities influence the behavior of materials under shock loading, eventually leading to spall formation and failure. Experimental, analytical, and computational techniques have matured to the point where systematic studies of materials with complex microstructures under shock loading and the associated failure mechanisms are feasible. This is enabled by more accurate diagnostics as well as characterization methods. As interest in complex materials grows, understanding and predicting the role of heterogeneities in determining the dynamic behavior becomes crucial. Early computational studies, hydrocodes, in particular, historically preclude any irregularities in the form of defects and impurities in the material microstructure for the sake of simplification and to retain the hydrodynamic conservation equations. Contemporary computational methods, notably molecular dynamics simulations, can overcome this limitation by incorporating inhomogeneities albeit at a much lower length and time scale. This review discusses literature that has focused on investigating the role of various imperfections in the shock and spall behavior, emphasizing mainly heterogeneities such as second-phase particles, inclusions, and voids under both shock compression and release. Pre-existing defects are found in most engineering materials, ranging from thermodynamically necessary vacancies, to interstitial and dislocation, to microstructural features such as inclusions, second phase particles, voids, grain boundaries, and triple junctions. This literature review explores the interaction of these heterogeneities under shock loading during compression and release. Systematic characterization of material heterogeneities before and after shock loading, along with direct measurements of Hugoniot elastic limit and spall strength, allows for more generalized theories to be formulated. Further, continuous improvement toward time-resolved, in situ experimental data strengthens the ability to elucidate upon results gathered from simulations and analytical models, thus improving the overall ability to understand and predict how materials behave under dynamic loading.

36 MATERIALS SCIENCE↗
DOE OSTI2023

MOOSE Reactor Module: An Open-Source Capability for Meshing Nuclear Reactor Geometries

The U.S. Department of Energy (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has developed numerous physics solvers utilizing the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) framework for multiphysics reactor analysis. These solvers require input finite element meshes representing the discretized spatial domain. Typically, reactor analysts turn to licensed tools for the creation of reactor geometry meshes. Recently, open-source functionality has been added to the MOOSE framework to mesh common reactor geometries and improve MOOSE-based nuclear reactor application user workflows. The new functionality is primarily contained in the new Reactor module of MOOSE and includes support for hexagonal pins, assemblies, and cores, extended Cartesian geometry support, options for modeling static and rotating control drums within a hexagonal assembly, core periphery triangulation, and automatic tagging of pin, assembly, plane, and depletion regions for easier post processing of physics results. A set of reactor geometry mesh builder objects further streamlines the construction of hexagonal and Cartesian cores and allows mapping of materials to regions during mesh generation. The meshes produced with the MOOSE Reactor module may be used directly within MOOSE-based applications or exported as Exodus II files for use in other finite element solvers. The tools have been demonstrated and verified using a variety of NEAMS physics solvers on a range of reactor applications, including a sodium-cooled fast reactor core analysis using Griffin, a fast reactor assembly thermal deformation analysis using MOOSE Tensor Mechanics, and a heat pipe–cooled microreactor coupled analysis using Griffin, Bison, and Sockeye. MOOSE’s Reactor module provides significant advantages compared to the use of external meshing tools when analyzing Cartesian and hexagonal reactor lattices using MOOSE-based applications: immediate accessibility (open-source) to the end user, low barrier to entry for new users, speed of mesh generation, volume preservation of meshed fuel pins, and simplification of analysis workflow when used in conjunction with MOOSE-based applications.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2023

Techno-economic analysis of advanced small modular nuclear reactors

Here, small modular nuclear reactors (SMRs) represent a robust opportunity to develop low-carbon and reliable power with the potential to meet cost parity with conventional power systems. This study presents a detailed, bottom-up economic evaluation of a 12 × 77 MW e (924 MW e total) light-water SMR (LW-SMR) plant, a 4 × 262 MW e (1,048 MW e ) gas-cooled SMR (GC-SMR) plant, and a 5 × 200 MW e (1,000 MW e total) molten salt SMR (MS-SMR) plant. Cost estimates are derived from equipment costs, labor hours, material inputs, and process-engineering models. The advanced SMRs are compared to natural gas combined cycle plants with and without post-combustion carbon capture and a conventional large nuclear reactor. Overnight capital cost (OCC) and levelized cost of energy (LCOE) estimates are developed. The OCC of the LW-SMR, GC-SMR, and MS-SMR are found to be $\$4,844$/kW, $\$4,355$/kW, and $\$3,985$/kW respectively. The LCOE of the LW-SMR, GC-SMR, and MS-SMR are found to be $\$89.6$/MWh, $\$81.5$/MWh, and $\$80.6$/MWh respectively. A Monte Carlo analysis is performed, for which the OCC and construction time of the LW-SMR is found to have a lower mean and standard deviation than a conventional large reactor. The LW-SMR OCC is found to have a mean of $\$5,233$/kW with a standard deviation of $\$658$/kW and a 90 % probability of remaining between $\$4,254$/kW and $\$6,399$/kW, while the construction duration is found to have a mean of 4.5 years with a standard deviation of 0.8 years and a 90 % probability of remaining between 3.4 and 6.0 years. The economic impact of economies of scale, simplification, modularization, and construction time for SMRs are discussed. Additionally, policy implications for direct SMR capital subsidies and the impact of a carbon tax on natural gas emissions are explored.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2024

Additively Manufactured Plastic Plasma Spectrometer

We report results in development and testing of a low resource top hat electrostatic analyzer (ESA) for space plasma measurements. This device has been additively manufactured (3D-printed) using fused deposition modeling. The classic top hat design is composed of four plastic pieces, without any surface coatings. The three conducting electrodes are printed from carbon nanotube infused polyether ether ketone (CNT-PEEK). The fourth piece, an insulating electrode support, uses pure PEEK. This ESA is designed to detect electrons in space from 10eV up to 30 keV. We demonstrate that the printed CNT-PEEK is sufficiently electrically conductive to support the fast high voltage slewing often required for high time resolution measurements. The plastic ESA has been successfully vibrated beyond standard pre-flight levels, tested under keV electron beam illumination over a wide range of temperatures, and tested under UV illumination, simulating the solar Ly-α flux. In comparison with an identical machined aluminum ESA, our CNT-PEEK ESA provides nominal energy/angle band passes, closely matching simulation. These band passes imply minimal impact from surface charging at beam energies of 2-3 keV, although more investigation is needed. We also find that the CNT-PEEK ESA provides far superior out-of-band electron rejection and UV photon rejection compared to the machined aluminum ESA. We do not detect any problems with trapped gases or outgassing. This development offers potential for significant mass savings, implementation of otherwise unattainable geometric configurations, and dramatic simplification in manufacturing and assembly processes required for development of space plasma instruments.

Quetzal Larrick↗
NASA NTRS1944

Experimental Verification of the Rudder-Free Stability Theory for an Airplane Model Equipped with Rudders Having Negative Floating Tendency and Negligible Friction

An investigation has been made in the Langley free-flight tunnel to obtain an experimental verification of the theoretical rudder-free stability characteristics of an airplane model equipped with conventional rudders having negative floating tendencies and negligible friction. The model used in the tests was equipped with a conventional single vertical tail having rudder area 40 percent of the vertical tail area. The model was tested both in free flight and mounted on a strut that allowed freedom only in yaw. Tests were made with three different amounts of rudder aerodynamic balance and with various values of mass, moment of inertia, and center-of-gravity location of the rudder. Most of the stability derivatives required for the theoretical calculations were determined from forced and free-oscillation tests of the particular model tested. The theoretical analysis showed that the rudder-free motions of an airplane consist largely of two oscillatory modes - a long-period oscillation somewhat similar to the normal rudder-fixed oscillation and a short-period oscillation introduced only when the rudder is set free. It was found possible in the tests to create lateral instability of the rudder-free short-period mode by large values of rudder mass parameters even though the rudder-fixed condition was highly stable. The results of the tests and calculation indicated that for most present-day airplanes having rudders of negative floating tendency, the rudder-free stability characteristics may be examined by simply considering the dynamic lateral stability using the value of the directional-stability parameter Cn(sub p) for the rudder-free condition in the conventional controls-fixed lateral-stability equations. For very large airplanes having relatively high values of the rudder mass parameters with respect to the rudder aerodynamic parameters, however, analysis of the rudder-free stability should be made with the complete equations of motion. Good agreement between calculated and measured rudder-free stability characteristics was obtained by use of the general rudder-free stability theory, in which four degrees of lateral freedom are considered. When this assumption is made that the rolling motions alone or the lateral and rolling motions may be neglected in the calculations of rudder-free stability, it is possible to predict satisfactorily the characteristics of the long-period (Dutch roll type) rudder-free oscillation for airplanes only when the effective-dihedral angle is small. With these simplifying assumptions, however, satisfactory prediction of the short-period oscillation may be obtained for any dihedral. Further simplification of the theory based on the assumption that the rudder moment of inertia might be disregarded was found to be invalid because this assumption made it impossible to calculate the characteristics of the short-period oscillations.

STABILITY- FAIRCHILD R2K-1↗