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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

TokaMaker: An open-source time-dependent Grad-Shafranov tool for the design and modeling of axisymmetric fusion devices

Here, in this paper, we present a new static and time-dependent MagnetoHydroDynamic (MHD) equilibrium code, TokaMaker, for axisymmetric configurations of magnetized plasmas, based on the well-known Grad-Shafranov equation. This code utilizes finite element methods on an unstructured triangular grid to enable capturing accurate machine geometry and simple mesh generation from engineering-like descriptions of present and future devices. The new code is designed for ease of use without sacrificing capability and speed through a combination of Python, Fortran, and C/C++ components. A detailed description of the numerical methods of the code, including a novel formulation of the boundary conditions for free-boundary equilibria, and validation of the implementation of those methods using both analytic test cases and cross-code validation is shown. Results show expected convergence across tested polynomial degree for analytic and cross-code test cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

Tension/torsion loading of composite laminates with free-edge boundary conditions

A finite element analysis of a family of quasi-isotropic graphite/epoxy laminates was conducted for tension, torsion, and combined tension/torsion loading. The model was employed to investigate the effects of delaminations on torsional stiffness and also to evaluate the total strain energy release rates associated with these delaminations.

Hooper, Steven J.↗
DOE OSTI2025

Physicochemical and Molecular Insights into the Boundary Layer and Free Troposphere Aerosol Interactions over the Southern Great Plains

Ambient aerosols’ vertical profiles are critical for evaluating the role of aerosols in atmospheric chemistry and radiative transfer, but limited data on these profiles hinders our ability to fully assess their impact on the Earth's radiative balance. Here, in this study, we investigated the size-, time-, and altitude resolved composition of individual particles and bulk molecular composition of particle samples collected by an uncrewed aerial system–ArcticShark over the Southern Great Plains. Single particle microanalysis shows that, the free tropospheric (FT) samples are dominated (56-66%) by carbonaceous sulfate particles, while boundary layer (BL) samples are dominated (57-74%) by carbonaceous particles. Back trajectory simulations suggest that FT particles are likely influenced by long-range transport and have undergone aqueous-phase processing. Conversely, in-situ size distribution data shows evidence of particle growth in the upper BL and just below the FT. This observation may indicate vertical transport of particles from an elevated aerosol layer in the FT, possibly linked to a new particle formation event. This observation is further supported by high resolution molecular composition data, which reveals particle volatility increasing with increasing size, which aligns with the growth event. This study aids in fundamental understanding of the compositional and molecular specificity of vertically resolved organic aerosols to provide insights into particle size evolution for future atmospheric models.

ArcticShark↗
NASA NTRS1989

Finite element analysis in fluids; Proceedings of the Seventh International Conference on Finite Element Methods in Flow Problems, University of Alabama, Huntsville, Apr. 3-7, 1989

Recent advances in computational fluid dynamics are examined in reviews and reports, with an emphasis on finite-element methods. Sections are devoted to adaptive meshes, atmospheric dynamics, combustion, compressible flows, control-volume finite elements, crystal growth, domain decomposition, EM-field problems, FDM/FEM, and fluid-structure interactions. Consideration is given to free-boundary problems with heat transfer, free surface flow, geophysical flow problems, heat and mass transfer, high-speed flow, incompressible flow, inverse design methods, MHD problems, the mathematics of finite elements, and mesh generation. Also discussed are mixed finite elements, multigrid methods, non-Newtonian fluids, numerical dissipation, parallel vector processing, reservoir simulation, seepage, shallow-water problems, spectral methods, supercomputer architectures, three-dimensional problems, and turbulent flows.

Chung, T. J.↗
NASA NTRS1996

Calculation of Turbulent Subsonic Diffuser Flows Using the NPARC Navier-Stokes Code

Axisymmetric subsonic diffuser flows were calculated with the NPARC Navier-Stokes code in order to determine the effects various code features have on the flow solutions. The code features examined in this work were turbulence models and boundary conditions. Four turbulence models available in NPARC were used: the Baldwin-Lomax algebraic model, the Baldwin-Barth one-equation model, and the Chien kappa-epsilon and Wilcox kappa-omega two-equation models. The three boundary conditions examined were the free boundary, the mass flux boundary and the subsonic outflow with variable static pressure. In addition to boundary condition type, the geometry downstream of the diffuser was varied to see if upstream influences were present. The NPARC results are compared with experimental data and recommendations are given for using NPARC to compute similar flows.

Dudek, J. C.↗
NASA NTRS1982

Exchange of sulfur between the free troposphere, marine boundary layer, and the sea surface

A steady state analysis has been made of the sources and sinks of sulfur in the remote marine atmosphere and of the steady state exhanges of SO2 and sulfate between the sea, boundary layer, and free troposphere. Two situations were examined: an area of high biological productivity (Gulf of Guinea) and an area of low productivity (tropical Atlantic). The results of the analysis suggest that about one-third of the background concentrations of SO2 and excess sulfate observed in the remote marine atmosphere mixes in from the free troposphere, with the balance, about 100 micrograms S/sq m per day, attributed to the local oxidation of reduced sulfur-bearing precursors in the boundary layer.

