SEARCH · Search NASA

Results for “3D equilibrium”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

The Thickness of Electric Current Sheets

The thickness of electric current sheets is extremely important, especially as it relates to the onset of fast magnetic reconnection. Onset determines how much magnetic free energy can build up in a field before it is explosively released. This has implications for coronal heating, CMEs, and many other phenomena on the Sun and throughout the universe. Significant effort has been devoted to the question of whether equilibrium current sheets in realistic geometries have finite or zero thickness. Using a simple force balance analysis, we show why current sheets without a guide field (2D) and with a guide field that is invariant in the guide field direction (2.5D) cannot be in equilibrium if they have both finite thickness and finite width. We then estimate the conditions under which the tension of a curved line-tied guide field can produce a 3D equilibrium sheet that is finite in all dimensions. Finally, we speculate that some quasi-statically evolving current sheets undergoing slow stressing may reach a critical point at which equilibrium is no longer possible. They then rapidly thin and reconnect.

James A Klimchuk↗

NASA NTRS1992

Bending response of thin-walled laminated composite cylinders

Closed-form linear response solutions are presented for the bending response of thin-walled laminated composite cylinders. Displacements and intralaminar stresses are computed using the Donnell theory, and interlaminar stresses are computed using 3D equilibrium equations of elasticity. Laminates with stacking sequencies and an aluminum cylinder are used to illustrate the displacement response, intralaminar, and interlaminar stresses for two length-to-radius ratios of 1 and 5. Data obtained indicate that if only fiber breakage causes ultimate failure, material failure will occur on the compression side of the cylinders, at the ends, due to fiber compression. The interlaminar normal stress is an order of magnitude smaller than the interlaminar shear stresses and it is not likely to initiate failure.

Fuchs, H. P.↗

NASA NTRS1975

Photoionization of barium clouds via the 3D metastable levels

The photoionization of optically thin barium clouds is analyzed and shown to occur primarily by a two-step process involving the 3D metastable term as the intermediate state. The equilibrium populations of the 1D and 3D metastable levels are calculated and found to differ significantly from the values now in the literature. These populations are combined with newly available photoionization data to calculate the resulting photoionization rate in the upper atmosphere.

Carlsten, J. L.↗

NASA NTRSDate not supplied

Turbo-Design: Open-Source Radial Equilibrium Turbomachinery Solver: Part I - Turbines

Advances in 3D Geometrical Designs and Cooling have played a significant role in improving the efficiency of turbomachinery. However, these advancements must be effectively translated back to the modeler. Machine learning can facilitate this transition. Specifically, machine learning–based loss models can bridge the gap between 3D and 1D designs, enabling modelers not only to predict velocity triangles but also to extract additional geometric features. Currently, the design tools used at NASA have not been updated to support such integration—until now. TurboDesign is an open-source, Python-based framework that replaces TD2 (LEW-11029-1) and AXOD2 (LEW-16323-1), both of which are radial equilibrium solvers for axial turbines. The goal of this update is to enable the integration of machine learning loss models into radial equilibrium equations. Additionally, TurboDesign is designed to support radial machines. This paper presents the governing equations, the assumptions underlying the code, the integration of legacy loss models, an example of machine learning model integration, and a validation comparison with CFD. All code, tutorials, and documentation are available at: https://www.github.com/nasa/turbo-design

Radial Equilibrium↗

NASA NTRS1989

Application of unstructured grid methods to steady and unsteady aerodynamic problems

The purpose is to describe the development of unstructured grid methods which have several advantages when compared to methods which make use of structured grids. Unstructured grids, for example, easily allow the treatment of complex geometries, allow for general mesh movement for realistic motions and structural deformations of complete aircraft configurations which is important for aeroelastic analysis, and enable adaptive mesh refinement to more accurately resolve the physics of the flow. Steady Euler calculations for a supersonic fighter configuration to demonstrate the complex geometry capability; unsteady Euler calculations for the supersonic fighter undergoing harmonic oscillations in a complete-vehicle bending mode to demonstrate the general mesh movement capability; and vortex-dominated conical-flow calculations for highly-swept delta wings to demonstrate the adaptive mesh refinement capability are discussed. The basic solution algorithm is a multi-stage Runge-Kutta time-stepping scheme with a finite-volume spatial discretization based on an unstructured grid of triangles in 2D or tetrahedra in 3D. The moving mesh capability is a general procedure which models each edge of each triangle (2D) or tetrahedra (3D) with a spring. The resulting static equilibrium equations which result from a summation of forces are then used to move the mesh to allow it to continuously conform to the instantaneous position or shape of the aircraft. The adaptive mesh refinement procedure enriches the unstructured mesh locally to more accurately resolve the vortical flow features. These capabilities are described in detail along with representative results which demonstrate several advantages of unstructured grid methods. The applicability of the unstructured grid methodology to steady and unsteady aerodynamic problems and directions for future work are discussed.

