Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “forward modeling”

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 1,747 records · Page 97

Comparison of wind tunnel and flight test afterbody and nozzle pressures for a twin-jet fighter aircraft at transonic speeds

Afterbody and nozzle pressures measured on a 1/12-scale model and in flight on a twin-jet fighter aircraft were compared as Mach number varied from 0.6 to 1.2, Reynolds number from 17.5 million to 302.5 million, and angle of attack from 1 to 7 deg. At Mach 0.6 and 0.8, nozzle pressure coefficient distributions and nozzle axial force coefficients agreed and showed good recompression. At Mach 0.9 and 1.2, flow complexity caused a loss in recompression for both flight and wind tunnel nozzle data. The flight data exhibited less negative values of pressure coefficient and lower axial force coefficients than did the wind tunnel data. Reynolds number effects were noted only at these Mach numbers. Jet temperature and mass flux ratio did not affect the comparisons of nozzle axial flow coefficient. At subsonic speeds, the levels of pressure coefficient distributions on the upper fuselage and lower nacelle surfaces for flight were less negative than those for the model. The model boundary layer thickness at the aft rake station exceeded that for the forward rake station and increased with increasing angle of attack. The flight boundary layer thickness at the aft rake station was less than that for the forward rake station and decreased with increasing angle of attack.

Nugent, Jack↗

The Communication Link and Error ANalysis (CLEAN) simulator

During the period July 1, 1993 through December 30, 1993, significant developments to the Communication Link and Error ANalysis (CLEAN) simulator were completed and include: (1) Soft decision Viterbi decoding; (2) node synchronization for the Soft decision Viterbi decoder; (3) insertion/deletion error programs; (4) convolutional encoder; (5) programs to investigate new convolutional codes; (6) pseudo-noise sequence generator; (7) soft decision data generator; (8) RICE compression/decompression (integration of RICE code generated by Pen-Shu Yeh at Goddard Space Flight Center); (9) Markov Chain channel modeling; (10) percent complete indicator when a program is executed; (11) header documentation; and (12) help utility. The CLEAN simulation tool is now capable of simulating a very wide variety of satellite communication links including the TDRSS downlink with RFI. The RICE compression/decompression schemes allow studies to be performed on error effects on RICE decompressed data. The Markov Chain modeling programs allow channels with memory to be simulated. Memory results from filtering, forward error correction encoding/decoding, differential encoding/decoding, channel RFI, nonlinear transponders and from many other satellite system processes. Besides the development of the simulation, a study was performed to determine whether the PCI provides a performance improvement for the TDRSS downlink. There exist RFI with several duty cycles for the TDRSS downlink. We conclude that the PCI does not improve performance for any of these interferers except possibly one which occurs for the TDRS East. Therefore, the usefulness of the PCI is a function of the time spent transmitting data to the WSGT through the TDRS East transponder.

Ebel, William J.↗

Low gravity transfer line chilldown

A code has been developed that solves for the transfer line chilldown time and flow and heat transfer characteristics in one-g environment. The code solves the transient, one dimensional, space averaged mass, momentum, and energy conservation equations for liquid-vapor two-phase flow in tubes. The physical configuration solved is that appropriate for bottom coolant injection in a vertically supported heated tube. Four distinct regions are considered consecutively: fully liquid, inverted annular, dispersed, and fully vapor flow. The conservation equations for both the liquid and the vapor are solved in each region separately. Also, in each region the mass and energy transport between each phase as well as the energy and momentum transport between the tube wall and the fluid are accounted for. A finite wall thickness is also considered. The model described above was solved numerically through a mixed finite difference scheme with forward time marching. The inverted annular regime was resolved using a semi-implicit finite differencing while the dispersed regime was solved explicitly. Also, a staggered mesh was used in which the velocity was resolved at mesh boundaries while all other field variables were resolved at the mesh centroids. Different mesh sizes were used depending on the region of solution. A coarse mesh was used in the dispersed flow region while a much finer mesh was used in both the inverted annular flow region and the tube walls.

