Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Algebraic models”

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 541 records · Page 30

Third-moment closure of turbulence for predictions of separating and reattaching shear flows: A study of Reynolds-stress closure model

A numerical study of computations in backward-facing steps with flow separation and reattachment, using the Reynolds stress closure is presented. The highlight of this study is the improvement of the Reynold-stress model (RSM) by modifying the diffusive transport of the Reynolds stresses through the formulation, solution and subsequent incorporation of the transport equations of the third moments, bar-u(i)u(j)u(k), into the turbulence model. The diffusive transport of the Reynolds stresses, represented by the gradients of the third moments, attains greater significance in recirculating flows. The third moments evaluated by the development and solution of the complete transport equations are superior to those obtained by existing algebraic correlations. A low-Reynolds number model for the transport equations of the third moments is developed and considerable improvement in the near-wall profiles of the third moments is observed. The values of the empirical constants utilized in the development of the model are recommended. The Reynolds-stress closure is consolidated by incorporating the equations of k and e, containing the modified diffusion coefficients, and the transport equations of the third moments into the Reynolds stress equations. Computational results obtained by the original k-e model, the original RSM and the consolidated and modified RSM are compared with experimental data. Overall improvement in the predictions is seen by consolidation of the RMS and a marked improvement in the profiles of bar-u(i)u(j)u(k) is obtained around the reattachment region.

Amano, R. S.↗

Numerical investigation of turbulence models for shock separated boundary-layer flows

Numerical solutions of the Navier-Stokes equations for shock separated turbulent boundary-layer flows are presented. Several turbulence models are investigated and assessed by their ability to predict the physical phenomena associated with two extensively documented experiments. The experimental flows consist of shock-wave boundary-layer interactions in axisymmetric internal and external geometries at Mach numbers of 1.5 and 7, respectively. Algebraic and one-equation eddy viscosity models are used to describe the Reynolds shear stress. Calculated values of skin friction, wall pressure distribution, kinetic energy of turbulence, and heat transfer are compared with measurements.

Viegas, J. R.↗

A moving-window parameter adaptive control system for the F8-DFBW aircraft

This paper describes the moving-window parameter adaptive control system developed for the NASA F8-DFBW aircraft. The control system employs a parameter identification process that, iteratively, adjusts parameters of a model of the aircraft motions in a batch processing manner so that responses generated from the model fit the outputs of sensors stored in a finite record referred to as the moving window. Tests are made on the validity of the parameter estimates before using the parameters in an on-line design process. The on-line design process is an algebraic mapping of the parameters of the model into primary control system feedback and feedforward gains. The mapping was selected to satisfy specific flying quality characteristics over the range of parameter variations expected. Results are presented from simulation studies on the identification algorithm, tests for parameter validity, and the on-line design process.

Montgomery, R. C.↗

A moving window parameter adaptive control system for the F8-DFBW aircraft

This paper describes the moving window parameter adaptive control system developed for the NASA F8-DFBW aircraft. The control system employs a parameter identification process that, iteratively, adjusts parameters of a model of the aircraft motions in a batch-processing manner so that responses generated from the model fit the outputs of sensors stored in a finite record referred to as the moving window. Tests are made on the validity of the parameter estimates before using the parameters in an on-line design process. The on-line design process is an algebraic mapping of the parameters of the model into primary control system feedback and feedforward gains. The mapping was selected to satisfy specific flying quality characteristics over the range of parameter variations expected. Results are presented from simulation studies on the identification algorithm made during the development of the system.

Dunn, H. J.↗

Identification of unsteady aerodynamics and aeroelastic integro-differential systems

The problem of estimating integro-differential models based on test or simulation data is dealt with. The identification techniques proposed for estimating parameters in models described by differential equations need to be considerably extended to deal with the integral terms. Conditions under which the integral terms may be approximated by algebraic values are discussed. The integro-differential models discussed are related to indicial models proposed by aerodynamicists to describe unsteady flow.

Gupta, N. K.↗

Application of Navier-Stokes code PAB3D with kappa-epsilon turbulence model to attached and separated flows

