Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear inverse modeling”

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 361 records · Page 20

Proton and helium cyclotron anisotropy instability thresholds in the magnetosheath

Both the protons and the helium ions of the terrestrial magnetosheath typically display T (sub perpendicular) greater than T (sub parallel), where perpendicular to and parallel to denote directions perpendicular and parallel to the background magnetic field. Observations of the highly compressed magnetosheath show an inverse correlation between these ion temperature anisotropies and the parallel proton beta. Computer simulations have demonstrated that these correlations are due to wave-particle scattering by electromagnetic ion cyclotron anisotropy instabilities. These correlations correspond to linear theory thresholds of the proton cyclotron and the helium cyclotron instabilities. This paper uses linear Vlasov theory and the assumption of a constant maximum growth rate to obtain closed-form expressions for these thresholds as a function of the relative helium density and the parallel proton beta in a parameter model of the magnetosheath.

Gary, S. Peter↗

Less Common Topological Phenomena in Bulk Materials

"Prediction and subsequent discovery of topological insulators is considered to be one of the main results in condensed matter physics in the last decades. Not surprisingly, it has received major attention of both researchers and funding agencies. This attention is well-deserved; yet, one cannot but note that nearly all this research revolves about essentially the same concept: electronic excitations with linear dispersion, covering, of course, such diverse and intriguing phenomena as topological insulator, bulk Dirac states (or Weyl, if not spin degenerate), Mayorana fermions. In this project, we will address, mainly, other topological phenomena, such as topologically nontrivial magnetic patterns (as, for instance, topological Hall and related phenomena). Specifically, we propose three interrelated trusts: (1) Time-reversal symmetry breaking nonrelativistic antiferromagnets, called altermagnets. These are materials that break Kramers degeneracy of electronic bands, despite having zero net magnetization by symmetry and being fully collinear, and not necessarily non-centrosymmetric. The corresponding band structure is very similar to the band structure in non-centrosymmetric spin-orbital materials, but materials that we propose to study are distinctly different, first and foremost in the sense that despite sharing many aspects of their electronic properties with the latter, they break the time-reversal symmetry without either spin-orbit coupling or lack of inversion symmetry. (2) Topologically nontrivial magnetic spirals. The PI has been engaged with the experimental group of Dr. Ghimire at GMU investigating Dirac materials with helical magnetism, based on stacked magnetic Kagome layers, with a generic formula of RMn6Sn6. In particular, Y Mn6Sn6 demonstrates a component of the Hall effect that is naturally interpreted in terms of a topological spin texture, as well as linear magnetoresistance. Our calculation identify Dirac states that are robust with respect to the spiral formation, and let us derive an advance mean-field model explaining the observed phase diagram. This model predicted four distinct phases, with very distinct properties, which have now been seen in neutron experiments. The same compound is known to demonstrate topological Hall effect in a particular magnetic phase, and only at elevated temperature. Based on our understanding of the phase diagram, we have worked out a phenomenological theory of a chiral (skyrmionic) response to an external magnetic field, similar to the nematic response to external strain in Fe-based superconductors, which is possible in a centrosymmetric lattice and without interplanar Dzyaloshinskii-Moriya interaction. This phenomenological theory agrees quantitatively with the experiment. It is in our plans to research other similar materials for this effect. (3) Search for 3D analogues of Fe-based superconductors. We want to investigate materials that can be viewed as 3D analogues of FeSe. Specifically, we want materials that are good metals and host antiferromagnetism, which can be suppressed by pressure and generate an s-wave superconductivity, as in Fe-based superconductors. We have in mind some candidates already. This work will proceed in close collaboration with the experimental group of Prof. Nirmal Ghimire in the same department, whose expertise lies in sample making, magnetometry and transport measurements of materials with complex magnetic structures."

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Radiative equilibrium model of Titan's atmosphere

The present global radiative equilibrium model for the Saturn satellite Titan is restricted to the two-stream approximation, is vertically homogeneous in its scattering properties, and is spectrally divided into one thermal and two solar channels. Between 13 and 33% of the total incident solar radiation is absorbed at the planetary surface, and the 30-60 ratio of violet to thermal IR absorption cross sections in the stratosphere leads to the large temperature inversion observed there. The spectrally integrated mass absorption coefficient at thermal wavelengths is approximately constant throughout the stratosphere, and approximately linear with pressure in the troposphere, implying the presence of a uniformly mixed aerosol in the stratosphere. There also appear to be two regions of enhanced opacity near 30 and 500 mbar.

Samuelson, R. E.↗

Local Excitations of a Charged Nitrogen Vacancy in Diamond with Multireference Density Matrix Embedding Theory

