Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mathematical modelling”

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,333 records · Page 74

Dynamic responses of railroad car models to vertical and lateral rail inputs

Simplified dynamic models were applied in a study of vibration in a high-speed railroad car. The mathematical models used were a four-degree-of-freedom model for vertical responses to vertical rail inputs and a ten-degree-of-freedom model for lateral response to lateral or rolling (cross-level) inputs from the rails. Elastic properties of the passenger car body were represented by bending and torsion of a uniform beam. Rail-to-car (truck) suspensions were modeled as spring-mass-dashpot oscillators. Lateral spring nonlinearities approximating certain complicated truck mechanisms were introduced. The models were excited by displacement and, in some cases, velocity inputs from the rails by both deterministic (including sinusoidal) and random input functions. Results were obtained both in the frequency and time domains. Solutions in the time domain for the lateral model were obtained for a wide variety of transient and random inputs generated on-line by an analog computer. Variations in one of the damping properties of the lateral car suspension gave large fluctuations in response over a range of car speeds for a given input. This damping coefficient was significant in reducing lateral car responses that were higher for nonlinear springs for three different inputs.

Sewall, J. L.↗

Prospects of a mathematical theory of human behavior in complex man-machine systems tasks

A hierarchy of human activities is derived by analyzing automobile driving in general terms. A structural description leads to a block diagram and a time-sharing computer analogy. The range of applicability of existing mathematical models is considered with respect to the hierarchy of human activities in actual complex tasks. Other mathematical tools so far not often applied to man machine systems are also discussed. The mathematical descriptions at least briefly considered here include utility, estimation, control, queueing, and fuzzy set theory as well as artificial intelligence techniques. Some thoughts are given as to how these methods might be integrated and how further work might be pursued.

Johannsen, G.↗

Preliminary results of experimental and analytical investigations of the tumbling phenomenon for an advanced configuration

A sustained autorotative pitching motion usually called 'tumbling' has been observed during dynamic model tests of the X-29A configuration. The X-29 is an advanced design incorporating forward-swept wings and canards in a highly relaxed static stability condition. Beginning with a historical review of the tumbling phenomenon, this paper discusses the current experimental results of dynamic model tumbling tests of the X-29 and the initial efforts to establish an aerodynamic and mathematical model for analysis.

Whipple, R. D.↗

Investigation of control law for active flutter suppression

Discussed is the technology of active flutter suppression research by the use of the wind-tunnel model of a delta wing with low aspect ratio. The research includes the establishment of a mathematical model for the servo-aeroelastic system, and the synthesis and analysis of sub-optimal output feedback control laws for active flutter suppression.

Wang, Weipeng↗

Model Of Orbital Density Of Air For Computing Drag

Simple, Orbital Density Model for Drag Equations program useful for computing effect of drag over one or more orbits. Mathematical model embodied in program incorporates major changes in density due to solar activity and magnetic activity of Earth. Diurnal (day/night) effects on orbit averaged out. Based on Jacchia daily-average density, evaluated at average time of year. Advantages, right ascension and declination of Sun not needed and computation time much reduced. Written in FORTRAN 77.

Lear, W. M.↗

Design-Tradeoff Model For Space Station

System Design Tradeoff Model (SDTM) computer program produces information which helps to enforce consistency of design objectives throughout system. Mathematical model of set of possible designs for Space Station Freedom. Program finds particular design enabling station to provide specified amounts of resources to users at lowest total (or life-cycle) cost. Compares alternative design concepts by changing set of possible designs, while holding specified services to users constant, and then comparing costs. Finally, both costs and services varied simultaneously when comparing different designs. Written in Turbo C 2.0.

Chamberlain, Robert G.↗

Finite element simulations of thermosolutal convection in vertical solidification of binary alloys

Dendritic vertical solidification of a binary alloy is modeled using the finite element method to assess the effect of thermosolutal convection in macrosegregation. The mathematical model assumes steady-state solidification with a planar, undeformable surface defined by the dendrite tips and the eutectic isotherm. The dendritic region is assumed to advance at a constant solidification velocity v. The stability of the modeled system has been investigated and nonlinear calculations performed that show finger-like convection when the system is unstable. Results for lead-tin alloys show that when the system is unstable, convection is only significant in the uppermost part of the mush and is entirely driven by convection in the bulk fluid.

Heinrich, J. C.↗

