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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

Hybrid NN/SVM Computational System for Optimizing Designs

A computational method and system based on a hybrid of an artificial neural network (NN) and a support vector machine (SVM) (see figure) has been conceived as a means of maximizing or minimizing an objective function, optionally subject to one or more constraints. Such maximization or minimization could be performed, for example, to optimize solve a data-regression or data-classification problem or to optimize a design associated with a response function. A response function can be considered as a subset of a response surface, which is a surface in a vector space of design and performance parameters. A typical example of a design problem that the method and system can be used to solve is that of an airfoil, for which a response function could be the spatial distribution of pressure over the airfoil. In this example, the response surface would describe the pressure distribution as a function of the operating conditions and the geometric parameters of the airfoil. The use of NNs to analyze physical objects in order to optimize their responses under specified physical conditions is well known. NN analysis is suitable for multidimensional interpolation of data that lack structure and enables the representation and optimization of a succession of numerical solutions of increasing complexity or increasing fidelity to the real world. NN analysis is especially useful in helping to satisfy multiple design objectives. Feedforward NNs can be used to make estimates based on nonlinear mathematical models. One difficulty associated with use of a feedforward NN arises from the need for nonlinear optimization to determine connection weights among input, intermediate, and output variables. It can be very expensive to train an NN in cases in which it is necessary to model large amounts of information. Less widely known (in comparison with NNs) are support vector machines (SVMs), which were originally applied in statistical learning theory. In terms that are necessarily oversimplified to fit the scope of this article, an SVM can be characterized as an algorithm that (1) effects a nonlinear mapping of input vectors into a higher-dimensional feature space and (2) involves a dual formulation of governing equations and constraints. One advantageous feature of the SVM approach is that an objective function (which one seeks to minimize to obtain coefficients that define an SVM mathematical model) is convex, so that unlike in the cases of many NN models, any local minimum of an SVM model is also a global minimum.

Rai, Man Mohan↗
NASA NTRS1991

Microwave Brightness Of Land Surfaces From Outer Space

Mathematical model approximates microwave radiation emitted by land surfaces traveling to microwave radiometer in outer space. Applied to measurements made by Scanning Multichannel Microwave Radiometer (SMMR). Developed for interpretation of microwave imagery of Earth to obtain distributions of various chemical, physical, and biological characteristics across its surface. Intended primarily for use in mapping moisture content of soil and fraction of Earth covered by vegetation. Advanced Very-High-Resolution Radiometer (AVHRR), provides additional information on vegetative cover, thereby making possible retrieval of soil-moisture values from SMMR measurements. Possible to monitor changes of land surface during intervals of 5 to 10 years, providing significant data for mathematical models of evolution of climate.

Kerr, Yann H.↗
NASA NTRS1976

Control of thermal balance by a liquid circulating garment based on a mathematical representation of the human thermoregulatory system

Test data and a mathematical model of the human thermoregulatory system were used to investigate control of thermal balance by means of a liquid circulating garment (LCG). The test data were derived from five series of experiments in which environmental and metabolic conditions were varied parametrically as a function of several independent variables, including LCG flowrate, LCG inlet temperature, net environmental heat exchange, surrounding gas ventilation rate, ambient pressure, metabolic rate, and subjective/obligatory cooling control. The resultant data were used to relate skin temperature to LCG water temperature and flowrate, to assess a thermal comfort band, to demonstrate the relationship between metabolic rate and LCG heat dissipation, and so forth. The usefulness of the mathematical model as a tool for data interpretation and for generation of trends and relationships among the various physiological parameters was also investigated and verified.

Kuznetz, L. H.↗
NASA NTRS1982

Fracture and crack growth in orthotropic laminates. Part 2: Experimental determination of internal damage growth in unidirectional boron/aluminum composite laminates

