Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Functional coupling”

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 829 records · Page 46

Coupled evolution of the atmospheres and interiors of planets and satellites

The evolution of a planetary atmosphere can be powerfully influenced by the planetary interior's function as both a source and a sink of atmospheric constituents; the interior can in turn be strongly influenced by the atmosphere because the mechanism of interior heat loss depends on a volatile content for which the atmosphere can serve both as sink and source. The dependence of mantle rheology on volatile content could furnish a feedback mechanism tending to keep regassing/degassing in balance, thereby maintaining a relatively constant atmospheric mass. Consideration of the abundances of radiogenic and nonradiogenic noble gases in the earth's atmosphere, and of the fluxes of these gases from the mantle, support a large degassing event early on, followed by a decrease in degassing efficiency with time and relatively inefficient outgassing over most of geologic time.

Schubert, G.↗

Stationary turbulent closure via the Hopf functional equation

Explicit solutions of the stationary Hopf equation are discussed and their computational possibilities are explored. The motivation is to circumvent the infinite hierarchy of coupled equations for the velocity moments and obtain an exact closure of the steady-state 3D Navier-Stokes equations, without modeling assumptions or truncation. The Hopf formulation of the Navier-Stokes equation is reviewed. A stationary homogeneous solution for 2D flow is displayed and discussed. It is shown how depletion of nonlinearity may arise for 3D forced homogeneous flow. The general 3D forced case is considered and a method for closing the 3D unforced equations with arbitrary boundary conditions is derived.

Shen, Hubert H.↗

Supersonic Testing of 0.8 m Disk Gap Band Parachutes in the Wake of a 70 Deg Sphere Cone Entry Vehicle

Supersonic wind tunnel testing of Viking-type 0.8 m Disk-Gap-Band (DGB) parachutes was conducted in the NASA Glenn Research Center 10'x10' wind-tunnel. The tests were conducted in support of the Mars Science Laboratory Parachute Decelerator System development and qualification program. The aerodynamic coupling of the entry-vehicle wake to parachute flow-field is under investigation to determine the cause and functional dependence of a supersonic canopy breathing phenomenon referred to as area oscillations, characteristic of DGB's above Mach 1.5 operation. Four percent of full-scale parachutes (0.8 m) were constructed similar to the flight-article in material and construction techniques. The parachutes were attached to a 70-deg sphere-cone entry-vehicle to simulate the Mars flight configuration. The parachutes were tested in the wind-tunnel from Mach 2 to 2.5 in a Reynolds number range of 2x105 to 1x106, representative of a Mars deployment. Three different test configurations were investigated. In the first two configurations, the parachutes were constrained horizontally through the vent region to measure canopy breathing and wake interaction for fixed trim angles of 0 and 10 degrees from the free-stream. In the third configuration the parachute was unconstrained, permitted to trim and cone, similar to free-flight (but capsule motion is constrained), varying its alignment relative to the entry-vehicle wake. Non-intrusive test diagnostics were chosen to quantify parachute performance and provide insight into the flow field structure. An in-line loadcell provided measurement of unsteady and mean drag. Shadowgraph of the upstream parachute flow field was used to capture bow-shock motion and wake coupling. Particle image velocimetry provided first and second order flow field statistics over a planar region of the flow field, just upstream of the parachute. A photogrammetric technique was used to quantify fabric motion using multiple high speed video cameras to record the location in time and space of reflective targets placed on the canopy interior. The experimental findings including an updated drag model and the physical basis of the area oscillation phenomenon will be discussed.

EDL↗

Machine-Learned Committor Functions for Reactive Molecular Dynamics

