Search NASA⌕ Search

SEARCH · Search NASA

Results for “Analytic linearization”

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 865 records · Page 48

The analytical representation of viscoelastic material properties using optimization techniques

This report presents a technique to model viscoelastic material properties with a function of the form of the Prony series. Generally, the method employed to determine the function constants requires assuming values for the exponential constants of the function and then resolving the remaining constants through linear least-squares techniques. The technique presented here allows all the constants to be analytically determined through optimization techniques. This technique is employed in a computer program named PRONY and makes use of commercially available optimization tool developed by VMA Engineering, Inc. The PRONY program was utilized to compare the technique against previously determined models for solid rocket motor TP-H1148 propellant and V747-75 Viton fluoroelastomer. In both cases, the optimization technique generated functions that modeled the test data with at least an order of magnitude better correlation. This technique has demonstrated the capability to use small or large data sets and to use data sets that have uniformly or nonuniformly spaced data pairs. The reduction of experimental data to accurate mathematical models is a vital part of most scientific and engineering research. This technique of regression through optimization can be applied to other mathematical models that are difficult to fit to experimental data through traditional regression techniques.

Hill, S. A.↗

A modal parameter extraction procedure applicable to linear time-invariant dynamic systems

Modal analysis has emerged as a valuable tool in many phases of the engineering design process. Complex vibration and acoustic problems in new designs can often be remedied through use of the method. Moreover, the technique has been used to enhance the conceptual understanding of structures by serving to verify analytical models. A new modal parameter estimation procedure is presented. The technique is applicable to linear, time-invariant systems and accommodates multiple input excitations. In order to provide a background for the derivation of the method, some modal parameter extraction procedures currently in use are described. Key features implemented in the new technique are elaborated upon.

Kurdila, A. J.↗

Rayleigh-Taylor instabilities with sheared magnetic fields

Magnetic Rayleigh-Taylor (MRT) instabilities may play a relevant role in many astrophysical problems. In this work the effect of magnetic shear on the growth rate of the MRT instability is investigated. The eigenmodes of an interface and a slab model under the presence of gravity are analytically calculated assuming that the orientation of the magnetic field changes in the equilibrium, i.e., there is magnetic shear. We solve the linearized magnetohydrodynamic equations in the incompressible regime. We find that the growth rate is bounded under the presence of magnetic shear. We have derived simple analytical expressions for the maximum growth rate, corresponding to the most unstable mode of the system. These expressions provide the explicit dependence of the growth rate on the various equilibrium parameters. For small angles the growth time is linearly proportional to the shear angle, and in this regime the single interface problem and the slab problem tend to the same result. On the contrary, in the limit of large angles and for the interface problem the growth time is essentially independent of the shear angle. In this regime we have also been able to calculate an approximate expression for the growth time for the slab configuration. Magnetic shear can have a strong effect on the growth rates of the instability. As an application of the results found in this paper we have indirectly determined the shear angle in solar prominence threads using their lifetimes and the estimation of the Alfvén speed of the structure.

79 ASTRONOMY AND ASTROPHYSICS↗

Shock compression modeling of metallic single crystals: comparison of finite difference, steady wave, and analytical solutions

Background: The shock response of metallic single crystals can be captured using a micro-mechanical description of the thermoelastic-viscoplastic material response; however, using a such a description within the context of traditional numerical methods may introduce a physical artifacts. Advantages and disadvantages of complex material descriptions, in particular the viscoplastic response, must be framed within approximations introduced by numerical methods. Methods: Three methods of modeling the shock response of metallic single crystals are summarized: finite difference simulations, steady wave simulations, and algebraic solutions of the Rankine-Hugoniot jump conditions. For the former two numerical techniques, a dislocation density based framework describes the rate- and temperature-dependent shear strength on each slip system. For the latter analytical technique, a simple (two-parameter) rate- and temperature-independent linear hardening description is necessarily invoked to enable simultaneous solution of the governing equations. For all models, the same nonlinear thermoelastic energy potential incorporating elastic constants of up to order 3 is applied. Results: Solutions are compared for plate impact of highly symmetric orientations (all three methods) and low symmetry orientations (numerical methods only) of aluminum single crystals shocked to 5 GPa (weak shock regime) and 25 GPa (overdriven regime). Conclusions: For weak shocks, results of the two numerical methods are very similar, regardless of crystallographic orientation. For strong shocks, artificial viscosity affects the finite difference solution, and effects of transverse waves for the lower symmetry orientations not captured by the steady wave method become important. The analytical solution, which can only be applied to highly symmetric orientations, provides reasonable accuracy with regards to prediction of most variables in the final shocked state but, by construction, does not provide insight into the shock structure afforded by the numerical methods.

