Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fictitious forces”

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.

36 records · Page 2

Free vibration analysis by BEM using particular integrals

A new method for the free-vibration analysis using the boundary element technique is presented. The method utilizes a fictitious vector function to approximate the inertia forces and then uses the well-known concept of complementary functions and particular integrals to solve the resulting governing differential equations. The necessary particular integrals are defined for the two and three-dimensional analyses, and the present formulation is applied to a number of two-dimensional problems to show its accuracy and efficiency in the solution of realistic engineering problems.

Ahmad, S.↗

Complexity of many-body interactions in transition metals via machine-learned force fields from the TM23 data set

Abstract This work examines challenges associated with the accuracy of machine-learned force fields (MLFFs) for bulk solid and liquid phases ofd-block elements. In exhaustive detail, we contrast the performance of force, energy, and stress predictions across the transition metals for two leading MLFF models: a kernel-based atomic cluster expansion method implemented using sparse Gaussian processes (FLARE), and an equivariant message-passing neural network (NequIP). Early transition metals present higher relative errors and are more difficult to learn relative to late platinum- and coinage-group elements, and this trend persists across model architectures. Trends in complexity of interatomic interactions for different metals are revealed via comparison of the performance of representations with different many-body order and angular resolution. Using arguments based on perturbation theory on the occupied and unoccupieddstates near the Fermi level, we determine that the large, sharpddensity of states both above and below the Fermi level in early transition metals leads to a more complex, harder-to-learn potential energy surface for these metals. Increasing the fictitious electronic temperature (smearing) modifies the angular sensitivity of forces and makes the early transition metal forces easier to learn. This work illustrates challenges in capturing intricate properties of metallic bonding with current leading MLFFs and provides a reference data set for transition metals, aimed at benchmarking the accuracy and improving the development of emerging machine-learned approximations.

Chemistry↗

Ephemeral penalty functions for contact-impact dynamics

The use of penalty functions to treat a class of structural contact-impact problems is investigated, with emphasis on ones in which the impact phenomena are primarily nondestructive in nature and in which only the gross characterization of the response is required. The dynamic equations of motion are integrated by the difference method. The penalty is represented as an ephemeral fictitious nonlinear spring that is inserted on anticipation of contact. The magnitude and variation of the penalty force is determined through energy balancing considerations. The 'bell shape' of the penalty force function for positive gap was found to be satisfactory, as it depends on only two parameters that can be directly assigned the physical meaning of force and distance. The determination of force law parameters by energy balance worked well. The incorporation of restitution coefficients by the area balancing method yielded excellent results, and no substantial modifications are anticipated. Extensional penalty springs are obviously sufficient for the simple examples treated.

De La Fuente, Horacio M.↗

Towards sharp error analysis of extended Lagrangian molecular dynamics

The extended Lagrangian molecular dynamics (XLMD) method provides a useful framework for reducing the computational cost of a class of molecular dynamics simulations with constrained latent variables. The XLMD method relaxes the constraints by introducing a fictitious mass ε for the latent variables and solving a set of singularly perturbed ordinary differential equations. While favorable numerical performance of XLMD has been demonstrated in several different contexts in the past decade, mathematical analysis of the method remains scarce. Here, we propose the first error analysis of the XLMD method in the context of a classical polarizable force field model. While the dynamics with respect to the atomic degrees of freedom are general and nonlinear, the key mathematical simplification of the polarizable force field model is that the constraints on the latent variables are given by a linear system of equations. We prove that when the initial value of the latent variables is compatible in a sense that we define, XLMD converges as the fictitious mass ε is made small with $\mathscr{O}$(ε) error for the atomic degrees of freedom and with $\mathscr{O}$($\sqrt{ε}$) error for the latent variables, when the dimension of the latent variable d' is 1. Furthermore, when the initial value of the latent variables is improved to be optimally compatible in a certain sense, we prove that the convergence rate can be improved to $\mathscr{O}$(ε) for the latent variables as well. Numerical results verify that both estimates are sharp not only for d'=1, but also for arbitrary d'. In the setting of general d', we do obtain convergence, but with the non-sharp rate of $\mathscr{O}$($\sqrt{ε}$) for both the atomic and latent variables.

74 ATOMIC AND MOLECULAR PHYSICS↗

Artificial dynamical effects in quantum field theory

