Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Maxwell equations”

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 325 records · Page 18

Computational fluid dynamics near the continuum limit

Two dimensional Navier-Stokes equations for a perfect gas are solved for hypersonic flow over a flat plate at Reynolds numbers ranging from continuum to the free molecule flow. In the transition flow regime, new slip boundary conditions are introduced, which reduce to the well-known slip conditions of Maxwell at small Knudsen numbers and yield the correct shear stress and heat transfer in the limiting case of free molecule flow. Comparison of the computed results with the existing experimental data and Monte Carlo calculations indicates that the continuum Navier-Stokes equations give surprisingly realistic results throughout the transition flow regime, when the new slip boundary conditions are used.

Maccormack, Robert W.↗

A Field-Particle Correlation Analysis of a Perpendicular Magnetized Collisionless Shock

Using the field-particle correlation technique, we examine the particle energization in a 1D-2V continuum Vlasov-Maxwell simulation of a perpendicular magnetized collisionless shock. The combination of the field-particle correlation technique with the high fidelity representation of the particle distribution function provided by a direct discretization of the Vlasov equation allows us to ascertain the details of the exchange of energy between the electromagnetic fields and the particles in phase space. We identify the velocity-space signatures of shock-drift acceleration of the ions and adiabatic heating of the electrons due to the perpendicular collisionless shock by constructing a simplified model with the minimum ingredients necessary to produce the observed energization signatures in the self-consistent Vlasov-Maxwell simulation. We are thus able to completely characterize the energy transfer in the perpendicular collisionless shock considered here and provide predictions for the application of the field-particle correlation technique to spacecraft measurements of collisionless shocks.

James Juno↗

Viscous Effects in the Elastodynamics of Thick Beams

A viscoelastic higher-order thick beam finite element formulation is extended to include elastodynamic deformations. The material constitutive law is a special differential form of the Maxwell solid. In the constitutive model, the elastic strains and the conjugate viscous strains are coupled through a system of first- order ordinary differential equations. The total time-dependent stress is the superposition of its elastic and viscous components. The elastodynamic equations of motion are derived from the virtual work principle. Computational examples are carried out for a thick orthotropic cantilevered beam. A quasi-static relaxation problem is employed as a validation test for the elastodynamic algorithm. The elastodynamic code is demonstrated by analyzing the damped vibrations of the beam which is deformed and then released to freely vibrate.

Johnson, A. R.↗

A viscoelastic higher-order beam finite element

A viscoelastic internal variable constitutive theory is applied to a higher-order elastic beam theory and finite element formulation. The behavior of the viscous material in the beam is approximately modeled as a Maxwell solid. The finite element formulation requires additional sets of nodal variables for each relaxation time constant needed by the Maxwell solid. Recent developments in modeling viscoelastic material behavior with strain variables that are conjugate to the elastic strain measures are combined with advances in modeling through-the-thickness stresses and strains in thick beams. The result is a viscous thick-beam finite element that possesses superior characteristics for transient analysis since its nodal viscous forces are not linearly dependent an the nodal velocities, which is the case when damping matrices are used. Instead, the nodal viscous forces are directly dependent on the material's relaxation spectrum and the history of the nodal variables through a differential form of the constitutive law for a Maxwell solid. The thick beam quasistatic analysis is explored herein as a first step towards developing more complex viscoelastic models for thick plates and shells, and for dynamic analyses. The internal variable constitutive theory is derived directly from the Boltzmann superposition theorem. The mechanical strains and the conjugate internal strains are shown to be related through a system of first-order, ordinary differential equations. The total time-dependent stress is the superposition of its elastic and viscous components. Equations of motion for the solid are derived from the virtual work principle using the total time-dependent stress. Numerical examples for the problems of relaxation, creep, and cyclic creep are carried out for a beam made from an orthotropic Maxwell solid.

Johnson, Arthur R.↗

Elasto visco-plastic flow with special attention to boundary conditions

A simple but nontrivial steady-state creeping elasto visco-plastic (Maxwell fluid) radial flow problem is analyzed, with special attention given to the effects of the boundary conditions. Solutions are obtained through integration of a governing equation on stress using the Runge-Kutta method for initial value problems and finite differences for boundary value problems. A more general approach through the finite element method, an approach that solves for the velocity field rather than the stress field and that is applicable to a wide range of problems, is presented and tested using the radial flow example. It is found that steady-state flows of elasto visco-plastic materials are strongly influenced by the state of stress of material as it enters the region of interest. The importance of this boundary or initial condition in analyses involving materials coming into control volumes from unusual stress environments is emphasized.

