Search NASA⌕ Search

SEARCH · Search NASA

Results for “CONJUGATED SYSTEM”

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

Perturbation analysis of the limit cycle of the free van der Pol equation

A power series expansion in the damping parameter, epsilon, of the limit cycle of the free van der Pol equation is constructed and analyzed. Coefficients in the expansion are computed in exact rational arithmetic using the symbolic manipulation system MACSYMA and using a FORTRAN program. The series is analyzed using Padeapproximants. The convergence of the series for the maximum amplitude of the limit cycle is limited by two pair of complex conjugate singularities in the complex epsilon-plane. A new expansion parameter is introduced which maps these singularities to infinity and leads to a new expansion for the amplitude which converges for all real values of epsilon. Amplitudes computed from this transformed series agree very well with reported numerical and asymptotic results. For the limit cycle itself, convergence of the series expansion is limited by three pair of complex conjugate branch point singularities. Two pair remain fixed throughout the cycle, and correspond to the singularities found in the maximum amplitude series, while the third pair moves in the epsilon-plane as a function of t from one of the fixed pairs to the other. The limit cycle series is transformed using a new expansion parameter, which leads to a new series that converges for larger values of epsilon.

Dadfar, M. B.↗

A numerically exact full wave packet approach to molecule-surface scattering

A numerically exact spectral method for solving the time-dependent Schroedinger equation in spherical coordinates is described. The angular dependence of the wave function is represented on a two-dimensional grid of evenly spaced points. The fast Fourier transform algorithm is used to transform between the angle space representation of the wave function and its conjugate representation in momentum space. The time propagation of the wave function is evaluated using an expansion of the time evolution operator as a series of Chebyshev polynomials. Calculations performed for a model system representing H2 scattering from a rectangular corrugated surface yield transition probabilities that are in excellent agreement with those obtained using the close-coupling wave packet (CCWP) method. However, the new method is found to require substantially more computation time than the CCWP method because of the large number of grid points needed to represent the angular dependence of the wave function and the variation in the number of terms required in the Chebyshev representation of the time evolution operator.

Mowrey, R. C.↗

Phase-Conjugate Mirror Removes Wave-Front Distortions

Phase-conjugate mirrors remove distortions or aberrations from wave fronts in laser beams. With help of mirrors, one obtains high-power laser beams of diffraction-limited quality. Potential applications include multistage laser systems.

Seery, Bernard D.↗

Biconjugate Impedance Matching Control of the HERO WEC: Preprint

Effective control strategies are essential for maximizing the power output of wave energy converters (WECs). In this study, we explore the methodology of port-based modelling which can be used for bi-conjugate impedance matching, a type of reactive control that involves maximizing power transfer by setting the load (PTO) impedance to be the complex conjugate of the source impedance (WEC). Despite its foundational role in electrical engineering, this method has seen limited application in wave energy and has not yet been demonstrated in wave-powered desalination systems. We represent our complex multidisciplinary system as an equivalent circuit to explore frequency-based relations. The HERO WEC, a small-scale, modular point absorber (PA) designed for reverse osmosis desalination is used as a case study to demonstrate port-based modelling for control applications. Improving the performance of wave-powered water desalination systems could play a critical role in addressing freshwater scarcity for coastal communities, especially those in remote or disaster-prone regions.

16 TIDAL AND WAVE POWER↗

A computational procedure for the dynamics of flexible beams within multibody systems

