Search NASA⌕ Search

SEARCH · Search NASA

Results for “Integral equations”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Deuteron-Deuteron Elastic and Three and Four-Body Breakup Scattering Using the Faddeev-Yakubovskii Equations

The deuteron-deuteron elastic and three and four-body breakup scattering cross section have been calculated using the Faddeev-Yakubovskii (FY) chain-of-partition momentum-space equations. In this calculation the initial two-cluster potential is split into separable and non-separable components, and the effective potential is reduced to elastic two-body and three and four-body breakup open-channels and a closed-channel many-body contribution. The closed-channel contribution is determined by minimizing a variational bound. The Coulomb interaction was included by expanding the initial and final Coulomb states in a Coulomb-Sturmian basis. The three sets of chain-of-partition integral equations were solved for the elastic and three and four-body breakup scattering amplitudes. The calculations were performed for the S = 2 spin/L = 0 angular-momentum state. The elastic and double-breakup calculations were performed for energies up to E = 5.48 MeV, while the single-breakup calculations were performed for energies up to E = 4.17 MeV. In the case of elastic scattering, the calculated scattering length of 5 a dd = 7.8 ± 0.3 fm is in good agreement with a FY cluster reduction calculation. The calculated phase shift is smaller than that predicted by the resonating group model and this difference is believed to be due to the differences in the potential and calculation methods. The breakup cross sections were calculated as a function of initial deuteron momentum and fragmented-deuteron momentum. Here, the d+d → d+n+p cross sections were compared with neutron yield measurements and, while the measurements also included the L > 0 components, the general features were consistent. Estimates of the calculational uncertainties/bias are provided.

42 ENGINEERING↗
NASA NTRS1992

Shape reanalysis and sensitivities utilizing preconditioned iterative boundary solvers

The computational advantages associated with the utilization of preconditined iterative equation solvers are quantified for the reanalysis of perturbed shapes using continuum structural boundary element analysis (BEA). Both single- and multi-zone three-dimensional problems are examined. Significant reductions in computer time are obtained by making use of previously computed solution vectors and preconditioners in subsequent analyses. The effectiveness of this technique is demonstrated for the computation of shape response sensitivities required in shape optimization. Computer times and accuracies achieved using the preconditioned iterative solvers are compared with those obtained via direct solvers and implicit differentiation of the boundary integral equations. It is concluded that this approach employing preconditioned iterative equation solvers in reanalysis and sensitivity analysis can be competitive with if not superior to those involving direct solvers.

Guru Prasad, K.↗
DOE OSTI2020

Coordinate parameterisation and spectral method optimisation for Beltrami field solver in stellarator geometry

The numerical solution of the stepped pressure equilibrium (Hudson et al 2012 Phys. Plasmas 19 112502) requires a fast and robust solver to obtain the Beltrami field in three-dimensional geometry such as stellarators. The spectral method implemented in the stepped pressure equilibrium code (SPEC) is efficient when the domain is a hollow torus, but ill-conditioning of the discretised linear equations occurs in the (solid) toroid due to the artificially singular coordinate parameterisation near the axis. Here, we propose an improved choice for the reference axis to prevent coordinates surfaces from overlapping. Then, we examine the parity and asymptotics of the magnetic vector potential near the axis and suggest the use of recombined and rescaled Zernike radial basis functions. The maximum relative error in the magnetic field of the Wendelstein 7-X geometry is shown to reach 10 –9 at high resolution in a series of convergence tests and benchmarks against the boundary integral equation solver for Taylor states. The new method is also reported to significantly improve the accuracy of multi-volume SPEC calculations. A comparison between free-boundary SPEC and the analytical Dommaschk potential is presented with higher-than-usual Fourier resolution. It is illustrated that we are able to resolve low amplitude current sheets when an interface is placed where there is no flux surface in the analytic solution. This was previously concealed because of insufficient numerical resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Multiscale modeling of HBV infection integrating intra- and intercellular viral propagation to analyze extracellular viral markers

