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 163 records · Page 9

DIRECT INTEGRATION OF THE COLLISIONLESS BOLTZMANN EQUATION IN SIX-DIMENSIONAL PHASE SPACE: SELF-GRAVITATING SYSTEMS

We present a scheme for numerical simulations of collisionless self-gravitating systems which directly integrates the Vlasov-Poisson equations in six-dimensional phase space. Using the results from a suite of large-scale numerical simulations, we demonstrate that the present scheme can simulate collisionless self-gravitating systems properly. The integration scheme is based on the positive flux conservation method recently developed in plasma physics. We test the accuracy of our code by performing several test calculations, including the stability of King spheres, the gravitational instability, and the Landau damping. We show that the mass and the energy are accurately conserved for all the test cases we study. The results are in good agreement with linear theory predictions and/or analytic solutions. The distribution function keeps the property of positivity and remains non-oscillatory. The largest simulations are run on 64{sup 6} grids. The computation speed scales well with the number of processors, and thus our code performs efficiently on massively parallel supercomputers.

79 ASTRONOMY AND ASTROPHYSICS↗

Effects of non-linearities on magnetic field generation

Magnetic fields are present on all scales in the Universe. While we understand the processes which amplify the fields fairly well, we do not have a ''natural'' mechanism to generate the small initial seed fields. By using fully relativistic cosmological perturbation theory and going beyond the usual confines of linear theory we show analytically how magnetic fields are generated. This is the first analytical calculation of the magnetic field at second order, using gauge-invariant cosmological perturbation theory, and including all the source terms. To this end, we have rederived the full set of governing equations independently. Our results suggest that magnetic fields of the order of 10{sup -30}- 10{sup -27} G can be generated (although this depends on the small scale cut-off of the integral), which is largely in agreement with previous results that relied upon numerical calculations. These fields are likely too small to act as the primordial seed fields for dynamo mechanisms.

79 ASTRONOMY AND ASTROPHYSICS↗

Elastically mediated interactions between grain boundaries and precipitates in two-phase coherent solids

In this work, we investigate analytically and numerically the interaction between grain boundaries and second phase precipitates in two-phase coherent solids in the presence of misfit strain. Our numerical study uses amplitude equations that describe the interaction of composition and stress [R. Spatschek and A. Karma, Phys. Rev. B 81, 214201 (2010)] and free-energies corresponding to two-dimensional hexagonal and three-dimensional BCC crystal structures that exhibit isotropic and anisotropic elastic properties, respectively. We consider two experimentally motivated geometries where (i) a lamellar precipitate nucleates along a planar grain boundary that is centered inside the precipitate, and (ii) a circular precipitate nucleates inside a grain at a finite distance to an initially planar grain boundary. For the first geometry, we find that the grain boundary becomes morphologically unstable due to the combination of long-range elastic interaction between the grain boundary and compositional domain boundaries, and shear-coupled grain boundary motion. We characterize this instability analytically by extending the linear stability analysis carried out recently [P.-A. Geslin, Y.-C. Xu, and A. Karma, Phys. Rev. Lett. 114, 105501 (2015)] to the more general case of elastic anisotropy. The analysis predicts that elastic anisotropy hinders but does not suppress the instability. Simulations also reveal that, in a well-developed nonlinear regime, this instability can lead to the break-up of low-angle grain boundaries when the misfit strain exceeds a threshold that depends on the grain boundary misorientation. For the second geometry, simulations show that the elastic interaction between an initially planar grain boundary and an adjacent circular precipitate causes the precipitate to migrate to and anchor at the grain boundary.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

From noise to information: The transfer function formalism for uncertainty quantification in reconstructing the nuclear density

The neutron distribution of neutron-rich nuclei provides critical information on the structure of finite nuclei and neutron stars. Parity violating experiments—such as PREX and CREX—provide a clean and largely model-independent determination of neutron densities. Such experiments, however, are challenging and expensive, which is why sound statistical arguments are required to maximize the information gained. We introduce a new framework, the “transfer function formalism,” aimed at uncertainty quantification, model selection, and experimental design in the context of neutron densities. The transfer functions (TFs) are built analytically by expressing the linear response of the objective function (e.g., χ2) to small perturbations of the data. Using the TF formalism, we are able to analyze the expected overall uncertainty—quantified in terms of bias and variance—of the mean square radius and interior density of 48 Ca and 208 Pb. Using relativistic mean field models as a proxy for the weak-charge density—and assuming that a total of five measurements could be performed on the weak form factor of 48 Ca and 208 Pb—we identify the optimal models and experimental locations that minimize the uncertainty in the extraction of the radius and interior density. We also explore the use of the TF formalism to understand the influence of prior distributions for the model parameters, as well as the optimization of model hyperparameters not constrained by the data. Here,we establish how the choice of experimental locations and the model that is used can have a significant impact on the final uncertainties of the extracted quantities of interest. For challenging experiments such as CREX and PREX, a proper quantification of such uncertainties is critical. We have demonstrated how the TF formalism provides several advantages for this type of analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Analytical solutions for transient and steady state beam loading in arbitrary traveling wave accelerating structures

