Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Analytical Balance”

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 289 records · Page 16

A modified Dodge algorithm for the parabolized Navier-Stokes equations and compressible duct flows

A revised version of Dodge's split-velocity method for numerical calculation of compressible duct flow was developed. The revision incorporates balancing of mass flow rates on each marching step in order to maintain front-to-back continuity during the calculation. The (checkerboard) zebra algorithm is applied to solution of the three dimensional continuity equation in conservative form. A second-order A-stable linear multistep method is employed in effecting a marching solution of the parabolized momentum equations. A checkerboard iteration is used to solve the resulting implicit nonlinear systems of finite-difference equations which govern stepwise transition. Qualitative agreement with analytical predictions and experimental results was obtained for some flows with well-known solutions. Previously announced in STAR as N82-16363

Cooke, C. H.↗

Coupled hydromagnetic wave excitation and ion acceleration upstream of the earth's bow shock

A self-consistent theory is set forth for the excitation of hydromagnetic waves and the acceleration of 'diffuse' ions upstream of the earth's bow shock in the quasi-equilibrium that results when the solar wind velocity and the interplanetary magnetic field are nearly parallel. For the waves, the quasi-equilibrium derives from a balance between excitation by the ions, which stream relative to the solar wind plasma, and convective loss to the magnetosheath. For the diffuse ions, the quasi-equilibrium derives from a balance between injection at the shock front, confinement to the foreshock by pitch angle scattering on the waves, and acceleration by compression at the shock front. It also results from loss to the magnetosheath, loss due to escape upstream of the foreshock, and loss via diffusion perpendicular to the average magnetic field onto field lines that do not connect to the shock front. Diffusion equations describing the ion transport and wave kinetic equations describing the hydromagnetic wave transport are solved self-consistently to give analytical expressions. These describe (1) the differential wave intensity spectrum as a function of frequency and distance from the bow shock and (2) the ion omnidirectional distribution functions and anisotropies as functions of energy and distance from the bow shock.

Lee, M. A.↗

The Martian climate and energy balance models with CO2/H2O atmospheres

The analysis begins with a seasonal energy balance model (EBM) for Mars. This is used to compute surface temperature versus x = sin(latitude) and time over the seasonal cycle. The core model also computes the evolving boundaries of the CO2 icecaps, net sublimational/condensation rates, and the resulting seasonal pressure wave. Model results are compared with surface temperature and pressure history data at Viking lander sites, indicating fairly good agreement when meridional heat transport is represented by a thermal diffusion coefficient D approx. 0.015 W/sq. m/K. Condensational wind distributions are also computed. An analytic model of Martian wind circulation is then proposed, as an extension of the EMB, which incorporates vertical wind profiles containing an x-dependent function evaluated by substitution in the equation defining the diffusion coefficient. This leads to a parameterization of D(x) and of the meridional circulation which recovers the high surface winds predicted by dynamic Mars atmosphere models (approx. 10 m/sec). Peak diffusion coefficients, D approx. 0.6 w/sq m/K, are found over strong Hadley zones - some 40 times larger than those of high-latitude baroclinic eddies. When the wind parameterization is used to find streamline patterns over Martian seasons, the resulting picture shows overturning hemispheric Hadley cells crossing the equator during solstices, and attaining peak intensities during the south summer dust storm season, while condensational winds are most important near the polar caps.

Hoffert, M. I.↗

Deviation of Long-Period Tides from Equilibrium: Kinematics and Geostrophy

New empirical estimates of the long-period fortnightly (Mf) tide obtained from TOPEX/Poseidon (T/P) altimeter data confirm significant basin-scale deviations from equilibrium. Elevations in the low-latitude Pacific have reduced amplitude and lag those in the Atlantic by 30 deg or more. These interbasin amplitude and phase variations are robust features that are reproduced by numerical solutions of the shallow-water equations, even for a constant-depth ocean with schematic interconnected rectangular basins. A simplified analytical model for cooscillating connected basins also reproduces the principal features observed in the empirical solutions. This simple model is largely kinematic. Zonally averaged elevations within a simple closed basin would be nearly in equilibrium with the gravitational potential, except for a constant offset required to conserve mass. With connected basins these offsets are mostly eliminated by interbasin mass flux. Because of rotation, this flux occurs mostly in a narrow boundary layer across the mouth and at the western edge of each basin, and geostrophic balance in this zone supports small residual offsets (and phase shifts) between basins. The simple model predicts that this effect should decrease roughly linearly with frequency, a result that is confirmed by numerical modeling and empirical T/P estimates of the monthly (Mm) tidal constituent. This model also explains some aspects of the anomalous nonisostatic response of the ocean to atmospheric pressure forcing at periods of around 5 days.