B. N. Antar↗

Recent National Transonic Facility Test Process Improvements (Invited)

This paper describes the results of two recent process improvements; drag feed-forward Mach number control and simultaneous force/moment and pressure testing, at the National Transonic Facility. These improvements have reduced the duration and cost of testing. The drag feed-forward Mach number control reduces the Mach number settling time by using measured model drag in the Mach number control algorithm. Simultaneous force/moment and pressure testing allows simultaneous collection of force/moment and pressure data without sacrificing data quality thereby reducing the overall testing time. Both improvements can be implemented at any wind tunnel. Additionally the NTF is working to develop and implement continuous pitch as a testing option as an additional method to reduce costs and maintain data quality.

Kilgore, W. A.↗

Failure Assessment

Three questions to which software developers want accurate, precise answers are "How can the software system fail?", "mat bad things will happen if the software fails?t', and "How many failures will the software experience?". Numerous techniques have been devised to answer these questions; three of the best known are: 1) Software Fault Tree Analysis (SFTA) 2) Software Failure Modes, Effects, and Criticality Analysis (SFMECA 3) Software Fault/Failure Modeling. SFTA and SFMECA have been successfully used to analyze the flight software for a number of robotic planetary exploration missions, including Galileo, Cassini, and Deep Space 1. Given the increasing interest in reusing software components from mission to mission, one of us has developed techniques for reusing the corresponding portions of the SFTA and SFMECA, reducing the effort required to conduct these analyses. SFTA has also been shown to be effective in analyzing the security aspects of software systems; intrusion mechanisms and effects can easily be modeled using these techniques. The Bi- Directional Safety Analysis (BDSA) method combines a forward search (similar to SFMECA) from potential failure modes to their effects, with a backward search (similar to SFTA) from feasible hazards to the contributing causes of each hazard. BDSA offers an efficient way to identify latent failures. Recent work has extended BDSA to product-line applications such as flight-instrumentation displays and developed tool support for the reuse of the failure-analysis artifacts within a product line. BDSA has also been streamlined to support those projects having tight cost and/or schedule constraints for their failure analysis efforts. We discuss lessons learned from practice, describe available tools, and identi@ some future directions for the topic. A substantial amount of research has been devoted to estimating the number of failures that a software system will experience during test and operations, as well as the number of faults that have been inserted into that system during its development. One of us has found that the amount of structural change to a system during its development is strongly related to the number of faults inserted into it. Using techniques requiring no additional effort on the part of the development organization, the required measurements of structural evolution can be easily obtained from a development effort's configuration management system and readily transformed into an estimate of fault content. So far, structure-fault relationships have been identified for source code; current work seeks to examine artifacts available earlier in the lifecycle to determine if similar relationships between structure and fault content can be found. In particular, relationships between requirements change requests and the number of faults inserted into the implemented system would provide a significant improvement in our ability to control software quality during the early development phases.

fault tree↗

Sensitivity of Space Launch System Buffet Forcing Functions to Buffet Mitigation Options

Time-varying buffet forcing functions arise from unsteady aerodynamic pressures and are one of many load environments, which contribute to the overall loading condition of a launch vehicle during ascent through the atmosphere. The buffet environment is typically highest at transonic conditions and can excite the vehicle dynamic modes of vibration. The vehicle response to these buffet forcing functions may cause high structural bending moments and vibratory environments, which can exceed the capabilities of the structure, or of vehicle components such as payloads and avionics. Vehicle configurations, protuberances, payload fairings, and large changes in stage diameter can trigger undesirable buffet environments. The Space Launch System (SLS) multi-body configuration and its structural dynamic characteristics presented challenges to the load cycle design process with respect to buffet-induced loads and responses. An initial wind-tunnel test of a 3-percent scale SLS rigid buffet model was conducted in 2012 and revealed high buffet environments behind the booster forward attachment protuberance, which contributed to reduced vehicle structural margins. Six buffet mitigation options were explored to alleviate the high buffet environments including modified booster nose cones and fences/strakes on the booster and core. These studies led to a second buffet test program that was conducted in 2014 to assess the ability of the buffet mitigation options to reduce buffet environments on the vehicle. This paper will present comparisons of buffet forcing functions from each of the buffet mitigation options tested, with a focus on sectional forcing function rms levels within regions of the vehicle prone to high buffet environments.