Shimazaki, Y.↗

Constitutive equations of ageing polymeric materials

The constitutive equation for the relaxation behavior of time-dependent, chemically unstable materials developed by Valanis and Peng (1983), which used the irreversible thermodynamics of internal variables in Eyring's absolute reaction theory and yielded a theoretical expression for the effect of chemical crosslink density on the relaxation rate, is presently applied to the creep behavior of a network polymer which is undergoing a scission process. In particular, two equations are derived which may for the first time show the relations between mechanical models and internal variables in the creep expressions, using a three-element model with a Maxwell element.

Peng, S. T. J.↗

Classical geometries defined by exterior differential systems on higher frame bundles

Exterior differential ideals are discussed, and sets of invariant generators presented, for Reimannian, conformal and projective geometries, and for specializations such as Ricci-flat, self-dual and Einstein-Maxwell theories. The Cartan characteristic integers are explicitly calculated, and involutory basis forms found, for each of these (specialized to four dimensions), exposing their algebraic structure and showing how they generate well-posed sets of partial differential equations.

Estabrook, Frank B.↗

Ionospheric plasma flow about a system of electrically biased flat plates

The steady state interaction of two electrically biased parallel plates immersed in a flowing plasma characteristic of low earth orbit is studied numerically. Fluid equations are developed to describe the motion of the cold positively charged plasma ions, and are solved using finite-differences in two dimensions on a Cartesian grid. The behavior of the plasma electrons is assumed to be described by the Maxwell-Boltzmann distribution. Results are compared to an analytical and a particle simulation technique for a simplified flow geometry consisting of a single semi-infinite negatively biased plate. Comparison of the extent of the electrical disturbance into the flowing plasma and the magnitude of the current collected by the plate is very good. The interaction of two equally biased parallel plates is studied as a function of applied potential. The separation distance at which the current collected by either plate decreases by five and twenty percent is determined as a function of applied potential. The percent decreases were based on a non-interacting case. The decrease in overall current is caused by a decrease in ionic density in the region between the plates. As the separation between the plates decreases, the plates collect the ions at a faster rate than they are supplied to the middle region by the oncoming plasma flow. The docking of spacecraft in orbit is simulated by moving two plates of unequal potential toward one another in a quasi-static manner. One plate is held at a large negative potential while the other floats electrically in the resulting potential field. It is found that the floating plate does not charge continuously negative as it approaches the other more negatively biased plate. Instead, it charges more and then less negative as ionic current decreases and then increases respectively upon approach. When the two plates come into contact, it is expected that the electrically floating plate will charge rapidly negative to a potential near that of the other plate.

Herr, Joel L.↗

The sun's magnetic sector structure

The synoptic appearance of solar magnetic sectors is studied using 454 sector boundaries observed at earth during 1959-1973. The sectors are clearly visible in the photospheric magnetic field. Sector boundaries can be clearly identified as north-south running demarcation lines between regions of persistent magnetic polarity imbalances. These regions extend up to about 35 deg of latitude on both sides of the equator. They generally do not extend into the polar caps. The polar cap boundary can be identified as an east-west demarcation line marking the poleward limit of the sectors. The typical flux imbalance for a magnetic sector is about 4 x 10 to the 21st power Maxwells.

Svalgaard, L.↗

High-order Accurate Finite-Volume Formulations for the Pressure Gradient Force in Layered Ocean Models

Discretisation of the horizontal pressure gradient force in layered ocean models is a challenging task, with non-trivial interactions between the thermodynamics of the fluid and the geometry of the layers often leading to numerical difficulties. We present two new finite-volume schemes for the pressure gradient operator designed to address these issues. In each case, the horizontal acceleration is computed as an integration of the contact pressure force that acts along the perimeter of an associated momentum control-volume. A pair of new schemes are developed by exploring different control-volume geometries. Non-linearities in the underlying equation-of-state definitions and thermodynamic profiles are treated using a high-order accurate numerical integration framework, designed to preserve hydrostatic balance in a non-linear manner. Numerical experiments show that the new methods achieve high levels of consistency, maintaining hydrostatic and thermobaric equilibrium in the presence of strongly-sloping layer geometries, non-linear equations-of-state and non-uniform vertical stratification profiles. These results suggest that the new pressure gradient formulations may be appropriate for general circulation models that employ hybrid vertical coordinates andor terrain-following representations.