Analytical solutions are derived for both transient and steady state gradient distributions in the traveling wave (TW) accelerating structures with arbitrary variation of parameters over the structure length. The results of the unloaded and beam loaded cases are presented. Finally, the exact analytical shape of the rf pulse waveform was found in order to apply the transient beam loading compensation scheme during the structure filling time. The obtained theoretical formulas were cross-checked by direct numerical simulations on the CLIC main linac accelerating structure and demonstrated a good agreement. The proposed methods provide a fast and reliable tool for the initial stage of the TW structure analysis.

43 PARTICLE ACCELERATORS↗

Impact of variable RBE on proton fractionation

Purpose: To explore the impact of variable proton relative biological effectiveness (RBE) on dose fractionation for clinically relevant situations. A generic RBE = 1.1 is generally used for isoeffect calculations, while experimental studies showed that proton RBE varies with tissue type, dose, and linear energy transfer (LET). Methods: An analytical expression for the LET and {alpha}/{beta} dependence of the linear-quadratic (LQ) model has been used for proton simulations in parallel with the assumption of a generic RBE = 1.1. Calculations have been performed for ranges of LET values and fractionation sensitivities to describe clinically relevant cases, such as the treatment of head and neck and prostate tumors. Isoeffect calculations were compared with predictions from a generic RBE value and reported clinical results. Results: The generic RBE = 1.1 appears to be a reasonable estimate for the proton RBE of rapidly growing tissues irradiated with low LET radiation. However, the use of a variable RBE predicts larger differences for tissues with low {alpha}/{beta} (both tumor and normal) and at low doses per fraction. In some situations these differences may appear in contrast to the findings from photon studies highlighting the importance of accurate accounting for the radiobiological effectiveness of protons. Furthermore, the use of variable RBE leads to closer predictions to clinical results. Conclusions: The LET dependence of the RBE has a strong impact on the predicted effectiveness of fractionated proton radiotherapy. The magnitude of the effect is modulated by the fractionation sensitivity and the fractional dose indicating the need for accurate analyses both in the target and around it. Care should therefore be employed for changing clinical fractionation patterns or when analyzing results from clinical studies for this type of radiation.

60 APPLIED LIFE SCIENCES↗

SELF-GRAVITATIONAL FORCE CALCULATION OF INFINITESIMALLY THIN GASEOUS DISKS ON NESTED GRIDS

We extend the work of Yen et al. and develop second-order formulae to accommodate a nested grid discretization for the direct self-gravitational force calculation for infinitesimally thin gaseous disks. This approach uses a two-dimensional kernel that is derived for infinitesimally thin disks and is free of artificial boundary conditions. The self-gravitational force calculation is presented in generalized convolution forms for a nested grid configuration. A numerical technique derived from a fast Fourier transform is employed to reduce the computational complexity to be nearly linear. By comparing with analytic potential–density pairs associated with the generalized Maclaurin disks, the extended approach is verified to be of second-order accuracy when using numerical simulations. The proposed method is accurate, computationally fast, and has the potential to be applied to studies of planetary migration and the gaseous morphology of disk galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

DNA Sequencing Using capillary Electrophoresis