Egbert, Gary D.↗

Space shuttle orbiter reaction control system jet interaction study

The space shuttle orbiter has forward mounted and rear mounted Reaction Control Systems (RCS) which are used for orbital maneuvering and also provide control during entry and abort maneuvers in the atmosphere. The effects of interaction between the RCS jets and the flow over the vehicle in the atmosphere are studied. Test data obtained in the NASA Langley Research Center 31 inch continuous flow hypersonic tunnel at a nominal Mach number of 10.3 is analyzed. The data were obtained with a 0.01 scale force model with aft mounted RCS nozzles mounted on the sting off of the force model balance. The plume simulations were accomplished primarily using air in a cold gas simulation through scaled nozzles, however, various cold gas mixtures of Helium and Argon were also tested. The effect of number of nozzles was tested as were limited tests of combined controls. The data show that RCS nozzle exit momentum ratio is the primary correlating parameter for effects where the plume impinges on an adjacent surface and mass flow ratio is the parameter where the plume interaction is primarily with the external stream. An analytic model of aft mounted RCS units was developed in which the total reaction control moments are the sum of thrust, impingement, interaction, and cross-coupling terms.

Rausch, J. R.↗

The Inner Workings of Magnetic Reconnection: Diffusion Region in the Balance

The question of whether the micro scale controls the macroscale or vice-versa remains one of the most challenging problems in plasmas. A particular topic of interest within this context is collisionless magnetic reconnection, where both points of views are espoused by different groups of researchers. This presentation will focus on this topic. We will begin by analyzing the properties of electron diffusion region dynamics both for guide field and anti-parallel reconnection, and how they can be scaled to different inflow conditions. As a next step, we will study typical temporal variations of the microscopic dynamics with the objective of understanding the potential for secular changes to the macroscopic system. The research will be based on a combination of analytical theory and numerical modeling.

Hesse, Michael↗

Friction damping studies in multiple turbine blade systems by lumped mass method

Analytical studies were conducted on multiple turbine blade systems using the lumped mass method. Each blade was idealized by a two mass-two spring model whose modal values were determined from the known frequencies corresponding to the first two bending modes and the frequency corresponding to the platform lockup condition. Two friction damping models were considered namely, the blade-to-blade and the blade-damper-blade. The equations of motion derived on the basis of these models were solved by a method of harmonic balance, assuming, in effect, that under cyclic excitation the blades will exhibit cyclic response at the same frequency. The solutions for 8 blade, 16 blade, and 4 blade systems were obtained using the computer program BLADE. The levels of damping produced by the two friction damping models were compared and evaluated. The optimal values of the friction force, for which the tip amplitude of the blades had a minimum value, were determined.

Raju, B. B.↗

A unidimensional model of comet ionosphere structure

The one-dimensional continuity and momentum equations for the plasma formed in the expanding coma of a comet near the sun are studied. An analytic expression for the magnetic field configuration in the presence of outflow, photoionization, dissociative recombination, plasma fluid pressure, and friction between the ions and neutrals is obtained. It is suggested that for a Halley-type comet there will be a region sunward of the nucleus from which the magnetic field is excluded, consistent with Giotto observations. Calculations have been performed for Halley and Giacobini-Zinner type comets, and in the field-free regions it is shown that the dominant terms in the momentum equation balancing the magnetic pressure gradient are the ion neutral friction and the net mass loading momentum gain.

Eviatar, Aharon↗

Algorithms for parallel flow solvers on message passing architectures

