Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conjugate gradient”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS2024

Achieving High Efficiency in Reduced Order Modeling for Large Scale Polycrystal Plasticity Simulations

Reduced order models for the nonlinear response of heterogeneous microstructures typically require a construction (or training) stage to build the reduced order basis. In this manuscript, an efficient model construction strategy for the eigenstrain homogenization method (EHM) is presented. The proposed strategy relies on a parallel, element-by-element, conjugate gradient solver. Near linear scaling has been achieved with respect to the number of degrees of freedom used to resolve the microstructure. Linear scaling with respect to the number of pre-analyses required to construct the reduced order model (ROM) follows from the EHM formulation. Furthermore, a parallel implementation for fast evaluation of the constructed ROM has been developed using shared memory parallelization. It has been shown that for large microstructures with ≈ 10,000 grains, the total computational cost of evaluating the nonlinear response of a polycrystal could be reduced by approximately an order of magnitude using 32 cores with respect to serial ROM simulation. The present methodology has been verified using an additively manufactured polycrystalline microstructure of a nickel-based superalloy, Inconel 625. The capability of the developed framework to construct a ROM for such large microstructures, as well as the ability of the ROM to predict average and local quantities of interest has been demonstrated.

microscale↗
DOE OSTI2024

An Efficient Numerical Algorithm for Solving Coupled Time-Dependent Ginzburg-Landau Equation for Superconductivity and Elasticity

A decoupled finite element algorithm is developed for simulating the vortex dynamics on an elastic superconductor which couples the time-dependent Ginzburg- Landau equation with the complex-valued superconducting order parameter and the vector-valued magnetic potential, and the elasticity equation. We present an iterative algorithm for the decoupled system arising from the time and spatial discretization using a combination of preconditioner, algebraic multigrid method (AMG) and preconditioned conjugate gradient method (PCG). The iterative algorithm allows us to perform large-scale three-dimensional simulations of mesoscale pattern formation during superconducting phase transitions with arbitrary elastic boundary conditions. Here, the performance and efficiency of the algorithm are numerically verified by several benchmark problems, exhibiting up to two orders of magnitude improvement depending on the scale of discrete system compared to the exact solver.

Efficiency↗
NASA NTRS1990

Domain decomposition methods in aerodynamics

Compressible Euler equations are solved for two-dimensional problems by a preconditioned conjugate gradient-like technique. An approximate Riemann solver is used to compute the numerical fluxes to second order accuracy in space. Two ways to achieve parallelism are tested, one which makes use of parallelism inherent in triangular solves and the other which employs domain decomposition techniques. The vectorization/parallelism in triangular solves is realized by the use of a recording technique called wavefront ordering. This process involves the interpretation of the triangular matrix as a directed graph and the analysis of the data dependencies. It is noted that the factorization can also be done in parallel with the wave front ordering. The performances of two ways of partitioning the domain, strips and slabs, are compared. Results on Cray YMP are reported for an inviscid transonic test case. The performances of linear algebra kernels are also reported.

Venkatakrishnan, V.↗
NASA NTRS1991

A Lagrange multiplier based divide and conquer finite element algorithm

A novel domain decomposition method based on a hybrid variational principle is presented. Prior to any computation, a given finite element mesh is torn into a set of totally disconnected submeshes. First, an incomplete solution is computed in each subdomain. Next, the compatibility of the displacement field at the interface nodes is enforced via discrete, polynomial and/or piecewise polynomial Lagrange multipliers. In the static case, each floating subdomain induces a local singularity that is resolved very efficiently. The interface problem associated with this domain decomposition method is, in general, indefinite and of variable size. A dedicated conjugate projected gradient algorithm is developed for solving the latter problem when it is not feasible to explicitly assemble the interface operator. When implemented on local memory multiprocessors, the proposed methodology requires less interprocessor communication than the classical method of substructuring. It is also suitable for parallel/vector computers with shared memory and compares favorably with factorization based parallel direct methods.

Farhat, C.↗
NASA NTRS1992

Iterative solution of linear systems

Recent advances in the field of iterative methods for solving large linear systems are reviewed. The main focus is on developments in the area of conjugate gradient-type algorithms and Krylov subspace methods for nonHermitian matrices.

Freund, Roland W.↗
NASA NTRS1995

Three-dimensional aerodynamic shape optimization using discrete sensitivity analysis