The overall goal of this program was to develop capillary electrophoresis as the tool to be used to sequence for the first time the Human Genome. Our program was part of the Human Genome Project. In this work, we were highly successful and the replaceable polymer we developed, linear polyacrylamide, was used by the DOE sequencing lab in California to sequence a significant portion of the human genome using the MegaBase multiple capillary array electrophoresis instrument. In this final report, we summarize our efforts and success. We began our work by separating by capillary electrophoresis double strand oligonucleotides using cross-linked polyacrylamide gels in fused silica capillaries. This work showed the potential of the methodology. However, preparation of such cross-linked gel capillaries was difficult with poor reproducibility, and even more important, the columns were not very stable. We improved stability by using non-cross linked linear polyacrylamide. Here, the entangled linear chains could move when osmotic pressure (e.g. sample injection) was imposed on the polymer matrix. This relaxation of the polymer dissipated the stress in the column. Our next advance was to use significantly lower concentrations of the linear polyacrylamide that the polymer could be automatically blown out after each run and replaced with fresh linear polymer solution. In this way, a new column was available for each analytical run. Finally, while testing many linear polymers, we selected linear polyacrylamide as the best matrix as it was the most hydrophilic polymer available. Under our DOE program, we demonstrated initially the success of the linear polyacrylamide to separate double strand DNA. We note that the method is used even today to assay purity of double stranded DNA fragments. Our focus, of course, was on the separation of single stranded DNA for sequencing purposes. In one paper, we demonstrated the success of our approach in sequencing up to 500 bases. Other application papers of sequencing up to this level were also published in the mid 1990's. A major interest of the sequencing community has always been read length. The longer the sequence read per run the more efficient the process as well as the ability to read repeat sequences. We therefore devoted a great deal of time to studying the factors influencing read length in capillary electrophoresis, including polymer type and molecule weight, capillary column temperature, applied electric field, etc. In our initial optimization, we were able to demonstrate, for the first time, the sequencing of over 1000 bases with 90% accuracy. The run required 80 minutes for separation. Sequencing of 1000 bases per column was next demonstrated on a multiple capillary instrument. Our studies revealed that linear polyacrylamide produced the longest read lengths because the hydrophilic single strand DNA had minimal interaction with the very hydrophilic linear polyacrylamide. Any interaction of the DNA with the polymer would lead to broader peaks and lower read length. Another important parameter was the molecular weight of the linear chains. High molecular weight (> 1 MDA) was important to allow the long single strand DNA to reptate through the entangled polymer matrix. In an important paper, we showed an inverse emulsion method to prepare reproducibility linear polyacrylamide polymer with an average MWT of 9MDa. This approach was used in the polymer for sequencing the human genome. Another critical factor in the successful use of capillary electrophoresis for sequencing was the sample preparation method. In the Sanger sequencing reaction, high concentration of salts and dideoxynucleotide remained. Since the sample was introduced to the capillary column by electrokinetic injection, these salt ions would be favorably injected into the column over the sequencing fragments, thus reducing the signal for longer fragments and hence reading read length. In two papers, we examined the role of individual components from the sequencing reaction and then developed a protocol to reduce the deleterious salts. We demonstrated a robust method for achieving long read length DNA sequencing. Continuing our advances, we next demonstrated the achievement of over 1000 bases in less than one hour with a base calling accuracy of between 98 and 99%. In this work, we implemented energy transfer dyes which allowed for cleaner differentiation of the 4 dye labeled terminal nucleotides. In addition, we developed improved base calling software to help read sequencing when the separation was only minimal as occurs at long read lengths. Another critical parameter we studied was column temperature. We demonstrated that read lengths improved as the column temperature was increased from room temperature to 60 C or 70 C. The higher temperature relaxed the DNA chains under the influence of the high electric field.

59 BASIC BIOLOGICAL SCIENCES↗

Gyro-induced acceleration of magnetic reconnection

The linear and nonlinear evolution of magnetic reconnection in collisionless high-temperature plasmas with a strong guide field is analyzed on the basis of a two-dimensional gyrofluid model. The linear growth rate of the reconnecting instability is compared to analytical calculations over the whole spectrum of linearly unstable wave numbers. In the strongly unstable regime (large Δ′), the nonlinear evolution of the reconnecting instability is found to undergo two distinctive acceleration phases separated by a stall phase in which the instantaneous growth rate decreases. The first acceleration phase is caused by the formation of strong electric fields close to the X-point due to ion gyration, while the second acceleration phase is driven by the development of an open Petschek-like configuration due to both ion and electron temperature effects. Furthermore, the maximum instantaneous growth rate is found to increase dramatically over its linear value for decreasing diffusion layers. This is a consequence of the fact that the peak instantaneous growth rate becomes weakly dependent on the microscopic plasma parameters if the diffusion region thickness is sufficiently smaller than the equilibrium magnetic field scale length. When this condition is satisfied, the peak reconnection rate asymptotes to a constant value.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Benchmarking a new closed-form thermal analysis technique against a traditional lumped parameter, finite-difference method