The purpose of this project has been to identify and test suitable technologies for implementation of fluid flow solvers -- possibly coupled with structures and heat equation solvers -- on MIMD parallel computers. In the course of this investigation much attention has been paid to efficient domain decomposition strategies for ADI-type algorithms. Multi-partitioning derives its efficiency from the assignment of several blocks of grid points to each processor in the parallel computer. A coarse-grain parallelism is obtained, and a near-perfect load balance results. In uni-partitioning every processor receives responsibility for exactly one block of grid points instead of several. This necessitates fine-grain pipelined program execution in order to obtain a reasonable load balance. Although fine-grain parallelism is less desirable on many systems, especially high-latency networks of workstations, uni-partition methods are still in wide use in production codes for flow problems. Consequently, it remains important to achieve good efficiency with this technique that has essentially been superseded by multi-partitioning for parallel ADI-type algorithms. Another reason for the concentration on improving the performance of pipeline methods is their applicability in other types of flow solver kernels with stronger implied data dependence. Analytical expressions can be derived for the size of the dynamic load imbalance incurred in traditional pipelines. From these it can be determined what is the optimal first-processor retardation that leads to the shortest total completion time for the pipeline process. Theoretical predictions of pipeline performance with and without optimization match experimental observations on the iPSC/860 very well. Analysis of pipeline performance also highlights the effect of uncareful grid partitioning in flow solvers that employ pipeline algorithms. If grid blocks at boundaries are not at least as large in the wall-normal direction as those immediately adjacent to them, then the first processor in the pipeline will receive a computational load that is less than that of subsequent processors, magnifying the pipeline slowdown effect. Extra compensation is needed for grid boundary effects, even if all grid blocks are equally sized.

Vanderwijngaart, Rob F.↗

Path Toward a Unifid Geometry for Radiation Transport

The Direct Accelerated Geometry for Radiation Analysis and Design (DAGRAD) element of the RadWorks Project under Advanced Exploration Systems (AES) within the Space Technology Mission Directorate (STMD) of NASA will enable new designs and concepts of operation for radiation risk assessment, mitigation and protection. This element is designed to produce a solution that will allow NASA to calculate the transport of space radiation through complex computer-aided design (CAD) models using the state-of-the-art analytic and Monte Carlo radiation transport codes. Due to the inherent hazard of astronaut and spacecraft exposure to ionizing radiation in low-Earth orbit (LEO) or in deep space, risk analyses must be performed for all crew vehicles and habitats. Incorporating these analyses into the design process can minimize the mass needed solely for radiation protection. Transport of the radiation fields as they pass through shielding and body materials can be simulated using Monte Carlo techniques or described by the Boltzmann equation, which is obtained by balancing changes in particle fluxes as they traverse a small volume of material with the gains and losses caused by atomic and nuclear collisions. Deterministic codes that solve the Boltzmann transport equation, such as HZETRN [high charge and energy transport code developed by NASA Langley Research Center (LaRC)], are generally computationally faster than Monte Carlo codes such as FLUKA, GEANT4, MCNP(X) or PHITS; however, they are currently limited to transport in one dimension, which poorly represents the secondary light ion and neutron radiation fields. NASA currently uses HZETRN space radiation transport software, both because it is computationally efficient and because proven methods have been developed for using this software to analyze complex geometries. Although Monte Carlo codes describe the relevant physics in a fully three-dimensional manner, their computational costs have thus far prevented their widespread use for analysis of complex CAD models, leading to the creation and maintenance of toolkit-specific simplistic geometry models. The work presented here builds on the Direct Accelerated Geometry Monte Carlo (DAGMC) toolkit developed for use with the Monte Carlo N-Particle (MCNP) transport code. The workflow for achieving radiation transport on CAD models using MCNP and FLUKA has been demonstrated and the results of analyses on realistic spacecraft/habitats will be presented. Future work is planned that will further automate this process and enable the use of multiple radiation transport codes on identical geometry models imported from CAD. This effort will enhance the modeling tools used by NASA to accurately evaluate the astronaut space radiation risk and accurately determine the protection provided by as-designed exploration mission vehicles and habitats

Lee, Kerry↗

Vehicle-Level System Impact of Boundary Layer Ingestion for the NASA D8 Concept Aircraft

The purpose of this study was to evaluate the vehicle-level impact of a boundary layer ingestion (BLI) propulsion system on a commercial transport aircraft concept. The NASA D8 (ND8) aircraft was chosen as the BLI concept aircraft to be studied. A power balance methodology developed by the Massachusetts Institute of Technology was adapted for use with the existing NASA sizing and performance tools to model the fuel consumption impact of BLI on the ND8. A key assumption for the BLI impact assessment was a 3.5% efficiency penalty associated with designing a fan for and operating in the distorted flow caused by BLI. The ND8 was compared to several other ND8-like aircraft that did not utilize BLI in order to determine the fuel consumption benefit attributable to BLI. Analytically “turning off” BLI on the ND8 without accounting for the physical requirements of redirecting the boundary layer or resizing the aircraft to meet the performance constraints resulted in a 2.8% increase in block fuel consumption to fly the design mission. When this non-physical aircraft was resized to meet the performance constraints, the block fuel consumption was 4.0% greater than the baseline ND8. The ND8 was also compared to an ND8-like aircraft with conventionally podded engines under the wing. This configuration had a 5.6% increase in block fuel consumption compared to the baseline ND8. This result is more reflective of the real world impact if BLI is not an available technology for the ND8 design. The BLI benefit results presented for this study should not be applied to other aircraft that have a propulsion-airframe integration design or BLI implementation different from the ND8.

