Numerical Issues for a Non-autonomous Logistic Model
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
We have recently conducted a controlled comparison of a number of approximations for gravitational clustering against the same n-body simulations. These include ordinary linear perturbation theory (Eulerian), the adhesion approximation, the frozen-flow approximation, the Zel'dovich approximation (describable as first-order Lagrangian perturbation theory), and its second-order generalization. In the last two cases we also created new versions of approximation by truncation, i.e., smoothing the initial conditions by various smoothing window shapes and varying their sizes. The primary tool for comparing simulations to approximation schemes was crosscorrelation of the evolved mass density fields, testing the extent to which mass was moved to the right place. The Zel'dovich approximation, with initial convolution with a Gaussian e(exp -k(exp 2)/k(exp 2, sub G)) where k(sub G) is adjusted to be just into the nonlinear regime of the evolved model (details in text) worked extremely well. Its second-order generalization worked slightly better. All other schemes, including those proposed as generalizations of the Zel'dovich approximation created by adding forces, were in fact generally worse by this measure. By explicitly checking, we verified that the success of our best-choice was a result of the best treatment of the phases of nonlinear Fourier components. Of all schemes tested, the adhesion approximation produced the most accurate nonlinear power spectrum and density distribution, but its phase errors suggest mass condensations were moved to slightly the wrong location. Due to its better reproduction of the mass density distribution function and power spectrum, it might be preferred for some uses. We recommend either n-body simulations or our modified versions of the Zel'dovich approximation, depending upon the purpose. The theoretical implication is that pancaking is implicit in all cosmological gravitational clustering, at least from Gaussian initial conditions, even when subcondensations are present.
The objective of this paper is to compare the results of several simulations performed to determine the worst-case location for a foam impact on the Space Shuttle wing leading edge. The simulations were performed using the commercial non-linear transient dynamic finite element code, LS-DYNA. These simulations represent the first in a series of parametric studies performed to support the selection of the worst-case impact scenario. Panel 9 was selected for this study to enable comparisons with previous simulations performed during the Columbia Accident Investigation. The projectile for this study is a 5.5-in cube of typical external tank foam weighing 0.23 lb. Seven locations spanning the panel surface were impacted with the foam cube. For each of these cases, the foam was traveling at 1000 ft/s directly aft, along the orbiter X-axis. Results compared from the parametric studies included strains, contact forces, and material energies for various simulations. The results show that the worst case impact location was on the top surface, near the apex.
The linear model of the Stirling system developed by NASA Glenn Research Center (GRC) has been extended to include a user-specified heat input. Previously developed linear models were limited to the Stirling convertor and electrical load. They represented the thermodynamic cycle with pressure factors that remained constant. The numerical values of the pressure factors were generated by linearizing GRC s non-linear System Dynamic Model (SDM) of the convertor at a chosen operating point. The pressure factors were fixed for that operating point, thus, the model lost accuracy if a transition to a different operating point were simulated. Although the previous linear model was used in developing controllers that manipulated current, voltage, and piston position, it could not be used in the development of control algorithms that regulated hot-end temperature. This basic model was extended to include the thermal dynamics associated with a hot-end temperature that varies over time in response to external changes as well as to changes in the Stirling cycle. The linear model described herein includes not only dynamics of the piston, displacer, gas, and electrical circuit, but also the transient effects of the heater head thermal inertia. The linear version algebraically couples two separate linear dynamic models, one model of the Stirling convertor and one model of the thermal system, through the pressure factors. The thermal system model includes heat flow of heat transfer fluid, insulation loss, and temperature drops from the heat source to the Stirling convertor expansion space. The linear model was compared to a nonlinear model, and performance was very similar. The resulting linear model can be implemented in a variety of computing environments, and is suitable for analysis with classical and state space controls analysis techniques.
The ESA-funded "Cross-Scale Technology Reference Study has been carried out with the primary aim to identify and analyse a mission concept for the investigation of fundamental space plasma processes that involve dynamical non-linear coupling across multiple length scales. To fulfill this scientific mission goal, a constellation of spacecraft is required, flying in loose formations around the Earth and sampling three characteristic plasma scale distances simultaneously, with at least two satellites per scale: electron kinetic (~10 km), ion kinetic (~100-2000 km), magnetospheric fluid (~3000-15000 km). The key Cross-Scale mission drivers identified are the number of S/C, the space segment configuration, the reference orbit design, the transfer and deployment strategy, the inter-satellite localization and synchronization process and the mission operations. This paper presents a comprehensive overview of the mission design and analysis for the Cross-Scale concept and outlines a technically feasible mission architecture for a multi-dimensional investigation of space plasma phenomena. The main effort has been devoted to apply a thorough mission-level trade-off approach and to accomplish an exhaustive analysis, so as to allow the characterization of a wide range of mission requirements and design solutions.
This report summarizes initial modeling of the local response of the Bigelow Expandable Activity Module (BEAM) to micrometeorite and orbital debris (MMOD) impacts using a structural, non-linear, transient dynamic finite element code. Complementary test results for a local BEAM structure are presented for both hammer and projectile impacts. Review of these data provided guidance for the transient dynamic model development. The local model is intended to support predictions using the global BEAM model, described in a companion report. Two types of local models were developed. One mimics the simplified Soft-Goods (fabric envelop) part of the BEAM NASTRAN model delivered by the project. The second investigates through-the-thickness modeling challenges for MMOD-type impacts. Both the testing and the analysis summaries contain lessons learned and areas for future efforts.
Future multi-vehicle operations will involve multiple scenarios that will require a planning tool for the rapid, interactive creation of fuel-efficient trajectories. The planning process must deal with higher-order, non-linear processes involving dynamics that are often counter-intuitive. The optimization of resulting trajectories can be difficult to envision. An interaction proximity operations planning system is being developed to provide the operator with easily interpreted visual feedback of trajectories and constraints. This system is hosted on an IRIS 4D graphics platform and utilizes the Clohessy-Wiltshire equations. An inverse dynamics algorithm is used to remove non-linearities while the trajectory maneuvers are decoupled and separated in a geometric spreadsheet. The operator has direct control of the position and time of trajectory waypoints to achieve the desired end conditions. Graphics provide the operator with visualization of satisfying operational constraints such as structural clearance, plume impingement, approach velocity limits, and arrival or departure corridors. Primer vector theory is combined with graphical presentation to improve operator understanding of suggested automated system solutions and to allow the operator to review, edit, or provide corrective action to the trajectory plan.
SAND2023-05107O JuliaOceanWaves is a collection of packages and functions for modeling wave energy converters. The open-source software is designed to be used by the wave energy community to model new concepts, perform design optimization, and more. Capabilities include: • A representation of ocean wave spectra and input/output (e.g. read from National Oceanic and Atmospheric Administration buoys) and post-processing functionalities; • A solver for the linear potential flow radiation and diffraction equations for wave forcing on floating structures. • Solving the floating body linear dynamics in the frequency domain. Future additions to the software will include solving non-linear floating body dynamics in the time-domain through time-marching Ordinary Differential Equation solvers and the pseudo-spectral optimization method. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.
For Titan our interest during the past few years was to explain the observed asymmetry in the albedo. We suggested earlier, from one-dimensional modeling studies, that vertical transport rates were comparable to particle fall speeds. Since heating of the upper atmosphere, which drives dynamical motions, is largely due to the aerosols, a non- linear interaction between dynamics, radiative heating and particle microphysics is possible. We pursued this interaction in a two-dimensional model. We showed that the observed variations in the albedo between the two hemispheres and over an orbital cycle, could be due to dynamical motions suspending particles so that particle sizes and optical depths vary across the planet. In the Cassini time frame, future studies of this interaction between dynamics, radiation and microphysics may be worthwhile using the strong modeling base that others, and we have developed. In our recently approved proposal, however, we plan to extend our modeling to hydrocarbon clouds that lie at lower levels. We know very little about such clouds, and numerical models for their properties are non-existent. These clouds may be observed by the Huygen's Probe, and by the Cassini orbiter, so predictions of their properties should help in the analysis of Cassini data. We have also developed a sophisticated model for the lower, condensational, clouds on Venus. In this model we explored the water vapor budget on Venus, and the properties of the clouds such as particle size distribution. Most researchers have investigated the upper clouds on Venus, which are essentially a photochemical aerosol with a long lifetime. The lower clouds, however, are similar to stratus clouds on Earth. These clouds have short lifetimes, and are tightly coupled to the dynamics at the base of the Venus cloud deck. We believe that the holes in these clouds seen at near infrared wavelengths by Galileo are related to some interaction between dynamics and cloud physics. One goal of our recently approved proposed work is to better understand this interaction, and to better understand how these clouds may vary if the climate of Venus were to vary. During the past few years we have applied our model of the water ice clouds on Mars to new data sets from Pathfinder, and Mars Global Surveyor. We have compared our predictions of cloud properties with those seen by Pathfinder, and found reasonable agreement. More recently we have improved our model by including a radiative transfer algorithm and a boundary layer transport scheme . The goal of this work is to understand the thermal inversions seen in Mars Global Surveyor data. We find that the inversions are created by clouds radiatively cooling the atmospheric layers in which the clouds form. We have also initiated laboratory work on the physics of carbon dioxide particle formation. In our recently approved work we propose to use these laboratory data, and to extend our modeling to carbon dioxide clouds in the Martian atmosphere.
Description of an efficient digital computer method for the determination of the propagation of elastic stresses and deformations in certain geometrically nonlinear structures subjected to high impact loading. The finite-element matrix displacement approach utilizing curved quadrilateral shell elements in conjunction with a nodewise predictor-corrector method employing Runge-Kutta extrapolation techniques has been adopted for the present solution. The related computer program written in FORTRAN V for the UNIVAC 1108 computer has proved to be effective for the solution of a range of practical problems including rectangular and cylindrical panels. Numerical results are presented for a relevant structure, the cell container, and the negative electrode of an impact-resistant battery subjected to high impact, simulating its free landing on a planetary surface.
In order to be cost effective, space structures must be extremely light weight, and subsequently, very flexible structures. The power system for Space Station Freedom is such a structure. Each array consists of a deployable truss mast and a split blanket of photovoltaic solar collectors. The solar arrays are deployed in orbit, and the blanket is stretched into position as the mast is extended. Geometric stiffness due to the preload make this an interesting non-linear problem. The space station will be subjected to various dynamic loads, during shuttle docking, solar tracking, attitude adjustment, etc. Accurate prediction of the natural frequencies and mode shapes of the space station components, including the solar arrays, is critical for determining the structural adequacy of the components, and for designing a dynamic controls system. The process used in developing and verifying the finite element dynamic model of the photo-voltaic arrays is documented. Various problems were identified, such as grounding effects due to geometric stiffness, large displacement effects, and pseudo-stiffness (grounding) due to lack of required rigid body modes. Analysis techniques, such as development of rigorous solutions using continuum mechanics, finite element solution sequence altering, equivalent systems using a curvature basis, Craig-Bampton superelement approach, and modal ordering schemes were utilized. The grounding problems associated with the geometric stiffness are emphasized.
In order to be cost effective, space structures must be extremely light weight, and subsequently, very flexible structures. The power system for Space Station Freedom is such a structure. Each array consists of a deployable truss mast and a split blanket of photo-voltaic solar collectors. The solar arrays are deployed in orbit, and the blanket is stretched into position as the mast is extended. Geometric stiffness due to the preload make this an interesting non-linear problem. The space station will be subjected to various dynamic loads, during shuttle docking, solar tracking, attitude adjustment, etc. Accurate prediction of the natural frequencies and mode shapes of the space station components, including the solar arrays, is critical for determining the structural adequacy of the components, and for designing a dynamic control system. The process used in developing and verifying the finite element dynamic model of the photo-voltaic arrays is documented. Various problems were identified, such as grounding effects due to geometric stiffness, large displacement effects, and pseudo-stiffness (grounding) due to lack of required rigid body modes. Analysis techniques, such as development of rigorous solutions using continuum mechanics, finite element solution sequence altering, equivalent systems using a curvature basis, Craig-Bampton superelement approach, and modal ordering schemes were utilized. The grounding problems associated with the geometric stiffness are emphasized.
The dynamic behavior of a reconfigurable adaptive truss structure with non-linear joints is investigated. The objective is to experimentally examine the effects of the local non-linearities on the global dynamics of the structure. Amplitude changes in the frequency response functions are measured at micron levels of motion. The amplitude and frequency variations of a number of modes indicate a non-linear Coulomb friction response. Hysteretic bifurcation behavior is also measured at an amplitude approximately equal to the specified free-play in the joint. Under the 1 g pre-load, however, the non-linearity was dominantly characteristic of Coulomb friction with little evidence of free-play stiffening.
The Final Technical Report on linear and non-linear studies of Alfven waves in space is presented. Areas of research included relaxation of magnetotail plasmas with field-aligned currents; the equilibrium dayside magnetosphere; macroscale particle simulation of kinetic Alfven wave physics; ballooning stability of plasmas with sheared equilibrium flows; theory of the drift-mirror instability; collisionless tearing instability in magnetotail plasmas; and nonadiabatic behavior of the magnetic moment of a charged particle in a dipole magnetic field and the development of stochastic webs.
In this work, we comprehensively analyse the cosmology dependence of counts-in-cells statistics. Furthermore, we focus on the shape of the one-point probability distribution function (PDF) of the matter density field at mildly non-linear scales. Based on large-deviation statistics, we parametrize the cosmology dependence of the matter PDF in terms of the linear power spectrum, the growth factor, the spherical collapse dynamics, and the non-linear variance. We extend our formalism to include massive neutrinos, finding that the total matter PDF is highly sensitive to the total neutrino mass Mν and can disentangle it from the clustering amplitude σ 8 . Additionally, using more than a million PDFs extracted from the Quijote simulations, we determine the response of the matter PDF to changing parameters in the νΛCDM model and successfully cross-validate the theoretical model and the simulation measurements. We present the first νΛCDM Fisher forecast for the matter PDF at multiple scales and redshifts, and its combination with the matter power spectrum. We establish that the matter PDF and the matter power spectrum are highly complementary at mildly non-linear scales. The matter PDF is particularly powerful for constraining the matter density Ωm, clustering amplitude σ 8 and the total neutrino mass Mν. Adding the mildly non-linear matter PDF to the mildly non-linear matter power spectrum improves constraints on Ω m by a factor of 5 and σ 8 by a factor of 2 when considering the three lowest redshifts. In our joint analysis of the matter PDF and matter power spectrum at three redshifts, the total neutrino mass is constrained to better than 0.01 eV with a total volume of 6 (Gpc h -1 ) 3 . We discuss how density-split statistics can be used to translate those encouraging results for the matter PDF into realistic observables in galaxy surveys.
Introduction: Physical and non-physical processes that occur in nature may influence biological processes, such as dissemination of infectious diseases. However, such processes may be hard to detect when they are complex systems. Because complexity is a dynamic and non-linear interaction among numerous elements and structural levels in which specific effects are not necessarily linked to any one specific element, cause-effect connections are rarely or poorly observed. Methods: To test this hypothesis, the complex and dynamic properties of geo-biological data were explored with high-resolution epidemiological data collected in the 2001 Uruguayan foot-and-mouth disease (FMD) epizootic that mainly affected cattle. County-level data on cases, farm density, road density, river density, and the ratio of road (or river) length/county perimeter were analyzed with an open-ended procedure that identified geographical clustering in the first 11 epidemic weeks. Two questions were asked: (i) do geo-referenced epidemiologic data display complex properties? and (ii) can such properties facilitate or prevent disease dissemination? Results: Emergent patterns were detected when complex data structures were analyzed, which were not observed when variables were assessed individually. Complex properties–including data circularity–were demonstrated. The emergent patterns helped identify 11 counties as ‘disseminators’ or ‘facilitators’ (F) and 264 counties as ‘barriers’ (B) of epidemic spread. In the early epidemic phase, F and B counties differed in terms of road density and FMD case density. Focusing on non-biological, geographical data, a second analysis indicated that complex relationships may identify B-like counties even before epidemics occur. Discussion: Geographical barriers and/or promoters of disease dispersal may precede the introduction of emerging pathogens. If corroborated, the analysis of geo-referenced complexity may support anticipatory epidemiological policies.
We demonstrate magnetic droplet soliton pairs in all-perpendicular spin-torque nano-oscillators (STNOs), where one droplet resides in the STNO free layer (FL) and the other in the reference layer (RL). Typically, theoretical, numerical, and experimental droplet studies have focused on the FL, with any additional dynamics in the RL entirely ignored. Here we show that there is not only significant magnetodynamics in the RL, but the RL itself can host a droplet driven by, and coexisting with, the FL droplet. Both single droplets and pairs are observed experimentally as stepwise changes and sharp peaks in the dc and differential resistance, respectively. While the single FL droplet is highly stable, the coexistence state exhibits high-power broadband microwave noise. Furthermore, micromagnetic simulations reveal that the pair dynamics display periodic, quasi-periodic, and chaotic signatures controlled by applied field and current. The strongly interacting and closely spaced droplet pair offers a unique platform for fundamental studies of highly non-linear soliton pair dynamics.
Systems and methods for high-speed compression of dynamic electrical signal waveforms to extend the measuring capabilities of conventional measuring devices such as oscilloscopes and high-speed data acquisition systems are discussed. Transfer function components and algorithmic transfer functions can be used to accurately measure signals that are within the frequency bandwidth but beyond the voltage range and voltage resolution capabilities of the measuring device.