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 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

TH-A-18C-07: Noise Suppression in Material Decomposition for Dual-Energy CT

Purpose: A general problem of dual-energy CT (DECT) is that the decomposition is sensitive to noise in the two sets of dual-energy projection data, resulting in severely degraded qualities of decomposed images. We have previously proposed an iterative denoising method for DECT. Using a linear decomposition function, the method does not gain the full benefits of DECT on beam-hardening correction. In this work, we expand the framework of our iterative method to include non-linear decomposition models for noise suppression in DECT. Methods: We first obtain decomposed projections, which are free of beam-hardening artifacts, using a lookup table pre-measured on a calibration phantom. First-pass material images with high noise are reconstructed from the decomposed projections using standard filter-backprojection reconstruction. Noise on the decomposed images is then suppressed by an iterative method, which is formulated in the form of least-square estimation with smoothness regularization. Based on the design principles of a best linear unbiased estimator, we include the inverse of the estimated variance-covariance matrix of the decomposed images as the penalty weight in the least-square term. Analytical formulae are derived to compute the variance-covariance matrix from the measured decomposition lookup table. Results: We have evaluated the proposed method via phantom studies. Using non-linear decomposition, our method effectively suppresses the streaking artifacts of beam-hardening and obtains more uniform images than our previous approach based on a linear model. The proposed method reduces the average noise standard deviation of two basis materials by one order of magnitude without sacrificing the spatial resolution. Conclusion: We propose a general framework of iterative denoising for material decomposition of DECT. Preliminary phantom studies have shown the proposed method improves the image uniformity and reduces noise level without resolution loss. In the future, we will perform more phantom studies to further validate the performance of the purposed method. This work is supported by a Varian MRA grant.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2017

Late-time growth rate, mixing, and anisotropy in the multimode narrowband Richtmyer–Meshkov instability: The θ-group collaboration

Turbulent Richtmyer–Meshkov instability (RMI) is investigated through a series of high resolution three-dimensional simulations of two initial conditions with eight independent codes. The simulations are initialised with a narrowband perturbation such that instability growth is due to non-linear coupling/backscatter from the energetic modes, thus generating the lowest expected growth rate from a pure RMI. By independently assessing the results from each algorithm and computing ensemble averages of multiple algorithms, the results allow a quantification of key flow properties as well as the uncertainty due to differing numerical approaches. A new analytical model predicting the initial layer growth for a multimode narrowband perturbation is presented, along with two models for the linear and non-linear regimes combined. Overall, the growth rate exponent is determined as θ=0.292±0.009, in good agreement with prior studies; however, the exponent is decaying slowly in time. Also, θ is shown to be relatively insensitive to the choice of mixing layer width measurements. The asymptotic integral molecular mixing measures Θ=0.792±0.014, Ξ=0.800±0.014, and Ψ=0.782±0.013 are lower than some experimental measurements but within the range of prior numerical studies. The flow field is shown to be persistently anisotropic for all algorithms, at the latest time having between 49% and 66% higher kinetic energy in the shock parallel direction compared to perpendicular and does not show any return to isotropy. The plane averaged volume fraction profiles at different time instants collapse reasonably well when scaled by the integral width, implying that the layer can be described by a single length scale and thus a single θ. In conclusion, quantitative data given for both ensemble averages and individual algorithms provide useful benchmark results for future research.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Statistical Significance Testing for Mixed Priors: A Combined Bayesian and Frequentist Analysis

In many hypothesis testing applications, we have mixed priors, with well-motivated informative priors for some parameters but not for others. The Bayesian methodology uses the Bayes factor and is helpful for the informative priors, as it incorporates Occam’s razor via the multiplicity or trials factor in the look-elsewhere effect. However, if the prior is not known completely, the frequentist hypothesis test via the false-positive rate is a better approach, as it is less sensitive to the prior choice. We argue that when only partial prior information is available, it is best to combine the two methodologies by using the Bayes factor as a test statistic in the frequentist analysis. We show that the standard frequentist maximum likelihood-ratio test statistic corresponds to the Bayes factor with a non-informative Jeffrey’s prior. We also show that mixed priors increase the statistical power in frequentist analyses over the maximum likelihood test statistic. We develop an analytic formalism that does not require expensive simulations and generalize Wilks’ theorem beyond its usual regime of validity. In specific limits, the formalism reproduces existing expressions, such as the p-value of linear models and periodograms. We apply the formalism to an example of exoplanet transits, where multiplicity can be more than 10 7 . We show that our analytic expressions reproduce the $p$-values derived from numerical simulations. We offer an interpretation of our formalism based on the statistical mechanics. We introduce the counting of states in a continuous parameter space using the uncertainty volume as the quantum of the state. We show that both the $p$-value and Bayes factor can be expressed as an energy versus entropy competition.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