The fracture behavior of unidirectional boron/aluminum composite laminates is investigated in order to verify the results of mathematical models. These models predict the stresses and displacements of fibers and the amount of damage growth in a center-notched lamina as a function of the applied remote stress and the matrix and fiber material properties. The damage may take the form of longitudinal yielding and splitting in the matrix as well as stable transverse damage consisting of broken fibers and matrix yielding ahead of the notch. A brittle lacquer coating is used to detect the yielding in the matrix while X-ray techniques are used to detemine the number of broken fibers in the laminate. The notched strengths and the amounts of damage found in the experimental specimens agree well with those predicted by the mathematical model.

Goree, J. G.↗
NASA NTRS1993

Use of high performance networks and supercomputers for real-time flight simulation

In order to meet the stringent time-critical requirements for real-time man-in-the-loop flight simulation, computer processing operations must be consistent in processing time and be completed in as short a time as possible. These operations include simulation mathematical model computation and data input/output to the simulators. In 1986, in response to increased demands for flight simulation performance, NASA's Langley Research Center (LaRC), working with the contractor, developed extensions to the Computer Automated Measurement and Control (CAMAC) technology which resulted in a factor of ten increase in the effective bandwidth and reduced latency of modules necessary for simulator communication. This technology extension is being used by more than 80 leading technological developers in the United States, Canada, and Europe. Included among the commercial applications are nuclear process control, power grid analysis, process monitoring, real-time simulation, and radar data acquisition. Personnel at LaRC are completing the development of the use of supercomputers for mathematical model computation to support real-time flight simulation. This includes the development of a real-time operating system and development of specialized software and hardware for the simulator network. This paper describes the data acquisition technology and the development of supercomputing for flight simulation.

Cleveland, Jeff I., II↗
NASA NTRS2000

Numerical Calculation of the Drag Force Acting on a Solid Particle Pushed by a Solid/Liquid Interface

The distribution of insoluble particles in a metal casting depends primarily on the interaction of the particles with the solid/liquid interface (SLI) during the solidification process. The balance of the forces acting on the particle essentially determines whether a particle will be engulfed or pushed by the SLI. An important component of this force balance is the drag force generated by the particle motion in front of the SLI. Previously developed mathematical models for particle/SLI interaction made use of steady-state solutions of this force provided by the lubrication theory. However, our numerical model based on the SLI tracking approach shows that the steady-state approach is inappropriate to model the interaction process and that at steady-state the theoretical solution underestimates the drag force. It was found that regression analysis of steady-state numerical solutions for cylindrical particles moving normal to a flat SLI gives a relationship of the form: Abstract The distribution of insoluble particles in a metal casting depends primarily on the interaction of the particles with the solid/liquid interface (SLI) during the solidification process. The balance of the forces acting on the particle essentially determines whether a particle will be engulfed or pushed by the SLI. An important component of this force balance is the drag force generated by the particle motion in front of the SLI. Previously developed mathematical models for particle/SLI interaction made use of steady-state solutions of this force provided by the lubrication theory. However, our numerical model based on the SLI tracking approach shows that the steady-state approach is inappropriate to model the interaction process and that at steady-state the theoretical solution underestimates the drag force. It was found that regression analysis of steady-state numerical solutions for cylindrical particles moving normal to a flat SLI gives a relationship of the form: F(sub D, sup num) =sqoare root of 3(pi)(eta)V(sub p)(R(sub p)/d)(sup 10(gamma)/3). This is to compared to the solution provided by the lubrication theory: F(sub D), sup theor) = 3 square root of 2(pi)(eta)V(sub p)(R(sub p)/d)(sub 3/2), where F(sub D) is the drag force, eta is the dynamic viscosity of the fluid, V(sub p) is the particle velocity, R(sub p) is the particle radius, gamma is Euler's constant (gamma approx. = 0.577), and d is the width of the gap between the particle and the SLI. The domain on which the equations proposed by the lubrication theory are relevant is clearly identified in the paper. The numerical model was then validated against the classical lubrication theory within the domain of its validity.

Catalina, Adrian V.↗
NASA NTRS1988

Vertically reciprocating auger