Chronic infection with hepatitis B virus (HBV) is caused by the persistence of closed circular DNA (cccDNA) in the nucleus of infected hepatocytes. Despite available therapeutic anti-HBV agents, eliminating the cccDNA remains challenging. Thus, quantifying and understanding the dynamics of cccDNA are essential for developing effective treatment strategies and new drugs. However, such study requires repeated liver biopsy to measure the intrahepatic cccDNA, which is basically not accepted because liver biopsy is potentially morbid and not common during hepatitis B treatment. We here aimed to develop a noninvasive method for quantifying cccDNA in the liver using surrogate markers in peripheral blood. We constructed a multiscale mathematical model that explicitly incorporates both intracellular and intercellular HBV infection processes. The model, based on age-structured partial differential equations, integrates experimental data from in vitro and in vivo investigations. By applying this model, we roughly predicted the amount and dynamics of intrahepatic cccDNA within a certain range using specific viral markers in serum samples, including HBV DNA, HBsAg, HBeAg, and HBcrAg. Our study represents a significant step towards advancing the understanding of chronic HBV infection. The noninvasive quantification of cccDNA using our proposed method holds promise for improving clinical analyses and treatment strategies. By comprehensively describing the interactions of all components involved in HBV infection, our multiscale mathematical model provides a valuable framework for further research and the development of targeted interventions.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1975

Theoretical investigations of plasma processes in the ion bombardment thruster

A physical model for a thruster discharge was developed, consisting of a spatially diverging plasma sustained electrically between a small ring cathode and a larger ring anode in a cylindrical chamber with an axial magnetic field. The associated boundary-value problem for the coupled partial differential equations with mixed boundary conditions, which describe the electric potential and the plasma velocity fields, was solved in closed form. By means of quantum-mechanical perturbation theory, a formula for the number S(E) of atoms sputtered on the average by an ion of energy E was derived from first principles. The boundary-value problem describing the diffusion of the sputtered atoms through the surrounding rarefied electron-ion plasma to the system surfaces of ion propulsion systems was formulated and treated analytically. It is shown that outer boundary-value problems of this type lead to a complex integral equation, which requires numerical resolution.

Wilhelm, H. E.↗
NASA NTRS1993

A class of hybrid finite element methods for electromagnetics: A review

Integral equation methods have generally been the workhorse for antenna and scattering computations. In the case of antennas, they continue to be the prominent computational approach, but for scattering applications the requirement for large-scale computations has turned researchers' attention to near neighbor methods such as the finite element method, which has low O(N) storage requirements and is readily adaptable in modeling complex geometrical features and material inhomogeneities. In this paper, we review three hybrid finite element methods for simulating composite scatterers, conformal microstrip antennas, and finite periodic arrays. Specifically, we discuss the finite element method and its application to electromagnetic problems when combined with the boundary integral, absorbing boundary conditions, and artificial absorbers for terminating the mesh. Particular attention is given to large-scale simulations, methods, and solvers for achieving low memory requirements and code performance on parallel computing architectures.

Volakis, J. L.↗
NASA NTRS1989

Structural zooming research and development of an interactive computer graphical interface for stress analysis of cracks

Engineering problems sometimes involve the numerical solution of boundary value problems over domains containing geometric feature with widely varying scales. Often, a detailed solution is required at one or more of these features. Small details in large structures may have profound effects upon global performance. Conversely, large-scale conditions may effect local performance. Many man-hours and CPU-hours are currently spent in modeling such problems. With the structural zooming technique, it is now possible to design an integrated program which allows the analyst to interactively focus upon a small region of interest, to modify the local geometry, and then to obtain highly accurate responses in that region which reflect both the properties of the overall structure and the local detail. A boundary integral equation analysis program, called BOAST, was recently developed for the stress analysis of cracks. This program can accurately analyze two-dimensional linear elastic fracture mechanics problems with far less computational effort than existing finite element codes. An interactive computer graphical interface to BOAST was written. The graphical interface would have several requirements: it would be menu-driven, with mouse input; all aspects of input would be entered graphically; the results of a BOAST analysis would be displayed pictorially but also the user would be able to probe interactively to get numerical values of displacement and stress at desired locations within the analysis domain; the entire procedure would be integrated into a single, easy to use package; and it would be written using calls to the graphic package called HOOPS. The program is nearing completion. All of the preprocessing features are working satisfactorily and were debugged. The postprocessing features are under development, and rudimentary postprocessing should be available by the end of the summer. The program was developed and run on a VAX workstation, and must be ported to the SUN workstation. This activity is currently underway.

Gerstle, Walter↗
NASA NTRS1978

Radiative transfer calculated from a Markov chain formalism

