Search NASA⌕ Search

SEARCH · Search NASA

Results for “Weak form”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

A Meshless Local Petrov-Galerkin Method for Euler-Bernoulli Beam Problems

An accurate and yet simple Meshless Local Petrov-Galerkin (MLPG) formulation for analyzing beam problems is presented. In the formulation, simple weight functions are chosen as test functions. The use of these functions shows that the weak form can be integrated with conventional Gaussian integration. The MLPG method was evaluated by applying the formulation to a variety of patch test and thin beam problems. The formulation successfully reproduced exact solutions to machine accuracy when test functions with C2 continuity and an appropriate order of basis functions are used.

Raju, I. S.↗
NASA NTRS2016

Simple Test Functions in Meshless Local Petrov-Galerkin Methods

Two meshless local Petrov-Galerkin (MLPG) methods based on two different trial functions but that use a simple linear test function were developed for beam and column problems. These methods used generalized moving least squares (GMLS) and radial basis (RB) interpolation functions as trial functions. These two methods were tested on various patch test problems. Both methods passed the patch tests successfully. Then the methods were applied to various beam vibration problems and problems involving Euler and Beck's columns. Both methods yielded accurate solutions for all problems studied. The simple linear test function offers considerable savings in computing efforts as the domain integrals involved in the weak form are avoided. The two methods based on this simple linear test function method produced accurate results for frequencies and buckling loads. Of the two methods studied, the method with radial basis trial functions is very attractive as the method is simple, accurate, and robust.

Raju, Ivatury S.↗
NASA NTRS2011

Electromagnetic Scattering Analysis of Large Size Asteroids/Comets for Reflection/Transmission Tomography (RTT)

A precise knowledge of the interior structure of asteroids, comets, and Near Earth Objects (NEO) is important to assess the consequences of their impacts with the Earth and develop efficient mitigation strategies. Knowledge of their interior structure also provides opportunities for extraction of raw materials for future space activities. Low frequency radio sounding is often proposed for investigating interior structures of asteroids and NEOs. For designing and optimizing radio sounding instrument it is advantageous to have an accurate and efficient numerical simulation model of radio reflection and transmission through large size bodies of asteroid shapes. In this presentation we will present electromagnetic (EM) scattering analysis of electrically large size asteroids using (1) a weak form formulation and (2) also a more accurate hybrid finite element method/method of moments (FEM/MOM) to help estimate their internal structures. Assuming the internal structure with known electrical properties of a sample asteroid, we first develop its forward EM scattering model. From the knowledge of EM scattering as a function of frequency and look angle we will then present the inverse scattering procedure to extract its interior structure image. Validity of the inverse scattering procedure will be presented through few simulation examples.

Deshpande, Manohar↗
NASA NTRS2019

Reliable Bonding of Composite Laminates Using Reflowable Epoxy Resins

Epoxy matrix composites assembled with adhesives maximize the performance of aerospace structures, but the possibility of forming weak bonds requires the installation of redundant fasteners, which add weight and manufacturing cost. Co-cured joints (e.g. unitized composite structures) are immune to weak bonds because the uncured resin undergoes diffusion and mixing through the joint. A means of co-curing complex structures may reduce the need for redundant fasteners in bondlines. To this end, NASA started the AERoBOND project to develop novel joining materials to enable a “secondary-co-cure” assembly process. Aerospace epoxy resin systems reformulated with offset stoichiometry prevented the resin from advancing beyond the gel point during a conventional autoclave cure cycle up to 180 °C. The offset resins were applied to the joining surfaces of laminate preforms as prepreg. Two surfaces with complimentary offset resins were joined using conventional secondary bonding techniques. Preliminary efforts have indicated that the resulting joint has no discernable interface and appears as a conventional co-cured laminate under optical magnification. This report will discuss the initial work performed regarding formulation of the epoxy resin system using calorimetry, rheology, and mechanical testing.

Palmieri, Frank L.↗
NASA NTRS1983

Applications of the contravariant form of the Navier-Stokes equations