36 MATERIALS SCIENCE↗

Effective Pressure Law for the Intrinsic Formation Factor in Low Permeability Sandstones: EPL law for the formation factor

The effective pressure law for the intrinsic formation factor of low permeability sandstones from the Ordos Basin (northwest of China) was studied experimentally by measuring the electrical resistivity of brine saturated rock samples while cycling the pore and confining pressures, p p and p c . We used a 100,000 ppm NaCl solution so that the results could be directly interpreted in terms of the intrinsic formation factor F. The Response Surface method was used to construct a quadratic function of pp and pc fitting the experimental data, from which the coefficient α of the effective pressure law, p eff = p c – αp p , was determined. We found that the coefficient α generally had low values, mainly from 0.2 to 0.4, thus contradicting the theoretical prediction that α should be equal to 1 for scale-invariant properties, such as F, in linear elastic materials having a homogeneous solid matrix. Since the two underlying assumptions, linear elasticity and homogeneity of the solid matrix, are likely violated in reservoir rocks, we tried to assess their effects on α using simple analytical models, based on idealized geometries of the pore and solid phases. Analysis of these models suggests that non-linear elasticity associated with the formation of solid–solid contacts in the compressed rocks (e.g., during closure of rough cracks or crack-like pores along the grain boundaries), may be more effective in producing the low α values observed than inhomogeneity of the solid matrix.

58 GEOSCIENCES↗

Gauge-invariant gravitational waves in matter beyond linearized gravity

Modeling the propagation of gravitational waves (GWs) in media other than vacuum is complicated by the gauge freedom of linearized gravity in that, once nonlinearities are taken into consideration, gauge artifacts can cause spurious acceleration of the matter. To eliminate these artifacts, we propose how to keep the theory of dispersive GWs gauge-invariant beyond the linear approximation and, in particular, obtain an unambiguous gauge-invariant expression for the energy–momentum of a GW in a dispersive medium. Using analytic tools from plasma physics, we propose an exactly gauge-invariant 'quasilinear' theory, in which GWs are governed by linear equations and also affect the background metric on scales large compared to their wavelength. As a corollary, the gauge-invariant geometrical optics of linear dispersive GWs in a general background is formulated. As an example, we show how the well-known properties of vacuum GWs are naturally and concisely yielded by our theory in a manifestly gauge-invariant form. We also show how the gauge invariance can be maintained within a given accuracy to an arbitrary order in the GW amplitude. These results are intended to form a physically meaningful framework for studying dispersive GWs in matter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A slip-spring framework to study relaxation dynamics of entangled wormlike micelles with kinetic Monte Carlo algorithm

Hypothesis: Wormlike micelles (WLMs) formed due to the self-assembly of amphiphiles in aqueous solution have similar viscoelastic properties as polymers. Owing to this similarity, in this work, it is postulated that kinetic Monte Carlo (kMC) sampling of slip-springs dynamics, which is able to model the rheology of polymers, can also be extended to capture the relaxation dynamics of WLMs. Theory: The proposed modeling framework considers the following relaxation mechanisms: reptation, union-scission, and constraint release. Specifically, each of these relaxation mechanisms is simulated as separate kMC events that capture the relaxation dynamics while considering the living nature of WLMs within the slip-spring framework. As a case study, the model is implemented to a system of sodium oleate and sodium chloride to predict the linear rheology and the characteristic relaxation times associated with the individual relaxation mechanisms at different pH and salt concentrations. Findings: Linear rheology predictions were found to be in good agreement with experimental data. Furthermore, the calculated relaxation times highlighted that reptation contributed to a continuous increase in viscosity while union-scission contributed to the decrease in viscosity of WLM solutions at a higher salinity and pH. Finally, this manifests the proposed model’s capability to provide insights into the key processes governing WLM’s rheology.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Two-dimensional models for the optical response of thin films