In Newtonian mechanics, studying a system in a non-Galilean reference frame can lead to inertial pseudoforces appearing, such as the centrifugal force that seems to arise in dynamics analysed in a rotating frame. Likewise, artificial effects may arise in relativistic quantum field theory (QFT) if a system is studied in a framework that violates Poincaré invariance. Here, we highlight how such issues complicate the traditional canonical quantization of QFTs and can lead to a subjective description of natural phenomena. By contrast, the treatment of the same problem using light-front quantization is free from spurious pseudoeffects because Poincaré invariance is effectively preserved for all practical intents and purposes. We illustrate these statements using several examples: the Gerasimov-Drell-Hearn (GDH) relation, a fundamental feature of QFT; the absence of any measurable impact of Lorentz contraction in high-energy collisions; and the fictitious character of vacuum fluctuation contributions to the cosmological constant.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dynamics and Control of a Tethered Enhanced Gravity Tractor Performing Asteroid Deflection

The dynamics and control of an Enhanced Gravity Tractor (EGT) augmented with a tether for deflecting an asteroid are studied. A conventional EGT consists of collected asteroidal mass collocated with the spacecraft. Because of the presence of a tether, the collected mass is placed where the EGT would have been without a tether, and the spacecraft is placed farther away from the asteroid. Doing so improves the fuel efficiency and safety margin of the EGT operation without significantly sacrificing the gravitational attraction between the asteroid and the EGT. The tether is modeled as a series of particles connected by spring-dashpot systems. Physical properties of the tether are selected to be similar to those of the SPECTRA-1000, Kevlar-29, and Kevlar-49 fibers. It is assumed that control is applied only to the spacecraft, and there is no active control associated with the collected mass. A Proportional-Derivative (PD) controller is employed to maintain the spacecraft and the collected mass at desired positions relative to the asteroid. Numerical simulations of tethered EGT operations at 2008 EV5, Itokawa, Apophis, and a fictitious ellipsoidal asteroid are performed. It is demonstrated that a PD controller is capable of accomplishing the control objectives. The gravity gradient and the control force keep the tether stretched throughout a normal tethered EGT operation, and the load on the tether is well within the design limit of the tether material. While including multiple particles in the tether model is essential in capturing details of tether vibration, the number of particles does not significantly affect the motions of the collected mass and the spacecraft. In addition, the distance from the asteroid mass center to the collected mass should be chosen judiciously in the case of a rotating slender asteroid; some distance ranges should be avoided as excessive lateral oscillations can be excited by resonance between the asteroid rotation and tether pendular motion.

Shen, Haijun↗

Statics of wrinkling films

An analytical investigation of the equilibrium of wrinkling films is conducted. Zak (1979) has shown that wrinkling occurs in connection with the instability of a smooth film having no resistance to bending in the case of compression. The governing equation for the equilibrium of a film with possible regions of wrinkling is considered. The introduction of fictitious stress reduces the governing equation to a form which formally coincides with the governing equation for a string. Equilibrium conditions in the case of an absence of external forces are explored, taking into account the stretching of a semispherical film, the torsion of a convex film of revolution, and stress singularities. A study is conducted of the equilibrium under conditions in which external forces normal to the surface of a film are present. Attention is also given to the equilibrium in a potential field.

Zak, M.↗

A Theoretical Investigation of the Rolling Oscillations of an Airplane with Ailerons Free

An analysis is made of the stability of an airplane with ailerons free, with particular attention to the motions when the ailerons have a tendency to float against the wind. The present analysis supersedes the aileron investigation contained in NACA Report No. 709. The equations of motion are first written to include yawing and sideslipping, and it is demonstrated that the principal effects of freeing the ailerons can be determined without regard to these motions. If the ailerons tend to float against the wind and have a high degree of aerodynamic balance, rolling oscillations, in addition to the normal lateral oscillations, are likely to occur. On the basis of the equations including only the rolling motion and the aileron deflection, formulas are derived for the stability and damping of the rolling oscillations in terms of the hinge moment derivatives and other characteristics of the ailerons and airplane. Charts are also presented showing the oscillatory regions and stability boundaries for a fictitious airplane of conventional proportions. The effects of friction in the control system are investigated and discussed. If the ailerons tend to trail with the wind, the condition for stable variation of stick force with aileron deflection is found to determine the amount of aerodynamic balance that may be used. If the ailerons tend to float against the wind, the period and damping of the rolling oscillations are found to be satisfactory (in a mass-balanced system) so long as the restoring moment is not completely balanced out. Unbalanced mass behind the hinge, however, has an unfavorable effect on the damping of the oscillations and so shifts the boundary that close aerodynamic balance may not be attainable.

Cohen, Doris↗

Time simulation of flutter with large stiffness changes

