Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Poisson equation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Efficient reactive Brownian dynamics

We develop a Split Reactive Brownian Dynamics (SRBD) algorithm for particle simulations of reaction-diffusion systems based on the Doi or volume reactivity model, in which pairs of particles react with a specified Poisson rate if they are closer than a chosen reactive distance. In our Doi model, we ensure that the microscopic reaction rules for various association and dissociation reactions are consistent with detailed balance (time reversibility) at thermodynamic equilibrium. The SRBD algorithm uses Strang splitting in time to separate reaction and diffusion and solves both the diffusion-only and reaction-only subproblems exactly, even at high packing densities. To efficiently process reactions without uncontrolled approximations, SRBD employs an event-driven algorithm that processes reactions in a time-ordered sequence over the duration of the time step. A grid of cells with size larger than all of the reactive distances is used to schedule and process the reactions, but unlike traditional grid-based methods such as reaction-diffusion master equation algorithms, the results of SRBD are statistically independent of the size of the grid used to accelerate the processing of reactions. We use the SRBD algorithm to compute the effective macroscopic reaction rate for both reaction-limited and diffusion-limited irreversible association in three dimensions and compare to existing theoretical predictions at low and moderate densities. We also study long-time tails in the time correlation functions for reversible association at thermodynamic equilibrium and compare to recent theoretical predictions. Finally, we compare different particle and continuum methods on a model exhibiting a Turing-like instability and pattern formation. Our studies reinforce the common finding that microscopic mechanisms and correlations matter for diffusion-limited systems, making continuum and even mesoscopic modeling of such systems difficult or impossible. We also find that for models in which particles diffuse off lattice, such as the Doi model, reactions lead to a spurious enhancement of the effective diffusion coefficients.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Hydrodynamic limit of Wigner-Poisson kinetic theory: Revisited

In this paper, we revisit the hydrodynamic limit of the Langmuir wave dispersion relation based on the Wigner-Poisson model in connection with that obtained directly from the original Lindhard dielectric function based on the random-phase-approximation. It is observed that the (fourth-order) expansion of the exact Lindhard dielectric constant correctly reduces to the hydrodynamic dispersion relation with an additional term of fourth-order, beside that caused by the quantum diffraction effect. It is also revealed that the generalized Lindhard dielectric theory accounts for the recently discovered Shukla-Eliasson attractive potential (SEAP). However, the expansion of the exact Lindhard static dielectric function leads to a k{sup 4} term of different magnitude than that obtained from the linearized quantum hydrodynamics model. It is shown that a correction factor of 1/9 should be included in the term arising from the quantum Bohm potential of the momentum balance equation in fluid model in order for a correct plasma dielectric response treatment. Finally, it is observed that the long-range oscillatory screening potential (Friedel oscillations) of type cos(2k{sub F}r)/r{sup 3}, which is a consequence of the divergence of the dielectric function at point k = 2k{sub F} in a quantum plasma, arises due to the finiteness of the Fermi-wavenumber and is smeared out in the limit of very high electron number-densities, typical of white dwarfs and neutron stars. In the very low electron number-density regime, typical of semiconductors and metals, where the Friedel oscillation wavelength becomes much larger compared to the interparticle distances, the SEAP appears with a much deeper potential valley. It is remarked that the fourth-order approximate Lindhard dielectric constant approaches that of the linearized quantum hydrodynamic in the limit if very high electron number-density. By evaluation of the imaginary part of the Lindhard dielectric function, it is shown that the Landau-damping region in ω-k plane increases dramatically by increase of the electron number-density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

GPU-BASED MONTE CARLO DUST RADIATIVE TRANSFER SCHEME APPLIED TO ACTIVE GALACTIC NUCLEI

