Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear mechanics”

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 1,369 records · Page 76

Temperature dependence of mechanical and thermal expansion properties of T300/5208 graphite/epoxy

The 116-394 K temperature range elastic constants, ultimate strengths and thermal expansion behavior of T300/5208 graphite/epoxy laminates of various different lay-up types have been determined. It is noted that most properties exhibit nearly linear temperature dependencies, although the transverse modulus and axial thermal expansion vary nonlinearly with temperature, and change in value by an order of magnitude over the present temperature range.

Hyer, M. W.↗

Effects of groundwater pumping on ground surface temperature: A regional modeling study in the North China Plain

Over-exploitation of groundwater (GW) in the North China Plain (NCP) since the 1960s has many environmental consequences. However, understanding of the dominant mechanisms remains limited, particularly at regional scale. In this study, the coupled ParFlow.CLM model representing subsurface and land-surface processes and their interactions was applied in the NCP at high spatio-temporal resolutions. The model was validated using the water and energy fluxes reported in previous studies and from the JRA-55 reanalysis. Numerical experiments were designed to examine the relative impacts of GW pumping and irrigation on the ground surface temperature (GST). Results show significant effects of GW pumping on GST in the NCP. Generally, the subsurface acts as a buffer to temporal variations in heat fluxes at the land-surface, but long-term pumping can gradually weaken this buffer, resulting in increases in the spatio-temporal variability of GST, as exemplified by hotter summers and colder winters. Considering that changes of water table depth (WTD) can significantly affect land surface heat fluxes when WTD ranges roughly between 0.01–10 m, the 0.5 m/year increase of WTD simulated by the model due to pumping can continue to increase the average WTD and hence, disturb GST, for about 20 years from the pre-pumping WTD in the NCP, before the WTD exceeds 10 m. The variations of GST are expected to be faster initially and gradually slow down due to the nonlinear behaviors of GST with WTD. The findings from this study in the NCP may also have implications for other regions with GW depletion.

Yang, Chen↗

Identifying Heterogeneous Micromechanical Properties of Biological Tissues via Physics–Informed Neural Networks

The heterogeneous micromechanical properties of biological tissues have profound implications across diverse medical and engineering domains. However, identifying full-field heterogeneous elastic properties of soft materials using traditional engineering approaches is fundamentally challenging due to difficulties in estimating local stress fields. Recently, there has been a growing interest in data-driven models for learning full-field mechanical responses, such as displacement and strain, from experimental or synthetic data. However, research studies on inferring full-field elastic properties of materials, a more challenging problem, are scarce, particularly for large deformation, hyperelastic materials. Here, a physics-informed machine learning approach is proposed to identify the elasticity map in nonlinear, large deformation hyperelastic materials. This study reports the prediction accuracies and computational efficiency of physics-informed neural networks (PINNs) in inferring the heterogeneous elasticity maps across materials with structural complexity that closely resemble real tissue microstructure, such as brain, tricuspid valve, and breast cancer tissues. Further, the improved architecture is applied to three hyperelastic constitutive models: Neo-Hookean, Mooney Rivlin, and Gent. Furthermore, the improved network architecture consistently produces accurate estimations of heterogeneous elasticity maps, even when there is up to 10% noise present in the training data.

59 BASIC BIOLOGICAL SCIENCES↗

Examination of stiff ion temperature gradient mode physics in simulations of DIII-D H-mode transport

