Nuclear fragmentation models and uncertainties in cosmic ray transport and radiobiology studies
No abstract available
NASA Discipline Radiation Health↗
SEARCH · Engineering Papers
Results for “linear response”
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.
Short latency compound action potentials from mammalian gravity receptor organs
Gravity receptor function was characterized in four mammalian species using far-field vestibular evoked potentials (VsEPs). VsEPs are compound action potentials of the vestibular nerve and central relays that are elicited by linear acceleration ramps applied to the cranium. Rats, mice, guinea pigs, and gerbils were studied. In all species, response onset occurred within 1.5 ms of the stimulus onset. Responses persisted during intense (116 dBSPL) wide-band (50 to 50 inverted question mark omitted inverted question mark000 Hz) forward masking, whereas auditory responses to intense clicks (112 dBpeSPL) were eliminated under the same conditions. VsEPs remained after cochlear extirpation but were eliminated following bilateral labyrinthectomy. Responses included a series of positive and negative peaks that occurred within 8 ms of stimulus onset (range of means at +6 dBre: 1.0 g/ms: P1=908 to 1062 micros, N1=1342 to 1475 micros, P2=1632 to 1952 micros, N2=2038 to 2387 micros). Mean response amplitudes at +6 dBre: 1.0 g/ms ranged from 0.14 to 0.99 microV. VsEP input/output functions revealed latency slopes that varied across peaks and species ranging from -19 to -51 micros/dB. Amplitude-intensity slopes also varied ranging from 0.04 to 0.08 microV/dB for rats and mice. Latency values were comparable to those of birds although amplitudes were substantially smaller in mammals. VsEP threshold values were considerably higher in mammals compared to birds and ranged from -8.1 to -10.5 dBre 1.0 g/ms across species. These results support the hypothesis that mammalian gravity receptors are less sensitive to dynamic stimuli than are those of birds.
Non-NASA Center↗
Generalized Jeans' Escape of Pick-Up Ions in Quasi-Linear Relaxation
Jeans escape is a well-validated formulation of upper atmospheric escape that we have generalized to estimate plasma escape from ionospheres. It involves the computation of the parts of particle velocity space that are unbound by the gravitational potential at the exobase, followed by a calculation of the flux carried by such unbound particles as they escape from the potential well. To generalize this approach for ions, we superposed an electrostatic ambipolar potential and a centrifugal potential, for motions across and along a divergent magnetic field. We then considered how the presence of superthermal electrons, produced by precipitating auroral primary electrons, controls the ambipolar potential. We also showed that the centrifugal potential plays a small role in controlling the mass escape flux from the terrestrial ionosphere. We then applied the transverse ion velocity distribution produced when ions, picked up by supersonic (i.e., auroral) ionospheric convection, relax via quasi-linear diffusion, as estimated for cometary comas [1]. The results provide a theoretical basis for observed ion escape response to electromagnetic and kinetic energy sources. They also suggest that super-sonic but sub-Alfvenic flow, with ion pick-up, is a unique and important regime of ion-neutral coupling, in which plasma wave-particle interactions are driven by ion-neutral collisions at densities for which the collision frequency falls near or below the gyro-frequency. As another possible illustration of this process, the heliopause ribbon discovered by the IBEX mission involves interactions between the solar wind ions and the interstellar neutral gas, in a regime that may be analogous [2].
Moore, T. E.↗
An Approach to Reasoning Service Migration in Data and Reasoning Fabric (DRF) Implementation
In this paper we consider service migration problem for Data and Reasoning Fabric (DRF) enabled airspace operations assuming a fixed cloud/edge infrastructure with allocated computing, storage and power resources, where cloud/edge servers, and communication stations are in a wired connected network, while vehicles use a wireless network for communication. The objective is to automatically select the best location for the requested service execution, which achieves minimum cost while satisfying the user quality of service (QoS) and available resources constraints. To this end, estimates of the response time, consumed energy and total cost are defined for each potential compute location. A mixed-integer linear program is then formulated and solved to identify optimal compute locations given QoS constraints, network infrastructure limitations, with worst-case vehicle positioning. The approach is applied to trajectory re-planning use case to avoid a collision with an emergency vehicle in real time.
Air mobility↗
Automated tuning of airframe vibration by structural optimization
Numerical optimization techniques are used to modify the dynamic response at a specified point(s) of a helicopter airframe structure due to a steady state narrow band excitation. Calculation of steady state vibration amplitude is reduced to the solution of linear equations with complex coefficients. The sensitivity of the dynamic amplitudes with respect to the structural parameter perturbations can be computed with the same technique as the one used in the static displacement sensitivity without requiring eigenvector sensitivity. Approximate models for critical structural responses are created based on the sensitivity data to reduce the amount of computational effort and to enable the design of structures of practical scale and complexity. This approach is general in that it accommodates static, dynamic, and frequency constraints simultaneously as long as their computational models are available. It can be used in optimizing mass distribution as well as stiffness modifications of practical structures.
Miura, H.↗
Differential equation based method for accurate approximations in optimization
A method to efficiently and accurately approximate the effect of design changes on structural response is described. The key to this method is to interpret sensitivity equations as differential equations that may be solved explicitly for closed form approximations, hence, the method is denoted the Differential Equation Based (DEB) method. Approximations were developed for vibration frequencies, mode shapes and static displacements. The DEB approximation method was applied to a cantilever beam and results compared with the commonly-used linear Taylor series approximations and exact solutions. The test calculations involved perturbing the height, width, cross-sectional area, tip mass, and bending inertia of the beam. The DEB method proved to be very accurate, and in most cases, was more accurate than the linear Taylor series approximation. The method is applicable to simultaneous perturbation of several design variables. Also, the approximations may be used to calculate other system response quantities. For example, the approximations for displacements are used to approximate bending stresses.
Pritchard, Jocelyn I.↗
Estimation of Sonic Fatigue by Reduced-Order Finite Element Based Analyses
A computationally efficient, reduced-order method is presented for prediction of sonic fatigue of structures exhibiting geometrically nonlinear response. A procedure to determine the nonlinear modal stiffness using commercial finite element codes allows the coupled nonlinear equations of motion in physical degrees of freedom to be transformed to a smaller coupled system of equations in modal coordinates. The nonlinear modal system is first solved using a computationally light equivalent linearization solution to determine if the structure responds to the applied loading in a nonlinear fashion. If so, a higher fidelity numerical simulation in modal coordinates is undertaken to more accurately determine the nonlinear response. Comparisons of displacement and stress response obtained from the reduced-order analyses are made with results obtained from numerical simulation in physical degrees-of-freedom. Fatigue life predictions from nonlinear modal and physical simulations are made using the rainflow cycle counting method in a linear cumulative damage analysis. Results computed for a simple beam structure under a random acoustic loading demonstrate the effectiveness of the approach and compare favorably with results obtained from the solution in physical degrees-of-freedom.
Rizzi, Stephen A.↗
Sustained linear acceleration
The subjective effects of sustained acceleration are discussed, including positive, negative, forward, backward, and lateral acceleration effects. Physiological effects, such as retinal and visual response, unconsciousness and cerebral function, pulmonary response, and renal output, are studied. Human tolerance and performance under sustained acceleration are ascertained.
Fraser, T. M.↗
Observed Seasonal to Decadal-Scale Responses in Mesospheric Water Vapor
The 14-yr (1991-2005) time series of mesospheric water vapor from the Halogen Occultation Experiment (HALOE) are analyzed using multiple linear regression (MLR) techniques for their6 seasonal and longer-period terms from 45S to 45N. The distribution of annual average water vapor shows a decrease from a maximum of 6.5 ppmv at 0.2 hPa to about 3.2 ppmv at 0.01 hPa, in accord with the effects of the photolysis of water vapor due to the Lyman-flux. The distribution of the semi-annual cycle amplitudes is nearly hemispherically symmetric at the low latitudes, while that of the annual cycles show larger amplitudes in the northern hemisphere. The diagnosed 11-yr, or solar cycle, max minus min, water vapor values are of the order of several percent at 0.2 hPa to about 23% at 0.01 hPa. The solar cycle terms have larger values in the northern than in the southern hemisphere, particularly in the middle mesosphere, and the associated linear trend terms are anomalously large in the same region. Those anomalies are due, at least in part, to the fact that the amplitudes of the seasonal cycles were varying at northern mid latitudes during 1991-2005, while the corresponding seasonal terms of the MLR model do not allow for that possibility. Although the 11-yr variation in water vapor is essentially hemispherically-symmetric and anti-phased with the solar cycle flux near 0.01 hPa, the concurrent temperature variations produce slightly colder conditions at the northern high latitudes at solar minimum. It is concluded that this temperature difference is most likely the reason for the greater occurrence of polar mesospheric clouds at the northern versus the southern high latitudes at solar minimum during the HALOE time period.
Remsberg, Ellis↗
Measured and predicted structural behavior of the HiMAT tailored composite wing
A series of load tests was conducted on the HiMAT tailored composite wing. Coupon tests were also run on a series of unbalanced laminates, including the ply configuration of the wing, the purpose of which was to compare the measured and predicted behavior of unbalanced laminates, including - in the case of the wing - a comparison between the behavior of the full scale structure and coupon tests. Both linear and nonlinear finite element (NASTRAN) analyses were carried out on the wing. Both linear and nonlinear point-stress analyses were performed on the coupons. All test articles were instrumented with strain gages, and wing deflections measured. The leading and trailing edges were found to have no effect on the response of the wing to applied loads. A decrease in the stiffness of the wing box was evident over the 27-test program. The measured load-strain behavior of the wing was found to be linear, in contrast to coupon tests of the same laminate, which were nonlinear. A linear NASTRAN analysis of the wing generally correlated more favorably with measurements than did a nonlinear analysis. An examination of the predicted deflections in the wing root region revealed an anomalous behavior of the structural model that cannot be explained. Both hysteresis and creep appear to be less significant in the wing tests than in the corresponding laminate coupon tests.
Nelson, Lawrence H.↗
Imaging a Free-Geometry Magma Reservoir Using FEM-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea
We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluidfilled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.
Inverse theory↗
Imaging the Complex Geometry of a Magma Reservoir Using Fem-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea
We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluid-filled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.
Inverse theory↗
Characterizing Observed Limit Cycles in the Cassini Main Engine Guidance Control System
The Cassini spacecraft dynamics-related telemetry during long Main Engine (ME) burns has indicated the presence of stable limit cycles between 0.03-0.04 Hz frequencies. These stable limit cycles cause the spacecraft to possess non-zero oscillating rates for extended periods of time. This indicates that the linear ME guidance control system does not model the complete dynamics of the spacecraft. In this study, we propose that the observed limit cycles in the spacecraft dynamics telemetry appear from a stable interaction between the unmodeled nonlinear elements in the ME guidance control system. Many nonlinearities in the control system emerge from translating the linear engine gimbal actuator (EGA) motion into a spacecraft rotation. One such nonlinearity comes from the gear backlash in the EGA system, which is the focus of this paper. The limit cycle characteristics and behavior can be predicted by modeling this gear backlash nonlinear element via a describing function and studying the interaction of this describing function with the overall dynamics of the spacecraft. The linear ME guidance controller and gear backlash nonlinearity are modeled analytically. The frequency, magnitude, and nature of the limit cycle are obtained from the frequency response of the ME guidance controller and nonlinear element. In addition, the ME guidance controller along with the nonlinearity is simulated. The simulation response contains a limit cycle with similar characterstics as predicted analytically: 0.03-0.04 Hz frequency and stable, sustained oscillations. The analytical and simulated limit cycle responses are compared to the flight telemetry for long burns such as the Saturn Orbit Insertion and Main Engine Orbit Trim Maneuvers. The analytical and simulated limit cycle characteristics compare well with the actual observed limit cycles in the flight telemetry. Both have frequencies between 0.03-0.04 Hz and stable oscillations. This work shows that the stable limit cycles occur due to the interaction between the unmodeled nonlinear elements and linear ME guidance controller.
stable limit cycle↗
Simultaneous analysis and design
Optimization techniques are increasingly being used for performing nonlinear structural analysis. The development of element by element (EBE) preconditioned conjugate gradient (CG) techniques is expected to extend this trend to linear analysis. Under these circumstances the structural design problem can be viewed as a nested optimization problem. There are computational benefits to treating this nested problem as a large single optimization problem. The response variables (such as displacements) and the structural parameters are all treated as design variables in a unified formulation which performs simultaneously the design and analysis. Two examples are used for demonstration. A seventy-two bar truss is optimized subject to linear stress constraints and a wing box structure is optimized subject to nonlinear collapse constraints. Both examples show substantial computational savings with the unified approach as compared to the traditional nested approach.
Haftka, R. T.↗
Simultaneous analysis and design
Optimization techniques are increasingly being used for performing nonlinear structural analysis. The development of element by element (EBE) preconditioned conjugate gradient (CG) techniques is expected to extend this trend to linear analysis. Under these circumstances the structural design problem can be viewed as a nested optimization problem. There are computational benefits to treating this nested problem as a large single optimization problem. The response variables (such as displacements) and the structural parameters are all treated as design variables in a unified formulation which performs simultaneously the design and analysis. Two examples are used for demonstration. A seventy-two bar truss is optimized subject to linear stress constraints and a wing box structure is optimized subject to nonlinear collapse constraints. Both examples show substantial computational savings with the unified approach as compared to the traditional nested approach.
Haftka, R. T.↗
Structural/control system optimization with variable actuator masses
A method is presented to integrate the design space for structural/control system optimization problems in the case of linear state feedback control. Nonstructural lumped masses and control system design variables as well as structural sizing variables are all treated equally as independent design variables in the optimization process. Structural and control design variable linking schemes are used in order to avoid a prohibitively large increase in the total number of independent design variables. When actuator masses are treated as nonstructural lumped mass design variables, special consideration is given to the relation between the transient peak responses and the required actuator masses which is formulated as a behavior constraint form. The original nonlinear mathematical programming problem based on a finite element formulation and linear state feedback is replaced by a sequence of explicit approximate problems exploiting various approximation concepts such as design variable linkings, temporary constraint deletion and first order Taylor series expansion of nonlinear behavior constraints in terms of intermediate design variables. Examples which involve a variety of dynamic behavior constraints (including constraints on closed-loop eigenvalues, peak transient displacements, peak actuator forces, and relations between the peak responses and the actuator masses) are effectively solved by using the method presented.
Jin, Ik M.↗
Seeing Double: ASASSN-18bt Exhibits a Two-component Rise in the Early-time K2 Light Curve
On 2018 February 4.41, the All-Sky Automated Survey for SuperNovae (ASAS-SN) discovered ASASSN-18bt in the K2 Campaign 16 field. With a redshift of z=0.01098 and a peak apparent magnitude of B(max)=14.31, ASASSN-18bt is the nearest and brightest SNe Ia yet observed by the Kepler spacecraft. Here we present the discovery of ASASSN-18bt, the K2 light curve, and prediscovery data from ASAS-SN and the Asteroid Terrestrial-impact Last Alert System. The K2 early-time light curve has an unprecedented 30-minute cadence and photometric precision for an SN Ia light curve, and it unambiguously shows a ∼4 day nearly linear phase followed by a steeper rise. Thus, ASASSN-18bt joins a growing list of SNe Ia whose early light curves are not well described by a single power law. We show that a double-power-law model fits the data reasonably well, hinting that two physical processes must be responsible for the observed rise. However, we find that current models of the interaction with a nondegenerate companion predict an abrupt rise and cannot adequately explain the initial, slower linear phase. Instead, we find that existing published models with shallow 56Ni are able to span the observed behavior and, with tuning, may be able to reproduce the ASASSN-18bt light curve. Regardless, more theoretical work is needed to satisfactorily model this and other early-time SNe Ia light curves. Finally, we use Swift X-ray nondetections to constrain the presence of circumstellar material (CSM) at much larger distances and lower densities than possible with the optical light curve. For a constant-density CSM, these nondetections constrain ρ<4.5×10(exp 5)per cu.cm at a radius of 4×10(exp 15) cm from the progenitor star. Assuming a wind-like environment, we place mass loss limits of M˙ < 8 x 10(exp -6) M(ʘ)per yr for v(w)=100 km/s, ruling out some symbiotic progenitor systems. This work highlights the power of well-sampled early-time data and the need for immediate multiband, high-cadence follow-up for progress in understanding SNe Ia.
B. J. Shappee↗
Stability and resonance in grooved-channel flows
Numerical simulation using the spectral element method is used to study the stability and resonant response of incompressible moderate-Reynolds-number flows in periodically-grooved channels. For Reynolds numbers below a critical value, the flow is shown to approach a stable steady-state, and the least stable modes resemble Tollmien-Schlichting channel waves. For Reynolds numbers greater than the critical value, self-sustained oscillations result. Oscillatory perturbation of the grooved-channel flow at the frequency of the least stable mode of the linearized system is found to result in subcritical resonant excitation as the critical Reynolds number is approached. The importance of nonhomogeneous geometry in the forced response of the flow is discussed.
Ghaddar, Nasreen K.↗