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,405 records · Page 78

Predicting Fractures In Fiber/Matrix Composites

NASA technical memorandum discusses status of emerging methodology of computational simulation of fractures in fiber/matrix composite structural components. Part of methodology consists of step-by-step procedures to simulate fractures in individual and mixed models in variety of generic composite-material components. Another part of methodology embodied in computer program CODSTRAN, which incorporates mathematical models of mechanical properties of composites and of fracture mechanics. Authors assert methodology ready for widespread application.

Chamis, C. C.↗

Inner loop flight control for the High-Speed Civil Transport

High-speed aerospace vehicles which employ high strength, light weight, yet deformable materials may exhibit significant interaction between the rigid-body and vibrational dynamics. Preliminary High-Speed Civil Transport (HSCT) configurations are a prime example. Traditionally, separate control systems have been used to augment the rigid-body and vibrational dynamics. In the HSCT arena, the highly coupled motions may not allow this design freedom. The research activity addresses two specific issues associated with the design and development of an integrated flight control system (FCS) for HSCT configurations, which are discussed next. The HSCT is expected to have a short period instability at subsonic speeds. Flight vehicles with this characteristic (i.e., F-16, F-22, X-29, Space Shuttle) are stabilized with what is called a superaugmented pitch rate loop. One concern is 'Will this stability augmentation logic work for a HSCT?' Studies show that an idealized pitch rate design would be acceptable, but is not realistic. Investigations using a contaminated pitch rate design reveal serious hurdles to overcome in the FCS design. Mounting location for the pitch rate sensor is critical. Results indicate a forward location leads to destabilizing pick-up of aeroelastic modes, while aft locations lead to undesirable coupling of the dominate pitch mode with the first aeroelastic mode. Intermediate locations for the sensor may not be acceptable. The source of the problem is the presence of low frequency aeroelastic modes in HSCT configurations, which are not present in vehicles currently using the superaugmented logic. To say the least, a conventional superaugmented pitch rate loop strategy may have undesirable characteristics. An unconventional strategy, which attempts to eliminate the above deficiencies by blending several pitch rate signals, indicates an improvement in the FCS architecture feasibility, but is still lacking in some respects. The HSCT configuration does not have aerodynamic surfaces in the vicinity of the nose (i.e., no canard or vane). A second concern is 'Can the fuselage bending/torsion aeroelastic modes be effectively augmented without sufficient control input near the vehicle nose?' The superaugmented FCS results above may be suggesting the necessity of a secondary feedback loop to achieve an acceptable integrated FCS. Preliminary analysis of HSCT aeroelastic mode shapes indicate the use of existing wing leading edge devices as a second control input may be lacking in control authority for the rigid-body attitude and aeroelastic modes. An effort is underway to incorporate generic wing leading edge devices and canards into a generic HSCT model for the purpose of assessing additional control authority and it's use in candidate FCS designs. A generic HSCT mathematical model was necessary for the studies above. A HSCT category model is available in NASA-CR-172201. This model describes the linear, longitudinal dynamics about the following flight condition: ascent, W = 730,000 lbs, h = 6,500 ft, M = 0.6. The model incorporates the full rigid-body variable set, as well as eighteen aeroelastic modes. Elevator deflection serves as the control input. Modifications to the model include the incorporation of relaxed static stability (i.e., static margin from -7.3% to +10%) and additional control inputs.

Newman, Brett A.↗

Robust Online Sequential RVFLNs for Data Modeling of Dynamic Time-Varying Systems with Application of an Ironmaking Blast Furnace

