Search NASA⌕ Search

SEARCH · Search NASA

Results for “FREE BOUNDARY”

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 127 records · Page 7

Toroidal modeling of plasma flow damping and density pump-out by RMP during ELM mitigation in HL-2A

Reduction of both the plasma density and toroidal flow speed, due to application of the predominantly n = 1 (n is the toroidal mode number) resonant magnetic perturbation (RMP) for controlling the edge localized mode in the HL-2A tokamak, is numerically investigated utilizing the quasi-linear initial-value code MARS-Q (Liu et al 2013 Phys. Plasmas 20 042503). Simulation results reveal that the neoclassical toroidal viscosity (NTV) due to three dimensional fields plays the key role in modifying the plasma momentum and particle transport in the HL-2A discharge. By comparing the modeling results with the measured density pump-out in the experiment, the electron NTV particle flux model, in combination with the free-boundary condition for the axisymmetric change of the density at the plasma edge, is found to yield the best agreement in terms of both the pump-out level and the overall time scale. Further sensitivity studies show that the simulated density pump-out level is reasonably robust against variations in the model assumptions, including the particle diffusion model and the non-ambipolar versus ambipolar NTV particle flux. The latter however affects the time scale for reaching the steady state solution. Finally, it is found that the plasma edge-peeling response, the NTV torque, as well as the plasma momentum and particle transport, all are sensitive to the toroidal phase difference between the upper and lower rows of the RMP coil currents in HL-2A, with the 30 ∘ coil phasing producing the minimal side effects on the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Current ramp-up studies for NSTX-U in support of the development of non-inductive current ramp-up scenarios in STs

The Spherical Tokamak (ST) is a low-aspect-ratio tokamak projected to enable simultaneous operation at high normalized toroidal beta and high bootstrap current fraction. However, a compact ST reactor cannot be built with a full-sized solenoid. The remaining current required for sustained operation must be ramped up non-inductively. Sustained 100% non-inductive current drive at the MA level has not yet been demonstrated on an ST. Non-inductive current ramp-up is therefore a more challenging objective that the NSTX-U research program aims to address. In support of this development, free-boundary TRANSP simulations have been used to develop a current ramp-up scenario in which a significant fraction of the current is driven non-inductively by a combination of neutral-beam and bootstrap current drive. The resulting discharge has properties favorable for non-inductive current drive, such as a low normalized internal plasma inductance of approximately 0.45, a minimum safety factor exceeding 2, and a slightly reversed-shear q-profile. In these discharges, the plasma current ramps to 900 kA within 250 ms in an L-mode configuration, with varying electron density and density profiles. These discharges should be suitable for next-step NSTX-U experiments aimed at generating a robust L-mode fiducial scenario with low central solenoid flux consumption, which could be progressively improved as new NSTX-U experimental discharges are developed.