The theory of Markov chains is used to formulate the radiative transport problem in a general way by modeling the successive interactions of a photon as a stochastic process. Under the minimal requirement that the stochastic process is a Markov chain, the determination of the diffuse reflection or transmission from a scattering atmosphere is equivalent to the solution of a system of linear equations. This treatment is mathematically equivalent to, and thus has many of the advantages of, Monte Carlo methods, but can be considerably more rapid than Monte Carlo algorithms for numerical calculations in particular applications. We have verified the speed and accuracy of this formalism for the standard problem of finding the intensity of scattered light from a homogeneous plane-parallel atmosphere with an arbitrary phase function for scattering. Accurate results over a wide range of parameters were obtained with computation times comparable to those of a standard 'doubling' routine. The generality of this formalism thus allows fast, direct solutions to problems that were previously soluble only by Monte Carlo methods. Some comparisons are made with respect to integral equation methods.

Esposito, L. W.↗
NASA NTRS1999

Bidirectional Reflectance of a Macroscopically Flat, High-Albedo Particulate Surface: An Efficient Radiative Transfer Solution and Applications to Regoliths

Many remote sensing applications rely on accurate knowledge of the bidirectional reflection function (BRF) of surfaces composed of discrete, randomly positioned scattering particles. Theoretical computations of BRFs for plane-parallel particulate layers are usually reduced to solving the radiative transfer equation (RTE) using one of existing exact or approximate techniques. Since semi-empirical approximate approaches are notorious for their low accuracy, violation of the energy conservation law, and ability to produce unphysical results, the use of numerically exact solutions of RTE has gained justified popularity. For example, the computation of BRFs for macroscopically flat particulate surfaces in many geophysical publications is based on the adding-doubling (AD) and discrete ordinate (DO) methods. A further saving of computer resources can be achieved by using a more efficient technique to solve the plane-parallel RTE than the AD and DO methods. Since many natural particulate surfaces can be well represented by the model of an optically semi-infinite, homogeneous scattering layer, one can find the BRF directly by solving the Ambartsumian's nonlinear integral equation using a simple iterative technique. In this way, the computation of the internal radiation field is avoided and the computer code becomes highly efficient and very accurate and compact. Furthermore, the BRF thus obtained fully obeys the fundamental physical laws of energy conservation and reciprocity. In this paper, we discuss numerical aspects and the computer implementation of this technique, examine the applicability of the Henyey-Greenstein phase function and the sigma-Eddington approximation in BRF and flux calculations, and describe sample applications demonstrating the potential effect of particle shape on the bidirectional reflectance of flat regolith surfaces. Although the effects of packing density and coherent backscattering are currently neglected, they can also be incorporated. The FORTRAN implementation of the technique is available on the World Wide Web, and can be applied to a wide range of remote sensing problems. BRF computations for undulated (macroscopically rough) surfaces are more complicated and often rely on time consuming Monte Carlo procedures. This approach is especially inefficient for optically thick, weakly absorbing media (e.g., snow and desert surfaces at visible wavelengths since a photon may undergo many internal scattering events before it exists the medium or is absorbed. However, undulated surfaces can often be represented as collections of locally flat tilted facets characterized by the BRF found from the traditional plane parallel RTE. In this way the MOnte Carlo procedure could be used only to evaluate the effects of surface shadowing and multiple surface reflections, thereby bypassing the time-consuming ray tracing inside the medium and providing a great savings of CPU time.

Mishchenko, Michael I.↗
NASA NTRS2011

Large Ka-Band Slot Array for Digital Beam-Forming Applications

This work describes the development of a large Ka Band Slot Array for the Glacier and Land Ice Surface Topography Interferometer (GLISTIN), a proposed spaceborne interferometric synthetic aperture radar for topographic mapping of ice sheets and glaciers. GLISTIN will collect ice topography measurement data over a wide swath with sub-seasonal repeat intervals using a Ka-band digitally beamformed antenna. For technology demonstration purpose a receive array of size 1x1 m, consisting of 160x160 radiating elements, was developed. The array is divided into 16 sticks, each stick consisting of 160x10 radiating elements, whose outputs are combined to produce 16 digital beams. A transmit array stick was also developed. The antenna arrays were designed using Elliott's design equations with the use of an infinite-array mutual-coupling model. A Floquet wave model was used to account for external coupling between radiating slots. Because of the use of uniform amplitude and phase distribution, the infinite array model yielded identical values for all radiating elements but for alternating offsets, and identical coupling elements but for alternating positive and negative tilts. Waveguide-fed slot arrays are finding many applications in radar, remote sensing, and communications applications because of their desirable properties such as low mass, low volume, and ease of design, manufacture, and deployability. Although waveguide-fed slot arrays have been designed, built, and tested in the past, this work represents several advances to the state of the art. The use of the infinite array model for the radiating slots yielded a simple design process for radiating and coupling slots. Method of moments solution to the integral equations for alternating offset radiating slots in an infinite array environment was developed and validated using the commercial finite element code HFSS. For the analysis purpose, a method of moments code was developed for an infinite array of subarrays. Overall the 1x1 m array was found to be successful in meeting the objectives of the GLISTIN demonstration antenna, especially with respect to the 0.042deg, 1/10th of the beamwidth of each stick, relative beam alignment between sticks.