Kritz, M. A.↗
NASA NTRS1983

A series expansion of the acoustic power radiated from planar sources

A series expansion in ascending powers of the wavenumber k is derived for the acoustic power delivered by baffled or unbaffled planar sources. This series provides a relatively simple means of derving expressions for the power radiated by a baffled source with a known velocity distribution and can be used for unbaffled plates when the velocity field outside the plate is also known. The terms in the series are calculated from the moments of this velocity distribution in the plane containing the source. If these moments are written as derivaties in wavenumber space, it is shown that a MacLaurin expansion of the Fourier transformed velocity provides an easy technique for computing the first few terms of the acoustic power. Examples are provided for baffled, rectangular plates with various boundary conditions. The arbirarily shaped plate with free boundaries is particularly interesting. It is proven that the volume flow across it surface must be zero and as a result corner and edge mode radiation cannot exist for this kind of source.

Willams, E. G.↗
NASA NTRS1993

Laplace's equation and the Dirichlet-Neumann map in multiply connected domains

A variety of problems in material science and fluid dynamics require the solution of Laplace's equation in multiply connected domains. Integral equation methods are natural candidates for such problems, since they discretize the boundary alone, require no special effort for free boundaries, and achieve superalgebraic convergence rates on sufficiently smooth domains in two space dimensions, regardless of shape. Current integral equation methods for the Dirichlet problem, however, require the solution of M independent problems of dimension N, where M is the number of boundary components and N is the total number of points in the discretization. In this paper, we present a new boundary integral equation approach, valid for both interior and exterior problems, which requires the solution of a single linear system of dimension N + M. We solve this system by making use of an iterative method (GMRES) combined with the last multipole method for the rapid calculation of the necessary matrix vector products. For a two-dimensional system with 200 components and 100 points on each boundary, we gain a speedup of a factor of 100 from the new analytic formulation and a factor of 50 from the fast multipole method. The resulting scheme brings large scale calculations in extremely complex domains within practical reach.

Greenbaum, A.↗
NASA NTRS1985

On the penetration of a hot diapir through a strongly temperature-dependent viscosity medium

The ascent of a hot spherical body through a fluid with a strongly temperature-dependent viscosity has been studied using an axisymmetric finite element method. Numerical solutions range over Peclet numbers of 0.1 - 1000 from constant viscosity up to viscosity variations of 100,000. Both rigid and stress-free boundary conditions were applied at the surface of the sphere. The dependence of drag on viscosity variation was shown to have no dependence on the stress boundary condition except for a Stokes flow scaling factor. A Nusselt number parameterization based on the stress-free constant viscosity functional dependence on the Peclet number scaled by a parameter depending on the viscosity structure fits both stress-free and rigid boundary condition data above viscosity variations of 100. The temperature scale height was determined as a function of sphere radius. For the simple physical model studied in this paper pre-heating is required to reduce the ambient viscosity of the country rock to less than 10 to the 22nd sq cm/s in order for a 10 km diapir to penetrate a distance of several radii.

Daly, S. F.↗
DOE OSTI2025

Feedforward equilibrium trajectory optimization with GSPulse

One of the common tasks required for designing new plasma scenarios or evaluating capabilities of a tokamak is to design the desired equilibria using a Grad-Shafranov (GS) equilibrium solver. However, most standard equilibrium solvers are time-independent and do not include dynamic effects such as plasma current flux consumption, induced vessel currents, or voltage constraints. Another class of tools, plasma equilibrium evolution simulators, do include time-dependent effects. These are generally structured to solve the forward problem of evolving the plasma equilibrium given feedback-controlled voltages. In this work, we introduce GSPulse, a novel algorithm for equilibrium trajectory optimization, that is more akin to a pulse planner than a pulse simulator. GSPulse includes time-dependent effects and solves the inverse problem: given a user-specified set of target equilibrium shapes, as well as limits on the coil currents and voltages, the optimizer returns trajectories of the voltages, currents, and achievable equilibria. This task is useful for scoping performance of a tokamak and exploring the space of achievable pulses. The computed equilibria satisfy both Grad-Shafranov force balance and axisymmetric circuit dynamics. The optimization is performed by restructuring the free-boundary equilibrium evolution equations into a form where it is computationally efficient to optimize the entire dynamic sequence. GSPulse can solve for hundreds of equilibria simultaneously within a few minutes. GSPulse has been validated against NSTX-U and MAST-U experiments and against SPARC feedback control simulations, and is being used to perform scenario design for SPARC. The computed trajectories can be used as feedforward inputs that are connected to the feedback controller to inform and improve feedback performance. The code for GSPulse is available open-source at github.com/jwai-cfs/GSPulse_public.