Time simulation of flutter, involving large local structural changes, is formulated with a state-space model that is based on a relatively small number of generalized coordinates. Free-free vibration modes are first calculated for a nominal finite-element model with relatively large fictitious masses located at the area of structural changes. A low-frequency subset of these modes is then transformed into a set of structural modal coordinates with which the entire simulation is performed. These generalized coordinates and the associated oscillatory aerodynamic force coefficient matrices are used to construct an efficient time-domain, state-space model for a basic aeroelastic case. The time simulation can then be performed by simply changing the mass, stiffness, and damping coupling terms when structural changes occur. It is shown that the size of the aeroelastic model required for time simulation with large structural changes at a few apriori known locations is similar to that required for direct analysis of a single structural case. The method is applied to the simulation of an aeroelastic wind-tunnel model. The diverging oscillations are followed by the activation of a tip-ballast decoupling mechanism that stabilizes the system but may cause significant transient overshoots.

Karpel, Mordechay↗

Time simulation of flutter with large stiffness changes

Time simulation of flutter, involving large local structural changes, is formulated with a state-space model that is based on a relatively small number of generalized coordinates. Free-free vibration modes are first calculated for a nominal finite-element model with relatively large fictitious masses located at the area of structural changes. A low-frequency subset of these modes is then transformed into a set of structural modal coordinates with which the entire simulation is performed. These generalized coordinates and the associated oscillatory aerodynamic force coefficient matrices are used to construct an efficient time-domain, state-space model for basic aeroelastic case. The time simulation can then be performed by simply changing the mass, stiffness and damping coupling terms when structural changes occur. It is shown that the size of the aeroelastic model required for time simulation with large structural changes at a few a priori known locations is similar to that required for direct analysis of a single structural case. The method is applied to the simulation of an aeroelastic wind-tunnel model. The diverging oscillations are followed by the activation of a tip-ballast decoupling mechanism that stabilizes the system but may cause significant transient overshoots.

Karpel, M.↗

Vibration of a liquid with a free surface in a spinning spherical tank

The problem of forced fluid vibrations in a partially filled spinning spherical tank is solved numerically by using the finite element method. The governing equations include Coriolis acceleration and spatially homogeneous vorticity. An exponential instability is detected in the present simulation for fill ratios below 0.5 and centrifugal acceleration to thrust ratios less than 1.7. This fictitious instability appears in the model as a result of the homogeneous vortex assumption since the free slosh equations are neutrally stable in the Liapunov sense.

El-Raheb, M.↗

Topological magnon band structure of emergent Landau levels in a skyrmion lattice

We report the motion of a spin excitation across topologically nontrivial magnetic order exhibits a deflection that is analogous to the effect of the Lorentz force on an electrically charged particle in an orbital magnetic field. We used polarized inelastic neutron scattering to investigate the propagation of magnons (i.e., bosonic collective spin excitations) in a lattice of skyrmion tubes in manganese silicide. For wave vectors perpendicular to the skyrmion tubes, the magnon spectra are consistent with the formation of finely spaced emergent Landau levels that are characteristic of the fictitious magnetic field used to account for the nontrivial topological winding of the skyrmion lattice. This provides evidence of a topological magnon band structure in reciprocal space, which is borne out of the nontrivial real-space topology of a magnetic order.

36 MATERIALS SCIENCE↗

Advanced Solar Electric Propulsion for Planetary Defense

Once the current survey to identify all near-Earth objects larger than 140-m diameter is complete, the major hazard from asteroid impacts may be from asteroids in the 50- to 100-m diameter size range. High-power solar electric propulsion systems employed in “slow-push” techniques may be the best way to deflect asteroids in this size range if deflection is warranted. The relative performance of three slow-push techniques—gravity tractor (GT), enhanced gravity tractor (EGT), and ion beam deflection (IBD)—are compared, assuming solar electric propulsion vehicles derived from NASA’s Asteroid Redirect Robotic Mission concept vehicle. Both the enhanced gravity tractor and ion beam deflection concepts are shown to be significantly better than the standard gravity tractor concept. The Hall-thruster based enhanced gravity tractor systems and ion beam deflection systems (based on the use of high-power gridded ion thrusters) are shown to provide comparable performance, i.e., similar deflection times and propellant required. Enhanced gravity tractor systems require the acquisition of material from the surface of the hazardous object in order to achieve the required “enhancement” of the gravitational coupling force. This makes EGT systems sensitive to the rotational state and surface properties of the unknown object, and potentially severely limits its applicability. Ion beam deflection is completely independent of the characteristics of the threat object. In fact, it is the only asteroid deflection technique, slow-push or otherwise, that can make this claim, thus potentially greatly increasing its applicability relative to other deflection techniques. Finally, high-power IBD systems are shown to be capable of deflecting the fictitious asteroid 2015 PDC used in the hypothetical asteroid impact exercise conducted at the 2015 Planetary Defense Conference.