The dynamic analysis of three dimensional elastic beams which experience large rotational and large deformational motions are examined. The beam motion is modeled using an inertial reference for the translational displacements and a body-fixed reference for the rotational quantities. Finite strain rod theories are then defined in conjunction with the beam kinematic description which accounts for the effects of stretching, bending, torsion, and transverse shear deformations. A convected coordinate representation of the Cauchy stress tensor and a conjugate strain definition is introduced to model the beam deformation. To treat the beam dynamics, a two-stage modification of the central difference algorithm is presented to integrate the translational coordinates and the angular velocity vector. The angular orientation is then obtained from the application of an implicit integration algorithm to the Euler parameter/angular velocity kinematical relation. The combined developments of the objective internal force computation with the dynamic solution procedures result in the computational preservation of total energy for undamped systems. The present methodology is also extended to model the dynamics of deployment/retrieval of the flexible members. A moving spatial grid corresponding to the configuration of a deployed rigid beam is employed as a reference for the dynamic variables. A transient integration scheme which accurately accounts for the deforming spatial grid is derived from a space-time finite element discretization of a Hamiltonian variational statement. The computational results of this general deforming finite element beam formulation are compared to reported results for a planar inverse-spaghetti problem.

Downer, Janice Diane↗

Perturbation analysis of the limit cycle of the free van der Pol equation

A power series expansion in the damping parameter, epsilon, of the limit cycle of the free van der Pol equation is constructed and analyzed. Coefficients in the expansion are computed in exact rational arithmetic using the symbolic manipulation system MACSYMA and using a FORTRAN program. The series is analyzed using Pade approximants. The convergence of the series for the maximum amplitude of the limit cycle is limited by two pair of complex conjugate singularities in the complex epsilon-plane. A new expansion parameter is introduced which maps these singularities to infinity and leads to a new expansion for the amplitude which converges for all real values of epsilon. Amplitudes computed from this transformed series agree very well with reported numerical and asymptotic results. For the limit cycle itself, convergence of the series expansion is limited by three pair of complex conjugate branch point singularities. Two pair remain fixed throughout the cycle, and correspond to the singularities found in the maximum amplitude series, while the third pair moves in the epsilon-plane as a function of t from one of the fixed pairs to the other. The limit cycle series is transformed using a new expansion parameter, which leads to a new series that converges for larger values of epsilon.

Dadfar, M. B.↗

An implicit numerical scheme for the simulation of internal viscous flows on unstructured grids

The Navier-Stokes equations are solved numerically for two-dimensional steady viscous laminar flows. The grids are generated based on the method of Delaunay triangulation. A finite-volume approach is used to discretize the conservation law form of the compressible flow equations written in terms of primitive variables. A preconditioning matrix is added to the equations so that low Mach number flows can be solved economically. The equations are time marched using either an implicit Gauss-Seidel iterative procedure or a solver based on a conjugate gradient like method. A four color scheme is employed to vectorize the block Gauss-Seidel relaxation procedure. This increases the memory requirements minimally and decreases the computer time spent solving the resulting system of equations substantially. A factor of 7.6 speed up in the matrix solver is typical for the viscous equations. Numerical results are obtained for inviscid flow over a bump in a channel at subsonic and transonic conditions for validation with structured solvers. Viscous results are computed for developing flow in a channel, a symmetric sudden expansion, periodic tandem cylinders in a cross-flow, and a four-port valve. Comparisons are made with available results obtained by other investigators.

Jorgenson, Philip C. E.↗

Domain decomposition methods for the parallel computation of reacting flows

Domain decomposition is a natural route to parallel computing for partial differential equation solvers. Subdomains of which the original domain of definition is comprised are assigned to independent processors at the price of periodic coordination between processors to compute global parameters and maintain the requisite degree of continuity of the solution at the subdomain interfaces. In the domain-decomposed solution of steady multidimensional systems of PDEs by finite difference methods using a pseudo-transient version of Newton iteration, the only portion of the computation which generally stands in the way of efficient parallelization is the solution of the large, sparse linear systems arising at each Newton step. For some Jacobian matrices drawn from an actual two-dimensional reacting flow problem, comparisons are made between relaxation-based linear solvers and also preconditioned iterative methods of Conjugate Gradient and Chebyshev type, focusing attention on both iteration count and global inner product count. The generalized minimum residual method with block-ILU preconditioning is judged the best serial method among those considered, and parallel numerical experiments on the Encore Multimax demonstrate for it approximately 10-fold speedup on 16 processors.

Keyes, David E.↗