equilibrium↗
NASA NTRS1987

The fully nonlinear development of Goertler vortices in growing boundary layers

The fully nonlinear development of small wavelength Goertler vortices in a growing boundary layer is investigated using a combination of asymptotic and numerical methods. The starting point for the analysis is the weakly nonlinear theory of Hall (1982b) who discussed the initial development of small amplitude vortices in a neighborhood of the location where they first become linearly unstable. That development is unusual in the context of nonlinear stability theory in that it is not described by the Stuart-Watson approach. In fact the development is governed by a pair of coupled nonlinear partial differential evolution equations for the vortex flow and the mean flow correction. Here the further development of this interaction is considered for vortices so large that the mean flow correction driven by them is as large as the basic state. Surprisingly it is found that such a nonlinear interaction can still be described by asymptotic means. It is shown that the vortices spread out across the boundary layer and effectively drive the boundary layer. In fact the system obtained by writing down the equations for the fundamental component of the vortex generate a differential equation for the basic state. Thus the mean flow adjusts so as to make these large amplitude vortices locally neutral. Moreover in the region where the vortices exist the mean flow has a square-root profile and the vortex velocity field can be written down in closed form. The upper and lower boundaries of the region of vortex activity are determined by a free-boundary problem involving the boundary layer equations. In general it is found that this region ultimately includes almost all of the original boundary layer and much of the free-stream. In this situation the mean flow has essentially no relationship to the flow which exists in the absence of the vortices.

Hall, Philip↗
NASA NTRS1983

On the maintenance of the Marine boundary layer-free tropospheric gradient of transient tracers

A simple model of the marine boundary layer is used to arrive at a steady-state budget for trace gases in this layer. Five processes maintain the tracer balance in the layer. The first is a large-scale horizontal convergence of the tracer; the second is the pumping by clouds of the tracer out of the boundary layer into the free troposphere; the third is penetration by the tracer from the free troposphere into the marine boundary from compensating subsidence in the cloud-free air surrounding cumulus clouds; the fourth is sea surface emission or deposition of the tracer; and the fifth is in situ photochemical production or loss of the tracer within the marine boundary layer. The budget equation is used in deriving an expression for the mixing ratio difference between the two regions. A scale analysis suggests that the mixing ratio difference between the marine boundary layer and the free troposphere can be an order of magnitude smaller in disturbed weather conditions than in suppressed conditions, showing that the difference in chemical concentrations between the free troposphere and the marine boundary layer cannot be rationalized solely on the basis of chemical measurements.

Gidel, L. T.↗
DOE OSTI2021

A weighted Shifted Boundary Method for free surface flow problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods and was recently introduced for the Poisson, linear advection/diffusion, Stokes, Navier-Stokes, acoustics, and shallow-water equations. By reformulating the original boundary value problem over a surrogate (approximate) computational domain, the SBM avoids integration over cut cells and the associated problematic issues regarding numerical stability and matrix conditioning. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we extend the SBM to the simulation of incompressible Navier-Stokes flows with moving free-surfaces, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach prevents spurious pressure oscillations in time, which would otherwise be produced if the total active fluid volume were to change abruptly over a time step. In fact, the proposed weighted SBM method induces small mass (i.e., volume) conservation errors, which converge quadratically in the case of piecewise-linear finite element interpolations, as the grid is refined. Finally, we present an extensive set of two- and three-dimensional tests to demonstrate the robustness and accuracy of the method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1985

Optimization studies for the development of adaptive wall wind tunnel

The development of adaptive wall wind tunnels makes it possible to establish interference-free boundaries in the region of the model. It is, therefore, feasible to simulate free air conditions. However, the amount of total testing time involved in establishing the required conditions and the complexity of the needed instrumentation prevent currently the use of such methods for routine production testing. Optimization studies have been proposed with the objective to find approaches which will make a more extended application of adaptive wall techniques possible. A new 25cm x 11cm adaptive wall test section was designed and installed to carry out some of these studies. Investigations envisaged are related to the development of a one step convergence method to reduce total testing time, and to measurement procedures which will permit a simplification of the instrumentation. Attention is given to the proposed concept, the test section, instrumentation, and the influence coefficients.

Bodapati, S.↗
NASA NTRS1989

The structure of sheared turbulence near a plane boundary

The effects of shear on the structure of turbulence near a plane boundary are examined. Inhomogeneous rapid distortion theory is applied to a uniform-shear boundary layer (USBL) with a slip velocity condition at the surface. The effects of blocking the vertical component of the turbulence by the boundary in a USBL are compared with the effects in a shear-free boundary layer. It is shown that the blocking is only slightly modified by the presence of shear. The shear does alter the way in which the boundary influences the horizontal components and the turbulent shear stress near the horizontal component.

Lee, Moon Joo↗