Rengarajan, Sembiam↗
NASA NTRS2011

Collocation and Galerkin Time-Stepping Methods

We study the numerical solutions of ordinary differential equations by one-step methods where the solution at tn is known and that at t(sub n+1) is to be calculated. The approaches employed are collocation, continuous Galerkin (CG) and discontinuous Galerkin (DG). Relations among these three approaches are established. A quadrature formula using s evaluation points is employed for the Galerkin formulations. We show that with such a quadrature, the CG method is identical to the collocation method using quadrature points as collocation points. Furthermore, if the quadrature formula is the right Radau one (including t(sub n+1)), then the DG and CG methods also become identical, and they reduce to the Radau IIA collocation method. In addition, we present a generalization of DG that yields a method identical to CG and collocation with arbitrary collocation points. Thus, the collocation, CG, and generalized DG methods are equivalent, and the latter two methods can be formulated using the differential instead of integral equation. Finally, all schemes discussed can be cast as s-stage implicit Runge-Kutta methods.

Huynh, H. T.↗
DOE OSTI2026

Finite-element boundary-integral simulation of thin wires and inhomogeneous penetrable bodies in subsurface multilayered anisotropic media

With the prevailing presence of drilling wells near the subsurface in mature oil and gas fields, the application of electromagnetic methods can be particularly challenging where the electromagnetic field is affected by the steel casing. In the past decades, borehole-to-surface and crosswell electromagnetic methods have been utilized for monitoring of reservoir and underground CO 2 storage. This paper presents a unified finite-element boundary-integral (FEBI) method capable of simultaneously modeling the complex electromagnetic interactions between thin metallic wires (representing steel casings) with 3D trajectory and arbitrary 3D inhomogeneous penetrable bodies (such as CO 2 plumes or hydrocarbon reservoirs) within anisotropic multilayered subsurface environments. Unlike existing approaches that treat these components separately or require dense discretization, or are limited to vertical wells, our unified formulation preserves flexible electromagnetic coupling while delivering improved computational efficiency. Assuming the background formation is multilayered anisotropic media, the surface integral equation method is applied to model the thin wires and boundaries of the inhomogeneous bodies. Meanwhile, the finite element method is applied to model the volume of inhomogeneous bodies. Here, the performance of the proposed FEBI method is assessed through comparison with reference numerical results and its practical significance is demonstrated through CO 2 plume monitoring scenarios.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2020

Simultaneous Retrieval of Surface Roughness Parameters for Bare Soils From Combined Active–Passive Microwave SMAP Observations

An active–passive microwave retrieval algorithm for simultaneous determination of soil surface roughness parameters [vertical root-mean-square (RMS) height ( s ) and horizontal correlation length ( l )] is presented for bare soils. The algorithm is based on active–passive microwave covariation, including the improved Integral Equation Method (I(exp 2) EM), and is tested with global soil moisture active passive (SMAP) observations. The estimated retrieval results for s and l are overall consistent with values in the literature, indicating the validity of the proposed algorithm. Sensitivity analyses showed that the developed roughness retrieval algorithm is independent of permittivity for E(sub s) > 10 [-]. Furthermore, the physical model basis of this approach (I(exp 2) EM) allows the application of different autocorrelation functions (ACF), such as Gaussian and exponential ACFs. Global roughness retrieval results confirm bare areas in deserts such as Sahara or Gobi. However, the type of ACF used within roughness parameter estimation is important. Retrieval results for the Gaussian ACF describe a rougher surface than retrieval results for the exponential ACF. No correlations were found between roughness results and the amount of precipitation or the soil texture, which could be due to the coarse spatial resolution of the SMAP data. The extension of this approach to vegetated soils is planned as an add-on study.

Correlation length↗
DOE OSTI2023

Experimental and theoretical investigations of rotating algae biofilm reactors (RABRs): Areal productivity, nutrient recovery, and energy efficiency