Batina, John T.↗

NASA NTRS1994

Three-dimensional radiative transfer on a massively parallel computer

We perform 3D radiative transfer calculations in non-local thermodynamic equilibrium (NLTE) in the simple two-level atom approximation on the Mas-Par MP-1, which contains 8192 processors and is a single instruction multiple data (SIMD) machine, an example of the new generation of massively parallel computers. On such a machine, all processors execute the same command at a given time, but on different data. To make radiative transfer calculations efficient, we must re-consider the numerical methods and storage of data. To solve the transfer equation, we adopt the short characteristic method and examine different acceleration methods to obtain the source function. We use the ALI method and test local and non-local operators. Furthermore, we compare the Ng and the orthomin methods of acceleration. We also investigate the use of multi-grid methods to get fast solutions for the NLTE case. In order to test these numerical methods, we apply them to two problems with and without periodic boundary conditions.

Vath, H. M.↗

NASA NTRS2012

Axisymmetric Plume Simulations with NASA's DSMC Analysis Code

A comparison of axisymmetric Direct Simulation Monte Carlo (DSMC) Analysis Code (DAC) results to analytic and Computational Fluid Dynamics (CFD) solutions in the near continuum regime and to 3D DAC solutions in the rarefied regime for expansion plumes into a vacuum is performed to investigate the validity of the newest DAC axisymmetric implementation. This new implementation, based on the standard DSMC axisymmetric approach where the representative molecules are allowed to move in all three dimensions but are rotated back to the plane of symmetry by the end of the move step, has been fully integrated into the 3D-based DAC code and therefore retains all of DAC s features, such as being able to compute flow over complex geometries and to model chemistry. Axisymmetric DAC results for a spherically symmetric isentropic expansion are in very good agreement with a source flow analytic solution in the continuum regime and show departure from equilibrium downstream of the estimated breakdown location. Axisymmetric density contours also compare favorably against CFD results for the R1E thruster while temperature contours depart from equilibrium very rapidly away from the estimated breakdown surface. Finally, axisymmetric and 3D DAC results are in very good agreement over the entire plume region and, as expected, this new axisymmetric implementation shows a significant reduction in computer resources required to achieve accurate simulations for this problem over the 3D simulations.

Stewart, B. D.↗

NASA NTRS1992

Three-dimensional structure of a curved wake

The effects of streamwise curvature on the 3D structure of a plane wake are studied. Preliminary results indicate that the wake is initially grossly 3D with large spanwise distortions in the mean velocity and Reynolds stress contours. Local extrema distributed in a quasi-periodic manner are seen. The spanwise variations are caused by the presence of relatively strong quadrupoles of mean streamwise vorticity. Further downstream the mean vorticity decays at a rate lower on the unstable side of the wake so that only the unstable side row of vortices persist to the far-field region. As a result, the spanwise distortions and local extrema seem to persist with a relatively slow decay. The effects of curvature are also apparent in the Reynolds stress results. In particular, the shear stress is significantly affected by curvature right from the start such that significantly higher magnitudes are observed on the unstable side. The curved wake clearly develops a very 3D asymmetric structure which shows no signs of reaching an equilibrium state.

Weygandt, J. H.↗

NASA NTRS1993

A comparison of the predictive capabilities of several turbulence models using upwind and central-difference computer codes

Four turbulence models are described and evaluated for transonic flows using the upwind code CFL3D and the central-difference code TLNS3D. In particular, the effects of recent modifications to the half-equation model of Johnson-King are explored in detail, and different versions of the model are compared. This model can obtain good results for both two-dimensional (2D) and three-dimensional (3D) separated flows. The one-equation models of Baldwin-Barth and Spalart-Allmaras perform well for separated airfoil flows, but can predict the shock too far forward at the outboard stations of a separated wing. The equilibrium model of Baldwin-Lomax predicts the shock location too far aft for both 2D and 3D separated flows, as expected. In general, all models perform well for attached or mildly separated flows.