The mathematical model and test results developed for the Vertically Reciprocating Auger (VRA) are summarized. The VRA is a device capable of transporting cuttings that result from below surface drilling. It was developed chiefly for the lunar surface, where conventional fluid flushing while drilling would not be practical. The VRA uses only reciprocating motion and transports material through reflections with the surface above. Particles are reflected forward and land ahead of radially placed fences, which prevent the particles from rolling back down the auger. Three input wave forms are considered to drive the auger. A modified sawtooth wave form was chosen for testing, over a modified square wave or sine wave, due to its simplicity and effectiveness. The three-dimensional mathematical model predicted a sand throughput rate of 0.2667 pounds/stroke, while the actual test setup transported 0.075 pounds/stroke. Based on this result, a correction factor of 0.281 is suggested for a modified sawtooth input.

Etheridge, Mark↗
DOE OSTI2023

Optimal Strategies for Hybrid Battery‐Storage Systems Design

As stationary hybrid energy‐storage systems (HESS) for power systems applications have recently drawn interest due to their enhanced performance and decreasing cost, developing systematic approaches for HESS design while considering controls is gaining traction. Herein, a method is presented to optimally design hybrid battery storage by proposing a mathematical modeling framework, formulated as a mixed integer linear programming model. The optimization is capable of handling multiple subsystems of batteries, considering their economic and technological performance. Decisions involve sizing of the batteries, optimal temporal and strategic dispatch to end uses, and energy sources for charging each battery. The applicability of the model is tested on four case studies for three battery chemistries representing distinct objectives: high‐power, high‐energy, and second life. Compared to traditionally designed battery storage with a homogeneous battery, optimally designed hybrid systems can save 12%–26% of system costs, depending on the nature of the dispatch profile. Findings point to design preference toward the second life battery supplemented with some high‐power or high‐energy battery capacity, or both. With the utilized electricity price structure, customers can experience approximately 10%–35% reduction in their bills.

Koleva, Mariya↗
NASA NTRS2013

Prediction of Ablation Rates from Solid Surfaces Exposed to High Temperature Gas Flow

A mathematical model and a solution algorithm is developed to study the physics of high temperature heat transfer and material ablation and identify the problems associated with the flow of hydrogen gas at very high temperatures and velocities through pipes and various components of Nuclear Thermal Rocket (NTR) motors. Ablation and melting can be experienced when the inner solid surface of the cooling channels and the diverging-converging nozzle of a Nuclear Thermal Rocket (NTR) motor is exposed to hydrogen gas flow at temperatures around 2500 degrees Kelvin and pressures around 3.4 MPa. In the experiments conducted on typical NTR motors developed in 1960s, degradation of the cooling channel material (cracking in the nuclear fuel element cladding) and in some instances melting of the core was observed. This paper presents the results of a preliminary study based on two types of physics based mathematical models that were developed to simulate the thermal-hydrodynamic conditions that lead to ablation of the solid surface of a stainless steel pipe exposed to high temperature hydrogen gas near sonic velocities. One of the proposed models is one-dimensional and assumes the gas flow to be unsteady, compressible and viscous. An in-house computer code was developed to solve the conservations equations of this model using a second-order accurate finite-difference technique. The second model assumes the flow to be three-dimensional, unsteady, compressible and viscous. A commercial CFD code (Fluent) was used to solve the later model equations. Both models assume the thermodynamic and transport properties of the hydrogen gas to be temperature dependent. In the solution algorithm developed for this study, the unsteady temperature of the pipe is determined from the heat equation for the solid. The solid-gas interface temperature is determined from an energy balance at the interface which includes heat transfer from or to the interface by conduction, convection, radiation, and ablation. Two different ablation models are proposed to determine the heat loss from the solid surface due to the ablation of the solid material. Both of them are physics based. Various numerical simulations were carried out using both models to predict the temperature distribution in the solid and in the gas flow, and then predict the ablation rates at a typical NTR motor hydrogen gas temperature and pressure. Solid mass loss rate per foot of a pipe was also calculated from these predictions. The results are presented for fully developed turbulent flow conditions in a sample SS pipe with a 6 inch diameter.