Brophy, John↗

High-performance parallel analysis of coupled problems for aircraft propulsion

This research program deals with the application of high-performance computing methods to the numerical simulation of complete jet engines. The program was initiated in 1993 by applying two-dimensional parallel aeroelastic codes to the interior gas flow problem of a by-pass jet engine. The fluid mesh generation, domain decomposition and solution capabilities were successfully tested. Attention was then focused on methodology for the partitioned analysis of the interaction of the gas flow with a flexible structure and with the fluid mesh motion driven by these structural displacements. The latter is treated by an ALE technique that models the fluid mesh motion as that of a fictitious mechanical network laid along the edges of near-field fluid elements. New partitioned analysis procedures to treat this coupled 3-component problem were developed in 1994. These procedures involved delayed corrections and subcycling, and have been successfully tested on several massively parallel computers. For the global steady-state axisymmetric analysis of a complete engine we have decided to use the NASA-sponsored ENG10 program, which uses a regular FV-multiblock-grid discretization in conjunction with circumferential averaging to include effects of blade forces, loss, combustor heat addition, blockage, bleeds and convective mixing. A load-balancing preprocessor for parallel versions of ENG10 has been developed. It is planned to use the steady-state global solution provided by ENG10 as input to a localized three-dimensional FSI analysis for engine regions where aeroelastic effects may be important.

Felippa, C. A.↗

High-Performance Parallel Analysis of Coupled Problems for Aircraft Propulsion

Applications are described of high-performance computing methods to the numerical simulation of complete jet engines. The methodology focuses on the partitioned analysis of the interaction of the gas flow with a flexible structure and with the fluid mesh motion driven by structural displacements. The latter is treated by a ALE technique that models the fluid mesh motion as that of a fictitious mechanical network laid along the edges of near-field elements. New partitioned analysis procedures to treat this coupled three-component problem were developed. These procedures involved delayed corrections and subcycling, and have been successfully tested on several massively parallel computers, including the iPSC-860, Paragon XP/S and the IBM SP2. The NASA-sponsored ENG10 program was used for the global steady state analysis of the whole engine. This program uses a regular FV-multiblock-grid discretization in conjunction with circumferential averaging to include effects of blade forces, loss, combustor heat addition, blockage, bleeds and convective mixing. A load-balancing preprocessor for parallel versions of ENG10 was developed as well as the capability for the first full 3D aeroelastic simulation of a multirow engine stage. This capability was tested on the IBM SP2 parallel supercomputer at NASA Ames.

Felippa, C. A.↗

Four-Dimensional Oceanic and Atmosperic Data Assimilation with Tropical Rainfall Measuring Mission Data