A three-dimensional Navier-Stokes solver was used to determine how accurately computations can predict local and average skin friction coefficients for attached and separated flows for simple experimental geometries. Algebraic and transport equation closures were used to model turbulence. To simulate anisotropic turbulence, the standard two-equation turbulence model was modified by adding nonlinear terms. The effects of both grid density and the turbulence model on the computed flow fields were also investigated and compared with available experimental data for subsonic and supersonic free-stream conditions.

Abdol-Hamid, Khaled S.↗

Predicting Buoyant Shear Flows Using Anisotropic Dissipation Rate Models

This paper examines the modeling of two-dimensional homogeneous stratified turbulent shear flows using the Reynolds-stress and Reynolds-heat-flux equations. Several closure models have been investigated-, the emphasis is placed on assessing the effect of modeling the dissipation rate tensor in the Reynolds-stress equation. Three different approaches are considered: one is an isotropic approach while the other two are anisotropic approaches. The isotropic approach is based on Kolmogorov's hypothesis and a dissipation rate equation modified to account for vortex stretching. One of the anisotropic approaches is based on an algebraic representation of the dissipation rate tensor, while another relies on solving a modeled transport equation for this tensor. In addition, within the former anisotropic approach, two different algebraic representations are examined one is a function of the Reynolds-stress anisotropy tensor, and the other is a function of' the mean velocity gradients. The performance of these closure models is evaluated against experimental and direct numerical simulation data of pure shear flows. pure buoyant flows and buoyant shear flows. Calculations have been carried out over a range of Richardson numbers (Ri) and two different Prandtl numbers (Pr); thus the effect of Pr on the development of counter-gradient heat flux in a stratified shear flow can be assessed. At low Ri, the isotropic model performs well in the predictions of stratified shear flows; however, its performance deteriorates as Ri increases. At high Ri, the transport equation model for the dissipation rate tensor gives the best result. Furthermore, the results also lend credence to the algebraic dissipation rate model based on the Reynolds stress anisotropy tensor. Finally, it is found that Pr has an effect on the development of counter-gradient heat flux. The calculations show that, under the action of shear, counter-gradient heat flux does not occur even at Ri = 1 in an air flow.

So, R. M. C.↗

Expedited Model Building Method for Rocket Nozzle Side Loading

Statistical analysis of test data is currently the preferred method of creating models for uncharacterized side loading in rocket nozzles. However, statistical analysis requires multiple data sets to be reliable; getting these sets requires building and testing multiple rocket motors, which is not feasible given the incredible time and cost required to do so. By using a combination of geometry, simplified finite element analysis (FEA) analogs, and algebra we can create a limited predictive model from only one test data set. Once created, the model requires only the input from one strain gauge located at each yaw and pitch attachment bracket for the actuation system. This process is not computationally intensive, therefore saving an abundance of time and budget over using full 3D FEA models to develop a similar predictive model. To develop the model, conversion factors between the strain at the gauge locations and corresponding actuator load is derived from the data. These conversion factors are used to create a “conversion ellipse” that bounds the constitutive properties of the nozzle as one lumped factor. Using a 2D cross-sectional analog of the nozzle exit plane, we create a function that allows us to determine the azimuthal location of the unexpected load based on the ratio between yaw and pitch strain readings. Combining this with the “conversion ellipse” yields a function that given the strain reading outputs an induced actuator load, regardless of the azimuthal location of the unexpected load. As previously mentioned, this predictive model is limited in its application as there is inaccuracy introduced by the use of data by only one test as well as simplifying assumptions. This approach makes a useful tool to quickly determine whether uncharacterized actuator loading events can damage actuation systems or induce loading that exceeds design limits.

Seetha A Kolli↗

Evaluation of turbulence models for three primary types of shock-separated boundary layers

Zero-equation (algebraic), one-equation (kinetic energy), and two-equation (kinetic energy plus length scale) turbulence eddy viscosity models were used in computing three basic types of shock-separated boundary-layer flows. The three basic types of shock boundary-layer interaction discussed are: (1) a normal shock wave at transonic speeds, (2) a compression corner shock at supersonic speeds, and (3) an incident oblique shock at hypersonic speeds. The models tested are simple, unmodified models used extensively for incompressible, unseparated flows. A comparison of computed and measured results for the compressible, separated flows described herein indicates that model performance is dependent on flow configuration with no distinct superiority of one model over the other for all three flow configurations.