Microalgae biofilms have been demonstrated to recover nutrients from wastewater and serve as biomass feedstock for bioproducts. However, there is a need to develop a platform to quantitatively describe microalgae biofilm production, which can provide guidance and insights for improving biomass areal productivity and nutrient uptake efficiency. Here, this paper proposes a unified experimental and theoretical framework to investigate algae biofilm growth on a rotating algae biofilm reactor (RABR). Experimental laboratory setups are used to conduct controlled experiments on testing environmental and operational factors for RABRs. We propose a differential–integral equation‐based mathematical model for microalgae biofilm cultivation guided by laboratory experimental findings. The predictive mathematical model development is coordinated with laboratory experiments of biofilm areal productivity associated with ammonia and inorganic phosphorus uptake by RABRs. The unified experimental and theoretical tool is used to investigate the effects of RABR rotating velocity, duty cycle (DC), and light intensity on algae biofilm growth, areal productivity, nutrient uptake efficiency, and energy efficiency in wastewater treatment. Our framework indicates that maintaining a reasonable light intensity range improves biomass areal productivity and nutrient uptake efficiency. Our framework also indicates that faster RABR rotation benefits biomass areal productivity. However, maximizing the nutrient uptake efficiency requires a reasonably low RABR rotating speed. Energy efficiency is strongly correlated with RABR rotating speed and DC.

09 BIOMASS FUELS↗
DOE OSTI2024

An integrated EOS, pore‐crush, strength and damage model framework for near‐field ground‐shock

Abstract An integrated Equation of State (EOS) and strength/pore‐crush/damage model framework is provided for modeling near to source (near‐field) ground‐shock response, where large deformations and pressures necessitate coupling EOS with pressure‐dependent plastic yield and damage. Nonlinear pressure‐dependence of strength up to high‐pressures is combined with a Modified Cam‐Clay‐like cap‐plasticity model in a way to allow degradation of strength from pore‐crush damage, what we call the “Yp‐Cap” model. Nonlinear hardening under compaction allows modeling the crush‐out of pores in combination with a fully saturated EOS, that is, for modeling partially saturated ground‐shock response, where air‐filled voids crush. Attention is given to algorithmic clarity and efficiency of the provided model, and the model is employed in example numerical simulations, including finite element simulations of underground explosions to exemplify its robustness and utility.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Finite- and infinite-volume study of DDπ scattering

We develop a comprehensive framework for extracting the pole position and properties of the doubly-charmed tetraquark $T$$^{+}_{cc}$(3875) from lattice QCD data using the relativistic three-particle formalism. This approach incorporates the effect of the one-pion exchange diagram in DDπ and DD ∗ scattering, making it applicable at energies coinciding with the left-hand cut in the partial-wave projected DD ∗ amplitude. We present an example application of this framework to existing lattice QCD data at mπ = 280 MeV. We solve the integral equations describing the DDπ reaction, use LSZ reduction to determine the corresponding DD ∗ amplitude, and find the values of the infinite-volume two- and three-body K matrices that lead to agreement with lattice DD ∗ phase shifts within their uncertainties. Using these K matrices in the three-particle quantization condition, we describe the finite- volume DD ∗ spectrum and find good agreement with the lattice QCD energies. Our results suggest that, at this pion mass, the tetraquark appears as a pair of subthreshold complex poles whose precise location strongly depends on the value of the DDπ three-particle K matrix.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

The reflected entanglement spectrum for free fermions

We consider the reflected entropy and the associated entanglement spectrum for free fermions reduced to two intervals in 1 + 1 dimensions. Working directly in the continuum theory the reflected entropy can be extracted from the spectrum of a singular integral equation whose kernel is determined by the known free fermion modular evolved correlation function. We find the spectrum numerically and analytically in certain limits. For intervals that almost touch the reflected entanglement spectrum approaches the spectrum of the thermal density matrix. This suggests that the reflected entanglement spectrum is well suited to the task of extracting physical data of the theory directly from the ground state wave function.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Three-body resonances in the φ 4 theory

We study the properties of three-body resonances using a lattice complex scalar φ4 theory with two scalars, with parameters chosen such that one heavy particle can decay into three light ones. We determine the two- and three-body spectra for several lattice volumes using variational techniques, and then analyze them with two versions of the three-particle finite-volume formalism: the Relativistic Field Theory approach and the Finite-Volume Unitarity approach. We find that both methods provide an equivalent description of the energy levels, and we are able to fit the spectra using simple parametrizations of the scattering quantities. By solving the integral equations of the corresponding three-particle formalisms, we determine the pole position of the resonance in the complex energy plane and thereby its mass and width. We find very good agreement between the two methods at different values of the coupling of the theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