In a world where the increasing complexity of modern industrial processes brings difficulties for accurate mathematical modeling, taking advantage of data has become an efficient solution to complex dynamic process modeling issue. In this paper, we develop a novel robust online sequential version of random vector functional-link networks (RVFLNs) for data-driven modeling of dynamic time-varying system and applied it in a blast furnace (BF) ironmaking process. First, to overcome the time-varying dynamics of process and to enable the RVFLNs to learn online with avoiding data saturation, an improved online sequential version of RVFLNs (OS-RFVLNs) is first presented by online sequential learning with forgetting factor. This improved OS-RVFLNs algorithm is not only suitable for the real-time and large data transfer situation, but also can adjust the sensitivity of the algorithm to different samples with the help of the introduced forgetting factor. Second, since the output weights of the improved OS-RVFLNs as well as other RVFLNs algorithms are obtained by the least squares approach, a robustness problem may occur when the training dataset is contaminated with various outliers. To solve this problem, a Cauchy distribution weighted M-estimator is introduced to improve the robustness of the improved OS- RVFLNs. For this proposed robust OS-RVFLNs (R-OS- RVFLNs), since the weights of different outlier data are properly determined by the Cauchy distribution function, their corresponding contribution on modeling can be properly distinguished. Thus robust and better modeling results can be achieved. Experiments using actual industrial data of BF ironmaking process and comparative studies have demonstrated that the proposed method produces a better estimation accuracy and stronger robustness than other methods.

Blast furnace (BF), Modelling, Dynamic systems↗

(Invited) Modeling Electrolyzers: Exploring the Applied Voltage Breakdown

To commercialize successfully polymer-electrolyte electrolyzers and optimize their performance, one requires a detailed understanding of the underlying physics and phenomena. Mathematical modeling is ideally suited to explore such intricacies. In this Tutorial, the modeling equations and approaches towards both proton- and hydroxide-conducting polymer electrolyzers will be detailed. This includes the introduction of the applied voltage breakdown that separates the overall polarization curve into its constitutive parts such that the limiting mechanisms can be ascertained. Throughout, different case studies will be explored including impacts of bubble generation, liquid versus vapor feed, and pH changes due to carbonate and hydroxide solution feeds. A key focus of the Tutorial will be in exploring and describing the various transport and competing phenomena within the electrolyzer cell.

42 ENGINEERING↗

Model verification of large structural systems

A methodology was formulated, and a general computer code implemented for processing sinusoidal vibration test data to simultaneously make adjustments to a prior mathematical model of a large structural system, and resolve measured response data to obtain a set of orthogonal modes representative of the test model. The derivation of estimator equations is shown along with example problems. A method for improving the prior analytic model is included.

Lee, L. T.↗

Low order equivalent models of highly augmented aircraft determined from flight data using maximum likelihood estimation

This paper presents the results of a study of the feasibility of using low order equivalent mathematical models of a highly augmented aircraft, the F-8 digital fly-by-wire (DFBW), for flying qualities research. Increasingly complex models were formulated and evaluated using flight data and maximum likelihood estimation techniques. The aircraft actuator was modeled alone first. Next the equivalent derivatives were used to model the longitudinal unaugmented F-8 DFBW aircraft dynamics. The most complex model incorporated a pure time shift of the pilot input, a first order lag, and the basic longitudinal airframe model. This same model was then implemented for the F-8 DFBW aircraft in a highly augmented mode. Excellent matching of the dynamics resulted for this model, indicating that low order equivalent models which are good representations of the highly augmented F-8 DFBW aircraft can be formulated with these methods.

Shafer, M. F.↗

Modeling Non-homologous End Joining