Preconditioned conjugate gradient methods for the compressible Navier-Stokes equations

The compressible Navier-Stokes equations are solved for a variety of two-dimensional inviscid and viscous problems by preconditioned conjugate gradient-like algorithms. Roe's flux difference splitting technique is used to discretize the inviscid fluxes. The viscous terms are discretized by using central differences. An algebraic turbulence model is also incorporated. The system of linear equations which arises out of the linearization of a fully implicit scheme is solved iteratively by the well known methods of GMRES (Generalized Minimum Residual technique) and Chebyschev iteration. Incomplete LU factorization and block diagonal factorization are used as preconditioners. The resulting algorithm is competitive with the best current schemes, but has wide applications in parallel computing and unstructured mesh computations.

Venkatakrishnan, V.↗

Parallelization of the Physical-Space Statistical Analysis System (PSAS)

Atmospheric data assimilation is a method of combining observations with model forecasts to produce a more accurate description of the atmosphere than the observations or forecast alone can provide. Data assimilation plays an increasingly important role in the study of climate and atmospheric chemistry. The NASA Data Assimilation Office (DAO) has developed the Goddard Earth Observing System Data Assimilation System (GEOS DAS) to create assimilated datasets. The core computational components of the GEOS DAS include the GEOS General Circulation Model (GCM) and the Physical-space Statistical Analysis System (PSAS). The need for timely validation of scientific enhancements to the data assimilation system poses computational demands that are best met by distributed parallel software. PSAS is implemented in Fortran 90 using object-based design principles. The analysis portions of the code solve two equations. The first of these is the "innovation" equation, which is solved on the unstructured observation grid using a preconditioned conjugate gradient (CG) method. The "analysis" equation is a transformation from the observation grid back to a structured grid, and is solved by a direct matrix-vector multiplication. Use of a factored-operator formulation reduces the computational complexity of both the CG solver and the matrix-vector multiplication, rendering the matrix-vector multiplications as a successive product of operators on a vector. Sparsity is introduced to these operators by partitioning the observations using an icosahedral decomposition scheme. PSAS builds a large (approx. 128MB) run-time database of parameters used in the calculation of these operators. Implementing a message passing parallel computing paradigm into an existing yet developing computational system as complex as PSAS is nontrivial. One of the technical challenges is balancing the requirements for computational reproducibility with the need for high performance. The problem of computational reproducibility is well known in the parallel computing community. It is a requirement that the parallel code perform calculations in a fashion that will yield identical results on different configurations of processing elements on the same platform. In some cases this problem can be solved by sacrificing performance. Meeting this requirement and still achieving high performance is very difficult. Topics to be discussed include: current PSAS design and parallelization strategy; reproducibility issues; load balance vs. database memory demands, possible solutions to these problems.

Larson, J. W.↗

Vacuum ultraviolet spectra of the late twilight airglow.

Evaluation of sounding rocket spectra of the late twilight (solar-zenith angle of 120 deg) ultraviolet airglow between 1260 and 1900 A. The only observed features are O I 1304 and 1356. When the instrument looked at an elevation of 17 deg above the western horizon, the brightnesses were 70 and 33 rayleighs, respectively. The upper limits on the total intensity of the Lyman-Birge-Hopfield and Vegard-Kaplan systems of N2 were 26 plus or minus 26 and 55 plus or minus 55 rayleighs, respectively. An estimate shows that a large part of the O I emissions may be due to excitation by conjugate-point electrons.

Buckley, J. L.↗

Non-grey radiation in a liquid rocket combustion chamber

The effects of radiation on droplets gasification in liquid rocket combustion chambers has been studied. The modeling includes a Legendre pseudo-spectral collocation approximation to solve the Radiative Transfer Equation (RTE). The band model has been utilized to account for non-grey emitting and absorbing gases present in the comustion chamber. The GEMCHIP II code has been utilized to study the fuel and oxidizer droplets combustion and interaction. The submodels within this code are capable of accounting for group combustion and conjugate effects between many droplets. The radiative model has been coupled with the GEMCHIP II code accounting for radial effects only, to provide the tool for studying the combustion-radiation coupling effects in a bipropellant system. The results indicate that the gasification/combustion process will be enhanced upstream of the chamber causing thicker flame sheet and associated higher combustion efficiency.