Parallel computation using boundary elements in solid mechanics

The inherent parallelism of the boundary element method is shown. The boundary element is formulated by assuming the linear variation of displacements and tractions within a line element. Moreover, MACSYMA symbolic program is employed to obtain the analytical results for influence coefficients. Three computational components are parallelized in this method to show the speedup and efficiency in computation. The global coefficient matrix is first formed concurrently. Then, the parallel Gaussian elimination solution scheme is applied to solve the resulting system of equations. Finally, and more importantly, the domain solutions of a given boundary value problem are calculated simultaneously. The linear speedups and high efficiencies are shown for solving a demonstrated problem on Sequent Symmetry S81 parallel computing system.

Chien, L. S.↗
DOE OSTI2021

Advances in constraining intrinsic alignment models with hydrodynamic simulations

We use galaxies from the illustristng, massiveblack-ii, and illustris-1 hydrodynamic simulations to investigate the behaviour of large scale galaxy intrinsic alignments. Our analysis spans four redshift slices over the approximate range of contemporary lensing surveys z = 0–1. We construct comparable weighted samples from the three simulations, which we then analyse using an alignment model that includes both linear and quadratic alignment contributions. Our data vector includes galaxy–galaxy, galaxy–shape, and shape–shape projected correlations, with the joint covariance matrix estimated analytically. In all of the simulations, we report non-zero IAs at the level of several σ. For a fixed lower mass threshold, we find a relatively strong redshift dependence in all three simulations, with the linear IA amplitude increasing by a factor of ~2 between redshifts z = 0 and z = 1. We report no significant evidence for non-zero values of the tidal torquing amplitude, A 2 , in TNG, above statistical uncertainties, although MBII favours a moderately negative A 2 ~ –2. Examining the properties of the TATT model as a function of colour, luminosity and galaxy type (satellite or central), our findings are consistent with the most recent measurements on real data. We also outline a novel method for constraining the TATT model parameters directly from the pixelized tidal field, alongside a proof-of-concept exercise using TNG. This technique is shown to be promising, although comparison with previous results obtained via other methods is non-trivial.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Ka-band linearizer for the Ultra-Compact X-ray free-electron laser at UCLA

There is a strong demand for accelerating structures able to achieve higher gradients and more compact dimensions for the next generation of linear accelerators for research, industrial and medical applications. Notably innovative technologies will permit compact and affordable advanced accelerators as the linear collider and X-ray free-electron lasers (XFELs) with accelerating gradients over twice the value achieved with current technologies. In particular XFELs are able to produce coherent X-ray pulses with peak brightness 10 orders of magnitude greater than preceding approaches, which has revolutionized numerous research fields through imaging of the nanoscopic world at the time and length scale of atom-based systems, that is of femtosecond and Angstrom. There is a strong interest for combining these two fields, to form a proper tool with the goal of producing a very compact XFEL in order to investigate multi-disciplinary topics in chemistry, biology, materials science, medicine and physics. In the framework of the Ultra-Compact XFEL project under study at the University of California, Los Angeles, a high gradient radio-frequency accelerating structure for the longitudinal phase-space linearization with an integrated voltage of at least 15 MV working on 6th harmonic of the main Linac frequency is required. We here present the electromagnetic design of a cryogenic normal-conducting 8 cm long Ka-band standing-wave linearizer working on π mode with a target accelerating gradient beyond 100 MV/m. The studies have been performed analytically and numerically to investigate the beam dynamics and electromagnetic issues.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

Spin drift in highly doped n-type Si