Here, we present a systematic study of 2D optical models for the response of thin layers of material under excitation by normally incident light. The treatment, within the framework of classical optics, analyzes a thin film supported by a semi-infinite substrate, with both the thin layer and the substrate assumed to exhibit local, isotropic linear response. Starting from the conventional three-dimensional (3D) slab model of the system, we derive a two-dimensional (2D) sheet model for the thin film in which the optical response is described by a sheet optical conductivity. We develop criteria for the applicability of this 2D sheet model for a layer with an optical thickness far smaller than the wavelength of the light. We examine in detail atomically thin semi-metallic and semiconductor van-der-Waals layers and ultrathin metal films as representative examples. Excellent agreement of the 2D sheet model with the 3D slab model is demonstrated over a broad spectral range from the radio frequency limit to the near ultraviolet. A linearized version of system response for the 2D model is also presented for the case where the influence of the optically thin layer is sufficiently weak. Analytical expressions for the applicability and accuracy of the different optical models are derived, and the appropriateness of the linearized treatment for the materials is considered. We discuss the advantages, as well as limitations, of these models for the purpose of deducing the optical response function of the thin layer from experiment. We generalize the theory to take into account in-plane anisotropy, layered thin film structures, and more general substrates. Implications of the 2D model for the transmission of light by the thin film and for the implementation of half- and totally absorbing layers are discussed.

2D materials↗

Analysis of Nonclassical Fracture Problems for Prestressed Bodies with Interacting Cracks

We consider two types of nonclassical fracture mechanisms, namely, the fracture of cracked bodies with initial (residual) stresses acting along the crack planes and the fracture of solids under compression along parallel cracks. To investigate nonaxisymmetric and axisymmetric problems for infinite solids containing two parallel coaxial cracks or a periodic set of coaxial parallel cracks, we use a combined analytic-numerical method within the framework of the three-dimensional linearized mechanics of solids. The analysis involves the representation of stresses and displacements in the linearized theory via harmonic potential functions. With the use of the integral Fourier–Hankel transformations, the problems are reduced to the solution of Fredholm integral equations of the second kind. This approach allows us to investigate problems in a unified general form for compressible and noncompressible homogeneous isotropic or transversely isotropic elastic bodies with an arbitrary structure of the elastic potential, and the material specification of the model is carried out only in the stage of numerical analysis of the resolving equations obtained in the general form. The effects of initial stresses on the stress intensity factors are analyzed for highly elastic materials and layered composites (modeled as transversely isotropic elastic bodies). The “resonance-like” effects are revealed when compressive initial stresses reach the values corresponding to the local loss of stability of the material in the vicinity of cracks, which, according to the indicated combined method, allows one to determine the critical (limiting) load parameters under the conditions of compression of the body along the cracks. The conclusions concerning the dependences of the stress intensity factors and critical (limiting) parameters of compression on the geometric parameters of the problems are formulated as well as concerning the dependences on physical and mechanical characteristics of the materials.

42 ENGINEERING↗

Bistability analysis in mathematical model of M-phase control in Xenopus oocyte extracts

To gain more insight into the molecular mechanism that regulates eukaryotic cell cycle, we propose a modification to an existing Mathematical model of mitosis-promoting factor control in Xenopus oocyte extract. To achieve this, we redirect how cell exits from mitosis using ubiquitin ligase enzyme APCP:Cdc20 which is acting as a modifier in cyclin degradation. We present a modified wiring map using combination of mass action kinetics and Michaelis–Menten kinetics to represent the biochemical processes. This results in a system of non-linear ordinary differential equations. The equilibria of the system were analyzed analytically. Using MPF as a bifurcation parameter, we analytically derived conditions under which bistability can occur in the MPF activation module. We find that for bistability to occur, it is necessary that the rate of phosphorylation of Wee1 must be greater than the rate of its dephosphorylation and that the concentration of MPF must lie in certain range of parameter values. Numerical simulations conducted show that the time taken to exit the cell cycle is shorter than those published in some reports and the sensitivity analysis results suggest that the rate of synthesis of cyclin is the most sensitive parameter in the dynamics of cell cycle. Through numerical experiments, we found new parameter sets that will give rise to bistability in the MPF activation module.