A systematic evaluation of gyrokinetic and gyrofluid model predictions of ion temperature gradient (ITG) stability and transport using parameters from DIII-D high confinement mode (H-mode) plasmas has been performed. The nonlinear CGYRO code is used to make the gyrokinetic predictions, and the quasilinear TGLF model for the corresponding gyrofluid predictions. The assessments are made at three radii (normalized toroidal flux ρtor = 0.4, 0.55, and 0.7) in three different plasma scenarios with varying levels of neutral beam heating and torque. For each of the nine cases (3 radii × 3 scenarios) considered, ITG turbulence is found to be the dominant long-wavelength instability and transport mechanism. The inclusions of both transverse magnetic fluctuations and dynamic fast beam ions are stabilizing for all cases considered, with strongest effects seen at ρor = 0.4 where the fast ion population and normalized plasma pressure β = 2μ0nT/B2 are highest. The further inclusion of parallel magnetic fluctuations does not have a meaningful impact on the ITG turbulence in these scenarios, but does destabilize (in combination with fast ions) new high-frequency instabilities at ρtor = 0.4 in the high power scenarios. In each case the linear and nonlinear ITG critical gradients are predicted to be lower than the measured ITG scale lengths and their associated uncertainties. Inclusion of equilibrium flow shear in the transport predictions generally leads to an upshift in effective critical gradient rather than a qualitative change in the predicted stiffness, with stronger responses typically seen in the gyrokinetic predictions than in the gyrofluid results. However, in most cases these upshifted gradients still remain below the measured values and their uncertainties. Although the predicted critical gradients are below the measured gradients, both models predicted flux-matching gradients consistent with measured values in six of the nine cases considered, with no clear systematic over- or underprediction. Thus, while the experimental ion temperature profiles do not appear to be closely pinned to the ITG critical gradient, both gyrokinetic and gyrofluid models are able to accurately match the measured gradients reasonably well in most cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analysis of shell-type structures subjected to time-dependent mechanical and thermal loading

A general mathematical model and solution methodologies for analyzing structural response of thin, metallic shell-type structures under large transient, cyclic, or static thermomechanical loads are developed. Among the system responses, which are associated with these load conditions, are thermal buckling, creep buckling and ratcheting. Thus, geometric as well as material type nonlinearities (of high order) can be anticipated and must be considered in the development of the mathematical model. Furthermore, this must also be accommodated in the solution procedures.

Simitses, G. J.↗

Mechanical behaviour and failure phenomenon of an in situ toughened silicon nitride

The Weibull modulus, fracture toughness and crack growth resistance of an in-situ toughened, silicon nitride material used to manufacture a turbine combustor were determined from room temperature to 1371 C. The material exhibited an elongated grain structure that resulted in improved fracture toughness, nonlinear crack growth resistance, and good elevated temperature strength. However, low temperature strength was limited by grains of excessive length (30 to 100 microns). These excessively long grains were surrounded by regions rich in sintering additives.

Salem, J. A.↗

Computational mechanics analysis tools for parallel-vector supercomputers

Computational algorithms for structural analysis on parallel-vector supercomputers are reviewed. These parallel algorithms, developed by the authors, are for the assembly of structural equations, 'out-of-core' strategies for linear equation solution, massively distributed-memory equation solution, unsymmetric equation solution, general eigen-solution, geometrically nonlinear finite element analysis, design sensitivity analysis for structural dynamics, optimization algorithm and domain decomposition. The source code for many of these algorithms is available from NASA Langley.

Storaasli, O. O.↗

Analytical WKB theory for high-harmonic generation and its application to massive Dirac electrons

Here we propose an analytical approach to high-harmonic generation (HHG) for nonperturbative low-frequency and high-intensity fields based on the (Jeffreys-)Wentzel-Kramers-Brillouin (WKB) approximation. By properly taking into account Stokes phenomena of WKB solutions, we obtain wave functions that systematically include the repetitive dynamics of production and acceleration of electron-hole pairs and quantum interference due to phase accumulation between different pair production times (Stückelberg phase). Using the obtained wave functions without relying on any phenomenological assumptions, we explicitly compute electric current (including intra- and interband contributions) as the source of HHG for a massive Dirac system in (1+1) dimensions under an ac electric field. We demonstrate that the WKB approximation agrees well with numerical results obtained by solving the time-dependent Schrödinger equation and point out that the quantum interference is important in HHG. We also predict in the deep nonperturbative regime that (1) harmonic intensities oscillate with respect to electric-field amplitude E 0 and frequency Ω with a period determined by the Stückelberg phase, (2) the cutoff order of HHG is determined by 2eE 0 ℏΩ 2 , with e being the electron charge, and that (3) noninteger harmonics, controlled by the Stückelberg phase, appear as a transient effect. Our WKB theory is particularly suited for a parameter regime, where the Keldysh parameter γ = (Δ/2)Ω/eE 0 ,with Δ being the gap size, is small. This parameter regime corresponds to intense lasers in the terahertz regime for realistic massive Dirac materials. Our analysis implies that the so-called HHG plateau can be observed at the terahertz frequency within the current technology.