A quantitative estimation of spin drift velocity in highly doped n-type silicon (Si) at 8 K is presented in this letter. A local two-terminal Hanle measurement enables the detection of a modulation of spin signals from the Si as a function of an external electric field, and this modulation is analyzed by using a spin drift-diffusion equation and an analytical solution of the Hanle-type spin precession. The analyses reveal that the spin drift velocity is linearly proportional to the electric field. The contribution of the spin drift effect to the spin signals is crosschecked by introducing a modified nonlocal four-terminal method.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2011

Exact analytical solutions to the master equation of quantum Brownian motion for a general environment

Research Highlights: > We study the model of a quantum oscillator linearly coupled to a bath of oscillators. > We derive the master equation and solutions for general spectra and temperatures. > We generalize to cases with an external force and arbitrary number of oscillators. > Other derivations have incorrect diffusion and force response for nonlocal damping. > We give exact results for ohmic, sub-ohmic and supra-ohmic environments. - Abstract: We revisit the model of a quantum Brownian oscillator linearly coupled to an environment of quantum oscillators at finite temperature. By introducing a compact and particularly well-suited formulation, we give a rather quick and direct derivation of the master equation and its solutions for general spectral functions and arbitrary temperatures. The flexibility of our approach allows for an immediate generalization to cases with an external force and with an arbitrary number of Brownian oscillators. More importantly, we point out an important mathematical subtlety concerning boundary-value problems for integro-differential equations which led to incorrect master equation coefficients and impacts on the description of nonlocal dissipation effects in all earlier derivations. Furthermore, we provide explicit, exact analytical results for the master equation coefficients and its solutions in a wide variety of cases, including ohmic, sub-ohmic and supra-ohmic environments with a finite cut-off.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1992

Linear stability theory and three-dimensional boundary layer transition

The viewgraphs and discussion of linear stability theory and three dimensional boundary layer transition are provided. The ability to predict, using analytical tools, the location of boundary layer transition over aircraft-type configurations is of great importance to designers interested in laminar flow control (LFC). The e(sup N) method has proven to be fairly effective in predicting, in a consistent manner, the location of the onset of transition for simple geometries in low disturbance environments. This method provides a correlation between the most amplified single normal mode and the experimental location of the onset of transition. Studies indicate that values of N between 8 and 10 correlate well with the onset of transition. For most previous calculations, the mean flows were restricted to two-dimensional or axisymmetric cases, or have employed simple three-dimensional mean flows (e.g., rotating disk, infinite swept wing, or tapered swept wing with straight isobars). Unfortunately, for flows over general wing configurations, and for nearly all flows over fuselage-type bodies at incidence, the analysis of fully three-dimensional flow fields is required. Results obtained for the linear stability of fully three-dimensional boundary layers formed over both wing and fuselage-type geometries, and for both high and low speed flows are discussed. When possible, transition estimates form the e(sup N) method are compared to experimentally determined locations. The stability calculations are made using a modified version of the linear stability code COSAL. Mean flows were computed using both Navier Stokes and boundary-layer codes.

Spall, Robert E.↗
DOE OSTI2014

Closure of the single fluid magnetohydrodynamic equations in presence of electron cyclotron current drive

In the presence of electron cyclotron current drive (ECCD), the Ohm's law of single fluid magnetohydrodynamics is modified as E + v × B = η(J – J{sub EC}). This paper presents a new closure relation for the EC driven current density appearing in this modified Ohm's law. The new relation faithfully represents the nonlocal character of the EC driven current and its main origin in the Fisch-Boozer effect. The closure relation is validated on both an analytical solution of an approximated Fokker-Planck equation as well as on full bounce-averaged, quasi-linear Fokker-Planck code simulations of ECCD inside rotating magnetic islands. The new model contains the model put forward by Giruzzi et al. [Nucl. Fusion 39, 107 (1999)] in one of its limits.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

A scalable framework for efficient coupling of thermal and microstructural simulations in additive manufacturing