Piatak, David J.↗

NASA ACES V&V Alignment Briefing

Regulations to establish operational and performance requirements for unmanned aircraft systems (UAS) are being developed by a consortium of government, industry and academic institutions. Those requirements will apply to the new detect and avoid (DAA) systems and other equipment necessary to integrate UAS with the National Airspace System (NAS) and are determined according to their contribution to the overall safety case for such an integration. This briefing focuses on providing an overview of the Airspace Concept Evaluation System (ACES) platform, review of detect-and-avoid models incorprated in ACES, sumamry of two planned ACES studies, and a way forward to impact the SC-228 VV plan.

Santiago, Confesor↗

Using Virtual Reality and Motion Capture as tools for Human Factors Engineering at NASA Marshall Space Flight Center

NASA Marshall Space Flight Center (MSFC) Human Factors Engineering (HFE) Team is implementing virtual reality (VR) and motion capture (MoCap) into HFE analyses of various projects through its Virtual Environments Lab (VEL). VR allows for multiple analyses early in the design process and more opportunities to give design feedback. This tool can be used by engineers in most disciplines to compare design alternatives and is particularly valuable to HFE to give early input during these evaluations. These techniques are being implemented for concept development of Deep Space Habitats (DSH), and work is being done to implement VR for design aspects of the Space Launch System (SLS). VR utilization in the VEL will push the design to be better formulated before mockups are constructed, saving budget and time. The MSFC VEL will continue forward leaning implementation with VR technologies in these and other projects for better models earlier in the design process.

Andrews, Tanya↗

Parallelized Quadrupole Simulations of Thermographic Responses of Composites

Thermography has been shown to be a viable technique for inspection of composites. Model inversion of the thermography data requires a fast method for performing the forward problem. Viable numerical methods for the thermal response forward problem are finite element, finite difference and the quadrupole method. Normally both the finite element and finite difference methods solve for the thermal response in the time domain which limits one’s ability to increase the speed of the simulation by parallelization. In contrast, the quadrupole method solves for the Laplace transform of the thermal response. One of the features of the Laplace transform methodology is the solution at any discrete time is independent of the solution at all other times. Therefore, it is easy to separate into a set of independent calculations with each of the times of interest being performed in parallel. Additionally, the numeric inversion of the Laplace transform typically involves numerically solving for the Laplace transform at multiple Laplace frequencies. Each of those solutions are also independent of solutions at other frequencies and can be calculated in parallel. By parallelization of this method, it is possible to perform the simulations of three-dimensional configurations in seconds. When the input stimulus for thermal response is a delta function heat flux (a reasonable approximation for flash heating), the thermal response is smooth. For this case, it is possible to accurately estimate the thermal response at any time within a given time interval from a set of simulations separated by exponentially increasing time steps. From these simulations, it is possible to accurately interpolate to find the response at intermediate times by a spline interpolation of the logarithm of time versus logarithm of temperature. The thermal response with exponential time stepping is shown to produce values for the thermal response which are within 1% of values within the time interval. The simulations are compared to finite element simulations of the same inspection configurations. The simulations are also compared to the thermographic measurements on composites where shape and depth of the delaminations are obtained from other inspection methods.

Thermography↗

Parallelized Quadrupole Simulations of Thermographic Responses of Composites