97 MATHEMATICS AND COMPUTING↗

Bayesian seismic inversion based on rock-physics prior modeling for the joint estimation of acoustic impedance, porosity and lithofacies

We propose a Bayesian approach for seismic inversion to estimate acoustic impedance, porosity and lithofacies within the reservoir conditioned to post-stack seismic and well data. The link between elastic and petrophysical properties is given by a joint prior distribution for the logarithm of impedance and porosity, based on a rock-physics model. The well conditioning is performed through a background model obtained by well log interpolation. Two different approaches are presented: in the first approach, the prior is defined by a single Gaussian distribution, whereas in the second approach it is defined by a Gaussian mixture to represent the well data multimodal distribution and link the Gaussian components to different geological lithofacies. The forward model is based on a linearized convolutional model. For the single Gaussian case, we obtain an analytical expression for the posterior distribution, resulting in a fast algorithm to compute the solution of the inverse problem, i.e. the posterior distribution of acoustic impedance and porosity as well as the facies probability given the observed data. For the Gaussian mixture prior, it is not possible to obtain the distributions analytically, hence we propose a Gibbs algorithm to perform the posterior sampling and obtain several reservoir model realizations, allowing an uncertainty analysis of the estimated properties and lithofacies. Both methodologies are applied to a real seismic dataset with three wells to obtain 3D models of acoustic impedance, porosity and lithofacies. The methodologies are validated through a blind well test and compared to a standard Bayesian inversion approach. Using the probability of the reservoir lithofacies, we also compute a 3D isosurface probability model of the main oil reservoir in the studied field.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Zonally dominated dynamics and Dimits threshold in curvature-driven ITG turbulence

The saturated state of turbulence driven by the ion-temperature-gradient instability is investigated using a two-dimensional long-wavelength fluid model that describes the perturbed electrostatic potential and perturbed ion temperature in a magnetic field with constant curvature (a $Z$ -pinch) and an equilibrium temperature gradient. Numerical simulations reveal a well-defined transition between a finite-amplitude saturated state dominated by strong zonal-flow and zonal temperature perturbations, and a blow-up state that fails to saturate on a box-independent scale. We argue that this transition is equivalent to the Dimits transition from a low-transport to a high-transport state seen in gyrokinetic numerical simulations (Dimits et al. , Phys. Plasmas , vol. 7, 2000, 969). A quasi-static staircase-like structure of the temperature gradient intertwined with zonal flows, which have patch-wise constant shear, emerges near the Dimits threshold. The turbulent heat flux in the low-collisionality near-marginal state is dominated by turbulent bursts, triggered by coherent long-lived structures closely resembling those found in gyrokinetic simulations with imposed equilibrium flow shear (van Wyk et al. , J. Plasma Phys. , vol. 82, 2016, 905820609). The breakup of the low-transport Dimits regime is linked to a competition between the two different sources of poloidal momentum in the system – the Reynolds stress and the advection of the diamagnetic flow by the $\boldsymbol {E}\times \boldsymbol {B}$ flow. By analysing the linear ion-temperature-gradient modes, we obtain a semi-analytic model for the Dimits threshold at large collisionality.

Physics↗

Scalings for radiation from plasma bubbles

In this paper, electron trajectories are studied analytically in the rest frame of a plasma bubble using nonevolving, linear, radial electric and azimuthal magnetic fields in a spherical structure. The electron motion is broken into two distinct periods; one where it orbits around the periphery to the rear of the bubble, and one where it performs oscillations within the bubble interior. By using the first period as an initial condition for the second, general scalings are developed for the x-ray radiation produced by the electron oscillations. The equations are also analyzed to give self-trapping conditions for the electron and to examine the sensitivity of the transverse momentum to small variations from an orbit that is a circular arc. The scalings are in reasonable agreement with recent experiments on x-ray generation and predict a peak spectral brightness of S{approx_equal}6x10{sup 27} photons/s mrad mm 0.1%BW of radiation with a critical energy of 300 MeV using a single stage accelerator driven by a 120 PW laser.

07 ISOTOPE AND RADIATION SOURCES↗

Global gyrokinetic particle simulation of toroidal Alfvén eigenmodes excited by antenna and fast ions