P-polarized reflectance spectroscopy: A high sensitive real-time monitoring technique to study surface kinetics under steady state epitaxial deposition conditions

This paper describes the results of real-time optical monitoring of epitaxial growth processes by p-polarized reflectance spectroscopy (PRS) using a single wavelength application under pulsed chemical beam epitaxy (PCBE) condition. The high surface sensitivity of PRS allows the monitoring of submonolayer precursors coverage on the surface as shown for GaP homoepitaxy and GaP on Si heteroepitaxy as examples. In the case of heteroepitaxy, the growth rate and optical properties are revealed by PRS using interference oscillations as they occur during growth. Super-imposed on these interference oscillations, the PRS signal exhibits a fine structure caused by the periodic alteration of the surface chemistry by the pulsed supply of chemical precursors. This fine structure is modeled under conditions where the surface chemistry cycles between phosphorus supersaturated and phosphorus depleted surfaces. The mathematical model describes the fine structure using a surface layer that increases during the tertiarybutyl phosphine (TBP) supply and decreases during and after the triethylgallium (TEG) pulse, which increases the growing GaP film thickness. The imaginary part of the dielectric function of the surface layer is revealed from the turning points in the fine structure, where the optical response to the first precursor pulse in the cycle sequence changes sign. The amplitude of the fine structure is determined by the surface layer thickness and the complex dielectric functions for the surface layer with the underlying bulk film. Surface kinetic data can be obtained by analyzing the rise and decay transients of the fine structure.

Dietz, Nikolaus↗

Design Considerations of Polishing Lap for Computer-Controlled Cylindrical Polishing Process

The future X-ray observatory missions, such as International X-ray Observatory, require grazing incidence replicated optics of extremely large collecting area (3 m2) in combination with angular resolution of less than 5 arcsec half-power diameter. The resolution of a mirror shell depends ultimately on the quality of the cylindrical mandrels from which they are being replicated. Mid-spatial-frequency axial figure error is a dominant contributor in the error budget of the mandrel. This paper presents our efforts to develop a deterministic cylindrical polishing process in order to keep the mid-spatial-frequency axial figure errors to a minimum. Simulation studies have been performed to optimize the operational parameters as well as the polishing lap configuration. Furthermore, depending upon the surface error profile, a model for localized polishing based on dwell time approach is developed. Using the inputs from the mathematical model, a mandrel, having conical approximated Wolter-1 geometry, has been polished on a newly developed computer-controlled cylindrical polishing machine. We report our first experimental results and discuss plans for further improvements in the polishing process.

Khan, Gufran S.↗

International Space Station Modal Correction Analysis

This paper summarizes the on-orbit modal test and the related modal analysis, model validation and correlation performed for the ISS Stage ULF4, DTF S4-1A, October 11,2010, GMT 284/06:13:00.00. The objective of this analysis is to validate and correlate analytical models with the intent to verify the ISS critical interface dynamic loads and improve fatigue life prediction. For the ISS configurations under consideration, on-orbit dynamic responses were collected with Russian vehicles attached and without the Orbiter attached to the ISS. ISS instrumentation systems that were used to collect the dynamic responses during the DTF S4-1A included the Internal Wireless Instrumentation System (IWIS), External Wireless Instrumentation System (EWIS), Structural Dynamic Measurement System (SDMS), Space Acceleration Measurement System (SAMS), Inertial Measurement Unit (IMU) and ISS External Cameras. Experimental modal analyses were performed on the measured data to extract modal parameters including frequency, damping and mode shape information. Correlation and comparisons between test and analytical modal parameters were performed to assess the accuracy of models for the ISS configuration under consideration. Based on the frequency comparisons, the accuracy of the mathematical models is assessed and model refinement recommendations are given. Section 2.0 of this report presents the math model used in the analysis. This section also describes the ISS configuration under consideration and summarizes the associated primary modes of interest along with the fundamental appendage modes. Section 3.0 discusses the details of the ISS Stage ULF4 DTF S4-1A test. Section 4.0 discusses the on-orbit instrumentation systems that were used in the collection of the data analyzed in this paper. The modal analysis approach and results used in the analysis of the collected data are summarized in Section 5.0. The model correlation and validation effort is reported in Section 6.0. Conclusions and recommendations drawn from this analysis are included in Section 7.0.