Thermography has been shown to be a viable technique for inspection of composites. Model inversion of the thermography data requires a fast method for performing the forward problem. Viable numerical methods for the thermal response forward problem are finite element, finite difference and the quadrupole method. Normally both the finite element and finite difference methods solve for the thermal response in the time domain which limits one’s ability to increase the speed of the simulation by parallelization. In contrast, the quadrupole method solves for the Laplace transform of the thermal response. One of the features of the Laplace transform methodology is the solution at any discrete time is independent of the solution at all other times. Therefore, it is easy to separate into a set of independent calculations with each of the times of interest being performed in parallel. Additionally, the numeric inversion of the Laplace transform typically involves numerically solving for the Laplace transform at multiple Laplace frequencies. Each of those solutions are also independent of solutions at other frequencies and can be calculated in parallel. By parallelization of this method, it is possible to perform the simulations of three-dimensional configurations in seconds. When the input stimulus for thermal response is a delta function heat flux (a reasonable approximation for flash heating), the thermal response is smooth. For this case, it is possible to accurately estimate the thermal response at any time within a given time interval from a set of simulations separated by exponentially increasing time steps. From these simulations, it is possible to accurately interpolate to find the response at intermediate times by a spline interpolation of the logarithm of time versus logarithm of temperature. The thermal response with exponential time stepping is shown to produce values for the thermal response which are within 1% of values within the time interval. The simulations are compared to finite element simulations of the same inspection configurations. The simulations are also compared to the thermographic measurements on composites where shape and depth of the delaminations are obtained from other inspection methods.

Thermography↗

Photoproduction of K0K0-bar with the GlueX Experiment

In this dissertation, we study photoproduction of K0 ?K 0 through the KSKLp and KSKSp ?nal states with the GlueX Phase-I (GlueX-I) data set. We measure the spin-density matrix elements (SDMEs) and di?erential cross section of ?(1020) ! KSKL to better understand photoproduction of light vector mesons. A mass independent Partial Wave Analysis (PWA) of the KSKL and KSKS meson spectrum below 2 GeV is also undertaken to study the meson spectrum with JPC = even++ and odd??. The ?(1020) di?erential cross section is measured at E = 8:2 ? 8:8 GeV and ?t = 0:15 ? 1:00 GeV2. The di?erential cross section is well described by an exponential decay with slope 4:44?0:01 GeV?2 and the integrated cross section is determined to be 295:7?0:4 nb, only statistical uncertainties are quoted. Both measurements are consistent and far more precise than the previous measurement by Ballam et al. [1]. The ?(1020) SDMEs were measured in nine bins of ?t in the same range. At low ?t, we ?nd the data were consistent with s-channel helicity conservation, SCHC, i. e. the only non-zero SDMEs were ?1 1?1 = ?Im(?2 1?1) = 1=2. At higher ?t, the SDMEs deviate from SCHC and the measurements indicate this is due to natural parity exchange since we observe that the SDMEs are consistent with zero unnatural exchange. We also ?find that the contribution from helicity double-flip amplitudes is consistent with zero. The measured SDMEs are in poor agreement with theoretical predictions put forward by JPAC [35]. These measurements will serve as input to re?ne models of production processes, which will be essential for the interpretation of possible signals of exotic mesons in GlueX. Analysis of the spectrum above the ?(1020) indicates the presence of at least three particles at ? 1:50, ? 1:75, and 2:20 GeV. We employ simple parametrizations of the resonance line shape based on relativistic Breit-Wigner functions considering cases with and without interference. The resonance at ? 1:50 GeV is not identifi?ed as any specifi?c resonance, but could be due to interference between the ?(1450) and !(1420). Our models favor the resonance at ~1:75 GeV as the X(1750) rather than the ?(1680). We ?find that introducing a third Breit-Winger into the models improves the ?t quality with the ?2/ndf going from 1.75 to 1.59 and 1.41 to 1.18 for models without and with interference respectively. A PWA of the KSKL spectrum indicates that the spectrum is consistent with exclusively spin-1 contribution up to ? 1:6 GeV. Above ? 1:6 GeV we ?find evidence for a small spin-3 contribution. We also ?nd that the spectrum predominantly positive reflectivity up to ? 1:6 GeV, after which the spectrum becomes a nearly equal mix of both. The KSKS system shows a rich spectrum with multiple resonance structures but is more statistically limited than the KSKL system. Modelling the line shape we ?find evidence at greater than 10? level in favor of a resonance at ? 1:75 GeV. The model parameters suggest this resonance is the f0(1710). The PWA of the KSKS system was inconclusive. However, using a small set of amplitudes suggests that the spin-2 contribution to the spectrum is primarily in the range 1:2?1:6 GeV, where f2(1270), a2(1320) and f02 (1525) are expected.