An aerodynamic shape optimization procedure based on discrete sensitivity analysis is extended to treat three-dimensional geometries. The function of sensitivity analysis is to directly couple computational fluid dynamics (CFD) with numerical optimization techniques, which facilitates the construction of efficient direct-design methods. The development of a practical three-dimensional design procedures entails many challenges, such as: (1) the demand for significant efficiency improvements over current design methods; (2) a general and flexible three-dimensional surface representation; and (3) the efficient solution of very large systems of linear algebraic equations. It is demonstrated that each of these challenges is overcome by: (1) employing fully implicit (Newton) methods for the CFD analyses; (2) adopting a Bezier-Bernstein polynomial parameterization of two- and three-dimensional surfaces; and (3) using preconditioned conjugate gradient-like linear system solvers. Whereas each of these extensions independently yields an improvement in computational efficiency, the combined effect of implementing all the extensions simultaneously results in a significant factor of 50 decrease in computational time and a factor of eight reduction in memory over the most efficient design strategies in current use. The new aerodynamic shape optimization procedure is demonstrated in the design of both two- and three-dimensional inviscid aerodynamic problems including a two-dimensional supersonic internal/external nozzle, two-dimensional transonic airfoils (resulting in supercritical shapes), three-dimensional transport wings, and three-dimensional supersonic delta wings. Each design application results in realistic and useful optimized shapes.

Burgreen, Gregory W.↗
NASA NTRS2010

Measurements and Simulations of Nadir-Viewing Radar Returns from the Melting Layer at X- and W-Bands

Simulated radar signatures within the melting layer in stratiform rain, namely the radar bright band, are checked by means of comparisons with simultaneous measurements of the bright band made by the EDOP (X-band) and CRS (W-band) airborne Doppler radars during the CRYSTAL-FACE campaign in 2002. A stratified-sphere model, allowing the fractional water content to vary along the radius of the particle, is used to compute the scattering properties of individual melting snowflakes. Using the effective dielectric constants computed by the conjugate gradient-fast Fourier transform (CGFFT) numerical method for X and W bands, and expressing the fractional water content of melting particle as an exponential function in particle radius, it is found that at X band the simulated radar bright-band profiles are in an excellent agreement with the measured profiles. It is also found that the simulated W-band profiles usually resemble the shapes of the measured bright-band profiles even though persistent offsets between them are present. These offsets, however, can be explained by the attenuation caused by cloud water and water vapor at W band. This is confirmed by the comparisons of the radar profiles made in the rain regions where the un-attenuated W-band reflectivity profiles can be estimated through the X- and W band Doppler velocity measurements. The bright-band model described in this paper has the potential to be used effectively for both radar and radiometer algorithms relevant to the TRMM and GPM satellite missions.

Liao, Liang↗
NASA NTRS1989

Phase conjugation by four-wave mixing in inhomogeneous plasmas

The effects of density, temperature, and velocity gradients on four-wave mixing (FWM) in a plasma are investigated. A fluid model is used in which the stimulated Brillouin terms are included, but pump depletion is neglected. The steady state phase conjugate reflectivity and signal transmission coefficients are calculated and discussed for both degenerate and resonant FWM. The substantial effects of inhomogeneity on the use of FWM as a plasma diagnostic are discussed.

Williams, Edward A.↗
DOE OSTI2023

Diffusion of Brønsted acidic dopants in conjugated polymers

Many semiconductor devices (e.g., light emitting diodes and photovoltaics) utilize heterojunctions of doped and undoped layers or depend on gradients of electronic doping to control charge transport. Understanding of the formation and stability of gradients in doping requires an understanding of diffusion of dopants and the complex changes in polymer properties that arise during doping. Conjugated polymers can be electrically doped by strong acids, but the details of the reaction mechanism and subsequent stability are not understood. Here, we show a clear kinetic isotope effect in the doping of thin films of poly(3-hexylthiophene) (P3HT) by bis(trifluoromethane)sulfonimide (HTFSI) from solution indicating that this doping process is limited by proton transfer to the polymer. Complementary X-ray photoelectron spectroscopy and dynamic secondary ion mass spectrometry (DSIMS) depth profiling of dopant concentrations show definitive evidence of dopant enrichment at the P3HT surface. These surface-limited concentration profiles suggest that diffusivity of dopants vary inversely with dopant concentration due to doping-induced changes to the structure of the conjugated polymer.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2004

Modeling of Field-Aligned Guided Echoes in the Plasmasphere