Akyuzlu, Kazim M.↗
NASA NTRS1994

High performance real-time flight simulation at NASA Langley

In order to meet the stringent time-critical requirements for real-time man-in-the-loop flight simulation, computer processing operations must be deterministic and be completed in as short a time as possible. This includes simulation mathematical model computational and data input/output to the simulators. In 1986, in response to increased demands for flight simulation performance, personnel at NASA's Langley Research Center (LaRC), working with the contractor, developed extensions to a standard input/output system to provide for high bandwidth, low latency data acquisition and distribution. The Computer Automated Measurement and Control technology (IEEE standard 595) was extended to meet the performance requirements for real-time simulation. This technology extension increased the effective bandwidth by a factor of ten and increased the performance of modules necessary for simulator communications. This technology is being used by more than 80 leading technological developers in the United States, Canada, and Europe. Included among the commercial applications of this technology are nuclear process control, power grid analysis, process monitoring, real-time simulation, and radar data acquisition. Personnel at LaRC have completed the development of the use of supercomputers for simulation mathematical model computational to support real-time flight simulation. This includes the development of a real-time operating system and the development of specialized software and hardware for the CAMAC simulator network. This work, coupled with the use of an open systems software architecture, has advanced the state of the art in real time flight simulation. The data acquisition technology innovation and experience with recent developments in this technology are described.

Cleveland, Jeff I., II↗
NASA NTRS1999

The Effects of Partial Mechanical Loading and Ibandronate on Skeletal Tissues in the Adult Rat Hindquarter Suspension Model for Microgravity

We report initial data from a suspended rat model that quantitatively relates chronic partial weightbearing to bone loss. Chronic partial weightbearing is our simulation of the effect of limited artificial gravity aboard spacecraft or reduced planetary gravity. Preliminary analysis of bone by PQCT, histomorphometry, mechanical testing and biochemistry suggest that chronic exposure to half of Earth gravity is insufficient to prevent severe bone loss. The effect of episodic full weightbearing activity (Earth Gravity) on rats otherwise at 50% weightbearing was also explored. This has similarity to treatment by an Earth G-rated centrifuge on a spacecraft that normally maintained artificial gravity at half of Earth G. Our preliminary evidence, using the above techniques to analyze bone, indicate that 2 hours daily of full weightbearing was insufficient to prevent the bone loss observed in 50% weightbearing animals. The effectiveness of partial weightbearing and episodic full weightbearing as potential countermeasures to bone loss in spaceflight was compared with treatment by ibandronate. Ibandronate, a long-acting potent bisphosphonate proved more effective in preventing bone loss and associated functionality based upon structure than our first efforts at mechanical countermeasures. The effectiveness of ibandronate was notable by each of the testing methods we used to study bone from gross structure and strength to tissue and biochemistry. These results appear to be independent of generalized systemic stress imposed by the suspension paradigm. Preliminary evidence does not suggest that blood levels of vitamin D were affected by our countermeasures. Despite the modest theraputic benefit of mechanical countermeasures of partial weightbearing and episodic full weightbearing, we know that some appropriate mechanical signal maintains bone mass in Earth gravity. Moreover, the only mechanism that correctly assigns bone mass and strength to oppose regionally specific force applied to bone is mechanical, a process based upon bone strain. Substantial evidence indicates that the specifics of dynamic loading i.e. time-varying forces are critical. Bone strain history is a predictor of the effect that mechanical conditions have on bone structure mass and strength. Using servo-controlled force plates on suspended rats with implanted strain gauges we manipulated impact forces of ambulation in the frequency (Fourier) domain. Our results indicate that high frequency components of impact forces are particularly potent in producing bone strain independent of the magnitude of the peak force or peak energy applied to the leg. Because a servo-system responds to forces produced by the rat's own muscle activity during ambulation, the direction of ground-reaction loads act on bone through the rat's own musculature. This is in distinction to passive vibration of the floor where forces reach bone through the natural filters of soft tissue and joints. Passive vibration may also be effective, but it may or may not increase bone in the appropriate architectural pattern to oppose the forces of normal ambulatory activity. Effectiveness of high frequency mechanical stimulation in producing regional (muscle directed) bone response will be limited by 1. the sensitivity of bone to a particular range of frequencies and 2. the inertia of the muscles, limiting their response to external forces by increasing tension along insertions. We have begun mathematical modeling of normal ambulatory activity. Effectiveness of high frequency mechanical stimulation in producing regional (muscle directed) bone response will be limited by 1. the sensitivity of bone to a particular range of frequencies and 2. the inertia of the muscles, limiting their response to external forces by increasing tension along insertions. We have begun mathematical modeling of the rat forelimb as a transfer function between impact force and bone strain to predict optimal dynamic loading conditions for this system. We plan additional studies of mechanical counter-measures that incorporate improved dynamic loading, features relevant to anticipated evaluation of artificial gravity, exercise regimens and exposure to Martian gravity, The combination of mechanical countermeasures with ibandronate will also be investigated for signs of synergy.