Here, we investigate the negatively charged nitrogen-vacancy center in diamond using periodic density matrix embedding theory (pDMET). To describe the strongly correlated excited states of this system, the complete active space self-consistent field (CASSCF) followed by n-electron valence state second-order perturbation theory (NEVPT2) was used as the impurity solver. Since the NEVPT2-DMET energies show a linear dependence on the inverse of the size of the embedding subspace, we performed an extrapolation of the excitation energies to the nonembedding limit using a linear regression. The extrapolated NEVPT2-DMET first triplet–triplet excitation energy is 2.31 eV and that for the optically inactive singlet–singlet transition is 1.02 eV, both in agreement with the experimentally observed vertical excitation energies of ∼2.18 eV and ∼1.26 eV, respectively. This is the first application of pDMET to a charged periodic system and the first investigation of the NV – defect using NEVPT2 for periodic supercell models.

36 MATERIALS SCIENCE↗

A Method to Solve Interior and Exterior Camera Calibration Parameters for Image Resection

An iterative method is presented to solve the internal and external camera calibration parameters, given model target points and their images from one or more camera locations. The direct linear transform formulation was used to obtain a guess for the iterative method, and herein lies one of the strengths of the present method. In all test cases, the method converged to the correct solution. In general, an overdetermined system of nonlinear equations is solved in the least-squares sense. The iterative method presented is based on Newton-Raphson for solving systems of nonlinear algebraic equations. The Jacobian is analytically derived and the pseudo-inverse of the Jacobian is obtained by singular value decomposition.

Samtaney, Ravi↗

Charge order textures induced by non-linear couplings in a half-doped manganite

Abstract The self-organization of strongly interacting electrons into superlattice structures underlies the properties of many quantum materials. How these electrons arrange within the superlattice dictates what symmetries are broken and what ground states are stabilized. Here we show that cryogenic scanning transmission electron microscopy (cryo-STEM) enables direct mapping of local symmetries and order at the intra-unit-cell level in the model charge-ordered system Nd 1/2 Sr 1/2 MnO 3 . In addition to imaging the prototypical site-centered charge order, we discover the nanoscale coexistence of an exotic intermediate state which mixes site and bond order and breaks inversion symmetry. We further show that nonlinear coupling of distinct lattice modes controls the selection between competing ground states. The results demonstrate the importance of lattice coupling for understanding and manipulating the character of electronic self-organization and that cryo-STEM can reveal local order in strongly correlated systems at the atomic scale.

42 ENGINEERING↗

Analytical and numerical modeling of optical second harmonic generation in anisotropic crystals using #SHAARP package