Reactive molecular dynamics (MD) is a powerful tool for atomistic-scale modeling of a diverse range of chemical processes. However, scaling these simulations to large systems and long times scales remains a challenge because of the complexity of the potential energy function required. The authors previously developed a heuristic approach, called REACTER, that incorporates reactivity in MD simulations in a less general but much more computationally efficient manner. REACTER uses standard, fixed valence force fields as the underlying potentialenergy surface for describing all interatomic interactions but adds a procedure for enforcing user-defined reactions that occur when certain geometric constraints on relative atomic positions are satisfied. Further, these bonding changes can be accepted or rejected with a probability related tothe local thermal energy. This work seeks to generalize this approach by replacing the set of user defined geometric constraints and energetic criteria with a committor function that specifies the probability of a reaction occurring on the basis of the local atomic configuration. The committor function is a useful mathematical tool for modeling rare events but, unfortunately, is very difficult to compute for realistic systems in a general way. This work describes a method for approximating the committor function using a machine learning approach, specifically a deep neural network trained with data from reactive MD and DFT-based dynamics simulations. This network is coupled to the existing REACTER protocol, as implemented in the LAMMPS MD package, and used to make on-the-fly predictions of reaction probabilities without the more extensive user input previously required. The new method is demonstrated using the polymerization of polystyrene as a case study. Although very dependent on the quality and quantity of training data, machine-learned committor functions show promise as a method for incorporating reaction probability from higher level calculations into highly scalable MD simulations.

polymer simulations↗

Radiation Losses Due to Tapering of a Double-Core Optical Waveguide

The theoretical model we designed parameterizes the power losses as a function of .the profile shape for a tapered, single mode, optical dielectric coupler. The focus of this project is to produce a working model that determines the power losses experienced by the fibers when light crosses a taper region. This phenomenon can be examined using coupled mode theory. The optical directional coupler consists of a parallel, dual-channel, waveguide with minimal spacing between the channels to permit energy exchange. Thus, power transfer is essentially a function of the taper profile. To find the fields in the fibers, the approach used was that of solving the Helmholtz equation in cylindrical coordinates involving Bessel and modified Bessel functions depending on the location.

Lyons, Donald R.↗

Atmospheric Depolarization Lidar Experimental Receiver: A Space Shuttle Hitchhiker Payload

Development work is underway at the Goddard Space Flight Center to construct a depolarization measuring atmospheric lidar receiver. The Atmospheric Lidar (AL) is tentatively scheduled to fly on the Space Shuttle in a late-1999 time frame. The AL will fly in conjunction with the Shuttle Laser Altimeter (SLA) and the Infrared Spectral Imaging Radiometer (ISIR) to provide a comprehensive package of atmospheric aerosol and cloud information. The AL operates in conjunction with the SLA laser transmitter and measures profiles of atmospheric backscatter at 532 nm. The receiver system discriminates between the parallel and perpendicular polarizations of the backscattered signal, thus providing depolarization ratios for scattering from clouds. The lidar receiver also provides cloud height and thickness measurements to complement the brightness temperature measurements generated by the ISIR thermal imager. The function of the AL is twofold. The primary function is to provide range-resolved measurements of atmospheric aerosol backscatter and depolarization ratio at 532 nm with 75 m vertical resolution. The scientific purpose of these measurements is to determine composition of clouds based on the depolarization ratio (i.e cloud content is water or ice), to determine cloud height and thickness, and to gain further understanding of the global distribution of aerosols. This information, when coupled with the cloud brightness measured by the ISIR thermal imager will provide a significant amount of information on cloud composition and radiative effects, particularly for cirrus and sub-visual cirrus clouds. A secondary function of the AL is to serve as an in-space test bed for lidar technology advancements, including a fully fiber-coupled receiver and photon counting from space. In addition, the data obtained by the AL will be used to develop software for the Geoscience Laser Altimeter System (GLAS) flight mission.

McGill, Matthew J.↗

Shield to Pin Coupling of Lightning-Like Transients on Payload Umbilical Cables on a Launch Pad