pressure gradient force↗

Evaluation of approximate relations for Delta /Q/ using a numerical solution of the Boltzmann equation

Data obtained from a numerical solution of the Boltzmann equation for shock-wave structure are used to test the accuracy of accepted approximate expressions for the two moments of the collision integral Delta (Q) for general intermolecular potentials in systems with a large translational nonequilibrium. The accuracy of the numerical scheme is established by comparison of the numerical results with exact expressions in the case of Maxwell molecules. They are then used in the case of hard-sphere molecules, which are the furthest-removed inverse power potential from the Maxwell molecule; and the accuracy of the approximate expressions in this domain is gauged. A number of approximate solutions are judged in this manner, and the general advantages of the numerical approach in itself are considered.

Nathenson, M.↗

Velocity lag of solid particles in oscillating gases and in gases passing through normal shock waves

The velocity lag of micrometer size spherical particles is theoretically determined for gas particle mixtures passing through a stationary normal shock wave and also for particles embedded in an oscillating gas flow. The particle sizes and densities chosen are those considered important for laser Doppler velocimeter applications. The governing equations for each flow system are formulated. The deviation from Stokes flow caused by inertial, compressibility, and rarefaction effects is accounted for in both flow systems by use of an empirical drag coefficient. Graphical results are presented which characterize particle tracking as a function of system parameters.

Maxwell, B. R.↗

Expendable launch vehicle studies

Analytical support studies of expendable launch vehicles concentrate on the stability of the dynamics during launch especially during or near the region of maximum dynamic pressure. The in-plane dynamic equations of a generic launch vehicle with multiple flexible bending and fuel sloshing modes are developed and linearized. The information from LeRC about the grids, masses, and modes is incorporated into the model. The eigenvalues of the plant are analyzed for several modeling factors: utilizing diagonal mass matrix, uniform beam assumption, inclusion of aerodynamics, and the interaction between the aerodynamics and the flexible bending motion. Preliminary PID, LQR, and LQG control designs with sensor and actuator dynamics for this system and simulations are also conducted. The initial analysis for comparison of PD (proportional-derivative) and full state feedback LQR Linear quadratic regulator) shows that the split weighted LQR controller has better performance than that of the PD. In order to meet both the performance and robustness requirements, the H(sub infinity) robust controller for the expendable launch vehicle is developed. The simulation indicates that both the performance and robustness of the H(sub infinity) controller are better than that for the PID and LQG controllers. The modelling and analysis support studies team has continued development of methodology, using eigensensitivity analysis, to solve three classes of discrete eigenvalue equations. In the first class, the matrix elements are non-linear functions of the eigenvector. All non-linear periodic motion can be cast in this form. Here the eigenvector is comprised of the coefficients of complete basis functions spanning the response space and the eigenvalue is the frequency. The second class of eigenvalue problems studied is the quadratic eigenvalue problem. Solutions for linear viscously damped structures or viscoelastic structures can be reduced to this form. Particular attention is paid to Maxwell and Kelvin models. The third class of problems consists of linear eigenvalue problems in which the elements of the mass and stiffness matrices are stochastic. dynamic structural response for which the parameters are given by probabilistic distribution functions, rather than deterministic values, can be cast in this form. Solutions for several problems in each class will be presented.

Bainum, Peter M.↗

Momentum and energy exchange collision terms for interpenetrating bi-Maxwellian gases

For application to aeronomy and space physics problems involving strongly magnetized plasma flows, we derived momentum and energy exchange collision terms for interpenetrating bi-Maxwellian gases. Collision terms were derived for Coulomb, Maxwell molecule, and constant collision cross-section interaction potentials. The collision terms are valid for arbitrary flow velocity differences and temperature differences between the interacting gases as well as for arbitrary temperature anisotropies. The collision terms had to be evaluated numerically and the appropriate coefficients are presented in tables. However, the collision terms were also fitted with simplified expressions, the accuracy of which depends on both the interaction potential and the temperature anisotropy. In addition, we derived the closed set of transport equations that are associated with the momentum and energy collision terms.

Barakat, A. R.↗