Marien, Ty V.↗

Accelerated Simulation of Air Pollution Using NVIDIA RAPIDS

Atmospheric chemistry models are a central tool to study and forecast the impact of air pollution on the environment, vegetation, and human health. However, the numerical simulation of chemical kinetics is computationally expensive due to the stiffness of the system of ordinary differential equations that describes atmospheric chemistry. Here we present an alternative approach to the computation of atmospheric chemistry based on machine learning. Our training data set is produced using the NASA Goddard Earth Observing System (GEOS) model with GEOS-Chem chemistry, run on the NASA Center for Climate Simulation (NCCS) Discover supercomputing cluster on 384 Intel Xeon Haswell cores. This model spends more than 50% of total run time on solving atmospheric chemistry. The data set contains as input features the air pollution concentrations before solving the differential equations, together with some key physical parameters such as temperature and sun intensity. As target variables we define the air pollution concentrations after solving the differential equations. Using Dask-cuDF and Dask-XGBoost on the NVIDIA RAPIDS platform on 8 Tesla V100 GPUs, we generate from this training set gradient boosted decision tree models that can reproduce the simulation of chemical kinetics. We do this on the NCCS Advanced Data Analytics Platform (ADAPT) science cloud environment. Our application takes full advantage of recent advances in Dask-XGBoost, such as multi-node and multi-GPU scaling for distributed training with large data sets. The increase in training data size enabled by this is critical to capture the full range of chemical environments encountered across the globe and all annual seasons.The boosted tree models offer good predictability and show many of the features of the full chemistry reference simulation. Further improvements can be achieved through mass balance considerations and by accounting for error correlations. We incorporate the boosted tree models into the GEOS reference model using XGBoost's C API. This enables a seamless integration of the GPU trained models into GEOS-Chem, which is written in Fortran and optimized for use in a massively parallel CPU environment. We show the benefits of this approach and discuss the potential speedup of this machine learning accelerated atmospheric chemistry model.

Keller, Christoph A.↗

Two-dimensional modelling of compressive failure in delaminated laminates

An analytical model is developed to assess the compressive strength criticality of near-surface interlaminar defects in laminated composites. The delaminated region is elliptic in shape, separating a thick isotropic plate from a thin orthotropic layer whose material axes coincide with the ellipse axes. The growth conditions and growth behavior of this defect are studied by breaking the overall problem into an elastic stability problem and a fracture problem. Post-buckling solution for the elliptic section is obtained using the Rayleigh-Ritz method while an energy balance criterion based on a self-similar disbond growth governs the fracture. The parameters controlling the growth or arrest of the delamination damage are identified as the fracture energy, disbond depth and elastic properties of the materials from both sides of the delaminating interface. By varying the degree of material anisotropy relative to the loading axis a range in growth behavior was found including stable or unstable crack growth parallel to or normal to the loading axis.

Chai, H.↗

Calculations and experiments concerning lifting force and power in TEMPUS

A critical comparison is reported between the theoretically predicted and experimentally measured values for the electromagnetic lifting force and the heating rates which may be achieved, under simulated microgravity conditions, using the TEMPUS electromagnetic levitation device. The experiments involved the suspending of a metallic sample from one arm of a recording balance, such that it was carefully positioned between the heating and the positioning coils of the levitation device. The net force exerted by the sample was measured as a function of position, the coil currents, and the nature of the sample. Some calculations are also reported regarding the power absorption by the sample. The theoretical predictions, based on the numerical solution of Maxwell's equations using the volume integral technique, were found to be in excellent agreement with the measurements. For the idealized case of a spherical sample, analytical solutions describing the lifting force were also found to agree very well with the computed results.

Zong, J. H.↗

Plans and Status of Wind-Tunnel Testing Employing an Aeroservoelastic Semispan Model