Linear properties of toroidal Alfvén eigenmode (TAE) is studied in global gyrokinetic particle simulations using both fast ion and antenna excitations. A synthetic antenna provides a precise measurement of the Alfvén continuum gap width and the TAE eigenmode frequency, damping rate, and mode structures. The measured gap width exhibits a linear dependence on the aspect ratio, in agreement to a local analytic theory. The TAE frequency and mode structure excited by fast ions show a significant radial symmetry breaking relative to the ideal magnetohydrodynamic theory due to the non-perturbative contributions from the fast ions. The electromagnetic capability of the global gyrokinetic toroidal code (GTC) is verified through these global gyrokinetic simulations of Alfvén eigenmode in cylindrical and toroidal geometries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonlinear relativistic single-electron Thomson scattering power spectrum for incoming laser of arbitrary intensity

The classical nonlinear incoherent Thomson scattering power spectrum from a single relativistic electron with incoming laser radiation of any intensity, investigated numerically by the present authors in a previous publication, displayed both an approximate quadratic behavior in frequency and a redshift of the power spectrum for high intensity incoming radiation. The present work is devoted to justify, in a more general setup, those numerical findings. Those justifications are reinforced by extending suitably analytical approaches, as developed by other authors. Moreover, our analytical treatment exhibits differences between the Doppler-like frequencies for linear and circular polarization of the incoming radiation. Those differences depend nonlinearly on the laser intensity and on the electron initial velocity and do not appear to have been displayed by previous authors. Those Doppler-like frequencies and their differences are validated by new Monte Carlo computations beyond our previuos ones and reported here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Preplasma effects on the generation of high-energy protons in ultraintense laser interaction with foil targets

It is shown that the intense quasistatic electric and magnetic fields self-generated near the axis of the laser-driven channel in an appropriately profiled preplasma during ultraintense laser interaction with a thin target can create dense relativistic electron bunches. The latter easily penetrate through the target and can greatly enhance the sheath field at the rear, resulting in significant increase in the laser-to-ion energy conversion efficiency and the maximum energy of the target normal sheath accelerated ions. Particle-in-cell simulations show that with a hydrogen targets a proton beam of peak energy ∼38 MeV and energy conversion efficiency ≥6.5% can be produced by a linearly polarized 5 × 10{sup 19} W/cm{sup 2} laser. An analytical model is also proposed and its results agree well with those of the simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Exact scattering matrix of graphs in magnetic field and quantum noise

We consider arbitrary quantum wire networks modelled by finite, noncompact, connected quantum graphs in the presence of an external magnetic field. We find a general formula for the total scattering matrix of the network in terms of its local scattering properties and its metric structure. This is applied to a quantum ring with N external edges. Connecting the external edges of the ring to heat reservoirs, we study the quantum transport on the graph in ambient magnetic field. We consider two types of dynamics on the ring: the free Schrödinger and the free massless Dirac equations. For each case, a detailed study of the thermal noise is performed analytically. Interestingly enough, in presence of a magnetic field, the standard linear Johnson-Nyquist law for the low temperature behaviour of the thermal noise becomes nonlinear. The precise regime of validity of this effect is discussed and a typical signature of the underlying dynamics is observed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Low frequency waves and gravitational instability in homogeneous magnetized gyrotropic quantum plasma

In this work, the low frequency waves and gravitational (Jeans) instability of a homogeneous gyrotropic, magnetized, quantum plasma is investigated using the quantum magnetohydrodynamic and Chew-Goldberger-Low fluid models. An analytical dispersion relation for the considered system is obtained solving the linearized perturbations equations employing the Fourier transformation. The onset criterion of the “firehose” instability is retained in parallel propagation, which is unaffected due to the presence of quantum corrections. The gravitational mode modified by the quantum corrections is obtained separately along with the “firehose” mode. In perpendicular propagation, the quantum diffraction term is coupled with the Jeans and Alfven modes whereas in parallel propagation, the Alfven mode does not contribute to the dispersion characteristics as it leads to the “firehose” instability criterion in terms of quantum pressure anisotropy. The stabilizing influences of the quantum diffraction parameter and magnetic field on the growth rates of Jeans instability are examined. It is observed that the growth rate stabilizes much faster in transverse mode due to Alfven stabilization as compared to the longitudinal mode of propagation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