Kehtarnavaz, H.↗

Confined and deconfined chaos in classical spin systems

Weakly perturbed integrable many-body systems are typically chaotic, and thermal at late times. However, there are distinct relationships between the timescales for thermalization and chaos. The typical relationship is confined chaos: when trajectories are still confined to regions in phase space with constant conserved quantities (actions), the conjugate angle variables are already unstable. Chaotic instabilities thus far precede thermalization. In a different relationship, which we term deconfined chaos, chaotic instabilities and thermalization occur on the same timescale. We investigate these two qualitatively distinct scenarios through numerical and analytical studies of two perturbed integrable classical spin models: the Ishimori spin chain (confined chaos), and the central spin model with XX interactions (deconfined chaos). We analytically establish (super)-integrability in the latter model in a microcanonical shell. Deconfined chaos emerges through the separation of phase space into large quasi-integrable regions and a thin chaotic manifold. The latter leads to chaos and thermalization on the fastest possible timescale, which is proportional to the inverse perturbation strength. This behavior is reminiscent of the quantum SYK models and strange metals.

Chaotic Dynamics (nlin.CD)↗

Parametric Analysis of Cyclic Phase Change and Energy Storage in Solar Heat Receivers

A parametric study on cyclic melting and freezing of an encapsulated phase change material (PCM), integrated into a solar heat receiver, has been performed. The cyclic nature of the present melt/freeze problem is relevant to latent heat thermal energy storage (LHTES) systems used to power solar Brayton engines in microgravity environments. Specifically, a physical and numerical model of the solar heat receiver component of NASA Lewis Research Center's Ground Test Demonstration (GTD) project was developed. Multi-conjugate effects such as the convective fluid flow of a low-Prandtl-number fluid, coupled with thermal conduction in the phase change material, containment tube and working fluid conduit were accounted for in the model. A single-band thermal radiation model was also included to quantify reradiative energy exchange inside the receiver and losses through the aperture. The eutectic LiF-CaF2 was used as the phase change material (PCM) and a mixture of He/Xe was used as the working fluid coolant. A modified version of the computer code HOTTube was used to generate results in the two-phase regime. Results indicate that parametric changes in receiver gas inlet temperature and receiver heat input effects higher sensitivity to changes in receiver gas exit temperatures.

Hall, Carsie A., III↗

Magnetic Field Solver

The Magnetic Field Solver computer program calculates the magnetic field generated by a group of collinear, cylindrical axisymmetric electromagnet coils. Given the current flowing in, and the number of turns, axial position, and axial and radial dimensions of each coil, the program calculates matrix coefficients for a finite-difference system of equations that approximates a two-dimensional partial differential equation for the magnetic potential contributed by the coil. The program iteratively solves these finite-difference equations by use of the modified incomplete Cholesky preconditioned-conjugate-gradient method. The total magnetic potential as a function of axial (z) and radial (r) position is then calculated as a sum of the magnetic potentials of the individual coils, using a high-accuracy interpolation scheme. Then the r and z components of the magnetic field as functions of r and z are calculated from the total magnetic potential by use of a high-accuracy finite-difference scheme. Notably, for the finite-difference calculations, the program generates nonuniform two-dimensional computational meshes from nonuniform one-dimensional meshes. Each mesh is generated in such a way as to minimize the numerical error for a benchmark one-dimensional magnetostatic problem.

Ilin, Andrew V.↗

Coupling Interfacial Redox‐Reactions with In Situ Proton Generation for the Photoelectrochemical Separation of Rare‐Earth Elements