Microphysical Retrievals Over Stratiform Rain Using Measurements from an Airborne Dual-Wavelength Radar-Radiometer

The need to understand the complementarity of the radar and radiometer is important not only to the Tropical Rain Measuring Mission (TRMM) program but to a growing number of multi-instrumented airborne experiment that combine single or dual-frequency radars with multichannel radiometers. The method of analysis used in this study begins with the derivation of dual-wavelength radar equations for the estimation of a two-parameter drop size distribution (DSD). Defining a "storm model" as the set of parameters that characterize snow density, cloud water, water vapor, and features of the melting layer, then to each storm model there will usually correspond a set of range-profiled drop size distributions that are approximate solutions of the radar equations. To test these solutions, a radiative transfer model is used to compute the brightness temperatures for the radiometric frequencies of interest. A storm model or class of storm models is considered optimum if it provides the best reproduction of the radar and radiometer measurements. Tests of the method are made for stratiform rain using simulated storm models as well as measured airborne data. Preliminary results show that the best correspondence between the measured and estimated radar profiles usually can be obtained by using a moderate snow density (0.1-0.2 g/cu cm), the Maxwell-Garnett mixing formula for partially melted hydrometeors (water matrix with snow inclusions), and low to moderate values of the integrated cloud liquid water (less than 1 kg/sq m). The storm-model parameters that yield the best reproductions of the measured radar reflectivity factors also provide brightness temperatures at 10 GHz that agree well with the measurements. On the other hand, the correspondence between the measured and modeled values usually worsens in going to the higher frequency channels at 19 and 34 GHz. In searching for possible reasons for the discrepancies, It is found that changes in the DSD parameter Mu, the radar constants, or the path-integrated attenuation can affect the high frequency channels significantly. In particular, parameters that cause only modest increases in the median mass diameter of the snow, and which have a minor effect on the radar returns or the low frequency brightness temperature, can produce a strong cooling of the 34 GHz brightness temperature.

Meneghini, Robert↗

Some Thermodynamic Considerations on the Physical and Quantum Nature of Space and Time

It is suggested that the Planck h = m(sub k)c Lambda(sub k) and the Boltzmann k = m(sub k)c nu(sub k)Constants have stochastic foundation. It is further suggested that a body of fluid at equilibrium is composed of a spectrum of molecular clusters (energy levels) the size of which are governed by the Maxwell-Boltzmann distribution function. Brownian motions are attributed to equilibrium between suspensions and molecular clusters. Atomic (molecular) transition between different size atomic- (molecular-) clusters (energy levels) is shown to result in emission/absorption of energy in accordance with Bohr's theory of atomic spectra. Physical space is identified as a tachyonic fluid that is Dirac's stochastic ether or de Broglie's hidden thermostat. Compressibility of physical space, in accordance with Planck's compressible ether, is shown to result in the Lorentz-Fitzgerald contraction, thus providing a causal explanation of relativistic effect in accordance with the perceptions of Poincare and Lorentz. The invariant Schrodinger equation is derived from the invariant Bernoulli equation for incompressible potential flow. Following Heisenberg a temporal uncertainty relation is introduced as Delta(nu(sub Beta)) Delta(Rho(sub Beta)) > = k.

Sohrab, Siavash H.↗

Two-Step Multi-Physics Analysis of an Annular Linear Induction Pump for Fission Power Systems

One of the key technologies associated with fission power systems (FPS) is the annular linear induction pump (ALIP). ALIPs are used to circulate liquid-metal fluid for transporting thermal energy from the nuclear reactor to the power conversion device. ALIPs designed and built to date for FPS project applications have not performed up to expectations. A unique, two-step approach was taken toward the multi-physics examination of an ALIP using ANSYS Maxwell 3D and Fluent. This multi-physics approach was developed so that engineers could investigate design variations that might improve pump performance. Of interest was to determine if simple geometric modifications could be made to the ALIP components with the goal of increasing the Lorentz forces acting on the liquid-metal fluid, which in turn would increase pumping capacity. The multi-physics model first calculates the Lorentz forces acting on the liquid metal fluid in the ALIP annulus. These forces are then used in a computational fluid dynamics simulation as (a) internal boundary conditions and (b) source functions in the momentum equations within the Navier-Stokes equations. The end result of the two-step analysis is a predicted pump pressure rise that can be compared with experimental data.

Fission Power Systems↗