Electric-dipole optical second harmonic generation (SHG) is a second-order nonlinear process that is widely used as a sensitive probe to detect broken inversion symmetry and local polar order. Analytical modeling of the SHG polarimetry of a nonlinear optical material is essential to extract its point group symmetry and the absolute nonlinear susceptibilities. Current literature on SHG analysis involves numerous approximations and a wide range of (in)accuracies. We have developed an open-source package called the Second Harmonic Analysis of Anisotropic Rotational Polarimetry (#SHAARP.si) which derives analytical and numerical solutions of reflection SHG polarimetry from a single interface (.si) for bulk homogeneous crystals with arbitrary symmetry group, arbitrary crystal orientation, complex and anisotropic linear dielectric tensor with frequency dispersion, a general SHG tensor and arbitrary light polarization. #SHAARP.si enables accurate modeling of polarimetry measurements in reflection geometry from highly absorbing crystals or wedge-shaped transparent crystals. The package is extendable to multiple interfaces.

36 MATERIALS SCIENCE↗

Theoretical and numerical studies of inverse source problem for the linear parabolic equation with sparse boundary measurements

We consider the inverse source problem in the parabolic equation, where the unknown source possesses the semi-discrete formulation. Theoretically, we prove that the flux data from any nonempty open subset of the boundary can uniquely determine the semi-discrete source. This means the observed area can be extremely small, and that is the reason we call it sparse boundary data. For the numerical reconstruction, we formulate the problem from the Bayesian sequential prediction perspective and conduct the numerical examples which estimate the space-time-dependent source state by state. To better demonstrate the method’s performance, we solve two common multiscale problems from two models with a long source sequence. The numerical results illustrate that the inversion is accurate and efficient.

97 MATHEMATICS AND COMPUTING↗

Neural Networks to Find the Optimal Forcing for Offsetting the Anthropogenic Climate Change Effects

Abstract Of great relevance to climate engineering is the systematic relationship between the radiative forcing to the climate system and the response of the system, a relationship often represented by the linear response function (LRF) of the system. However, estimating the LRF often becomes an ill-posed inverse problem due to high-dimensionality and nonunique relationships between the forcing and response. Recent advances in machine learning make it possible to address the ill-posed inverse problem through regularization and sparse system fitting. Here, we develop a convolutional neural network (CNN) for regularized inversion. The CNN is trained using the surface temperature responses from a set of Green’s function perturbation experiments as imagery input data together with data sample densification. The resulting CNN model can infer the forcing pattern responsible for the temperature response from out-of-sample forcing scenarios. This promising proof of concept suggests a possible strategy for estimating the optimal forcing to negate certain undesirable effects of climate change. The limited success of this effort underscores the challenges of solving an inverse problem for a climate system with inherent nonlinearity. Significance Statement Predicting the climate response for a given climate forcing is a direct problem, while inferring the forcing for a given desired climate response is often an inverse, ill-posed, problem, posing a new challenge to the climate community. This study makes the first attempt to infer the radiative forcing for a given target pattern of global surface temperature response using a deep learning approach. The resulting deeply trained convolutional neural network inversion model shows promise in capturing the forcing pattern corresponding to a given surface temperature response, with a significant implication on the design of an optimal solar radiation management strategy for curbing global warming. This study also highlights the technical challenges that future research should prioritize in seeking feasible solutions to the inverse climate problem.

Ren, Huiying↗

Laser ultrasonic imaging of subsurface defects with the linear sampling method

Laser ultrasonics is a remote nondestructive evaluation technique suitable for real-time monitoring of fabrication processes in semiconductor metrology, advanced manufacturing, and other applications where non-contact, high fidelity measurements are required. Here we investigate laser ultrasonic data processing approaches to reconstruct images of subsurface side drilled holes in aluminum alloy specimens. We demonstrate through simulation that the model-based linear sampling method (LSM) can perform accurate shape reconstruction of single and multiple holes and produce images with well-defined boundaries. We experimentally confirm that LSM produces images that represent the internal geometric features of an object, some of which may be missed by conventional imaging.

42 ENGINEERING↗

The inverse problem of constructing a gravimetric geoid

Computation of a single geoidal height from gravity acceleration data formally requires that the latter be known everywhere on the earth. A computational procedure based on linear inverse theory for estimating geoidal heights from incomplete sets of data is presented. The same scheme can be used to estimate gravity accelerations from altimetry-derived geoids. The systematic error owing to lack of data and the choice of a particular inverse operator is described by using resolution functions and their spherical harmonic expansions. An rms value of this error is also estimated by assuming a spectrum for the unknown geoid. The influence of the size of the data region, the spacing between data, the filtering applied to the data, and the model weighting function chosen are all quantified in a spherical geometry. The examples presented show that when low degree spherical harmonic coefficients are available - from satellite orbit analysis - a band-passed version of the geoid can be constructed from local gravity data, even with a relatively restricted data set.

Zlotnicki, V.↗

Estimation of optimal Kalman filter gain from non-optimal filter residuals

This paper presents a novel method of estimating the optimal steady state Kalman filter gain of a linear discrete time-invariant system from a non-optimal Kalman filter residual sequence. The relation between the optimal residual sequence and a signal derived from the non-optimal residual sequence is described by a Moving Average (MA) model whose coefficients are expressed in terms of the state space parameters and the optimal steady state Kalman filter gain. In order to identify the MA model, a whitening filter of the derived signal, which corresponds to an AutoRegressive (AR) model of the signal, is first identified using the least-squares method. Then the inverse filter of the whitening filter, which corresponds to the MA model, is calculated. From the coefficients of the identified MA model, the optimal steady state Kalman filter gain can be obtained. Numerical example is provided to illustrate the feasibility of this approach.

Chen, Chung-Wen↗

A full potential inverse method based on a density linearization scheme for wing design

A mixed analysis inverse procedure based on the full potential equation in conservation form was developed to recontour a given base wing to produce density linearization scheme in applying the pressure boundary condition in terms of the velocity potential. The FL030 finite volume analysis code was modified to include the inverse option. The new surface shape information, associated with the modified pressure boundary condition, is calculated at a constant span station based on a mass flux integration. The inverse method is shown to recover the original shape when the analysis pressure is not altered. Inverse calculations for weakening of a strong shock system and for a laminar flow control (LFC) pressure distribution are presented. Two methods for a trailing edge closure model are proposed for further study.

Shankar, V.↗

Bending of Light in Modified Gravity at Large Distances

We discuss the bending of light in a recent model for gravity at large distances containing a Rindler type acceleration proposed by Grumiller. We consider the static, spherically symmetric metric with cosmological constant and Rindler-like term 2ar presented in this model, and we use the procedure by Rindler and Ishak. to obtain the bending angle of light in this metric. Earlier work on light bending in this model by Carloni, Grumiller, and Preis, using the method normally employed for asymptotically flat space-times, led to a conflicting result (caused by the Rindler-like term in the metric) of a bending angle that increases with the distance of closest approach r(sub 0) of the light ray from the centrally concentrated spherically symmetric matter distribution. However, when using the alternative approach for light bending in nonasymptotically flat space-times, we show that the linear Rindler-like term produces a small correction to the general relativistic result that is inversely proportional to r(sub 0). This will in turn affect the bounds on Rindler acceleration obtained earlier from light bending and casts doubts on the nature of the linear term 2ar in the metric

Grumiller↗

Eruption of a Multiple-Turn Helical Magnetic Flux Tube in a Large Flare: Evidence for External and Internal Reconnection that Fits the Breakout Model of Solar Magnetic Eruptions

We present observations and an interpretation of a unique multiple-turn spiral flux tube eruption from active region 10030 on 2002 July 15. The TRACE C IV observations clearly show a flux tube that is helical and erupting from within a sheared magnetic field. These observations are interpreted in the context of the breakout model for magnetic field explosions. The initiation of the helix eruption. as determined by a linear backward extrapolation, starts 25 s after the peak of the flare's strongest impulsive spike of microwave gyrosynchrotron radiation early in the flare s explosive phase, implying that the sheared core field is not the site of the initial reconnection. Within the quadrupolar configuration of the active region, the external and internal reconnection sites are identified in each of two consecutive eruptive flares that produce a double coronal mass ejection (CME). The first external breakout reconnection apparently releases an underlying sheared core field and allows it to erupt, leading to internal reconnection in the wake of the erupting helix. This internal reconnection releases the helix and heats the two-ribbon flare. These events lead to the first CME and are followed by a second breakout that initiates a second and larger halo CME. The strong magnetic shear in the region is compatible with the observed rapid proper motion and evolution of the active region. The multiple-turn helix originates from above a sheared-field magnetic inversion line within a filament channel. and starts to erupt only after fast breakout reconnection has started. These observations are counter to the standard flare model and support the breakout model for eruptive flare initiation.

Gary, G. Allen↗

Safety First: Stability and Dissipation of Line-tied Force-free Flux Tubes in Magnetized Coronae

Abstract Magnetized plasma columns and extended magnetic structures with both footpoints anchored to a surface layer are an important building block of astrophysical dissipation models. Current loops shining in X-rays during the growth of plasma instabilities are observed in the corona of the Sun and are expected to exist in highly magnetized neutron star magnetospheres and accretion disk coronae. For varying twist and system sizes, we investigate the stability of line-tied force-free flux tubes and the dissipation of twist energy during instabilities using linear analysis and time-dependent force-free electrodynamics simulations. Kink modes ( m = 1) and efficient magnetic energy dissipation develop for plasma safety factors q ≲ 1, where q is the inverse of the number of magnetic field line windings per column length. Higher-order fluting modes ( m > 1) can distort equilibrium flux tubes for q > 1 but induce significantly less dissipation. In our analysis, the characteristic pitch μ ˜ 0 of flux-tube field lines determines the growth rate ( ∝ μ ˜ 0 3 ) and minimum wavelength of the kink instability ( ∝ μ ˜ 0 − 1 ). We use these scalings to determine a minimum flux tube length for the growth of the kink instability for any given μ ˜ 0 . By drawing analogies to idealized magnetar magnetospheres with varying regimes of boundary shearing rates, we discuss the expected impact of the pitch-dependent growth rates for magnetospheric dissipation in magnetar conditions.

79 ASTRONOMY AND ASTROPHYSICS↗

Unsteady aerodynamic modeling and active aeroelastic control

Unsteady aerodynamic modeling techniques are developed and applied to the study of active control of elastic vehicles. The problem of active control of a supercritical flutter mode poses a definite design goal stability, and is treated in detail. The transfer functions relating the arbitrary airfoil motions to the airloads are derived from the Laplace transforms of the linearized airload expressions for incompressible two dimensional flow. The transfer function relating the motions to the circulatory part of these loads is recognized as the Theodorsen function extended to complex values of reduced frequency, and is termed the generalized Theodorsen function. Inversion of the Laplace transforms yields exact transient airloads and airfoil motions. Exact root loci of aeroelastic modes are calculated, providing quantitative information regarding subcritical and supercritical flutter conditions.

Edwards, J. W.↗

Mass and energy flow near sunspots

Sunspots block the flow of energy to the solar surface. The blocked energy heats the volume beneath the spot, producing a pressure excess which drives an outflow of mass. Linear numerical models of the mass and energy flow around spots were constructed to estimate the predictions of this physical picture against the observed properties of sunspot bright rings and moat flows. The width of the bright ring and moat are predicted to be proportional to the depth of the spot penumbra, in conflict with the observed proportionality of the moat width to the spot diameter. Postulating that spot depths are proportional to spot diameters would bury the moat flow too deeply to be observed, because the radial velocity at the surface is found to be inversely proportional to the depth of the spot penumbra. The radial velocity at the surface is of order a few hundred meters per second after 1 day, in agreement with the observed excess of moat velocities over supergranule velocities.

Nye, Alan↗