Non-homologous end joining (NHEJ) is the dominant DNA double strand break (DSB) repair pathway and involves several NHEJ proteins such as Ku, DNA-PKcs, XRCC4, Ligase IV and so on. Once DSBs are generated, Ku is first recruited to the DNA end, followed by other NHEJ proteins for DNA end processing and ligation. Because of the direct ligation of break ends without the need for a homologous template, NHEJ turns out to be an error-prone but efficient repair pathway. Some mechanisms have been proposed of how the efficiency of NHEJ repair is affected. The type of DNA damage is an important factor of NHEJ repair. For instance, the length of DNA fragment may determine the recruitment efficiency of NHEJ protein such as Ku [1], or the complexity of the DNA breaks [2] is accounted for the choice of NHEJ proteins and subpathway of NHEJ repair. On the other hand, the chromatin structure also plays a role of the accessibility of NHEJ protein to the DNA damage site. In this talk, some mathematical models of NHEJ, that consist of series of biochemical reactions complying with the laws of chemical reaction (e.g. mass action, etc.), will be introduced. By mathematical and numerical analysis and parameter estimation, the models are able to capture the qualitative biological features and show good agreement with experimental data. As conclusions, from the viewpoint of modeling, how the NHEJ proteins are recruited will be first discussed for connection between the classical sequential model [4] and recently proposed two-phase model [5]. Then how the NHEJ repair pathway is affected, by the length of DNA fragment [6], the complexity of DNA damage [7] and the chromatin structure [8], will be addressed

Li, Yongfeng↗

Component-wise reduced order model lattice-type structure design

Lattice-type structures can provide a combination of stiffness with light weight that is desirable in a variety of applications. Design optimization of these structures must rely on approximations of the governing physics to render solution of a mathematical model feasible. In this paper, we propose a topology optimization (TO) formulation that approximates the governing physics using component-wise reduced order modeling as introduced in Huynh et al. (2013); Eftang and Patera (2013), which can reduce solution time by multiple orders of magnitude over a full-order finite element model while providing a relative error in the solution of 1%. In addition, the offline training data set from such component-wise models is reusable, allowing its application to many design problems for only the cost of a single offline training phase, and the component-wise method is nearly embarrassingly parallel. We also show how the parameterization chosen in our optimization allows a simplification of the component-wise reduced order model (CWROM) not noted in previous literature, for further speedup of the optimization process. Furthermore, the sensitivity of the compliance with respect to the particular parameterization is derived solely at the component level. In numerical examples, we demonstrate a 1000x speedup over a full-order FEM model with relative error of 1% and show minimum compliance designs for two different cantilever beam examples, one smaller and one larger. Finally, error bounds for the displacement field, compliance, and compliance sensitivity of the CWROM are derived.

97 MATHEMATICS AND COMPUTING↗

Estimating reaction parameters in mechanism-enabled population balance models of nanoparticle size distributions: A Bayesian inverse problem approach

In order to quantitatively predict nano- as well as other particle-size distributions, one needs to have both a mathematical model and estimates of the parameters that appear in these models. Here, we show how one can use Bayesian inversion to obtain statistical estimates for the parameters that appear in recently derived mechanism-enabled population balance models (ME-PBM) of nanoparticle growth. The Bayesian approach addresses the question of “how well do we know our parameters, along with their uncertainties?.” The results reveal that Bayesian inversion statistical analysis on an example, prototype $\mathrm{lr(0)_n}$ nanoparticle formation system allows one to estimate not just the most likely rate constants and other parameter values, but also their SDs, confidence intervals, and other statistical information. Moreover, knowing the reliability of the mechanistic model's parameters in turn helps inform one about the reliability of the proposed mechanism, as well as the reliability of its predictions. Importantly, the paper can also be seen as a tutorial with the additional goal of achieving a “Gold Standard” Bayesian inversion ME-PBM benchmark that others can use as a control to check their own use of this methodology for other systems of interest throughout nature. Overall, the results provide strong support for the hypothesis that there is substantial value in using a Bayesian inversion methodology for parameter estimation in particle formation systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Integrative Circuit-Host Modeling of a Genetic Switch in Varying Environments