In this paper we describe in-situ testing of a long payload umbilical, on a launch site, injected with “lightning- like” transients and describe resulting pin-to-pin voltages. Injections and voltage measurements near the ground support equipment room, as well as at a location near the payload junction box, are made. The umbilical cables tested include an outer over-braid and the inner conductor coupling is examined for open circuit, short-circuit and various loads representative of spacecraft input impedances. This testing is important because the Kennedy Space Center (KSC) where the lightning occurrence is the highest in the United States, is the primary launch site for Launch Services Program spacecraft customers. Lightning planning is essential but developing a lightning plan is often overlooked or not adequately analyzed leaving the spacecraft vulnerable to time delays or even damage when lightning occurs. At other popular launch sites like Vandenberg Air Force Base (VAFB) where lightning occurs less often, although at the same or greater intensity when it does occur, lightning planning is often completely ignored by the spacecraft. The two major questions to be addressed in the lightning plan are what retesting should be done to establish a “goodness” level and what is the trigger criteria for this testing? The spacecraft will typically use a standard spacecraft check-out procedure to address the necessary retesting, but determining the trigger criteria is often an issue. For instance, a spacecraft needs to understand what their immunity is to a certain lightning magnitude and location. Determining the amount of current that can be coupled onto a spacecraft umbilical can be calculated by using worst case assumptions or measured with current probes and current measurement devices. Spacecraft can also determine what pin-to-pin voltages they are sensitive to, however pin-to-pin voltage measurements are not typically taken during the strike due to the invasive nature of this measurement. In this paper, we present detailed data on the shield to pin voltage transfer functions to provide insight to the spacecraft developers for lightning retest criteria planning. The results from this unique testing opportunity provide essential details on specific coupling mechanisms affecting spacecraft hardware that interfaces with the ground support equipment. This missing link between cable shield currents and payload susceptibility voltages has been methodically tested and representative data presented.

Trout, Dawn↗

Implicit Coupling Approach for Simulation of Charring Carbon Ablators

This study demonstrates that coupling of a material thermal response code and a flow solver with nonequilibrium gas/surface interaction for simulation of charring carbon ablators can be performed using an implicit approach. The material thermal response code used in this study is the three-dimensional version of Fully Implicit Ablation and Thermal response program, which predicts charring material thermal response and shape change on hypersonic space vehicles. The flow code solves the reacting Navier-Stokes equations using Data Parallel Line Relaxation method. Coupling between the material response and flow codes is performed by solving the surface mass balance in flow solver and the surface energy balance in material response code. Thus, the material surface recession is predicted in flow code, and the surface temperature and pyrolysis gas injection rate are computed in material response code. It is demonstrated that the time-lagged explicit approach is sufficient for simulations at low surface heating conditions, in which the surface ablation rate is not a strong function of the surface temperature. At elevated surface heating conditions, the implicit approach has to be taken, because the carbon ablation rate becomes a stiff function of the surface temperature, and thus the explicit approach appears to be inappropriate resulting in severe numerical oscillations of predicted surface temperature. Implicit coupling for simulation of arc-jet models is performed, and the predictions are compared with measured data. Implicit coupling for trajectory based simulation of Stardust fore-body heat shield is also conducted. The predicted stagnation point total recession is compared with that predicted using the chemical equilibrium surface assumption

Ablation↗

Transverse Mode Dynamics of VCSELs Through Space-Time Domain Simulation

Modeling and simulation are important to understand laser operation and to optimize and design device functions. Numerical simulation of VCSEL (Vertical Cavity Surface Emitting Lasers) has been largely based on solving time-independent Helmholtz equation or time dependent coupled mode equations. There are various advantages for choosing these approaches. However, the disadvantages are also apparent. The former cannot handle dynamical mode competition seen in VCSELs, while the latter assumes a given type and number of modes a priori. Furthermore, the microscopic physics of heterstructures and electron-hole plasma is very often represented by a few parameters such as linear gain coefficients and the linewidth enhancement factor. These are over simplification of space and frequency (wavelength) dependent gain and refractive index functions. When the space-time dynamical operation of VCSELs becomes important, these simple approximations become questionable. In this paper, we apply a recently developed model for edge-emitting lasers to a gain guided VCSEL for space-time domain simulation. This model takes into account the actual nonlinear dependence of gain and refractive index on frequency and carrier density within the frame work of the effective Bloch equations. The corresponding partial differential equations are solved directly by finite difference methods. Laser behavior with increasing pumping current is investigated in detail. Special attention is paid to the dynamical competition of the transverse modes.