Schultheis, Lester W.↗
DOE OSTI2024

Diffusioosmotic flow reversals due to ion–ion electrostatic correlations

Existing theories of diffusioosmosis have neglected ion–ion electrostatic correlations, which are important in concentrated electrolytes. Here, we develop a mathematical model to numerically compute the diffusioosmotic mobilities of binary symmetric electrolytes across low to high concentrations in a charged parallel-plate channel. We use the modified Poisson equation to model the ion–ion electrostatic correlations and the Bikerman model to account for the finite size of ions. We report two key findings. First, ion–ion electrostatic correlations can cause a unique reversal in the direction of diffusioosmosis. Such a reversal is not captured by existing theories, occurs at ≈ 0.4 M for a monovalent electrolyte, and at a much lower concentration of ≈ 0.003 M for a divalent electrolyte in a channel with the same surface charge. This highlights that diffusioosmosis of a concentrated electrolyte can be qualitatively different from that of a dilute electrolyte, not just in its magnitude but also its direction. Second, we predict a separate diffusioosmotic flow reversal, which is not due to electrostatic correlations but the competition between the underlying chemiosmosis and electroosmosis. This reversal can be achieved by varying the magnitude of the channel surface charge without changing its sign. However, electrostatic correlations can radically change how this flow reversal depends on the channel surface charge and ion diffusivity between a concentrated and a dilute electrolyte. Furthermore, the mathematical model developed here can be used to design diffusioosmosis of dilute and concentrated electrolytes, which is central to applications such as species mixing and separation, enhanced oil recovery, and reverse electrodialysis.

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

An Analysis of the Impact of Particle Growth on Transport Losses in Polymer-Electrolyte Fuel Cells

Voltage cycling causes catalyst nanoparticles in polymer-electrolyte fuel cells to grow. The concomitant loss of interfacial area results in larger kinetic, transport, and possibly ohmic overpotentials. This paper uses recently published experimental data and mathematical models to investigate the evolutions of transport and ohmic resistances to platinum nanoparticles located on the surface and inside the micropores of carbon black supports. Resistance to oxygen transport rises as surface area declines primarily because the flux to each remaining larger platinum particle increases. The path lengths governing oxygen diffusion to surface and buried platinum sites also increase as nanoparticles grow. Platinum nanoparticles on the surface become relatively less favorable as voltage cycling proceeds because they grow faster than platinum in micropores. Because voltage cycling causes total interfacial area to decline, and the fraction located inside micropores to increase, nanoscale ohmic losses increase as a catalyst layer decays. The practical importance of this effect is difficult to discern because proton conductivity in carbon micropores is not well characterized.

07 ISOTOPE AND RADIATION SOURCES↗
NASA NTRS1995

Actively Controlled Shaft Seals for Aerospace Applications