The conditions under which high frequency (f>>f(sub uh)) long-range extraordinary-mode discrete field-aligned echoes observed by the Radio Plasma Imager (RPI) on board the Imager for Magnetopause-to-Aurora Global Exploration (IMAGE) satellite in the plasmasphere are investigated by ray tracing modeling. Field-aligned discrete echoes are most commonly observed by RPI in the plasmasphere although they are also observed over the polar cap region. The plasmasphere field-aligned echoes appearing as multiple echo traces at different virtual ranges are attributed to signals reflected successively between conjugate hemispheres that propagate along or nearly along closed geomagnetic field lines. The ray tracing simulations show that field-aligned ducts with as little as 1% density perturbations (depletions) and less than 10 wavelengths wide can guide nearly field-aligned propagating high frequency X mode waves. Effective guidance of wave at a given frequency and wave normal angle (Psi) depends on the cross-field density scale of the duct, such that ducts with stronger density depletions need to be wider in order to maintain the same gradient of refractive index across the magnetic field. While signal guidance by field aligned density gradient without ducting is possible only over the polar region, conjugate field-aligned echoes that have traversed through the equatorial region are most likely guided by ducting.

Fung, Shing F.↗
DOE OSTI2020

The Generality of the GUGA MRCI Approach in COLUMBUS for Treating Complex Quantum Chemistry

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Synthesis of poly(3,4-ethylenedioxythuophene) derivatives using three-armed conjugated cross-linker and its thermoelectric properties

Conducting polymer based thermoelectric systems are very effective at harvesting electricity from waste heat with low-temperature gradients relative to environmental temperature. However, although various studies concerning the doping effect exist, there are insufficient reports on the effect of the change in polymer chain structure on the thermoelectric properties. Here we demonstrate different poly(3,4-ethylenedioxythiophene) (PEDOT)-derivative chain structures (such as linear, two-dimensional, co-polymer), by varying the amount of three-armed conjugated cross-linker (1,3,5-tri(2-thienyl)-benzene (TTB)). Further, at a small amount of TTB addition (until 0.5 mM), the electrochemical charge transport of the PEDOT films increased without shift of the UV–Vis absorbing property. However, from higher than 0.5 mM TTB, electrochemical charge transport decreased with blue shift of the UV–Vis absorbing curve. Moreover, the power factor of the PEDOT-derivative films without TTB to a small amount of TTB (0.5 mM) was improved from (27.03 to 49.43) µWm -1 K -2 , and at higher than 0.5 mM TTB concentration, was diminished, respectively, because a small amount of TTB enhanced the π–π interaction of the polymer chains without disrupting the conjugation, while excess TTB molecules blocked conjugation of the polymer.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Megahertz-rate shock-wave distortion cancellation via phase conjugate digital in-line holography

Abstract Holography is a powerful tool for three-dimensional imaging. However, in explosive, supersonic, hypersonic, cavitating, or ionizing environments, shock-waves and density gradients impart phase distortions that obscure objects in the field-of-view. Capturing time-resolved information in these environments also requires ultra-high-speed acquisition. To reduce phase distortions and increase imaging rates, we introduce an ultra-high-speed phase conjugate digital in-line holography (PCDIH) technique. In this concept, a coherent beam passes through the shock-wave distortion, reflects off a phase conjugate mirror, and propagates back through the shock-wave, thereby minimizing imaging distortions from phase delays. By implementing the method using a pulse-burst laser setup at up to 5 million-frames-per-second, time-resolved holograms of ultra-fast events are now possible. This technique is applied for holographic imaging through laser-spark plasma-generated shock-waves and to enable three-dimensional tracking of explosively generated hypersonic fragments. Simulations further advance our understanding of physical processes and experiments demonstrate ultra-high-speed PCDIH techniques for capturing dynamics.

Science & Technology - Other Topics↗
DOE OSTI2021

The impact of magnetic fields on momentum transport and saturation of shear-flow instability by stable modes

The Kelvin–Helmholtz (KH) instability of a shear layer with an initially uniform magnetic field in the direction of flow is studied in the framework of 2D incompressible magnetohydrodynamics with finite resistivity and viscosity using direct numerical simulations. The shear layer evolves freely, with no external forcing, and thus broadens in time as turbulent stresses transport momentum across it. As with hydrodynamic KH, the instability here features a conjugate stable mode for every unstable mode in the absence of dissipation. Stable modes are shown to transport momentum up its gradient, shrinking the layer width whenever they exceed unstable modes in amplitude. In simulations with weak magnetic fields, the linear instability is minimally affected by the field, but enhanced small-scale fluctuations relative to the hydrodynamic case are observed. These enhanced fluctuations coincide with increased energy dissipation and faster layer broadening, with these features more pronounced in simulations with stronger fields. These trends result from the magnetic field reducing the effects of stable modes relative to the transfer of energy to small scales. As field strength increases, stable modes become less excited, thus transporting less momentum against its gradient. Furthermore, the energy that would otherwise transfer back to the driving shear because of the stable modes is instead allowed to cascade to small scales, where it is lost to dissipation. Approximations of the turbulent state in terms of a reduced set of modes are explored. While the Reynolds stress is well-described using just two modes per wavenumber at large scales, the Maxwell stress is not.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Saturation physics of threshold heat-flux reduction