An oceanic data assimilation system which allows to utilize the forthcoming Tropical Rainfall Measuring Mission (TRMM) data has been developed and applied to the Pacific Ocean to produce the velocity field. The assimilated data will be indispensable to examine the effects of rainfall and its variability on the structure and circulation of the tropical oceans and to assess the impact of global warming due to the increase of carbon dioxide on the ocean circulation system and the marine pollution caused by oil spill and ocean damping of radionuclide. The data will also provide the verification for the oceanic and ocean-atmosphere coupled General Circulation Models (GCM's). The system consists of oceanic GCM, analysis scheme and data. In the system the flow field has been determined to be physically consistent with the observed density field and the sea surface winds derived from the Special Sensor Microwave Imagery (SSM/I) data which drive the ocean current. The time integration has been performed for five years until the flow field near the surface attained the steady state starting from the rest ocean with observed temperature and salinity fields, and the SSM/I surface wind velocity. The resultant flow field showed high producibility of the system. Especially the flow near the ocean surface agreed well with available observed data. The system, for the first time, succeeded to produce the eastward subtropical current which has been discovered in the joint investigation on Kuroshio current (CSK) in the 1960s. To verify the quality of the flow field a trajectory analysis has been carried out and compared with the Algos buoy data. BRIEF DESCRIPTION OF THE DATA ASSIMILATION SYSTEM ## Oceanic GCM and analysis scheme--The basic equations are much the same as used for the GCM's, except for the Newtonian damping terms introduced into the prediction equations for the potential temperature and salinity to maintain these fields as observed. The C grid of 2'lat. by 2'long. in horizontal and the 11 vertical levels are applied to the entire Pacific Ocean. At the east and west ocean boundaries the periodic boundary conditions are applied creating fictitious ocean there. The SMAC Method is used to increase the accuracy of mass conservation. * Data--The JODC temperature and salinity data obtained from 1906 to 1988 are used in the system between Long.100'E. and 60'W. The surface wind data are derived from the SSM/I data by Dr-R. Atlas of NASA/GSFC. The data set contains every 6 hours data from July 1987 to June 1989 on the grid of 2'lat. by 2.5'long. The averaged for the whole period and then interpolated into the 2'lat. by 2'long. grid data are used to force the system. The sea bottom topography data was based on the General Bathymetric Chart of the Ocean (GEBCO) supplied by the Canadian Hydrographic Service under contract with the International Hydrographic Organization and International Oceanographic Commission of UNESCO.

Takano, Kenji↗

High-Performance Parallel Analysis of Coupled Problems for Aircraft Propulsion

This research program dealt with the application of high-performance computing methods to the numerical simulation of complete jet engines. The program was initiated in January 1993 by applying two-dimensional parallel aeroelastic codes to the interior gas flow problem of a bypass jet engine. The fluid mesh generation, domain decomposition and solution capabilities were successfully tested. Attention was then focused on methodology for the partitioned analysis of the interaction of the gas flow with a flexible structure and with the fluid mesh motion driven by these structural displacements. The latter is treated by a ALE technique that models the fluid mesh motion as that of a fictitious mechanical network laid along the edges of near-field fluid elements. New partitioned analysis procedures to treat this coupled three-component problem were developed during 1994 and 1995. These procedures involved delayed corrections and subcycling, and have been successfully tested on several massively parallel computers, including the iPSC-860, Paragon XP/S and the IBM SP2. For the global steady-state axisymmetric analysis of a complete engine we have decided to use the NASA-sponsored ENG10 program, which uses a regular FV-multiblock-grid discretization in conjunction with circumferential averaging to include effects of blade forces, loss, combustor heat addition, blockage, bleeds and convective mixing. A load-balancing preprocessor tor parallel versions of ENG10 was developed. During 1995 and 1996 we developed the capability tor the first full 3D aeroelastic simulation of a multirow engine stage. This capability was tested on the IBM SP2 parallel supercomputer at NASA Ames. Benchmark results were presented at the 1196 Computational Aeroscience meeting.

Felippa, C. A.↗

Sea level variation

Published values for the long-term, global mean sea level rise determined from tide gauge records range from about one to three mm per year. The scatter of the estimates appears to arise largely from the use of data from gauges located at convergent tectonic plate boundaries where changes of land elevation give fictitious sea level trends, and the effects of large interdecadal and longer sea level variations on short (less than 50+ years) or sappy records. In addition, virtually all gauges undergo subsidence or uplift due to isostatic rebound from the last deglaciation at a rate comparable to or greater than the secular rise of sea level. Modeling rebound by the ICE-3G model of Tushingham and Peltier (1990) and avoiding tide gauge records in areas of converging tectonic plates produces a highly consistent set of long sea level records. A global set of 21 such stations in nine oceanic regions with an average record length of 76 years during the period 1880-1980 yields the global sea level rise value 1.8 mm/year +/- 0.1. Greenhouse warming scenarios commonly forecast an additional acceleration of global sea level in the next 5 or 6+ decades in the range 0.1-0.2 mm/yr2. Because of the large power at low frequencies in the sea level spectrum, very long tide gauge records (75 years minimum) have been examined for past apparent sea level acceleration. For the 80-year period 1905-1985, 23 essentially complete tide gauge records in 10 geographic groups are available for analysis. These yielded the apparent global acceleration -0.011 (+/- 0.012) mm/yr2. A larger, less uniform set of 37 records in the same 10 groups with 92 years average length covering the 141 years from 1850-1991 gave 0.001 (+/- 0.008) mm/yr2. Thus there is no evidence for an apparent acceleration in the past 100+ years that is significant either statistically, or in comparison to values associated with global warming. Unfortunately, the large interdecadal fluctuations of sea level severely affect estimates of global sea level acceleration for time spans of less than about 50 years. This means that tide gauges alone cannot serve as a reliable leading indicator of climate change in less than many decades. This time required can be significantly reduced if the interdecadal fluctuations of sea level can be understood in terms of their forcing mechanisms, and then removed from the tide gauge records.

Douglas, Bruce C.↗