This paper presents the research objectives, summarizes the pre-wind-tunnel-test experimental results to date, summarizes the analytical predictions to date, and outlines the wind-tunnel-test plans for an aeroservoelastic semispan wind-tunnel model. The model is referred to as the Supersonic Semispan Transport (S4T) Active Controls Testbed (ACT) and is based on a supersonic cruise configuration. The model has three hydraulically-actuated surfaces (all-movable horizontal tail, all-movable ride control vane, and aileron) for active controls. The model is instrumented with accelerometers, unsteady pressure transducers, and strain gages and will be mounted on a 5-component sidewall balance. The model will be tested twice in the Langley Transonic Dynamics Tunnel (TDT). The first entry will be an "open-loop" model-characterization test; the second entry will be a "closed-loop" test during which active flutter suppression, gust load alleviation and ride quality control experiments will be conducted.

Perry, Boyd, III↗

Sensitivity analysis for the dynamic response of thermoviscoplastic shells of revolution

A computational procedure is presented for evaluating the sensitivity coefficients of the dynamic axisymmetric, fully-coupled, thermoviscoplastic response of shells of revolution. The analytical formulation is based on Reissner's large deformation shell theory with the effects of large-strain, transverse shear deformation, rotatory inertia and moments turning around the normal to the middle surface included. The material model is chosen to be viscoplasticity with strain hardening and thermal hardening, and an associated flow rule is used with a von Mises effective stress. A mixed formulation is used for the shell equations with the fundamental unknowns consisting of six stress resultants, three generalized displacements and three velocity components. The energy-balance equation is solved using a Galerkin procedure, with the temperature as the fundamental unknown. Spatial discretization is performed in one dimension (meridional direction) for the momentum and constitutive equations of the shell, and in two dimensions (meridional and thickness directions) for the energy-balance equation. The temporal integration is performed by using an explicit central difference scheme (leap-frog method) for the momentum equation; a predictor-corrector version of the trapezoidal rule is used for the energy-balance equation; and an explicit scheme consistent with the central difference method is used to integrate the constitutive equations. The sensitivity coefficients are evaluated by using a direct differentiation approach. Numerical results are presented for a spherical cap subjected to step loading.

Kulkarni, Makarand↗

Analyzing Responses of Chemical Sensor Arrays

NASA is developing a third-generation electronic nose (ENose) capable of continuous monitoring of the International Space Station s cabin atmosphere for specific, harmful airborne contaminants. Previous generations of the ENose have been described in prior NASA Tech Briefs issues. Sensor selection is critical in both (prefabrication) sensor material selection and (post-fabrication) data analysis of the ENose, which detects several analytes that are difficult to detect, or that are at very low concentration ranges. Existing sensor selection approaches usually include limited statistical measures, where selectivity is more important but reliability and sensitivity are not of concern. When reliability and sensitivity can be major limiting factors in detecting target compounds reliably, the existing approach is not able to provide meaningful selection that will actually improve data analysis results. The approach and software reported here consider more statistical measures (factors) than existing approaches for a similar purpose. The result is a more balanced and robust sensor selection from a less than ideal sensor array. The software offers quick, flexible, optimal sensor selection and weighting for a variety of purposes without a time-consuming, iterative search by performing sensor calibrations to a known linear or nonlinear model, evaluating the individual sensor s statistics, scoring the individual sensor s overall performance, finding the best sensor array size to maximize class separation, finding optimal weights for the remaining sensor array, estimating limits of detection for the target compounds, evaluating fingerprint distance between group pairs, and finding the best event-detecting sensors.

Zhou, Hanying↗

Comparison of big event with calculations of the air shower development

The incidence of high energy hadrons and electron-photons in air showers at various stages of development is calculated. Numerical calculation is used to solve the diffusion equation for a nuclear cascade and analytical calculation for cascade shower induced gamma rays. From these calculations, one can get the longitudinal development of the high energy hadron and electron-photon components, and the energy spectra of these components at various depths of air shower development. The total number of hadrons (N sub H) and electron-photon components (N sub gamma) are related according to stages of the air shower development and primary energy. The relation of the total energy of hadron and electron-photon component above the threshold energy is given. The energy balance between both components is also a useful parameter to study high energy events accompanying air showers. The relation of N sub H and fractional hadronic energy E (sum E sub H sup gamma/sum E sub H sup gamma + Sum E sub gamma) is calculated. This relation is helpful to understand the stage of air shower development(t) and primary energy (E sub p).

Niwa, M.↗