Raman, Roger [Univ. of Washington, Seattle, WA (Un↗

A Virtual Element Method for the Full MHD system

In this manuscript we present a novel discretization for the incompressible MHD system. Our approach follows the framework of the Virtual Element Method and offers two main advantages. The method can be implemented in unstructured meshes making it highly versatile and capable of handling a wide array of problems involving interfaces, free-boundaries or adaptive refinements on the mesh. The second advantage involves the divergence of the magnetic field, our approach guarantees that it remains solenoidal. We include a theoretical proof of the condition on the magnetic field as well as energy estimates and a well-posedness study. The latter sheds light as to the stability properties of the method.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

FES 2023 Theory & Simulation Performance Target

The Report for the 2023 Theory & Simulation Performance Target (TSPT) is presented herein. The 1st, 2nd and 3rd Quarter Reports are given in Sec. I, Sec. II and Sec. III. These are unchanged from those previously submitted, except that in this final document the equations, figures and citations have been automatically renumbered. The 4th Quarter Report is given in Sec. IV. The following material presents verification calculations of several widely used three-dimensional (3D) magnetohydrodynamic (MHD) equilibrium codes in strongly shaped configurations; for both fixed and free-boundary equilibria; for the so-called vacuum case and for the case of high plasma pressure; and, for magnetic fields that have nested magnetic surfaces and for fields that have a magnetic island. The configurations studied includes the recently discovered “precise-QA” case [Phys. Rev. Lett., 128, 035001 (2002)] and other highly optimized configurations. Because of the strong (extreme) shaping, the latter of these pose challenges to the 3D MHD codes. A new method is introduced for constructing coordinates in arbitrarily shaped toroidal domains, and this can improve the robustness of 3D equilibrium calculations. Beyond the equilibrium calculation, this report demonstrates the capability of performing non-linear extended-MHD modelling of optimized stellarator equilibria. This capability allows the assessment not only of linear stability thresholds but also of the nonlinear evolution of unstable equilibria to new configurations that might be far from the original configuration, with significantly different transport properties and magnetic topology. They also enable the self-consistent determination of pressure profiles in non-integrable magnetic topology through the use of highly anisotropic thermal-transport equations. The calculations presented here represent the first demonstration of extended-MHD calculations in strongly-shaped, optimized quasi-symmetric configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

INVESTIGATION OF PROBLEMS IN THERMAL CONVECTION

Numerical solution of thermal convection equations for velocity and temperature fluctuations as well as the stability analysis and heat transport for free boundary conditions

THERMAL CONVECTION↗

A half plane and a strip with an arbitrarily located crack

A technique is presented for dealing with the problem of an elastic domain containing an arbitrarily oriented internal crack. The problem is formulated as a system of integral equations for a fictitious layer of body forces imbedded in the plane along a closed smooth curve encircling the original domain. The problems of a half plane with a crack in the neighborhood of its free boundary and of an infinite strip containing a symmetrically located internal crack with an arbitrary orientation are considered as examples. In each case the stress intensity factors are computed and are given as functions of the crack angle.

Erdogan, F.↗

Results of ground vibration tests on a YF-12 airplane

Ground vibration tests were conducted on a YF-12 airplane. To approximate a structural free-free boundary condition during the tests, each of the landing gears was supported on a support system designed to have a low natural frequency. The test equipment and the procedures used for the ground vibration tests are described. The results are presented in the form of frequency response data, measured mode lines, and elastic mode shapes for the wing/body, rudder, and fuselage ventral fin. In the frequency range between 3.4 cps and 28.8 cps, nine symmetrical wing/body modes, six antisymmetrical wing/body modes, two rudder modes, and one ventral fin mode were measured.

Wilson, R. J.↗

Vibrations of the earth's magnetic tail

The natural vibration period of the tail (regarded as bundle of magnetic line-of-force strings) and the characteristic time of recurrence of magnetospheric substorms are of the same order, suggesting the possibility of a relation between these phenomena. A model of the tail in the form of a plasma cylinder whose free boundary is immersed in the interplanetary plasma is examined. The nature of the natural oscillations of the tail excited by variations of the solar wind pressure is studied in this model. These vibrations have a total energy of approximately 10 to the 22nd power erg, which is sufficient to generate magnetospheric substorms.

Yershkovich, A. I.↗

A half plane and a strip with an arbitrarily located crack

The paper introduces a technique to deal with the problem of an elastic domain containing an arbitrarily oriented internal crack. The problem is formulated as a system of integral equations for a fictitious layer of body forces imbedded in the plane along a closed smooth curve encircling the original domain. The problems of a half plane with a crack in the neighborhood of its free boundary and of an infinite strip containing a symmetrically located internal crack with an arbitrary orientation are considered as examples. In each case the stress intensity factors are computed and are given as functions of the crack angle.

Erdogan, F.↗

Transient wave propagation in composite plates due to impact

Stress-wave propagation in a composite plate due to an impact has been examined using anisotropic elasticity theory. The plate is modelled as a number of identical anisotropic layers whose elastic moduli are given. Mindlin's (1959) approximate theory of plates is then applied to each layer to obtain a set of difference-differential equations of motion by using the interlaminar stresses and displacements as explicit variables. Dispersion relationships for harmonic waves are found when traction-free boundary conditions are applied to both surfaces of the plate. The difference-differential equations are reduced to difference equations via integral transforms. With given impact boundary conditions these equations are solved for an arbitrary number of layers in the plate, and the transient propagation of a stress wave is calculated by means of a fast Fourier transform algorithm.

Sung Kim, B.↗

Stress-free end problem in layered materials

In this paper the plane elastostatic problem for a medium which consists of periodically arranged two sets of bonded dissimilar layers or strips is considered. First it is assumed that one set of strips contains a crack which crosses the bimaterial interfaces. Then, by letting the collinear cracks join, the stress-free end problem is formulated. The singular behavior of the solutions at the point on intersection of the stress-free boundary and the interfaces is examined and appropriate stress intensity factors are defined. The results of some numerical examples are then presented which include the cases of both plane stress and plane strain.

Erdogan, F.↗

Approximate solutions to the Navier-Stokes initial value problem

A Galerkin-Ritz procedure for any arbitrary system of field equations is shown to follow generically from 'two-functional' variational conditions. The initial-value problem for boundary-free incompressible Navier-Stokes flow is solved analytically in the two-parameter-function approximation.

Rosen, G.↗

Stress intensity factors in two bonded elastic layers containing cracks perpendicular to and on the interface. Part 1: Analysis

The basic crack problem which is essential for the study of subcritical crack propagation and fracture of layered structural materials is considered. Because of the apparent analytical difficulties, the problem is idealized as one of plane strain or plane stress. An additional simplifying assumption is made by restricting the formulation of the problem to crack geometries and loading conditions which have a plane of symmetry perpendicular to the interface. The general problem is formulated in terms of a coupled system of four integral equations. For each relevant crack configuration of practical interest, the singular behavior of the solution near and at the ends and points of intersection of the cracks is investigated and the related characteristic equations are obtained. The edge crack terminating at and crossing the interface, the T-shaped crack consisting of a broken layer and a delamination crack, the cross-shaped crack which consists of a delamination crack intersecting a crack which is perpendicular to the interface, and a delamination crack initiating from a stress-free boundary of the bonded layers are some of the practical crack geometries considered.

Lu, M. C.↗

The heating of a thermally conducting stratified medium. I - Self-gravitating gas threaded by magnetic fields

Two theoretical problems are presented to illustrate the equilibrium structure of a gravitating gas threaded by magnetic fields. The gas is heated everywhere at a rate proportional to the local density and cooled by thermal conduction channeled along the magnetic field. Emphasis is placed on the influence of anisotropic thermal conduction on the distribution of the gas. In the first problem, the magnetic field has straight field lines that are nearly parallel. The simple field line geometry allows us to calculate, in closed form, the equilibria configurations of a gas slab which is cooled nonuniformly by thermal conduction along field lines which leave the gas slab at varying angles. The problem involves the gas having a free boundary, which is one of the unknowns to be determined. A discussion is given of the set of all possible equilibria so constructed. In the second problem there is no electric current flowing in the gas but there is a potential magnetic field having curved field lines. Exact solutions of equilibriums are presented to illustrate the heating of stellar and galactic coronas.

Lerche, I.↗

Profiles of electrified drops and bubbles

Axisymmetric equilibrium shapes of conducting drops and bubbles, (1) pendant or sessile on one face of a circular parallel-plate capacitor or (2) free and surface-charged, are found by solving simultaneously the free boundary problem consisting of the augmented Young-Laplace equation for surface shape and the Laplace equation for electrostatic field, given the surface potential. The problem is nonlinear and the method is a finite element algorithm employing Newton iteration, a modified frontal solver, and triangular as well as quadrilateral tessellations of the domain exterior to the drop in order to facilitate refined analysis of sharply curved drop tips seen in experiments. The stability limit predicted by this computer-aided theoretical analysis agrees well with experiments.

Basaran, O. A.↗

Eruptive solar magnetic fields

The quasi-steady evolution of solar magnetic fields in response to gradual photospheric changes is considered, with particular attention given to the threshold of a sudden eruption in the solar atmosphere. The formal model of an evolving, force-free field dependent on two Cartesian coordinates is extended to a field which is not force free but in static equilibrium with plasma pressure and gravity. The basic physics is illustrated through the evolution of a loop-shaped electric current sheet enclosing a potential bipolar field with footpoints rooted in the photosphere. A free-boundary problem is posed and then solved for the equilibrium configuration of the current sheet in a hydrostatically supported isothermal atmosphere. As the footpoints move apart to spread a constant photospheric magnetic flux over a larger region, the equilibria available extend the field to increasing heights.

Low, B. C.↗

Experimental and theoretical studies of the laws governing condensate deposition from combustion gases

A description is presented of the results of a research program directed at an improved understanding of condensate deposition rate phenomena in combustion systems. The conducted experiments make use of real-time optical laser reflectance-interference-polarization techniques in flame environments. The obtained new data and the results of previous gravimetric experiments are employed as a basis for the development of a comprehensive convective diffusion deposition theory, taking into account the assumption of a multicomponent vapor or multisize class particles 'source-free' boundary layer. The theory makes it possible to provide self-consistent salt/ash/soot deposition rate predictions over a wide variety of environmental conditions.

Rosner, D. E.↗