The contravariant Navier-Stokes equations in weak conservation form are well suited to certain fluid flow analysis problems. Three dimensional contravariant momentum equations may be used to obtain Navier-Stokes equations in weak conservation form on a nonplanar two dimensional surface with varying streamsheet thickness. Thus a three dimensional flow can be simulated with two dimensional equations to obtain a quasi-three dimensional solution for viscous flow. When the Navier-Stokes equations on the two dimensional nonplanar surface are transformed to a generalized body fitted mesh coordinate system, the resulting equations are similar to the equations for a body fitted mesh coordinate system on the Euclidean plane. Contravariant momentum components are also useful for analyzing compressible, three dimensional viscous flow through an internal duct by parabolic marching. This type of flow is efficiently analyzed by parabolic marching methods, where the streamwise momentum equation is uncoupled from the two crossflow momentum equations. This can be done, even for ducts with a large amount of turning, if the Navier-Stokes equations are written with contravariant components.

Katsanis, T.↗
NASA NTRS2012

The Mechanisms for the Onset and Explosive Eruption of Coronal Mass Ejections and Eruptive Flares

We have investigated the onset and acceleration of coronal mass ejections (CMEs) and eruptive flares. To isolate the eruption physics, our study uses the breakout model, which is insensitive to the energy buildup process leading to the eruption. We performed 2.5D simulations with adaptive mesh refinement that achieved the highest overall spatial resolution to date in a CME/eruptive flare simulation. The ultra-high resolution allows us to separate clearly the timing of the various phases of the eruption. Using new computational tools, we have determined the number and evolution of all X- and O-type nulls in the system, thereby tracking both the progress and the products of reconnection throughout the computational domain. Our results show definitively that CME onset is due to the start of fast reconnection at the breakout current sheet. Once this reconnection begins, eruption is inevitable; if this is the only reconnection in the system, however, the eruption will be slow. The explosive CME acceleration is triggered by fast reconnection at the flare current sheet. Our results indicate that the explosive eruption is caused by a resistive instability, not an ideal process. Moreover, both breakout and flare reconnections begin first as a form of weak tearing characterized by a slowly evolving plasmoids, but eventually transition to a fast form with well-defined Alfvenic reconnection jets and rapid flux transfer. This transition to fast reconnection is required for both CME onset and explosive acceleration. We discuss the key implications of our results for CME/flare observations and for theories of magnetic reconnection.

Karpen, Judith T.↗
NASA NTRS1987

Is there a weak mixed polarity background field? Theoretical arguments

A number of processes associated with the formation of active regions produce 'U-loops': fluxtubes having two ends at the photosphere but otherwise still embedded in the convection zone. The mass trapped on the field lines of such loops makes them behave in a qualitatively different way from the 'omega-loops' that form active regions. It is shown that U-loops will disperse though the convection zone and form a weak (down to a few gauss) field that covers a significant fraction of the solar surface. This field is tentatively identified with the inner-network fields observed at Kitt Peak and Big Bear. The process by which these fields escape through the surface is described; a remarkable property is that it can make active region fields apparently disappear in situ. The mixed polarity moving magnetic features near sunspots are interpreted as a locally intense form of this disappearance by escape of U-loops.

Spruit, H. C.↗
NASA NTRS1990

Finite element solution of optimal control problems with inequality constraints

A finite-element method based on a weak Hamiltonian form of the necessary conditions is summarized for optimal control problems. Very crude shape functions (so simple that element numerical quadrature is not necessary) can be used to develop an efficient procedure for obtaining candidate solutions (i.e., those which satisfy all the necessary conditions) even for highly nonlinear problems. An extension of the formulation allowing for discontinuities in the states and derivatives of the states is given. A theory that includes control inequality constraints is fully developed. An advanced launch vehicle (ALV) model is presented. The model involves staging and control constraints, thus demonstrating the full power of the weak formulation to date. Numerical results are presented along with total elapsed computer time required to obtain the results. The speed and accuracy in obtaining the results make this method a strong candidate for a real-time guidance algorithm.

Bless, Robert R.↗
NASA NTRS1983

A viscous shock-layer flowfield analysis by an explicit-implicit method

An implicit analogue of a widely used explicit method to external axisymmetric laminar flows with strong entropy gradients is extended. The details of the "numerics" of the implicit part are provided in a body oriented coordinate system with a moving outer (shock) boundary during the transient part of the solutions. The limiting values of the Courant number are obtained when the shock boundary is treated explicitly. The solution algorithm outlined includes the treatment of the source term associated with the equations in weak conservation form. From the results obtained for two sample problems, it becomes clear that accuracy of predictions is, indeed, very good at higher values of the Courant number. There is a significant saving in overall computing time, depending on the Courant number used and the flow Reynolds number. These properties combined with the simplicity of programming the implicit analog may appeal to researchers for using it in the analysis of 3-D flow problems.