A three-dimensional parallel Monte Carlo (MC) dust radiative transfer code is presented. To overcome the huge computing-time requirements of MC treatments, the computational power of vectorized hardware is used, utilizing either multi-core computer power or graphics processing units. The approach is a self-consistent way to solve the radiative transfer equation in arbitrary dust configurations. The code calculates the equilibrium temperatures of two populations of large grains and stochastic heated polycyclic aromatic hydrocarbons. Anisotropic scattering is treated applying the Heney-Greenstein phase function. The spectral energy distribution (SED) of the object is derived at low spatial resolution by a photon counting procedure and at high spatial resolution by a vectorized ray tracer. The latter allows computation of high signal-to-noise images of the objects at any frequencies and arbitrary viewing angles. We test the robustness of our approach against other radiative transfer codes. The SED and dust temperatures of one- and two-dimensional benchmarks are reproduced at high precision. The parallelization capability of various MC algorithms is analyzed and included in our treatment. We utilize the Lucy algorithm for the optical thin case where the Poisson noise is high, the iteration-free Bjorkman and Wood method to reduce the calculation time, and the Fleck and Canfield diffusion approximation for extreme optical thick cells. The code is applied to model the appearance of active galactic nuclei (AGNs) at optical and infrared wavelengths. The AGN torus is clumpy and includes fluffy composite grains of various sizes made up of silicates and carbon. The dependence of the SED on the number of clumps in the torus and the viewing angle is studied. The appearance of the 10 {mu}m silicate features in absorption or emission is discussed. The SED of the radio-loud quasar 3C 249.1 is fit by the AGN model and a cirrus component to account for the far-infrared emission.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A collision operator for describing dissipation in noncanonical phase space

The phase space of a noncanonical Hamiltonian system is partially inaccessible due to dynamical constraints (Casimir invariants) arising from the kernel of the Poisson tensor. When an ensemble of noncanonical Hamiltonian systems is allowed to interact, dissipative processes eventually break the phase space constraints, resulting in a thermodynamic equilibrium described by a Maxwell–Boltzmann distribution. However, the time scale required to reach Maxwell–Boltzmann statistics is often much longer than the time scale over which a given system achieves a state of thermal equilibrium. Examples include diffusion in rigid mechanical systems, as well as collisionless relaxation in magnetized plasmas and stellar systems, where the interval between binary Coulomb or gravitational collisions can be longer than the time scale over which stable structures are self-organized. Here, we focus on self-organizing phenomena over spacetime scales such that particle interactions respect the noncanonical Hamiltonian structure, but yet act to create a state of thermodynamic equilibrium. We derive a collision operator for general noncanonical Hamiltonian systems, applicable to fast, localized interactions. This collision operator depends on the interaction exchanged by colliding particles and on the Poisson tensor encoding the noncanonical phase space structure, is consistent with entropy growth and conservation of particle number and energy, preserves the interior Casimir invariants, reduces to the Landau collision operator in the limit of grazing binary Coulomb collisions in canonical phase space, and exhibits a metriplectic structure. We further show how thermodynamic equilibria depart from Maxwell–Boltzmann statistics due to the noncanonical phase space structure, and how self-organization and collisionless relaxation in magnetized plasmas and stellar systems can be described through the derived collision operator.

Boltzmann equation↗
DOE OSTI2024

Pion gravitational form factors in the QCD instanton vacuum. I

The pion form factors of the QCD energy-momentum tensor (EMT) are studied in the instanton liquid model (ILM) of the QCD vacuum. In this approach, the breaking of conformal symmetry is encoded in the form of stronger-than-Poisson fluctuations in the number of instantons. For the trace of the EMT, it is shown that the gluonic trace anomaly term contributes half the pion mass, with the other half coming from the quark-mass-dependent σ term. The Q 2 dependence of the form factors is governed by glueball and scalar meson exchanges. For the EMT, the spin-0 (trace) and spin-2 (traceless rank-2 tensor) form factors are computed at next-to-leading order in the instanton density using effective quark operators. Relations between the gluon and quark contributions to the EMT form factors are derived. The form factors are also expressed in terms of the pion light-front wave functions in the ILM. The results at the low resolution scale of the inverse instanton size are evolved to higher scales using the renormalization group equation. The ILM results compare well with those of recent lattice QCD calculations. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Hybrid Finite Element--Spectral Method for the Fractional Laplacian: Approximation Theory and Efficient Solver