Predicting microstructure evolution in metal additive manufacturing (AM) is important for process optimization, but spatiotemporal scale disparities between thermal transport and microstructure evolution create significant challenges for efficient data transfer between simulation codes. To address this, we present Stork, a scalable framework for coupling thermal and microstructural simulations. Stork uses a sparse data representation to identify and store active solidification sub-volumes, enabling highly parallel quad-linear interpolation from coarse thermal grids to fine microstructure grids without large intermediate storage. We demonstrate the framework by coupling the semi-analytic heat transfer code 3DThesis with the time-parallel cellular automata code Toucan. This approach achieves over two orders of magnitude reduction in data generation time and file size compared to prior workflows. Numerical studies show that quad-linear interpolation preserves grain morphology and crystallographic texture in laser powder bed fusion (LPBF) simulations for coarsening ratios up to 16. Overall, Stork provides a scalable pathway for high-throughput, component-scale AM simulations on modern high-performance computing systems.

36 MATERIALS SCIENCE↗
DOE OSTI2013

VECTOR TOMOGRAPHY FOR THE CORONAL MAGNETIC FIELD. II. HANLE EFFECT MEASUREMENTS

In this paper, we investigate the feasibility of saturated coronal Hanle effect vector tomography or the application of vector tomographic inversion techniques to reconstruct the three-dimensional magnetic field configuration of the solar corona using linear polarization measurements of coronal emission lines. We applied Hanle effect vector tomographic inversion to artificial data produced from analytical coronal magnetic field models with equatorial and meridional currents and global coronal magnetic field models constructed by extrapolation of real photospheric magnetic field measurements. We tested tomographic inversion with only Stokes Q, U, electron density, and temperature inputs to simulate observations over large limb distances where the Stokes I parameters are difficult to obtain with ground-based coronagraphs. We synthesized the coronal linear polarization maps by inputting realistic noise appropriate for ground-based observations over a period of two weeks into the inversion algorithm. We found that our Hanle effect vector tomographic inversion can partially recover the coronal field with a poloidal field configuration, but that it is insensitive to a corona with a toroidal field. This result demonstrates that Hanle effect vector tomography is an effective tool for studying the solar corona and that it is complementary to Zeeman effect vector tomography for the reconstruction of the coronal magnetic field.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1974

Mission Analysis Program for Solar Electric Propulsion (MAPSEP). Volume 1: Analytical manual

The mission analysis program for solar electric propulsion (MAPSEP) is comprised of the basic modes: TOPSEP (trajectory generation), GODSEP (linear error analysis), and SIMSEP (simulation). The program is designed to analyze any low thrust mission with respect to trajectory performance, guidance and navigation, and to provide system related requirements for the purpose of vehicle design. The MAPSEP organization is described along with all models and algorithms. Topics discussed include: trajectory and error covariance propagation methods, orbit determination processes, thrust modeling, and trajectory correction (guidance) schemes.

Hong, P. E.↗
NASA NTRS1977

The linear Fresnel lens - Solar optical analysis of tracking error effects

Real sun-tracking solar concentrators imperfectly follow the solar disk, operationally sustaining both transverse and axial misalignments. This paper describes an analysis of the solar concentration performance of a line-focusing flat-base Fresnel lens in the presence of small transverse tracking errors. Simple optics and ray-tracing techniques are used to evaluate the lens solar transmittance and focal-plane imaging characteristics. Computer-generated example data for an f/1.0 lens indicate that less than a 1% transmittance degradation occurs for transverse errors up to 2.5 deg. In this range, solar-image profiles shift laterally in the focal plane, the peak concentration ratio drops, and profile asymmetry increases with tracking error. With profile shift as the primary factor, the ninety-percent target-intercept width increases rapidly for small misalignments, e.g., almost threefold for a 1-deg error. The analytical model and computational results provide a design base for tracking and absorber systems for the linear-Fresnel-lens solar concentrator.

Cosby, R. M.↗
NASA NTRS1972

The prediction of nonlinear longitudinal combustion instability in liquid propellant rockets

An analytical technique was developed to solve nonlinear longitudinal combustion instability problems. The analysis yields the transient and limit cycle behavior of unstable motors and the threshold amplitude required to trigger a linearly stable motor into unstable operation. The limit cycle waveforms were found to exhibit shock wave characteristics for most unstable engine operating conditions. A method of correlating the analytical solutions with experimental data was developed. Calculated results indicate that a second-order solution adequately describes the behavior of combustion instability oscillations over a broad range of engine operating conditions, but that higher order effects must be accounted for in order to investigate engine triggering.