Coakley, T. J.↗

A unified mathematical theory of electrophoretic processes

A mathematical theory is presented which shows that each of the four classical electrophoretic modes (zone electrophoresis, moving boundary electrophoresis, isotachophoresis, and isoelectric focusing) is based on the same general principles and can collectively be described in terms of a single set of equations. This model can predict the evolution of the four electrophoretic modes as a function of time. The model system is one-dimensional, neglecting the effects of electroosmosis, temperature gradients, and any bulk flows of liquid. The model is based on equations which express the components' dissociation equilibria, the mass transport due to electromigration and diffusion, electroneutrality, and the conservation of mass and charge. The model consists of a system of coupled partial differential and nonlinear algebraic equations which can be solved numerically by use of a computer. The versatility of this model was verified using an example of a three-component system containing cacodylate, tris hydroxylmethylaminomethane, and histidine. Results show that this model not only correctly predicts the characteristic features of each electrophoretic mode, but also gives details of the concentration, pH, and conductivity profiles not easily amenable to direct experimental measurement.

Bier, M.↗

A model for rotorcraft flying qualities studies

This paper outlines the development of a mathematical model that is expected to be useful for rotorcraft flying qualities research. A computer model is presented that can be applied to a range of different rotorcraft configurations. The algorithm computes vehicle trim and a linear state-space model of the aircraft. The trim algorithm uses non linear optimization theory to solve the nonlinear algebraic trim equations. The linear aircraft equations consist of an airframe model and a flight control system dynamic model. The airframe model includes coupled rotor and fuselage rigid body dynamics and aerodynamics. The aerodynamic model for the rotors utilizes blade element theory and a three state dynamic inflow model. Aerodynamics of the fuselage and fuselage empennages are included. The linear state-space description for the flight control system is developed using standard block diagram data.

Mittal, Manoj↗

Projection-Based Reduced Order Modeling for Spacecraft Thermal Analysis

This paper presents a mathematically rigorous, subspace projection-based reduced order modeling (ROM) methodology and an integrated framework to automatically generate reduced order models for spacecraft thermal analysis. Two key steps in the reduced order modeling procedure are described: (1) the acquisition of a full-scale spacecraft model in the ordinary differential equation (ODE) and differential algebraic equation (DAE) form to resolve its dynamic thermal behavior; and (2) the ROM to markedly reduce the dimension of the full-scale model. Specifically, proper orthogonal decomposition (POD) in conjunction with discrete empirical interpolation method (DEIM) and trajectory piece-wise linear (TPWL) methods are developed to address the strong nonlinear thermal effects due to coupled conductive and radiative heat transfer in the spacecraft environment. Case studies using NASA-relevant satellite models are undertaken to verify the capability and to assess the computational performance of the ROM technique in terms of speed-up and error relative to the full-scale model. ROM exhibits excellent agreement in spatiotemporal thermal profiles (<0.5% relative error in pertinent time scales) along with salient computational acceleration (up to two orders of magnitude speed-up) over the full-scale analysis. These findings establish the feasibility of ROM to perform rational and computationally affordable thermal analysis, develop reliable thermal control strategies for spacecraft, and greatly reduce the development cycle times and costs.

0000↗

Post-glacial relaxation of a viscously stratified compressible mantle

The postglacial relaxation phenomenon of a viscously compressible nonself-gravitating spherical-shell model is investigated. Analytical solutions to an exponentially depth-dependent viscosity are developed with various types of compressible density models such as those with an exponential dependence of the radius and an algebraic-root dependence of the radius. The solutions of a viscously compressible multilayered model with constant thermodynamic properties and viscosity are developed by the propagator matrix method. The results show that inferences of deep mantle viscosity from postglacial rebound would be hampered by mantle compressibility for long-wavelength harmonics because of the smaller excitation of compressible eigenfunctions in the lower mantle. The relaxation times and velocity fields are more sensitive to higher viscosity contrast for the exponentially varying viscosity than for models with discrete jumps in the viscosity structure.