Here, a numerical scheme is presented for approximating fractional order Poisson problems in two and three dimensions. The scheme is based on reformulating the original problem posed over $\Omega$ on the extruded domain $\mathcal{C}=\Omega\times[0,\infty)$ following. The resulting degenerate elliptic integer order PDE is then approximated using a hybrid FEM-spectral scheme. Finite elements are used in the direction parallel to the problem domain $\Omega$, and an appropriate spectral method is used in the extruded direction. The spectral part of the scheme requires that we approximate the true eigenvalues of the integer order Laplacian over $\Omega$. We derive an a priori error estimate which takes account of the error arising from using an approximation in place of the true eigenvalues. We further present a strategy for choosing approximations of the eigenvalues based on Weyl's law and finite element discretizations of the eigenvalue problem. The system of linear algebraic equations arising from the hybrid FEM-spectral scheme is decomposed into blocks which can be solved effectively using standard iterative solvers such as multigrid and conjugate gradient. Numerical examples in two and three dimensions suggest that the approach is quasi-optimal in terms of complexity.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2011

Molecular Dynamics Simulation of the Antiamoebin Ion Channel: Linking Structure and Conductance

Molecular dynamics simulations were carried out in order to ascertain which of the potential multimeric forms of the transmembrane peptaibol channel, antiamoebin, is consistant with its measured conductance. Estimates of the conductance obtained through counting ions that cross the channel and by solving the Nernst-Planck equation yield consistent results, indicating that the motion of ions inside the channel can be satisfactorily described as diffusive.The calculated conductance of octameric channels is markedly higher than the conductance measured in single channel recordings, whereas the tetramer appears to be non-conducting. The conductance of the hexamer was estimated to be 115+/-34 pS and 74+/-20 pS, at 150 mV and 75 mV, respectively, in satisfactory agreement with the value of 90 pS measured at 75 mV. On this basis we propose that the antiamoebin channel consists of six monomers. Its pore is large enough to accommodate K(+) and Cl(-) with their first solvation shells intact. The free energy barrier encountered by K(+) is only 2.2 kcal/mol whereas Cl(-) encounters a substantially higher barrier of nearly 5 kcal/mol. This difference makes the channel selective for cations. Ion crossing events are shown to be uncorrelated and follow Poisson statistics. keywords: ion channels, peptaibols, channel conductance, molecular dynamics

Wilson, Michael A.↗
NASA NTRS2016

Estimating the Rate of Occurrence of Renal Stones in Astronauts

Changes in urine chemistry, during and post flight, potentially increases the risk of renal stones in astronauts. Although much is known about the effects of space flight on urine chemistry, no inflight incidence of renal stones in US astronauts exists and the question "How much does this risk change with space flight?" remains difficult to accurately quantify. In this discussion, we tackle this question utilizing a combination of deterministic and probabilistic modeling that implements the physics behind free stone growth and agglomeration, speciation of urine chemistry and published observations of population renal stone incidences to estimate changes in the rate of renal stone presentation. The modeling process utilizes a Population Balance Equation based model developed in the companion IWS abstract by Kassemi et al. (2016) to evaluate the maximum growth and agglomeration potential from a specified set of urine chemistry values. Changes in renal stone occurrence rates are obtained from this model in a probabilistic simulation that interrogates the range of possible urine chemistries using Monte Carlo techniques. Subsequently, each randomly sampled urine chemistry undergoes speciation analysis using the well-established Joint Expert Speciation System (JESS) code to calculate critical values, such as ionic strength and relative supersaturation. The Kassemi model utilizes this information to predict the mean and maximum stone size. We close the assessment loop by using a transfer function that estimates the rate of stone formation from combining the relative supersaturation and both the mean and maximum free stone growth sizes. The transfer function is established by a simulation analysis which combines population stone formation rates and Poisson regression. Training this transfer function requires using the output of the aforementioned assessment steps with inputs from known non-stone-former and known stone-former urine chemistries. Established in a Monte Carlo system, the entire renal stone analysis model produces a probability distribution of the stone formation rate and an expected uncertainty in the estimate. The utility of this analysis will be demonstrated by showing the change in renal stone occurrence predicted by this method using urine chemistry distributions published in Whitson et al. 2009. A comparison to the model predictions to previous assessments of renal stone risk will be used to illustrate initial validation of the model.

gravitational physiology↗
DOE OSTI2018

Material Model Testing in the DYNA3D Code: Material Models 1 through 5