Linera, Gabriel Rodriguez↗

Measurement of the double-differential cross section of muon-neutrino charged-current interactions with low hadronic energy in the NOvA Near Detector

The NOvA collaboration reports cross-section measurements for $\nu_{\mu}$ charged-current interactions with low hadronic energy (maximum kinetic energy of 250 MeV for protons and 175 MeV for pions) in the NOvA Near Detector. The results are presented as a double-differential cross section as a function of the direct observables of the final-state muon kinematics. Results are also presented as a single-differential cross section as a function of the derived square of the four-momentum transfer, $Q^{2}$, and as a function of the derived neutrino energy. The data correspond to an accumulated 8.09$\times10^{20}$ protons-on-target (POT) in the neutrino mode of the NuMI beam, with a narrow band of neutrino energies peaked at 1.8 GeV. The analysis provides a sample of neutrino-nucleus interactions with an enhanced fraction of quasi-elastic and two-particle-two-hole (2p2h) interactions. This enhancement allows quantitative comparisons with various nuclear models. We find strong disagreement between data and theory-based models in various regions of the muon kinematic phase space, especially in the forward muon direction.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Reverse-thrust technology for variable-pitch fan propulsion systems

Tests conducted to develop the technology necessary to meet the unique reverse-thrust performance requirements of a variable pitch fan propulsion system are discussed. The losses and distortion associated with the air entering the fan and core compressor from the rear of the engine, the direction of fan blade pitch rotation for best reverse-thrust aeroacoustic performance, and engine response and operating characteristics during forward- to reverse-thrust transients are among the factors studied. The test results of several scale fan models as well as a full-size variable pitch fan engine are summarized. Results show the following: a flared exhaust nozzle makes a good reverse-thrust inlet; acceptable core inlet duct recovery and distortion levels in reverse flow were demonstrated; adequate thrust levels were achieved; forward- to reverse-thrust response time achieved was better than the goal; thrust and noise levels strongly favor reverse through feather pitch; and finally, flight-type inlets make the establishment of reverse flow more difficult.

Sagerser, D. A.↗

Formation of slow shock pairs associated with coronal mass ejections

The formation of a forward-reverse slow shock pair in the solar corona is presently simulated by an MHD model that uses the Rankine-Hugoniot solution to calculate the flow-property jumps at all shock crossings. The shocks divide the solution-domain into several continuous flow regions whose respective governing characteristics are solved by the method of characteristics. The plasma impact compresses the plasma near the front of the coronal mass ejection (CME); as the CME-associated slow shock pair moves outwards in interplanetary space, it evolves into a pair of fast shocks. All three phenomena are eventually manifested in interplanetary space as a magnetic cloud accompanied by a fast shock pair, with a forward shock preceding the cloud and a reverse shock appearing either within or behind the cloud.

Whang, Y. V.↗

Impact studies of a 1/3-scale model of an air cushion vehicle

An experimental investigation was conducted to determine the effects of various parameters of the impact performance of a 1/3-scale dynamic model of an air cushion vehicle. Impact response was determined by measuring the maximum values of variables, including sidelobe, front lobe, and cavity pressures, normal acceleration, pitch and roll angles, and vertical displacement during impact, for various combinations of drop height, initial pitch and roll angles, and forward speed. Increasing initial pitch angle increased the maximum values of the front lobe pressure, normal acceleration, nose down pitch angle, and to some extent, vertical displacement, but it inversely affected the maximum cavity pressure. Increasing the drop height of the model increased the potential energy of the system and generally produced larger responses over the entire range of variables measured, except for the roll angle after impact, which remained constant. Forward speed had no effect on the impact performance of the model, except for essentially doubling the maximum nose down pitch angle after impact at the maximum speed tested.