Abstract A dual‐function photoelectrochemical (PEC) separation system is demonstrated for rare‐earth element (REE) recovery. The sustainable release of the captured REEs is promoted through the synergistic integration of a redox‐reaction for electrostatic repulsion, and in situ proton generation for ion‐exchange, all driven by photoelectrochemistry. The platform consists of a redox‐copolymer, poly(ferrocenylpropyl methacrylamide‐ co ‐methacrylic acid) (P(FPMAm‐ co ‐MAA)) (PFM), conjugated with carbon nanotubes (CNTs) and coated onto titanium dioxide nanorods (TNRs). The (PFM‐CNT)/TNR spontaneously adsorbs up to 214.2 mg of Yttrium/g PFM by ion‐exchange, and demonstrates broad applicability for other REEs. The adsorbed REEs are released through the PEC oxidation of ferrocene (Fc) to ferrocenium (Fc + ), and the simultaneous PEC water splitting reactions at the TNRs that protonate the carboxylate binding groups. This dual photoelectrochemically‐driven mechanism for REE release is investigated by in situ pH measurements, as well as vibrational and X‐ray photoelectron spectroscopy. Through PEC approaches, a 68.8% reduction in energy consumption during REE recovery has been achieved compared to purely electrochemical systems, with a regeneration efficiency close to 100%. For NdFeB magnets from waste hard disk drives, Nd and Dy recovery efficiencies of 59.2 and 61.1% are achieved. The dual‐functionality of these copolymer PEC systems offers a sustainable platform for modulating critical element recovery.

Cho, Ki‐Hyun [Department of Chemical and Biomolecu↗

Sensitivity calculations for iteratively solved problems

The calculation of sensitivity derivatives of solutions of iteratively solved systems of algebraic equations is investigated. A modified finite difference procedure is presented which improves the accuracy of the calculated derivatives. The procedure is demonstrated for a simple algebraic example as well as an element-by-element preconditioned conjugate gradient iterative solution technique applied to truss examples.

Haftka, R. T.↗

Phosphate binding protein as the biorecognition element in a biosensor for phosphate

This work explores the potential use of a member of the periplasmic family of binding proteins, the phosphate binding protein (PBP), as the biorecognition element in a sensing scheme for the detection of inorganic phosphate (Pi). The selectivity of this protein originates from its natural role which, in Escherichia coli, is to serve as the initial receptor for the highly specific translocation of Pi to the cytoplasm. The single polypeptide chain of PBP is folded into two similar domains connected by three short peptide linkages that serve as a hinge. The Pi binding site is located deep within the cleft between the two domains. In the presence of the ligand, the two globular domains engulf the former in a hinge-like manner. The resultant conformational change constitutes the basis of the sensor development. A mutant of PBP (MPBP), where an alanine was replaced by a cysteine residue, was prepared by site-directed mutagenesis using the polymerase chain reaction (PCR). The mutant was expressed, from plasmid pSD501, in the periplasmic space of E. coli and purified in a single chromatographic step on a perfusion anion-exchange column. Site-specific labeling was achieved by attaching the fluorophore, N-[2-(1-maleimidyl)ethyl]-7-(diethylamino)coumarin-3-carboxamide (MDCC), to the protein through the sulfhydryl group of the cysteine moiety. Steady-state fluorescence studies of the MPBP-MDCC conjugate showed a change in the intensity of the signal upon addition of Pi. Calibration curves for Pi were constructed by relating the intensity of the fluorescence signal with the amount of analyte present in the sample. The sensing system was first developed and optimized on a spectrofluorometer using ml volumes of sample. It was then adapted to be used on a microtiter plate arrangement with microliter sample volumes. The system's versatility was finally proven by developing a fiber optic fluorescence-based sensor for monitoring Pi. In all three cases the detection limits for the analyte were in the sub-microMolar range. It was also demonstrated that the sensing system was selective for phosphate over other structurally-similar anions, paving the way for the design and development of a new family of biosensors utilizing the specific binding properties of periplasmic proteins. c2003 Elsevier B.V. All rights reserved.

NASA Program Fundamental Space Biology↗