The saturation physics of ion-temperature-gradient-driven turbulence is examined in relation to the temperature-gradient variation of the heat flux, which can exhibit an upshift of the critical gradient for significant flux relative to the linear instability threshold. Gyrokinetic measurements of saturation properties and spectral energy transfer, which will be defined in Sec. II, are presented, indicating that the physics of saturation is fundamentally unchanged on either side of the upshifted gradient. To analyze heat transport below and above the upshifted critical gradient, a fluid model for toroidal ion-temperature-gradient turbulence is modified to include the kinetic instability threshold. The model and the heat flux are rendered in the eigenmode decomposition to track the dominant mode-coupling channel of zonal-flow-catalyzed transfer to a conjugate stable mode. Given linear and nonlinear symmetries, the stable mode level and the cross-correlation of the unstable and stable mode amplitudes are related to the unstable mode level via linear physics. Here, the heat flux can then be written in terms of the unstable-mode level, which through a nonlinear balance depends on the eigenmode-dependent coupling coefficients and the triplet correlation time of the dominant coupled modes. Resonance in these quantities leads to suppressed heat flux above the linear threshold, with a nonlinear upshift of the critical gradient set by the resonance broadening of a finite perpendicular wavenumber and collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2007

Self-consistent Model of Magnetospheric Ring Current and Propagating Electromagnetic Ion Cyclotron Waves. 2. Wave Induced Ring Current Precipitation and Thermal Electron Heating

This paper continues presentation and discussion of the results from our new global self-consistent theoretical model of interacting ring current ions and propagating electromagnetic ion cyclotron waves [Khazanov et al., 2006]. To study the effects of electromagnetic ion cyclotron wave propagation and refraction on the wave induced ring current precipitation and heating of the thermal plasmaspheric electrons, we simulate the May 1998 storm. The main findings after a simulation can be summarized as follows. Firstly, the wave induced ring current precipitation exhibits quite a lot of fine structure, and is highly organized by location of the plasmapause gradient. The strongest fluxes of about 4 x 10(exp 6) (cm(raised dot) s(raised dot) sr(raised dot) (sup -1)) are observed during the maill and early recovery phases of the storm. The very interesting and probably more important finding is that in a number of cases the most intense precipitating fluxes are not connected to the most intense waves in simple manner. The characteristics of the wave power spectral density distribution over the wave normal angle are extremely crucial for the effectiveness of the ring current ion scattering. Secondly, comparison of the global proton precipitating patterns with the results from RAM [Kozyra et al., 1997a] reveals that although we observe a qualitative agreement between the localizations of the wave induced precipitations in the models, there is no quantitative agreement between the magnitudes of the fluxes. The quantitative differences are mainly due to a qualitative difference between the characteristics of the wave power spectral density distributions over the wave normal angle in RAM and in our model. Thirdly, the heat fluxes to plasmaspheric electrons caused by Landau resonate energy absorption from electromagnetic ion cyclotron waves are observed in the postnoon-premidnight MLT sector, and can reach the magnitude of 10(exp 11) eV/(cm(sup 2)(raised dot)s). The Coulomb energy degradation of the RC H(+) and O(+) ions maximizes at about 10(exp 11) (eV/(cm(sup 2) (raised dot) s), and typically leads to electron energy deposition rates of about 2(raised dot) 10(exp 10) (eV/(cm(sup 2)(raised dot)s) which are observed during two periods; 32-48 hours, and 76-86 hours after 1 May, 0000 UT. The theoretically derived spatial structure of the thermal electron heating caused by interaction of the ring current with the plasmasphere is strongly supported by concurrent and conjugate plasma measurements from the plasmasphere, ring current, and topside ionosphere [Gurgiolo et al., 2005]. Finally, the wave induced intense electron heating has a structure of the spot-like patches along the most enhanced density gradients in the plasmasphere boundary layer and can be a possible driver to the observed but still not explained small-scale structures of enhanced emissions in the stable auroral red arcs.

Khazanov, G. V.↗