Daugherty, R. H.↗

Hadronic interaction model sibyll 2.3d and extensive air showers

We present a new version of the hadron interaction event generator SIBYLL. While the core ideas of the model have been preserved, the new version handles the production of baryon pairs and leading particles in a new way. In addition, production of charmed hadrons is included. Updates to the model are informed by high-precision measurements of the total and inelastic cross sections with the forward detectors at the LHC that constrain the extrapolation to ultrahigh energy. Minimum-bias measurements of particle spectra and multiplicities support the tuning of fragmentation parameters. This paper demonstrates the impact of these changes on air-shower observables such as X max and N μ , drawing comparisons with other contemporary cosmic-ray interaction models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Machine Learning Self-Diffusion Prediction for Lennard-Jones Fluids in Pores

Predicting the diffusion coefficient of fluids under nanoconfinement is important for many applications including the extraction of shale gas from kerogen and product turnover in porous catalysts. Due to the large number of important variables, including pore shape and size, fluid temperature and density, and the fluid–wall interaction strength, simulating diffusion coefficients using molecular dynamics (MD) in a systematic study could prove to be prohibitively expensive. Here, we use machine learning models trained on a subset of MD data to predict the self-diffusion coefficients of Lennard-Jones fluids in pores. Our MD data set contains 2280 simulations of ideal slit pore, cylindrical pore, and hexagonal pore geometries. We use the forward feature selection method to determine the most useful features (i.e., descriptors) for developing an artificial neutral network (ANN) model with an emphasis on easily acquired features. Our model shows good predictive ability with a coefficient of determination (i.e., R 2 ) of ~0.99 and a mean squared error of ~2.9 × 10 –5 . Finally, we propose an alteration to our feature set that will allow the ANN model to be applied to nonideal pore geometries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Parameter estimation with maximal updated densities

A recently developed measure-theoretic framework solves a stochastic inverse problem (SIP) for models where uncertainties in model output data are predominantly due to aleatoric (i.e., irreducible) uncertainties in model inputs (i.e., parameters). The subsequent inferential target is a distribution on parameters. Another type of inverse problem is to quantify uncertainties in estimates of “true” parameter values under the assumption that such uncertainties should be reduced as more data are incorporated into the problem, i.e., the uncertainty is considered epistemic. A major contribution of this work is the formulation and solution of such a parameter identification problem (PIP) within the measure-theoretic framework developed for the SIP. The approach is novel in that it utilizes a solution to a stochastic forward problem (SFP) to update an initial density only in the parameter directions informed by the model output data. In other words, this method performs “selective regularization” only in the parameter directions not informed by data. Further, the solution is defined by a maximal updated density (MUD) point where the updated density defines the measure-theoretic solution to the PIP. Another significant contribution of this work is the full theory of existence and uniqueness of MUD points for linear maps with Gaussian distributions. Data-constructed Quantity of Interest (QoI) maps are also presented and analyzed for solving the PIP within this measure-theoretic framework as a means of reducing uncertainties in the MUD estimate. We conclude with a demonstration of the general applicability of the method on two problems involving either spatial or temporal data for estimating uncertain model parameters. The first problem utilizes spatial data from a stationary partial differential equation to produce a MUD estimate of an uncertain boundary condition. The second problem utilizes temporal data obtained from the state-of-the-art ADvanced CIRCulation (ADCIRC) model to obtain a MUD estimate of uncertain wind drag coefficients for a simulated extreme weather event near the Shinnecock Inlet located in the Outer Barrier of Long Island, NY, USA.

42 ENGINEERING↗