Wu, J.↗

Investigation of third-order closure model of turbulence for the computation of incompressible flows in a channel with a backward-facing step

To predict the diffusion process of the Reynolds stresses in reattaching shear flows, the transport model for the triple-velocity products has been developed and tested for the computation of the flow in a channel with a backward-facing step. Upon comparison of the results of uuv, uvv, and vvv with those obtained by using existing algebraic correlations, it was shown that the present model improved the prediction of the triple-velocity products.

Amano, R. S.↗

Microwave Scattering Model for Grass Blade Structures

In this paper, the electromagnetic scattering solution for a grass blade with complex cross-section geometry is considered. It is assumed that the blade cross section is electrically small, but its length is large compared to the incident wavelength. In a recent study it has been shown that the scattering solution for such problems, in the form of a polarizability tensor, can be obtained using the low-frequency approximation in conjunction with the method of moments. In addition, the study shows that the relationship between the polarizability tensor of a dielectric cylinder and its dielectric constant can be approximated by a simple algebraic expression. The results of this study are used to show that this algebraic approximation is valid also for cylinders with cross sections the shape of grass blades, providing that proper values am selected for each of three constants appearing in the expression. These constants are dependent on cylinder shape, and if the relationship between the constants and the three parameters describing a grass blade shape can be determined, an algebraic approximation relating polarizability tensor to blade shape, as well as dielectric constant, can be formed. Since the elements of the polarizability tensor are dependent on only these parameters, this algebraic approximation can replace the cumbersome method of moments model. A conjugate gradient method is then implemented to correctly determine the three constants of the algebraic approximation for each blade shape. A third-order polynomial fit to the data is then determined for each constant, thus providing a complete analytic replacement to the numerical (moment method) scattering model. Comparisons of this approximation to the numerical model show an average error of less than 3%.

Stiles, James M.↗

Numerical solution for propagation of coupled longitudinal and transverse waves normal to the applied magnetic field in a three-fluid medium.

The solution was obtained without approximating the dispersion relation. This was made possible by the use of a computer method for carrying out the extremely large number of algebraic manipulations involved. The three-fluid model employed accounts for the interaction of the electron, ion, and neutral species. The solution for the e-folding distance describes the damping characteristics of each wave mode. Each of the e-folding solutions corresponds to one of the phase velocity solutions.

Mcclendon, R. W.↗

Simulation of a steady-state integrated human thermal system.

The mathematical model of an integrated human thermal system is formulated. The system consists of an external thermal regulation device on the human body. The purpose of the device (a network of cooling tubes held in contact with the surface of the skin) is to maintain the human body in a state of thermoneutrality. The device is controlled by varying the inlet coolant temperature and coolant mass flow rate. The differential equations of the model are approximated by a set of algebraic equations which result from the application of the explicit forward finite difference method to the differential equations. The integrated human thermal system is simulated for a variety of combinations of the inlet coolant temperature, coolant mass flow rate, and metabolic rates. Two specific cases are considered: (1) the external thermal regulation device is placed only on the head and (2) the devices are placed on the head and the torso. The results of the simulation indicate that when the human body is exposed to hot environment, thermoneutrality can be attained by localized cooling if the operating variables of the external regulation device(s) are properly controlled.

Hsu, F. T.↗

Simulation of a steady-state integrated human thermal system.

The mathematical model of an integrated human thermal system is formulated. The system consists of an external thermal regulation device on the human body. The purpose of the device (a network of cooling tubes held in contact with the surface of the skin) is to maintain the human body in a state of thermoneutrality. The device is controlled by varying the inlet coolant temperature and coolant mass flow rate. The differential equations of the model are approximated by a set of algebraic equations which result from the application of the explicit forward finite difference method to the differential equations. The integrated human thermal system is simulated for a variety of combinations of the inlet coolant temperature, coolant mass flow rate, and metabolic rates.

Hsu, F. T.↗