A benchmarking effort was conducted to determine the accuracy of a new analytic generic geology thermal repository model developed at LLNL relative to a more traditional, numerical, lumped parameter technique. The fast-running analytical thermal transport model assumes uniform thermal properties throughout a homogenous storage medium. Arrays of time-dependent heat sources are included geometrically as arrays of line segments and points. The solver uses a source-based linear superposition of closed form analytical functions from each contributing point or line to arrive at an estimate of the thermal evolution of a generic geologic repository. Temperature rise throughout the storage medium is computed as a linear superposition of temperature rises. It is modeled using the MathCAD mathematical engine and is parameterized to allow myriad gridded repository geometries and geologic characteristics [4]. It was anticipated that the accuracy and utility of the temperature field calculated with the LLNL analytical model would provide an accurate 'birds-eye' view in regions that are many tunnel radii away from actual storage units; i.e., at distances where tunnels and individual storage units could realistically be approximated as physical lines or points. However, geometrically explicit storage units, waste packages, tunnel walls and close-in rock are not included in the MathCAD model.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

CO-occurring exposure to perchlorate, nitrate and thiocyanate alters thyroid function in healthy pregnant women

Background: Adequate maternal thyroid function during pregnancy is necessary for normal fetal brain development, making pregnancy a critical window of vulnerability to thyroid disrupting insults. Sodium/iodide symporter (NIS) inhibitors, namely perchlorate, nitrate, and thiocyanate, have been shown individually to competitively inhibit uptake of iodine by the thyroid. Several epidemiologic studies examined the association between these individual exposures and thyroid function. Few studies have examined the effect of this chemical mixture on thyroid function during pregnancy Objectives: We examined the cross sectional association between urinary perchlorate, thiocyanate and nitrate concentrations and thyroid function among healthy pregnant women living in New York City using weighted quantile sum (WQS) regression. Methods: We measured thyroid stimulating hormone (TSH) and free thyroxine (FreeT4) in blood samples; perchlorate, thiocyanate, nitrate and iodide in urine samples collected from 284 pregnant women at 12 (±2.8) weeks gestation. We examined associations between urinary analyte concentrations and TSH or FreeT4 using linear regression or WQS adjusting for gestational age, urinary iodide and creatinine. Results: Individual analyte concentrations in urine were significantly correlated (Spearman's r 0.4–0.5, p<0.001). Linear regression analyses did not suggest associations between individual concentrations and thyroid function. The WQS revealed a significant positive association between the weighted sum of urinary concentrations of the three analytes and increased TSH. Perchlorate had the largest weight in the index, indicating the largest contribution to the WQS. Conclusions: Co-exposure to perchlorate, nitrate and thiocyanate may alter maternal thyroid function, specifically TSH, during pregnancy. - Highlights: • Perchlorate, nitrate, thiocyanate and iodide measured in maternal urine. • Thyroid function (TSH and Free T4) measured in maternal blood. • Weighted quantile sum (WQS) regression examined complex mixture effect. • WQS identified an inverse association between the exposure mixture and maternal TSH. • Perchlorate indicated as the ‘bad actor’ of the mixture.

54 ENVIRONMENTAL SCIENCES↗

Analytical admittance characterization of high mobility channel

In this contribution, we investigate the small-signal admittance of the high electron mobility transistors field-effect channels under a continuation branching of the current between channel and gate by using an analytical model. The analytical approach takes into account the linearization of the 2D Poisson equation and the drift current along the channel. The analytical equations discuss the frequency dependence of the admittance at source and drain terminals on the geometrical transistor parameters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Development of harmonic balance capability for Charon

We report on the development of a frequency domain method of analysis in the Panzer foundation of Charon. We first present a harmonic balance approach for calculating the frequency-domain response (in its weak form) of a non-linear system of partial differential equations (PDEs). Our approach is anemable to adaptation of Charon's transient PDE models for frequency domain analysis. We make an observation that allows us to analyze either small-signal or large-signal responses with minimal specialization of the algorithm. We conclude by confirming our small- and large-signal analyses of a transient, linear Helmholtz equation by comparing its analytic solution to our results. We include figures from a sequence of non-linear perturbations of this system, showcasing the fact that, when the non-linearities are insignificant, the small- and large-signal analyses obtain similar solutions. On the other hand, we depict the inadequacy of a small-signal analysis to accurately capture the response in the presence of a large non-linearity, and underscore the requirement to employ a large-signal analysis for modelling highly non-linear systems.

97 MATHEMATICS AND COMPUTING↗

Combined preheating on the lattice with applications to Higgs inflation

We use classical lattice simulations in 3+1 dimensions to study the interplay between the resonant production of particles during preheating and the subsequent decay of these into a set of secondary species. We choose to work in a simplified version of Higgs inflation in which the Higgs field non-minimally coupled to gravity plays the role of the inflaton. Our numerical results extend the analytical estimates in the literature beyond the linear regime and shed some light on the limitations of the analytical techniques. The inclusion of fast and inefficient decays postpones the onset of parametric resonance by depleting the particles produced at the bottom of the potential. In spite of this delay, fermions are shown to play an important role on the destruction of the inflaton field. The limitations of our approach and its applications to a realistic Higgs inflation scenario are also discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Flow adjustment inside homogeneous canopies after a leading edge – An analytical approach backed by LES

A two-dimensional analytical model for describing the mean flow behavior inside a vegetation canopy after a leading edge in neutral conditions was developed and tested by means of large eddy simulations (LES) employing the LES code PALM. The analytical model is developed for the region directly after the canopy edge, the adjustment region, where one-dimensional canopy models fail due to the sharp change in roughness. The derivation of this adjustment region model is based on an analytic solution of the two-dimensional Reynolds averaged Navier–Stokes equation in neutral conditions for a canopy with constant plant area density (PAD). The main assumptions for solving the governing equations are separability of the velocity components concerning the spatial variables and the neglection of the Reynolds stress gradients. These two assumptions are verified by means of LES. To determine the emerging model parameters, a simultaneous fitting scheme was applied to the velocity and pressure data of a reference LES simulation. Furthermore a sensitivity analysis of the adjustment region model, equipped with the previously calculated parameters, was performed varying the three relevant length, the canopy height ( h ), the canopy length and the adjustment length ( L c ), in additional LES. Even if the model parameters are, in general, functions of h / L c , it was found out that the model is capable of predicting the flow quantities in various cases, when using constant parameters. Subsequently the adjustment region model is combined with the one-dimensional model of Massman, which is applicable for the interior of the canopy, to attain an analytical model capable of describing the mean flow for the full canopy domain. As a result, the model is tested against an analytical model based on a linearization approach.

54 ENVIRONMENTAL SCIENCES↗

Hydrodynamic description of long-distance spin transport through noncollinear magnetization states: Role of dispersion, nonlinearity, and damping

Nonlocal compensation of magnetic damping by spin injection has been theoretically shown to display dynamic, noncollinear magnetization states that carry spin currents over micrometer distances. Such states can be generically referred to as dissipative exchange flows (DEFs) because spatially diffusing spin currents are established by the mutual exchange torque exerted by neighboring spins. Analytical studies to date have been limited to the weak spin injection assumption whereby the equation of motion for the magnetization is mapped to hydrodynamic equations describing spin flow and then linearized. In this work, we analytically and numerically study easy-plane ferromagnetic channels subject to spin injection of arbitrary strength at one extremum under a unified hydrodynamic framework. We propose that DEFs generally exhibit a nonlinear profile along the channel accompanied by a nonlinear frequency tunability. At large injection strengths, we fully characterize a magnetization state we call a contact-soliton DEF (CS-DEF) composed of a stationary soliton at the injection site, which smoothly transitions into a DEF and exhibits a negative frequency tunability. The transition between a DEF and a CS-DEF occurs at the maximum precessional frequency and coincides with the Landau criterion: a subsonic to supersonic flow transition. Leveraging the hydraulic-electrical analogy, the current-voltage characteristics of a nonlinear DEF circuit are presented. Micromagnetic simulations of nanowires that include magnetocrystalline anisotropy and nonlocal dipole fields are in qualitative agreement with the analytical results. The magnetization states found here along with their characteristic profile and spectral features provide quantitative guidelines to pursue an experimental demonstration of DEFs in ferromagnetic materials and establish a unified description for long-distance spin transport.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Fault tolerance in parity-state linear optical quantum computing

We use a combination of analytical and numerical techniques to calculate the noise threshold and resource requirements for a linear optical quantum computing scheme based on parity-state encoding. Parity-state encoding is used at the lowest level of code concatenation in order to efficiently correct errors arising from the inherent nondeterminism of two-qubit linear-optical gates. When combined with teleported error-correction (using either a Steane or Golay code) at higher levels of concatenation, the parity-state scheme is found to achieve a saving of approximately three orders of magnitude in resources when compared to the cluster state scheme, at a cost of a somewhat reduced noise threshold.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