Gupta, R. M.↗
NASA NTRS1983

Viscous shock-layer flowfield analysis by an explicit-implicit method

The present work extends the recently reported implicit analogue of MacCormack's earlier widely-used explicit method to external axisymmetric laminar flows with strong entropy gradients. The details of the 'numerics' of the implicit part are provided in a body-oriented coordinate system with a moving outer (shock) boundary during the transient part of the solutions. The limiting values of the Courant number are obtained when the shock boundary is treated explicitly. The solution algorithm outlined includes the treatment of the source term associated with the equations in weak conservation form. From the results obtained for two sample problems, it becomes clear that accuracy of predictions is, indeed, very good at higher values of the Courant number. There is a significant saving in overall computing time, depending on the Courant number used and the flow Reynolds number. These properties combined with the simplicity of programming the implicit analogue may appeal to researchers for using it in the analysis of three-dimensional flow problems.

Gupta, R. N.↗
NASA NTRS1987

Particle simulation of auroral double layers

Work on the simulation of auroral double layers (DLs) with realistic particle-in-cell models is presented. An early model simulated weak DLs formed in a self-consistent circuit but under conditions subject to the ion-acoustic instability. Recent work has focused on strong DLs formed when currentless jets are injected into a dipole magnetic field.

Smith, Bruce L.↗
NASA NTRS1989

Active vibration mitigation of distributed parameter, smart-type structures using Pseudo-Feedback Optimal Control (PFOC)

A new, near-optimal feedback control technique is introduced that is shown to provide excellent vibration attenuation for those distributed parameter systems that are often encountered in the areas of aeroservoelasticity and large space systems. The technique relies on a novel solution methodology for the classical optimal control problem. Specifically, the quadratic regulator control problem for a flexible vibrating structure is first cast in a weak functional form that admits an approximate solution. The necessary conditions (first-order) are then solved via a time finite-element method. The procedure produces a low dimensional, algebraic parameterization of the optimal control problem that provides a rigorous basis for a discrete controller with a first-order like hold output. Simulation has shown that the algorithm can successfully control a wide variety of plant forms including multi-input/multi-output systems and systems exhibiting significant nonlinearities. In order to firmly establish the efficacy of the algorithm, a laboratory control experiment was implemented to provide planar (bending) vibration attenuation of a highly flexible beam (with a first clamped-free mode of approximately 0.5 Hz).

Patten, W. N.↗
NASA NTRS1994

Gravity survey of the Mount Toondina impact structure, South Australia

Gravity and seismic reflection data, together with geologic mapping, indicate that the Mount Toondina feature in South Australia is best interpreted as an eroded 4-km-diameter impact structure consisting of a ring structural depression surrounding a pronounced central uplift. Beds at the center of the structure within the central uplift have been raised as much as 200 m from depth and deformed by convergent flow. Seismic reflection data indicate that deformation extends to depths of only approximately = 800 m; at greater depths the reflectors are nearly flat lying, indicating little or no deformation. Gravity data show residual anomalies of +1.0 mGal coincident with the central uplift and a -0.5 Mgal low associated with the ring structural depression. Modeling of the gravity data indicates that relatively high-density material occurs within the central uplift, whereas the ring depression is filled with low-density material. The deformation at Mount Toondina is typical of a complex impact crater; the 4-km diameter is consistent with the expected threshold size for complex craters formed in weak to moderate strength sedimentary rocks.

Plescia, J. B.↗
NASA NTRS1994

Gravity Survey of the Mt. Toondina Impact Structure South Australia