Lores, M. E.↗
DOE OSTI2017

Dose distribution characteristic study of the Rotating Gamma Knife (RGK) by using the geometry analytic method

The purpose of this study was to realize the processing of dose distribution of RGK at the treatment isocenter at any gantry rotational angle by using an analytic geometry method to avoid inadequate arc therapy angles when implementing the treatment plan. Gaf chromic film was used for dose evaluation. A calibration curve was first obtained using linear accelerator irradiation. The 50% dose relative to the central axis at fixed gantry angles at the x-, y- and z-axes was obtained using Gaf chromic film and was compared to the analytic geometry method. The full width half maximum (FWHM) on the x-, y- and z-axes was predicted for the RGK dose distribution characteristic analysis. The FWHM on the x-, y-, and z-axes varies with different gantry and rotational plate angles. The most dramatic intersection variation appeared at a static gantry angle of 25°. The ratio of the FWHM of the y- and z-axes to that of the x-axis was up to 9 and 10. The geometric analytic method can be used for an accurate analysis of dose distribution in RGK replacing the actual film exposure experiment. It is essential to select the best arc irradiation angle to prescribe the dose to avoid excess irradiation of normal tissue. The geometric method used in this study can also be applied for rotational arc therapy dose analysis such as tomotherapy, linear-based stereotactic radiotherapy, or volume matrices arc therapy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2022

Model‐Based Approach to Selecting Pegfilgrastim Dose for Pharmacokinetic and Pharmacodynamic Similarity Studies in Biosimilar Development

This study applied modeling and simulation (M&S) approaches to evaluate the sensitivity of pegfilgrastim pharmacokinetics (PKs) and pharmacodynamics (PDs) to changes in dose amount, and linear or nonlinear clearance (CL) over pegfilgrastim subcutaneous dose of 2–6 mg. A previously published model was adapted to better describe pegfilgrastim PK and PD data in healthy subjects and used in simulation. Nonlinear CL accounts for 98% and 77%, respectively, of the total CL of pegfilgrastim at 2 and 6 mg. The sensitivity analyses showed: (i) PK of 2 and 6 mg doses are similarly sensitive to detect differences for a 5% change in dose; (ii) PK of 2 mg dose is more sensitive to changes in receptor binding affinity, a model parameter for nonlinear CL, and a product quality attribute characterized with orthogonal methods as part of demonstrating analytical similarity between products; (iii) PK of approved 6 mg dose is more sensitive to changes in linear CL, which has not been associated with any specific product quality attributes, and (iv) the PDs are not sensitive to changes in linear or nonlinear CL. Taken together, our analyses support that the approved pegfilgrastim dose of 6 mg is appropriate for detecting differences between a biosimilar and the reference products in pegfilgrastim PK and PD similarity studies. The described M&S approaches can be adopted to support dose selection for biosimilars with nonlinear PK and complex PK–PD interplay.

Pharmacology & Pharmacy↗
NASA NTRS2000

Nonlinearity and Strain-Rate Dependence in the Deformation Response of Polymer Matrix Composites Modeled

There has been no accurate procedure for modeling the high-speed impact of composite materials, but such an analytical capability will be required in designing reliable lightweight engine-containment systems. The majority of the models in use assume a linear elastic material response that does not vary with strain rate. However, for containment systems, polymer matrix composites incorporating ductile polymers are likely to be used. For such a material, the deformation response is likely to be nonlinear and to vary with strain rate. An analytical model has been developed at the NASA Glenn Research Center at Lewis Field that incorporates both of these features. A set of constitutive equations that was originally developed to analyze the viscoplastic deformation of metals (Ramaswamy-Stouffer equations) was modified to simulate the nonlinear, rate-dependent deformation of polymers. Specifically, the effects of hydrostatic stresses on the inelastic response, which can be significant in polymers, were accounted for by a modification of the definition of the effective stress. The constitutive equations were then incorporated into a composite micromechanics model based on the mechanics of materials theory. This theory predicts the deformation response of a composite material from the properties and behavior of the individual constituents. In this manner, the nonlinear, rate-dependent deformation response of a polymer matrix composite can be predicted.

Goldberg, Robert K.↗