Fotz[atrocl. Lrostom↗

Dynamic modeling of structures from measured complex modes

A technique is presented to use a set of identified complex modes together with an analytical mathematical model of a structure under test to compute improved mass, stiffness and damping matrices. A set of identified normal modes, computed from the measured complex modes, is used in the mass orthogonality equation to compute an improved mass matrix. This eliminates possible errors that may result from using approximated complex modes as normal modes. The improved mass matrix, the measured complex modes and the higher analytical modes are then used to compute the improved stiffness and damping matrices. The number of degrees-of-freedom of the improved model is limited to equal the number of elements in the measured modal vectors. A simulated experiment shows considerable improvements, in the system's analytical dynamic model, over the frequency range of the given measured modal information.

Ibrahim, s. R.↗

Dynamic Modeling of Structures from Measured Complex Modes

A technique is presented to use a set of identified complex modes together with an analytical mathematical model of a structure under test to compute improved mass, stiffness and damping matrices. A set of identified normal modes, computed from the measured complex modes, is used in the mass orthogonality equation to compute an improved mass matrix. This eliminates possible errors that may result from using approximated complex modes as normal modes. The improved mass matrix, the measured complex modes and the higher analytical modes are then used to compute the improved stiffness and damping matrices. The number of degrees-of-freedom of the improved model is limited to equal the number of elements in the measured modal vectors. A simulated experiment shows considerable improvements, in the system's analytical dynamic model, over the frequency range of the given measured modal information.

Ibrahim, S. R.↗

A review of modelling techniques for the open and closed-loop dynamics of large space systems

This paper reviews the steps in the development of mathematical models that can be used to simulate the in-orbit dynamic behavior of large flexible systems. A general continuum formulation is compared with the hybrid coordinate formulation and also a finite element representation of the total system. A review of structural analysis routines emphasizes the use of computer generated graphics to help understand the different modal elastic shape functions of complex systems. Numerical techniques employed to synthesize shape and attitude control laws are summarized. Finally, the modeling of environmental disturbance torques due to the interaction of solar radiation pressure on vibrating and thermally deflected systems is discussed.

Bainum, Peter M.↗

Ocean Wave Studies with Applications to Ocean Modeling and Improvement of Satellite Altimeter Measurements

Combining analysis of satellite data (altimeter, scatterometer, high-resolution visible and infrared images, etc.) with mathematical modeling of non-linear wave processes, we investigate various ocean wave fields (on scales from capillary to planetary), their role in ocean dynamics and turbulent transport (of heat and biogeochemical quantities), and their effects on satellite altimeter measuring accuracy. In 1998 my attention was focused on long internal gravity waves (10 to 1000 km), known also as baroclinic inertia-gravity (BIG) waves. We found these waves to be a major factor of altimeter measurements "noise," resulting in a greater uncertainty [up to 10 cm in terms of sea surface height (SSH) amplitude] in the measured SSH signal than that caused by the sea state bias variations (up to 5 cm or so). This effect still remains largely overlooked by the satellite altimeter community. Our studies of BIG waves address not only their influence on altimeter measurements but also their role in global ocean dynamics and in transport and turbulent diffusion of biogeochemical quantities. In particular, in collaboration with Prof Peter Weichman, Caltech, we developed a theory of turbulent diffusion caused by wave motions of most general nature. Applied to the problem of horizontal turbulent diffusion in the ocean, the theory yielded the effective diffusion coefficient as a function of BIG wave parameters obtainable from satellite altimeter data. This effort, begun in 1997, has been successfully completed in 1998. We also developed a theory that relates spatial fluctuations of scalar fields (such as sea surface temperature, chlorophyll concentration, drifting ice concentration, etc.) to statistical characteristics of BIG waves obtainable from altimeter measurements. A manuscript is in the final stages of preparation. In order to verify the theoretical predictions and apply them to observations, we are now analyzing Sea-viewing Wide Field of view Sensor (SeaWiFS) and Field of view Sensor (SeaWiFS) and Advanced Very High-Resolution Radiometer (AVHRR) data on sea surface temperature (SST) and chlorophyll concentration jointly with TOPEX/POSEIDON data on SSH variations.

Glazman, Roman E.↗

Modeling and Analysis of the Reverse Water Gas Shift Process for In-Situ Propellant Production

This report focuses on the development of mathematical models and simulation tools developed for the Reverse Water Gas Shift (RWGS) process. This process is a candidate technology for oxygen production on Mars under the In-Situ Propellant Production (ISPP) project. An analysis of the RWGS process was performed using a material balance for the system. The material balance is very complex due to the downstream separations and subsequent recycle inherent with the process. A numerical simulation was developed for the RWGS process to provide a tool for analysis and optimization of experimental hardware, which will be constructed later this year at Kennedy Space Center (KSC). Attempts to solve the material balance for the system, which can be defined by 27 nonlinear equations, initially failed. A convergence scheme was developed which led to successful solution of the material balance, however the simplified equations used for the gas separation membrane were found insufficient. Additional more rigorous models were successfully developed and solved for the membrane separation. Sample results from these models are included in this report, with recommendations for experimental work needed for model validation.

Whitlow, Jonathan E.↗

Qualitative and quantitative enhancement of parameter estimation for model-based diagnostics using automatic differentiation with an application to inertial fusion

Parameter estimation using observables is a fundamental concept in the experimental sciences. Mathematical models that represent the physical processes can enable reconstructions of the experimental observables and greatly assist in parameter estimation by turning it into an optimization problem which can be solved by gradient-free or gradient-based methods. In this work, the recent rise in flexible frameworks for developing differentiable scientific computing programs is leveraged in order to dramatically accelerate data analysis of a common experimental diagnostic relevant to laser–plasma and inertial fusion experiments, Thomson scattering. A differentiable Thomson-scattering data analysis tool is developed that uses reverse-mode automatic differentiation (AD) to calculate gradients. By switching from finite differencing to reverse-mode AD, three distinct outcomes are achieved. First, gradient descent is accelerated dramatically to the extent that it enables near real-time usage in laser–plasma experiments. Second, qualitatively novel quantities which require $\mathcal{O}(10^3)$ parameters can now be included in the analysis of data which enables unprecedented measurements of small-scale laser–plasma phenomena. Third, uncertainty estimation approaches that leverage the value of the Hessian become accurate and efficient because reverse-mode AD can be used for calculating the Hessian.

97 MATHEMATICS AND COMPUTING↗

Neural message-passing for objective-based uncertainty quantification and optimal experimental design

Various real-world scientific applications involve the mathematical modeling of complex uncertain systems with numerous unknown parameters. Accurate parameter estimation is often practically infeasible in such systems, as the available training data may be insufficient and the cost of acquiring additional data may be high. In such cases, based on a Bayesian paradigm, we can design robust operators retaining the best overall performance across all possible models and design optimal experiments that can effectively reduce uncertainty to enhance the performance of such operators maximally. While objective-based uncertainty quantification (objective-UQ) based on MOCU (mean objective cost of uncertainty) provides an effective means for quantifying uncertainty in complex systems, the high computational cost of estimating MOCU has been a challenge in applying it to real-world scientific/engineering problems. In this work, we propose a novel scheme to reduce the computational cost for objective-UQ via MOCU based on a data-driven approach. We adopt a neural message-passing model for surrogate modeling, incorporating a novel axiomatic constraint loss that penalizes an increase in the estimated system uncertainty. As an illustrative example, we consider the optimal experimental design (OED) problem for uncertain Kuramoto models, where the goal is to predict the experiments that can most effectively enhance robust synchronization performance through uncertainty reduction. We show that our proposed approach can accelerate MOCU-based OED by four to five orders of magnitude, without any visible performance loss compared to the state-of-the-art. The proposed approach applies to general OED tasks, beyond the Kuramoto model.

97 MATHEMATICS AND COMPUTING↗

Measure this, not that: Optimizing the cost and model-based information content of measurements

Model-based design of experiments (MBDoE) is a powerful framework for selecting and calibrating science-based mathematical models from data. Here, this work extends popular MBDoE workflows by proposing a convex mixed integer (non)linear programming (MINLP) to optimize the selection of measurements. The solver MindtPy is modified to support calculating the D-optimality objective and its gradient via an external package, scipy, using the grey-box module in Pyomo. The new approach is demonstrated in two case studies: estimating highly correlated kinetics from a batch reactor and estimating transport parameters in a large-scale rotary packed bed for CO 2 capture. Both case studies show how examining the Pareto optimal trade-offs between information content measured by A- and D-optimality versus measurement budget offers practical guidance for selecting measurements for scientific experiments.

97 MATHEMATICS AND COMPUTING↗