This report is the first volume in an effort to study the material model implementations in DYNA3D, using the latest developmental version of the code. Input files were derived from a template supplied by Jerry Lin to perform an equation-of-state study. Most of the simulations use 10 cubes, each 1 cm on a side with a density of 8.03 (steel). Each of the cubes tests a different aspect of the material model. The first three cubes are testing normal stresses and normal strains (uniaxial stresses and strains) placed on the cube parallel to the x-axis, y-axis and z-axis, respectively. The x (or y or z) = 1 face of the cube is moved along the x-axis (or y-axis or z-axis) to test the model, while the opposite side remains in the x (or y or z) = 0 plane. For some of the tests, the nodes are free to move perpendicular to the applied strain to allow the cross-sectional area of the cube to change in response to being placed under pressure or tension (to test Poisson’s Ratio), while in others the nodes are frozen to eliminate dynamic effects resulting from the code. The cube is placed in tension along the axis being tested and then released, followed by a cycle of compression and release. This tension/compression cycle is repeated 3 more times. The additional cycles allow for the model to be tested under different loads. For these calculations, the amount by which the cube is stretched and compressed is varied from cycle to cycle to allow testing of model features. Cubes number 4 and 5 are intended to test xy shear stresses/strains. This is done by sliding one face of the cube in its plane parallel to the axis of interest. For example, for cube 4, the upper face of the cube is slid in the y=1 plane in a direction parallel to the x-axis (so there are no changes in the y and z coordinates of the points in the plane) and then returned to its starting position. This preserves the volume of the cube, so there is no change in density or pressure. The cycle is repeated 8 times. For cube 5, the plane of the cube satisfying x=1 is slid in this plane in a direction parallel to the y-axis and then returned to its starting position. This preserves the volume of the cube, so there is no change in density or pressure. The cycle is repeated 8 times. Cubes number 6 through 9 repeat these tests for yz shear stresses/strains (Cubes 6 and 7) and zx shear stresses and strains (Cubes 8 and 9). For reasons discussed in the next paragraph, it is unclear if these tests represent a state of pure shear stress and shear strain in the code. Application of shear stresses and strains should not, in general, result in the appearance of normal stresses and strains in Cubes 4 - 9. In the infinitesimal limit, this result might be obtained. However, DYNA3D is a finite deformation code and in the real deformations used in the code, this result is not achieved and the individual components of the normal stresses and strains are non-zero during the cycling. However, the sums of the normal strains are zero in these tests and, except in the case of thermally dependent models, so are the sums of the normal stresses. Cube 10 tests the response of the material model to a hydrostatic loading of the cube. The load curve used for the first three tests is used again and applied to each of the 8 nodes in a manner causing the cube to be either compressed or expanded. This compression/tension cycle is repeated 4 times. The additional cycles allow for the model to be tested under different loads. The amount by which the cube is stretched and compressed is varied from cycle to cycle to allow testing of model features. This document will discuss the testing of Material Models 1 through 5. Relevant results from the test cubes will be discussed. A discussion of the results from all 10 test cubes will be presented only if necessary.

36 MATERIALS SCIENCE↗
DOE OSTI2016

An extended micromechanics method for probing interphase properties in polymer nanocomposites [An extended micromechanics method for overlapping geometries with application to polymer nanocomposites]

Inclusions comprised on filler particles and interphase regions commonly form complex morphologies in polymer nanocomposites. Addressing these morphologies as systems of overlapping simple shapes allows for the study of dilute particles, clustered particles, and interacting interphases all in one general modeling framework. To account for the material properties in these overlapping geometries, weighted-mean and additive overlapping conditions are introduced and the corresponding inclusion-wise integral equations are formulated. An extended micromechanics method based on these overlapping conditions for linear elastic and viscoelastic heterogeneous material is then developed. An important feature of the proposed approach is that the effect of both the geometric overlapping (clustered particles) and physical overlapping (interacting interphases) on the effective properties can be distinguished. Lastly, we apply the extended micromechanics method to a viscoelastic polymer nanocomposite with interphase regions, and estimate the properties and thickness of the interphase region based on experimental data for carbon-black filled styrene butadiene rubbers.

42 ENGINEERING↗
DOE OSTI2026

A projection method for particle resampling