Synthetic biology is advancing into a new phase where real-world applications are emphasized. There is hence an urgent need for mathematical modeling that can quantitatively describe the behaviors of genetic devices in natural, fluctuating environments. We utilize an integrative circuit-host modeling framework to examine the dynamics of a genetic switch and its host cell in varying environments. For both steady-state and transient cases, we find increasing nutrient reduces the bistability region of the phase space and eventually drives the switch from bistability to monostability. In response, cellular growth and proteome partitioning experience the same transition. Antibiotic perturbations cause the similar circuit and host responses as nutrient variations. However, one difference is the trend of growth rate, which augments with nutrient but declines with antibiotic levels. The framework provides a mechanistic scheme to account for both the dynamic and static characteristics of the circuit-host system upon environmental perturbations, underscoring the intimacy of gene circuits and their hosts and elucidating the complexity of circuit behaviors arising from environmental variations.

59 BASIC BIOLOGICAL SCIENCES↗

Symmetric linear systems - An application of algebraic systems theory

Dynamical systems which contain several identical subsystems occur in a variety of applications ranging from command and control systems and discretization of partial differential equations, to the stability augmentation of pairs of helicopters lifting a large mass. Linear models for such systems display certain obvious symmetries. In this paper, we discuss how these symmetries can be incorporated into a mathematical model that utilizes the modern theory of algebraic systems. Such systems are inherently related to the representation theory of algebras over fields. We will show that any control scheme which respects the dynamical structure either implicitly or explicitly uses the underlying algebra.

Hazewinkel, M.↗

Study Of Compressibility Corrections To Turbulence Models

Effects on shear layers in simulated confined and unconfined flows studied. Report presents comparative study of some terms that correct for effects of compressibility in standard k-epsilon mathematical model of turbulence where k denotes turbulence kinetic energy and epsilon denotes rate of dissipation of turbulence kenetic energy. Involved simulation of flows by numerical solution of Reynolds-averaged Navier-Stokes equations.

Viegas, J. R.↗

Generating A Strobed Laser Light Sheet

An optoelectronic system generating synchronous, strobed sheet of laser light developed for use in making visible flow of air about model helicopter rotor. Used in wind-tunnel tests to determine actual locations of vortices for comparison with locations predicted by mathematical models to validate models. Each blade tip produces vortex. By establishing successive vortex locations, researcher determines trajectory of vortex pattern. Light-sheet strobe circuits provide selection of blade positions, strobe-pulse durations, and multiple pulses per revolution for rotors having two to nine blades. To make flow visible, vaporizing propylene glycol injected upstream of model. System also provides calibrated trigger delay of strobe pulses, adjustable strobe-pulse durations, selectable number of blades, and slip-sync mode to make flow visible as though in slow motion.

Leighty, Bradley D.↗

Marangoni instability induced convection in an evaporating liquid droplet

The processes occurring when a liquid drop undergoes evaporation are described. When a liquid drop undergoes evaporation, its surface temperature decreases. If the droplet is free floating in a microgravity environment, the heat transfer process inside the droplet is initially condition controlled. As the process continues, a radial temperature gradient builds up at the free surface until the critical Marangoni number is exceeded. Then the onset of instability induces thermocapillary convective flows, which in turn speed up the evaporation. The convective flows will subside when the interior of the droplet reaches a certain equilibrium temperature, and the process will return to the diffusion controlled mode. Both preliminary modeling and recent laboratory data have confirmed that Marangoni instability induced convection can and does occur in the droplet evaporation process. Mathematical models representing Marangoni instability in an evaporating liquid drop are presented. An ideal space experiment to study and characterize the onset of Marangoni instability in an evaporating liquid droplet and to establish the effect of Marangoni instability induced convection on the droplet evaporation rate is outlined and the need for conducting such experiments in space is highlighted. However, before an opportunity to conduct experiments in space arises, ground based experiments have to be conducted to study feasibility issues and proof of concept. A ground based experiment of this type is outlined.

Chai, An-Ti↗

Constraints on galactic distributions of gamma-ray burst sources from BATSE observations