36 MATERIALS SCIENCE↗

Extracting dark-matter velocities from halo masses: A reconstruction conjecture

Increasing attention has recently focused on nontraditional dark-matter production mechanisms which result in primordial dark-matter velocity distributions with highly nonthermal shapes. Here, in this paper, we undertake an assessment of how the detailed shape of a general dark-matter velocity distribution impacts structure formation in the nonlinear regime. In particular, we investigate the impact on the halo-mass and subhalo-mass functions, as well as on astrophysical observables such as satellite and cluster-number counts. We find that many of the standard expectations no longer hold in situations in which this velocity distribution takes a highly nontrivial, even multimodal shape. For example, we find that the nominal free-streaming scale alone becomes insufficient to characterize the effect of free-streaming on structure formation. In addition, we propose a simple one-line conjecture which can be used to “reconstruct” the primordial dark-matter velocity distribution directly from the shape of the halo-mass function. Although our conjecture is completely heuristic, we show that it successfully reproduces the salient features of the underlying dark-matter velocity distribution even for nontrivial distributions which are highly nonthermal and/or multimodal, such as might occur for nonminimal dark sectors. Moreover, since our approach relies only on the halo-mass function, our conjecture provides a method of probing dark-matter properties even for scenarios in which the dark and visible sectors interact only gravitationally.

79 ASTRONOMY AND ASTROPHYSICS↗

Ion-rich Acceleration during an Eruptive Flux Rope Event in a Multiple Null-point Configuration

We report on the γ-ray emission above 100 MeV from the GOES M3.3 flare SOL2012-06-03. The hard X-ray (HXR) and microwave emissions have typical time profiles with a fast rise to a well-defined peak followed by a slower decay. The >100 MeV emission during the prompt phase displayed a double-peaked temporal structure with the first peak following the HXR and microwaves, and the second one, about 3 times stronger, occurring 17 ± 2 s later. The time profiles seem to indicate two separate acceleration mechanisms at work, where the second γ-ray peak reveals a potentially pure or at least largely dominant ion acceleration. The Atmospheric Imaging Assembly imaging shows a bright elliptical ribbon and a transient brightening in the northwestern (NW) region. Nonlinear force-free extrapolations at the time of the impulsive peaks show closed field lines connecting the NW region to the southeastern part of the ribbon, and the magnetic topology revealed clusters of nulls. These observations suggest a spine-and-fan geometry, and based on these observations, we interpret the second γ-ray peak as being due to the predominant acceleration of ions in a region with multiple null points. The >100 MeV emission from this flare also exhibits a delayed phase with an exponential decay of roughly 350 s. We find that the delayed emission is consistent with ions being trapped in a closed flux tube with gradual escape via their loss cone to the chromosphere.