Particle discretizations of partial differential equations are advantageous for high-dimensional kinetic models in phase-space due to their better scalability than continuum approaches with respect to dimension. Complex processes collectively referred to as particle noise hamper long time simulations with particle methods. One approach to address this problem is particle mesh adaptivity, or remapping, known as particle resampling and remeshing. Here, this work introduces a resampling method that projects particles to and from a (finite element) function space. The method is simple, using standard sparse linear algebra and finite element techniques, and it preserves all moments up to the order of a polynomial represented exactly by the continuum function space. It is distinguished from most other mesh-based methods in that new particle positions and number are decoupled from the mesh, allowing particle and continuum meshes to be adapted relatively independently. While this work is developed with structured particle and continuum phase-space grids on 1X + 1V Vlasov-Poisson models of Landau damping and two-stream instability, the method is well-suited to unstructured grids. Stable long time dynamics are demonstrated up to time T = 500. Reproducibility artifacts and data are publicly available.

Kinetic methods↗
DOE OSTI2012

Effect of a short weak prepulse on laser-triggered front-surface heavy-ion acceleration

A suppression of light-ion acceleration (from surface water contaminants) was observed when a moderate-intensity subpicosecond laser pulse was focused on a thick metal target. Simultaneously, an effective generation of high-energy multicharge ions of the target material (Fe) was experimentally observed. A numerical simulation based on the Boltzmann-Vlasov-Poisson model revealed that this is due to the very specific regime of cleaning contaminants from the target surface by the short weak prepulse preceding the main pulse by more than 10 ns and having an intensity below the surface breakdown threshold. Because this prepulse causes the contaminant layer to boil explosively, a low-density gap forms above the target surface. These conditions are consequently favorable for boosting the energy of heavy ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Complete quasilinear model for the acceleration-driven lower hybrid drift instability and a computational assessment of its validity

A complete quasilinear model is derived for the electrostatic acceleration-driven lower hybrid drift instability in a uniform two-species low-beta plasma in which current is perpendicular to the background magnetic field. The model consists of coupled nonlinear velocity space diffusion equations for the volume-averaged ion and electron distribution functions. Each species' diffusion coefficient depends on a time-evolving spectral density of the electric-field energy per unit volume and a time-evolving dispersion relation. The dispersion relation is expressed analytically in integral form without the use of asymptotic limits and applies to arbitrary distribution functions, so long as they can be expressed as a function of one velocity coordinate, e.g., f⁡(vy) or f⁡(v⊥). The quasilinear model conserves energy and is complete in that it fully describes the evolution of the distribution functions, including resonant and nonresonant particle-wave interactions, while accounting for distribution-function-dependent mixed-complex frequencies. Further, the quasilinear diffusion model is solved numerically and self-consistently using a Crank-Nicolson temporal discretization and a second-order finite-volume velocity-space discretization. Numerical solutions are compared to nonlinear fourth-order accurate continuum kinetic Vlasov-Poisson simulations. Evolution of electric-field energy, growth rates, distribution functions, and diffusion coefficients are shown to be in agreement with Vlasov simulations. The quasilinear model is shown to predict anomalous transport terms, like resistivity and heating, to within a factor of order unity. Discrepancies between the quasilinear model and Vlasov simulations are assessed and attributed primarily to lack of damping in the quasilinear description and to the use of unperturbed-orbit susceptibilities in the linear theory dispersion relation. The results illuminate the predictive accuracy of the quasilinear model, place approximate bounds on its validity, and provide much needed vetting of quasilinear theory's ability to predict the nonlinear state of a microturbulent plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Scalable Matrix-Free Adaptive Product-Convolution Approximation for Locally Translation-Invariant Operators

We present an adaptive grid matrix-free operator approximation scheme based on a "product-convolution" interpolation of convolution operators. This scheme is appropriate for operators that are locally translation-invariant, even if these operators are high rank or full rank. Such operators arise in Schur complement methods for solving partial differential equations (PDEs), as Hessians in PDE-constrained optimization and inverse problems, as integral operators, as covariance operators, and as Dirichlet-to-Neumann maps. Constructing the approximation requires computing the impulse responses of the operator to point sources centered on nodes in an adaptively refined grid of sample points. A randomized a posteriori error estimator drives the adaptivity. Once constructed, the approximation can be efficiently applied to vectors using the fast Fourier transform. The approximation can be efficiently converted to hierarchical matrix (H-matrix) format, then inverted or factorized using scalable H-matrix arithmetic. The quality of the approximation degrades gracefully as fewer sample points are used, allowing cheap lower quality approximations to be used as preconditioners. Furthermore, this yields an automated method to construct preconditioners for locally translation-invariant Schur complements. We directly address issues related to boundaries and prove that our scheme eliminates boundary artifacts. We test the scheme on a spatially varying blurring kernel, on the nonlocal component of an interface Schur complement for the Poisson operator, and on the data misfit Hessian for an advection dominated advection-diffusion inverse problem. Numerical results show that the scheme outperforms existing methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Constant residual electrostatic electron plasma mode in Vlasov-Ampere system