This study experimentally investigates an actively controlled mechanical seal for aerospace applications. The seal of interest is a gas seal, which is considerably more compact than previous actively controlled mechanical seals that were developed for industrial use. In a mechanical seal, the radial convergence of the seal interface has a primary effect on the film thickness. Active control of the film thickness is established by controlling the radial convergence of the seal interface with a piezoelectric actuator. An actively controlled mechanical seal was initially designed and evaluated using a mathematical model. Based on these results, a seal was fabricated and tested under laboratory conditions. The seal was tested with both helium and air, at rotational speeds up to 3770 rad/sec, and at sealed pressures as high as 1.48 x 10(exp 6) Pa. The seal was operated with both manual control and with a closed-loop control system that used either the leakage rate or face temperature as the feedback. The output of the controller was the voltage applied to the piezoelectric actuator. The seal operated successfully for both short term tests (less than one hour) and for longer term tests (four hours) with a closed-loop control system. The leakage rates were typically 5-15 slm (standard liters per minute), and the face temperatures were generally maintained below 100C. When leakage rate was used as the feedback signal, the setpoint leakage rate was typically maintained within 1 slm. However, larger deviations occurred during sudden changes in sealed pressure. When face temperature was used as the feedback signal, the setpoint face temperature was generally maintained within 3 C, with larger deviations occurring when the sealed pressure changes suddenly. the experimental results were compared to the predictions from the mathematical model. The model was successful in predicting the trends in leakage rate that occurred as the balance ratio and sealed pressure changed, although the leakage rates were not quantitatively predicted with a high degree of accuracy. This model could be useful in providing valuable design information for future actively controlled mechanical seals.

Salant, Richard F.↗
DOE OSTI2021

Melting rate correlation with batch properties and melter operating conditions during conversion of nuclear waste melter feeds to glasses

The rate of conversion of nuclear waste melter feed to glass is affected by the selection of melter feed materials and by melter design and operation. The melting rate correlation (MRC) is an equation that relates the glass production rate with two types of variables: (1) feed and melt properties: conversion heat, cold-cap bottom temperature, and glass melt viscosity; and (2) melter design and operation parameters: melter geometry, melter operating temperature, and gas bubbling rate. The MRC shows good agreement for an extended melting-rate data set of high-level waste (HLW) melter feeds and a data set generated for low-activity waste (LAW) melter feeds. Laboratory observation of heated melter feed samples is often used to assess the cold-cap bottom temperature of HLW melter feeds (moderately foaming feeds), but this technique appears inadequate for LAW melter feeds (vigorously foaming feeds). For LAW feeds, an adequate assessment of the cold-cap bottom temperature was achieved using evolved gas analysis, which allows identification of the collapse of primary foam for oxidized feeds. This assessment shows that the cold-cap bottom temperature for vigorously foaming LAW feeds is higher than that for moderately foaming HLW feeds. When the results of MRC are compared, LAW feeds are generally less sensitive to the bubbling rate and melt viscosity, and more sensitive to the cold-cap bottom temperature than HLW feeds. The MRC qualifies as a promising tool to support the selection of melter feed materials and melter operating conditions, which is determined from expensive independent scaled melter experiments, and sophisticated mathematical models.

Lee, Seung Min↗
DOE OSTI2023

Modification of the MTS model for high strain-rate behavior of TI-6AL-4V

The Mechanical Threshold Stress (MTS) model provides excellent predictive capabilities for the material constitutive response for a wide range of temperatures and strain rates. However, the MTS model fails to capture the rapidly increasing yield stress at high strain rate behavior as the deformation controlling mechanism transitions from thermal activation to drag mechanisms, only capturing the linear behavior. Further, the model typically over predicts the flow stress behavior at yield and post yield due to its use of a constant work hardening rate parameter derived from the stress–strain response at constant saturation stress. An alternative approach to fitting portions of the MTS model is investigated and mathematical models are developed to address these issues. The results show that with appropriate experimental data, the mechanical threshold stress and work hardening rate parameters within the MTS model can quite easily and accurately be modified to extend applicability to high strain rate behavior and more accurately model the initial flow stress behavior at early work hardening rates without modification of the functions core to the MTS model itself.

36 MATERIALS SCIENCE↗