Goorjian, Peter M.↗

Oxidation Behavior of Carbon Fiber-Reinforced Composites

OXIMAP is a numerical (FEA-based) solution tool capable of calculating the carbon fiber and fiber coating oxidation patterns within any arbitrarily shaped carbon silicon carbide composite structure as a function of time, temperature, and the environmental oxygen partial pressure. The mathematical formulation is derived from the mechanics of the flow of ideal gases through a chemically reacting, porous solid. The result of the formulation is a set of two coupled, non-linear differential equations written in terms of the oxidant and oxide partial pressures. The differential equations are solved simultaneously to obtain the partial vapor pressures of the oxidant and oxides as a function of the spatial location and time. The local rate of carbon oxidation is determined at each time step using the map of the local oxidant partial vapor pressure along with the Arrhenius rate equation. The non-linear differential equations are cast into matrix equations by applying the Bubnov-Galerkin weighted residual finite element method, allowing for the solution of the differential equations numerically.

Sullivan, Roy M.↗

Photosynthesis-dependent Isoprene Emission from Leaf to Planet in a Global Carbon-chemistry-climate Model

We describe the implementation of a biochemical model of isoprene emission that depends on the electron requirement for isoprene synthesis into the FarquharBallBerry leaf model of photosynthesis and stomatal conductance that is embedded within a global chemistry-climate simulation framework. The isoprene production is calculated as a function of electron transport-limited photosynthesis, intercellular and atmospheric carbon dioxide concentration, and canopy temperature. The vegetation biophysics module computes the photosynthetic uptake of carbon dioxide coupled with the transpiration of water vapor and the isoprene emission rate at the 30 min physical integration time step of the global chemistry-climate model. In the model, the rate of carbon assimilation provides the dominant control on isoprene emission variability over canopy temperature. A control simulation representative of the present-day climatic state that uses 8 plant functional types (PFTs), prescribed phenology and generic PFT-specific isoprene emission potentials (fraction of electrons available for isoprene synthesis) reproduces 50 of the variability across different ecosystems and seasons in a global database of 28 measured campaign-average fluxes. Compared to time-varying isoprene flux measurements at 9 select sites, the model authentically captures the observed variability in the 30 min average diurnal cycle (R2 6496) and simulates the flux magnitude to within a factor of 2. The control run yields a global isoprene source strength of 451 TgC yr1 that increases by 30 in the artificial absence of plant water stress and by 55 for potential natural vegetation.

photosynthesis↗

A double-loop tracking system.