In a collisionless Vlasov-Poisson (V-P) electron plasma system, two types of modes for electric field perturbation exist: the exponentially Landau damped electron plasma waves and the initial-value sensitive ballistic modes. Here, the V-P system is modified slightly to a Vlasov-Ampere (V-A) system. A new constant residual mode is revealed. Mathematically, this mode comes from the Laplace transform of an initial electric field perturbation, and physically represents that an initial perturbation (e.g., external electric field perturbation) would not be damped away. Thus, this residual mode is more difficult to be damped than the ballistic mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Parameterized anomalous transport model for current-carrying collisionless plasmas in pulsed power inertial confinement fusion

Current delivery in pulsed power inertial confinement fusion is influenced by collisionless current-carrying microturbulent plasmas, which are sourced from electrode surfaces. In this setting, the lower hybrid drift instability—triggered by plasma acceleration—is a leading candidate driver of difficult-to-predict momentum and energy transport. To characterize the nonlinear state of the microturbulent plasma, a parameterized anomalous transport model is developed for the instability, with analytic formulas for anomalous collision frequency, resistivity, and species heating rates. The formulas are expressed in terms of linear-theory variables and four dimensionless parameters that characterize the macroscopic plasma state. The model is built on linear theory analysis, power law analysis, and quasilinear theory analysis, and is validated using a series of nonlinear continuum kinetic Vlasov–Poisson simulations. The theoretical and computational investigation demonstrates that the anomalous collision frequency associated with the instability can be reliably approximated, within about a factor of five or better, by the unscaled linear theory growth rate of the fastest-growing wavenumber mode. This finding enables efficient calculation of anomalous resistivity and species heating rates over a wide range of plasma conditions, resulting in improved predictive capabilities.

Complex functions↗
DOE OSTI2021

Properties of giant dipole resonances within an extended pairing model with a focus on spectral statistics

In this paper, we report on a study of the spectral features associated with dipole resonances in medium mass nuclei (E ≤ 12 MeV), as revealed in the framework of the spd-interacting boson model. The effect of pairing correlations on the theory follows from solutions obtained through an application of the Bethe Ansatz Equation. In general, calculated spectra around the critical point of the vibrational to γ-soft transitions appears to approach that of a Gaussian Orthogonal Ensemble, while near the rotational and vibrational limits of the theory the spectra show more regular behavior. Specifically, the results reveal that the statistical features of the spectra are sensitive to the vector boson pairing strength, c p , in the transition region; that is, when c p is zero, or when the system approaches one of its dynamical symmetries limits the spectrum display regular features, while for stronger c p values, or when near to the critical phase transition region, the spectral feature show more chaotic behavior. Overall, our results indicate that the statistical features are governed by the interplay between dipole resonant energies, pairing correlations, and interactions between and among the single and vector bosons modes of the theory. As part of this work we also found out that chaoticity occurs when results were fit to a Berry-Robnik distribution. Throughout our analyses, we used experimentally known information about both positive and negative parity states. Our findings suggest that dipole resonances appear to be best-described by Poisson statistics for A ≈ 32-138 nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

An asymptotic preserving automatic domain decomposition method for the Vlasov–Poisson–BGK system with applications to plasmas

In this work we present an efficient strategy to deal with plasma physics simulations in which localized departures from thermodynamical equilibrium are present. The method is based on the introduction of intermediate regions which allows smooth transitions between kinetic and fluid zones. In this paper we extend Domain Decomposition techniques, obtained through dynamic coupling and buffer zones, to the study of plasmas and, moreover, we combine them with Asymptotic Preserving and Asymptotically Accurate strategies for the time integration. We use a hybrid scheme in which both kinetic and fluid descriptions are considered and coupled together while the kinetic model is solved by asymptotic preserving and accurate methods, in order to guarantee high efficiency and accuracy in all regimes. The numerical scheme is validated and its performances are analyzed by numerical simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