The paradigm that gamma-ray bursts originate from Galactic sources is studied in detail using the angular and intensity distributions observed by the Burst and Transient Source Experiment (BATSE) on NASA's Compton Gamma Ray Observatory (CGRO). Monte Carlo models of gamma-ray burst spatial distributions and luminosity functions are used to simulate bursts, which are then folded through mathematical models of BATSE selection effects. The observed and computed angular intensity distributions are analyzed using modifications of standard statistical homogeneity and isotropy studies. Analysis of the BATSE angular and intensity distributions greatly constrains the origins and luminosities of burst sources. In particular, it appears that no single population of sources confined to a Galactic disk, halo, or localized spiral arm satisfactorily explains BATSE observations and that effects of the burst luminosity function are secondary when considering such models. One family of models that still satisfies BATSE observations comprises sources located in an extended spherical Galactic corona. Coronal models are limited to small ranges of burst luminosity and core radius, and the allowed parameter space for such models shrinks with each new burst BATSE observes. Multiple-population models of bursts are found to work only if (1) the primary population accounts for the general isotropy and inhomogeneity seen in the BATSE observations and (2) secondary populations either have characteristics similar to the primary population or contain numbers that are small relative to the primary population.

Hakkila, Jon↗

Simulation of Propellant Loading System Senior Design Implement in Computer Algorithm

Propellant loading from the Storage Tank to the External Tank is one of the very important and time consuming pre-launch ground operations for the launch vehicle. The propellant loading system is a complex integrated system involving many physical components such as the storage tank filled with cryogenic fluid at a very low temperature, the long pipe line connecting the storage tank with the external tank, the external tank along with the flare stack, and vent systems for releasing the excess fuel. Some of the very important parameters useful for design purpose are the prediction of pre-chill time, loading time, amount of fuel lost, the maximum pressure rise etc. The physics involved for mathematical modeling is quite complex due to the fact the process is unsteady, there is phase change as some of the fuel changes from liquid to gas state, then conjugate heat transfer in the pipe walls as well as between solid-to-fluid region. The simulation is very tedious and time consuming too. So overall, this is a complex system and the objective of the work is student's involvement and work in the parametric study and optimization of numerical modeling towards the design of such system. The students have to first become familiar and understand the physical process, the related mathematics and the numerical algorithm. The work involves exploring (i) improved algorithm to make the transient simulation computationally effective (reduced CPU time) and (ii) Parametric study to evaluate design parameters by changing the operational conditions

Bandyopadhyay, Alak↗

A Computer Model for Analyzing Volatile Removal Assembly

A computer model simulates reactional gas/liquid two-phase flow processes in porous media. A typical process is the oxygen/wastewater flow in the Volatile Removal Assembly (VRA) in the Closed Environment Life Support System (CELSS) installed in the International Space Station (ISS). The volatile organics in the wastewater are combusted by oxygen gas to form clean water and carbon dioxide, which is solved in the water phase. The model predicts the oxygen gas concentration profile in the reactor, which is an indicator of reactor performance. In this innovation, a mathematical model is included in the computer model for calculating the mass transfer from the gas phase to the liquid phase. The amount of mass transfer depends on several factors, including gas-phase concentration, distribution, and reaction rate. For a given reactor dimension, these factors depend on pressure and temperature in the reactor and composition and flow rate of the influent.

Guo, Boyun↗

External field effects on solidification - Macroscopic and microscopic models

Comparison of theoretical effects of magnetic and gravitational fields on melts of metals and semiconductors provides two kinds of insight into space processing materials problems. First, the task of comparing effects of the two fields yields mathematical models which may be used to estimate material properties and responses to processing under various field conditions. Then, by using these models, microgravity effects can be inferred from magnetic field data already obtained in laboratories on earth. Magnetohydrodynamics (static and oscillating fields) and the Free Volume Model are used to study macroscopic effects, and Sekerka's interface stability theory is used to show how external fields can affect microsegregation. It is demonstrated that the mathematical formalism is essentially the same for describing macroscopic effects of both microgravity and magnetic fields.

Miller, R. I.↗