A nonlinear analysis which can be used to assess certain statistical characteristics of double-loop tracking systems is presented. It takes into account the mutual coupling effects of the loops in the system. Two approaches are taken to obtain steady-state probability density functions (pdf's) of the system phase errors. From these pdf's, important system performance statistics, e.g., the phase-error variances, can be calculated, thus illustrating the application and usefulness of the analysis. The analysis is applied to a satellite transponder as an example.

Yuen, J. H.↗

Two-dimensional turbulence in inviscid fluids or guiding center plasmas

Numerical tests of the inviscid equilibrium theory of Kraichnan (1975) are described. The mathematical description applies equally well to the two-dimensional electrostatic guiding center plasma and to the two-dimensional inviscid Navier-Stokes fluid. The predictions of this analytic theory are discussed. A pair of coupled equations are derived for the two-time vorticity autocorrelation in Fourier space and the infinitesimal unit response function G of Kraichnan in the so-called Eulerian direct interaction approximation. Kraichnan's assertion that thermodynamic limits exist for the negative temperature states is questioned.

Seyler, C. E., Jr.↗

Meteorological radar facility (MRF) slot conductance investigations

A preliminary meteorological radar facility (MRF) array design was completed in support of the slot conductance measurement program. Three different slot measurement techniques were evaluated. The selection of the probe comparison measurement technique was selected as the principal experimental method with the impedance measurement technique chosen to measure a few higher conductance slots to be used as reference slots. The impedance of 43 slots in 0.9 x 0.4 inch standard waveguide and of 40 slots in 0.835 x 0.4 inch waveguide was measured. Also, impedance measurements were made of a few slots using image planes to simulate mutual coupling effects. The measured and theoretical conductance, susceptance, and radiation phase data are presented in graphic form as a function of slot displacement for constant slot length, and of slot length for constant slot displacement. It is concluded that the proposed MRF array design approach is a feasible one.

Ratkevich, A. E.↗

Edge softening of the Shuttle TPS strain isolation pad

Tensile tests and an analytical investigation were performed to characterize the edge softening behavior of the strain isolation pad (SIP) between the Orbiter skin and thermal protection system. The tensile tests were carried out with varying sizes of disk-shaped specimens bonded between aluminum disks. The specimens strength and stiffness were determined on the basis of specimen size, and an analytical model of the microstructural stress-strain characteristics was developed. Strength and stiffness were found to decrease near the free edges because through-the-thickness fibers located there were not anchored. No size dependence at maximum load was observed in specimens between 0.75-4.0 in. thick. In-plane and out-of-plane coupling in deformation was detected. The model gave accurate predictions of the tensile behavior of the SIP as a function of distance to a free edge.

Ransone, P. O.↗

A Galerkin method for the estimation of parameters in hybrid systems governing the vibration of flexible beams with tip bodies

An approximation scheme is developed for the identification of hybrid systems describing the transverse vibrations of flexible beams with attached tip bodies. In particular, problems involving the estimation of functional parameters are considered. The identification problem is formulated as a least squares fit to data subject to the coupled system of partial and ordinary differential equations describing the transverse displacement of the beam and the motion of the tip bodies respectively. A cubic spline-based Galerkin method applied to the state equations in weak form and the discretization of the admissible parameter space yield a sequence of approximating finite dimensional identification problems. It is shown that each of the approximating problems admits a solution and that from the resulting sequence of optimal solutions a convergent subsequence can be extracted, the limit of which is a solution to the original identification problem. The approximating identification problems can be solved using standard techniques and readily available software.

Banks, H. T.↗

Theoretical study of the X1Sigma(+) states of the alkali hydrides NaH-CsH

By means of near Hartree-Fock quality Slater basis sets, and the incorporation of electron correlation through the coupled-pair formalism, theoretical potentials are obtained for the X1Sigma(+) states of NaH, KH, and RbH. Electric dipole moment functions are given for NaH-RbH, as well as vibrationally averaged dipole moments, Einstein coefficients, and radiative lifetimes for the first 10 vibrational levels; an extensive study is made of the computational requirements for an accurate permanent dipole moment of KH.

Langhoff, Stephen R.↗

Neural network based decomposition in optimal structural synthesis

The present paper describes potential applications of neural networks in the multilevel decomposition based optimal design of structural systems. The generic structural optimization problem of interest, if handled as a single problem, results in a large dimensionality problem. Decomposition strategies allow for this problem to be represented by a set of smaller, decoupled problems, for which solutions may either be obtained with greater ease or may be obtained in parallel. Neural network models derived through supervised training, are used in two distinct modes in this work. The first uses neural networks to make available efficient analysis models for use in repetitive function evaluations as required by the optimization algorithm. In the second mode, neural networks are used to represent the coupling that exists between the decomposed subproblems. The approach is illustrated by application to the multilevel decomposition-based synthesis of representative truss and frame structures.

Hajela, P.↗