Gravity and seismic reflection data, together with geologic mapping, indicate that the Mt. Toondina feature in South Australia is best interpreted as an eroded 4-km-diameter impact structure consisting of a ring structural depression surrounding a pronounced central uplift. Beds at the center of the structure within the central uplift have been raised as much as 200 m from depth and deformed by convergent flow. Seismic reflection data indicate that deformation extends to depths of only about 800 m; at greater depths the reflectors are nearly flat lying indicating little or no deformation. Gravity data show residual anomalies of +1.0 mGal coincident with the central uplift and a -0.5 mGal low associated with the ring structural depression. Modeling of the gravity data indicates that relatively high-density material occurs within the central uplift, whereas the ring depression is filled with low-density material. The deformation at Mt. Toondina is typical of a complex impact crater; the 4-km diameter is consistent with the expected threshold size for complex craters formed in weak to moderate strength sedimentary rocks.

Mt. Toondina Australia↗
NASA NTRS2004

Solution-Assisted Optical Contacting

A modified version of a conventional optical-contact procedure has been found to facilitate alignment of optical components. The optical-contact procedure (called simply optical contacting in the art) is a standard means of bonding two highly polished and cleaned glass optical components without using epoxies or other adhesives. In its unmodified form, the procedure does not involve the use of any foreign substances at all: components to be optically contacted are dry. The main disadvantage of conventional optical contacting is that it is difficult or impossible to adjust the alignment of the components once they have become bonded. In the modified version of the procedure, a drop of an alcohol-based optical cleaning solution (isopropyl alcohol or similar) is placed at the interface between two components immediately before putting the components together. The solution forms a weak bond that gradually strengthens during a time interval of the order of tens of seconds as the alcohol evaporates. While the solution is present, the components can be slid, without loss of contact, to perform fine adjustments of their relative positions. After about a minute, most of the alcohol has evaporated and the optical components are rigidly attached to each other. If necessary, more solution can be added to enable resumption or repetition of the adjustment until the components are aligned to the required precision.

Shaddock, Daniel↗
NASA NTRS2019

Detection of Nonthermal Hard X-Ray Emission from the “Fermi Bubble” in an External Galaxy

We report new Chandra hard X-ray (>2 keV) and Jansky Very Large Array C-band observations of the nuclear super bubble of NGC 3079, an analog of the “Fermi bubble” in our Milky Way. We detect extended hard X-ray emission on the southwest (SW)side of the galactic nucleus with coherent multiwavelength features in radio, Hα, and soft X-ray. The hard X-ray feature has a cone shape with possibly a weak cap, forming a bubble-like structure with a diameter of∼1.1 kpc. A similar extended feature, however, is not detected on the northeast (NE)side, which is brighter in all other wavelengths such as radio, Hα, and soft X-ray. Scattered photons from the nuclear region or other nearby point-like X-ray bright sources, inverse Compton emission from cosmic-ray(CR)electrons via interaction with the cosmic microwave background, or any individually faint stellar X-ray source populations, cannot explain the extended hard X-ray emission on the SW side and the strongly NE/SW asymmetry. A synchrotron emission model, plus a thermal component accounting for the excess at∼1 keV, can well characterize the broadband radio/hard X-ray spectra. The broadband synchrotron spectra do not show any significant cutoff, and even possibly slightly flatten at higher energy. This rule out a loss-limited scenario in the acceleration of the CR electrons in or around this superbubble. As the first detection of kiloparsec-scale extended hard X-ray emission associated with a galactic nuclear superbubble, the spatial and spectral properties of the multiwavelength emissions indicate that the CR leptons responsible for the broadband synchrotron emission from the SW bubble must be accelerated in situ, instead of transported from the nuclear region of the galaxy.

Jiang-Tao Li↗
NASA NTRS1989

A weak Hamiltonian finite element method for optimal control problems

A temporal finite element method based on a mixed form of the Hamiltonian weak principle is developed for dynamics and optimal control problems. The mixed form of Hamilton's weak principle contains both displacements and momenta as primary variables that are expanded in terms of nodal values and simple polynomial shape functions. Unlike other forms of Hamilton's principle, however, time derivatives of the momenta and displacements do not appear therein; instead, only the virtual momenta and virtual displacements are differentiated with respect to time. Based on the duality that is observed to exist between the mixed form of Hamilton's weak principle and variational principles governing classical optimal control problems, a temporal finite element formulation of the latter can be developed in a rather straightforward manner. Several well-known problems in dynamics and optimal control are illustrated. The example dynamics problem involves a time-marching problem. As optimal control examples, elementary trajectory optimization problems are treated.

Hodges, Dewey H.↗