Pesce-Rollins, Melissa [Istituto Nazionale di Fisi↗

Microcomputer based controller for the Langley 0.3-meter Transonic Cryogenic Tunnel

Flow control of the Langley 0.3-meter Transonic Cryogenic Tunnel (TCT) is a multivariable nonlinear control problem. Globally stable control laws were generated to hold tunnel conditions in the presence of geometrical disturbances in the test section and precisely control the tunnel states for small and large set point changes. The control laws are mechanized as four inner control loops for tunnel pressure, temperature, fan speed, and liquid nitrogen supply pressure, and two outer loops for Mach number and Reynolds number. These integrated control laws have been mechanized on a 16-bit microcomputer working on DOS. This document details the model of the 0.3-m TCT, control laws, microcomputer realization, and its performance. The tunnel closed loop responses to small and large set point changes were presented. The controller incorporates safe thermal management of the tunnel cooldown based on thermal restrictions. The controller was shown to provide control of temperature to + or - 0.2K, pressure to + or - 0.07 psia, and Mach number to + or - 0.002 of a given set point during aerodynamic data acquisition in the presence of intrusive geometrical changes like flexwall movement, angle-of-attack changes, and drag rake traverse. The controller also provides a new feature of Reynolds number control. The controller provides a safe, reliable, and economical control of the 0.3-m TCT.

Balakrishna, S.↗

Evidence of chaotic pattern in solar flux through a reproducible sequence of period-doubling-type bifurcations

Presented here is a preliminary study of the limits to solar flux intensity prediction, and of whether the general lack of predictability in the solar flux arises from the nonlinear chaotic nature of the Sun's physical activity. Statistical analysis of a chaotic signal can extract only its most gross features, and detailed physical models fail, since even the simplest equations of motion for a nonlinear system can exhibit chaotic behavior. A recent theory by Feigenbaum suggests that nonlinear systems that can be led into chaotic behavior through a sequence of period-doubling bifurcations will exhibit a universal behavior. As the control parameter is increased, the bifurcation points occur in such a way that a proper ratio of these will approach the universal Feigenbaum number. Experimental evidence supporting the applicability of the Feigenbaum scenario to solar flux data is sparse. However, given the hypothesis that the Sun's convection zones are similar to a Rayleigh-Bernard mechanism, we can learn a great deal from the remarkable agreement observed between the prediction by theory (period doubling - a universal route to chaos) and the amplitude decrease of the signal's regular subharmonics. The authors show that period-doubling-type bifurcation is a possible route to a chaotic pattern of solar flux that is distinguishable from the logarithm of its power spectral density. This conclusion is the first positive step toward a reformulation of solar flux by a nonlinear chaotic approach. The ultimate goal of this research is to be able to predict an estimate of the upper and lower bounds for solar flux within its predictable zones. Naturally, it is an important task to identify the time horizons beyond which predictability becomes incompatible with computability.

Ashrafi, S.↗

Evidence of chaotic pattern in solar flux through a reproducible sequence of period-doubling-type bifurcations

A preliminary study of the limits to solar flux intensity prediction, and of whether the general lack of predictability in the solar flux arises from the nonlinear chaotic nature of the Sun's physical activity is presented. Statistical analysis of a chaotic signal can extract only its most gross features, and detailed physical models fail, since even the simplest equations of motion for a nonlinear system can exhibit chaotic behavior. A recent theory by Feigenbaum suggests that nonlinear systems that can be led into chaotic behavior through a sequence of period-doubling bifurcations will exhibit a universal behavior. As the control parameter is increased, the bifurcation points occur in such a way that a proper ratio of these will approach the universal Feigenbaum number. Experimental evidence supporting the applicability of the Feigenbaum scenario to solar flux data is sparse. However, given the hypothesis that the Sun's convection zones are similar to a Rayleigh-Bernard mechanism, we can learn a great deal from the remarkable agreement observed between the prediction by theory (period doubling - a universal route to chaos) and the amplitude decrease of the signal's regular subharmonics. It is shown that period-doubling-type bifurcation is a possible route to a chaotic pattern of solar flux that is distinguishable from the logarithm of its power spectral density. This conclusion is the first positive step toward a reformulation of solar flux by a nonlinear chaotic approach. The ultimate goal of this research is to be able to predict an estimate of the upper and lower bounds for solar flux within its predictable zones. Naturally, it is an important task to identify the time horizons beyond which predictability becomes incompatible with computability.

Ashrafi, S.↗

Investigation of the flight mechanics simulation of a hovering helicopter

The flight mechanics simulation of a hovering helicopter is investigated by comparing the results of two different numerical models with flight test data for a hovering AH-64 Apache. The two models are the U.S. Army BEMAP and the Technion model. These nonlinear models are linearized by applying a numerical linearization procedure. The results of the linear models are compared with identification results in terms of eigenvalues, stability and control derivatives, and frequency responses. Detailed time histories of the responses of the complete nonlinear models, as a result of various pilots' inputs, are compared with flight test results. In addition the sensitivity of the models to various effects are also investigated. The results are discussed and problematic aspects of the simulation are identified.

Chaimovich, M.↗

Shock tube experiments on the three-layer Richtmyer–Meshkov instability

Here, a vertical shock tube is used for experiments on the three-layer Richtmyer–Meshkov instability. Two closely spaced membrane-less interfaces are formed by the flow of two different sects of three gases: one with air above CO 2 above SF 6 and the other with helium above air above SF 6 . The lightest of the three gases enters the shock tube at the top of the driven section and flows downward. Conversely, the heaviest gas enters at the bottom of the shock tube and flows upward while the intermediate density gas enters at the middle through porous plates. All three gases are allowed to escape through holes at the layer location, leaving an approximately 30-mm layer of intermediate-density gas suspended between the lightest gas from above and the heaviest gas from below. A single-mode, two-dimensional initial perturbation is then imposed on the lower interface by oscillating the shock tube in the horizontal direction. The flow is visualized by seeding the intermediate gas with particles and illuminating it with a pulsed laser. Image sequences are then captured using high-speed video cameras. Perturbation amplitude measurements are made from the three-layer system and compared with measurements from 2, two-layer systems. It is observed that the presence of the upper, initially flat interface produces a decrease in growth of instability amplitude in the nonlinear phase over an equivalent single-interface configuration.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Neural units with time-dependent functionality

We show that the time-resolved dynamics of an underdamped harmonic oscillator can be used to do multifunctional computation, performing distinct computations at distinct times within a single dynamical trajectory. We consider the amplitude of an oscillator whose inputs influence its frequency. The activity of the oscillator at fixed times is a nonmonotonic function of its inputs, so it can solve problems such as XOR that are not linearly separable. The activity of the oscillator at fixed input is a nonmonotonic function of time, so it is multifunctional in a temporal sense, and able to carry out distinct nonlinear computations at distinct times within the same dynamical trajectory. We show that a single oscillator, observed at different times, can act as all of the elementary logic gates and perform binary addition, the latter usually implemented in hardware using five logic gates. We show that a set of n oscillators, observed at different times, can perform an arbitrary number of analog-to-n-bit digital conversions. We also show that oscillators can be trained by gradient descent to perform distinct classification tasks at distinct times. Computing with time-dependent functionality can be done in or out of equilibrium, and suggests a way of reducing the number of parameters or devices required to do nonlinear computations.

97 MATHEMATICS AND COMPUTING↗

Adaptive Shape Functions and Internal Mesh Adaptation for Modelling Progressive Failure in Adhesively Bonded Joints

Enhanced finite elements are elements with an embedded analytical solution that can capture detailed local fields, enabling more efficient, mesh independent finite element analysis. The shape functions are determined based on the analytical model rather than prescribed. This method was applied to adhesively bonded joints to model joint behavior with one element through the thickness. This study demonstrates two methods of maintaining the fidelity of such elements during adhesive non-linearity and cracking without increasing the mesh needed for an accurate solution. The first method uses adaptive shape functions, where the shape functions are recalculated at each load step based on the softening of the adhesive. The second method is internal mesh adaption, where cracking of the adhesive within an element is captured by further discretizing the element internally to represent the partially cracked geometry. By keeping mesh adaptations within an element, a finer mesh can be used during the analysis without affecting the global finite element model mesh. Examples are shown which highlight when each method is most effective in reducing the number of elements needed to capture adhesive nonlinearity and cracking. These methods are validated against analogous finite element models utilizing cohesive zone elements.

computational mechanics↗

Analysis of shell type structures subjected to time dependent mechanical and thermal loading

A general mathematical model and solution methodologies for analyzing structural response of thin, metallic shell-type structures under large transient, cyclic or static thermomechanical loads is considered. Among the system responses, which are associated with these load conditions, are thermal buckling, creep buckling and ratchetting. Thus, geometric as well as material-type nonlinearities (of high order) can be anticipated and must be considered in the development of the mathematical model.

Simitses, G. J.↗