Rumsey, Christopher L.↗

NASA NTRS1993

A Comparison of the Predictive Capabilities of Several Turbulence Models using Upwind and Central-Difference Computer Codes

Rumsey, Christopher L.↗

NASA NTRS2014

Thermal Non-equilibrium Consistent with Widespread Cooling

Time correlation analysis has been used to show widespread cooling in the solar corona; this cooling has been interpreted as a result of impulsive (nanoflare) heating. In this work, we investigate wide-spread cooling using a 3D model for a solar active region which has been heated with highly stratified heating. This type of heating drives thermal non-equilibrium solutions, meaning that though the heating is effectively steady, the density and temperature in the solution are not. We simulate the expected observations in narrowband EUV images and apply the time correlation analysis. We find that the results of this analysis are qualitatively similar to the observed data. We discuss additional diagnostics that may be applied to differentiate between these two heating scenarios.

Winebarger, A.↗

NASA NTRS2014

Thermal Non-equilibrium Consistent with Widespread Cooling

Winebarger, A.↗

NASA NTRS1989

Size distribution of dust grains: A problem of self-similarity

Distribution functions describing the results of natural processes frequently show the shape of power laws, e.g., mass functions of stars and molecular clouds, velocity spectrum of turbulence, size distributions of asteroids, micrometeorites and also interstellar dust grains. It is an open question whether this behavior is a result simply coming about by the chosen mathematical representation of the observational data or reflects a deep-seated principle of nature. The authors suppose the latter being the case. Using a dust model consisting of silicate and graphite grains Mathis et al. (1977) showed that the interstellar extinction curve can be represented by taking a grain radii distribution of power law type n(a) varies as a(exp -p) with 3.3 less than or equal to p less than or equal to 3.6 (example 1) as a basis. A different approach to understanding power laws like that in example 1 becomes possible by the theory of self-similar processes (scale invariance). The beta model of turbulence (Frisch et al., 1978) leads in an elementary way to the concept of the self-similarity dimension D, a special case of Mandelbrot's (1977) fractal dimension. In the frame of this beta model, it is supposed that on each stage of a cascade the system decays to N clumps and that only the portion beta N remains active further on. An important feature of this model is that the active eddies become less and less space-filling. In the following, the authors assume that grain-grain collisions are such a scale-invarient process and that the remaining grains are the inactive (frozen) clumps of the cascade. In this way, a size distribution n(a) da varies as a(exp -(D+1))da (example 2) results. It seems to be highly probable that the power law character of the size distribution of interstellar dust grains is the result of a self-similarity process. We can, however, not exclude that the process leading to the interstellar grain size distribution is not fragmentation at all. It could be, e.g., diffusion-limited growth discussed by Sander (1986), who applied the theory of fractal geometry to the classification of non-equilibrium growth processes. He received D=2.4 for diffusion-limited aggregation in 3d-space.

Henning, TH.↗

NASA NTRS1991

Comparison of hypersonic experiments and PNS predictions. I - Aerothermodynamics. II - Aerodynamics

A 3D parabolized Navier-Stokes (PNS) technique is examined in terms of applications to equilibrium-air and perfect-gas hypersonic flows at different Mach numbers and angles of attack. The method renders sublayer approximation unnecessary, and treats axis-normal coupling effects first followed by a shock solution and a corrector step. Predictions of hypersonic flowfields over several blunt conical configurations are found to agree with established wind-tunnel and flight data, although those for the Dhawan-Narasimha (1958) transition model vary from the experimental data. The PNS technique is also applied to predictions of wall-pressure distributions and force and moment data, and the results agree with experimental values. The PNS method allows the calculation of a wide variety of configurations with only a slight increase over traditional computing times.

Bhutta, Bilal A.↗

NASA NTRS1992

A three-dimensional upwind PNS code for chemically reacting scramjet flowfields

A new upwind, parabolized Navier-Stokes (PNS) code has been developed to compute the three-dimensional (3D) chemically reacting flow in scramjet (supersonic combustion ramjet) engines. The code is a modification of the 3D upwind PNS (UPS) airflow code which has been extended in the present study to permit internal flow calculations with hydrogen-air chemistry. With these additions, the new code has the capability of computing aerodynamic and propulsive flowfields simultaneously. The algorithm solves the PNS equations using a finite-volume, upwind TVD method based on Roe's approximate Riemann solver that has been modified to account for 'real gas' effects. The fluid medium is assumed to be a chemically reacting mixture of thermally perfect (but calorically imperfect) gases in thermal equilibrium. The new code has been applied to two test cases. These include the Burrows-Kurkov supersonic combustion experiment and a generic 3D scramjet flowfield. The computed results compare favorably with the available experimental data.

Wadawadigi, Ganesh↗

NASA NTRS1989

Progress in computing nozzle/plume flow fields

The long-term goal is to develop the capability to predict chemically-reacting, multi-stream nozzle and plume flow fields. Two basic Navier-Stokes solvers, including the widely used F-3D code, are upgraded to include several upwind difference schemes and portable chemistry packages. Current computational capabilities for solving equilibrium single-stream and multi-stream, frozen gas, and finite rate chemistry problems are described. A variety of complex nozzle and plume flows were computed. Solutions presented include axisymmetric plume flow for ideal and equilibrium air, 3-D NASP nozzle/afterbody flow, and an internal nozzle calculation comparing various finite-rate chemistry packages.

Ruffin, Stephen M.↗

NASA NTRS1993

A global three-dimensional model of the stratospheric sulfuric acid layer

A 3D model which encompasses SO2 production from OCS, followed by its oxidation to gaseous H2SO4, the condensation-evaporation equilibrium of gaseous and particulate H2SO4, and finally particle condensation and rainout, is presently used to study processes maintaining the nonvolcanically-perturbed stratosphere's sulfuric acid layer. A comparison of the results thus obtained with remotely sensed stratospheric aerosol extinction data shows the model to simulate the general behavior of stratospheric aerosol extinction.

Golombek, Amram↗

NASA NTRS2008

The CME Flare Arcade and the Width of the CME in the Outer Corona

Moore, Sterling, & Suess (2007, ApJ, 668, 1221) present evidence that (1) a CME is typically a magnetic bubble, a low-beta gplasmoid with legs h having roughly the 3D shape of a light bulb, and (2) in the outer corona the CME plasmoid is in lateral pressure equilibrium with the ambient magnetic field. They present three CMEs observed by SOHO/LASCO, each from a very different source located near the limb. One of these CMEs came from a compact ejective eruption from a small part of a sunspot active region, another came from a large quiet-region filament eruption, and the third CME, an extremely large and fast one, was produced in tandem with an X20 flare arcade that was centered on a huge delta sunspot. Each of these CMEs had more or less the classic lightbulb silhouette and attained a constant heliocentric angular width in the outer corona. This indicates that the CME plasmoid attained lateral magnetic pressure balance with the ambient radial magnetic field in the outer corona. This lateral pressure balance, together with the standard scenario for CME production by the eruption of a sheared-core magnetic arcade, yields the following simple estimate of the strength B(sub Flare) of the magnetic field in the flare arcade produced together with the CME: B(sub Flare) 1.4(theta CME/theta Flare)sup 2 G, where theta (sub CME) is the heliocentric angular width of the CME plasmoid in the outer corona and theta (sub Flare) is the heliocentric angular width of the full-grown flare arcade. Conversely, theta (sub CME) approximately equal to (R(sub Sun)sup -1(phi(sub Flare)/1.4)sup 1/2 radians, where Flare is the magnetic flux covered by the full-grown flare arcade. In addition to presenting the three CMEs of Moore, Sterling, & Suess (2007) and their agreement with this relation between CME and Flare, we present a further empirical test of this relation. For CMEs that erupt from active regions, the co-produced flare arcade seldom if ever covers the entire active region: if AR is the total magnetic flux of the active region, Flare . AR, and we predict that CME. (R(sub Sun))sup -1(theta AR/1.4)sup 1/2 radians. For a random sample of 31 CMEs that erupted from active regions within 30 of the limb, for each CME we have measured CME from LASCO/C3 and have measured AR from a SOHO/MDI magnetogram of the source active region when it was within 30 of disk center. We find that each CME obeys the above predicted inequality, none having width greater than half of the upper bound given by theta(sub AR). Thus, an active region's magnetic flux content, together with its location on the solar disk, largely determines whether the active region can possibly produce a CME that is wide enough to intercept the Earth.

Moore, Ron↗