Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solution regularity”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Physically-based parameterization of spatially variable soil and vegetation using satellite multispectral data

A stochastic-geometric landsurface reflectance model is formulated and tested for the parameterization of spatially variable vegetation and soil at subpixel scales using satellite multispectral images without ground truth. Landscapes are conceptualized as 3-D Lambertian reflecting surfaces consisting of plant canopies, represented by solid geometric figures, superposed on a flat soil background. A computer simulation program is developed to investigate image characteristics at various spatial aggregations representative of satellite observational scales, or pixels. The evolution of the shape and structure of the red-infrared space, or scattergram, of typical semivegetated scenes is investigated by sequentially introducing model variables into the simulation. The analytical moments of the total pixel reflectance, including the mean, variance, spatial covariance, and cross-spectral covariance, are derived in terms of the moments of the individual fractional cover and reflectance components. The moments are applied to the solution of the inverse problem: The estimation of subpixel landscape properties on a pixel-by-pixel basis, given only one multispectral image and limited assumptions on the structure of the landscape. The landsurface reflectance model and inversion technique are tested using actual aerial radiometric data collected over regularly spaced pecan trees, and using both aerial and LANDSAT Thematic Mapper data obtained over discontinuous, randomly spaced conifer canopies in a natural forested watershed. Different amounts of solar backscattered diffuse radiation are assumed and the sensitivity of the estimated landsurface parameters to those amounts is examined.

Jasinski, Michael F.↗
NASA NTRS2008

Interpolating Fields of Carbon Monoxide Data Using a Hybrid Statistical-Physical Model

Atmospheric Carbon Monoxide (CO) is a pollutant gas of which the US congress has mandated regular monitoring, and satellite sensors can be used to retrieve regional concentrations of CO over several vertical layers. However, CO at cloudy locations cannot be observed and have to be estimated from the observed data set, resulting in an interpolation problem. The current state-of-the-art solution is to combine prior information, computed by a deterministic physical model, with observations. However, the deterministic model may introduce uncertainties that do not derive from the data. While sharing certain features with the physical model, this paper presents a Bayesian hierarchical model for interpolating CO on a 3-dimensional spatial grid, across time. To our knowledge such a model has not been considered before. The model is applied to a hypothetical air-quality monitoring scenario, and is compared to existing interpolation methods. The results provide motivation for the use of the statistical model for regional to local applications.

Arellano, A. A.↗
NASA NTRSDate not supplied

Navigation for the ACS3 Solar Sail Mission

NASA’s Advanced Composite Solar Sail System (ACS3) mission consists of a spacecraft that plans to be launched in early 2024. The spacecraft carries an 80 m2 solar sail that can produce effective thrust to alter the initial 1000 km sun-synchronous orbit. The main objective of the mission is to demonstrate the capabilities of the solar sail to effectively change various orbital parameters such as the semi-major axis and the inclination. Various composite materials were used to produce the sail, together with lightweight booms that will deploy from a 12U CubeSat bus. The astrodynamics team at NASA Ames Research Center has built a Flight Dynamics System (FDS) to provide mission navigation and to produce regular ephemeris once in orbit. The FDS can compute the trajectories that the spacecraft will perform once the sail is deployed. To achieve that, GPS data is obtained from the spacecraft telemetry and then is used with a Kalman filter and a smoother to obtain an orbit determination solution. The outcome of this process reduces the position and velocity uncertainty in a daily cadence. After that, the state vector output is used to feed a propagation model that includes the updated attitude and orbit of the spacecraft at that given moment. The trajectory model considers the updated attitude plan of the spacecraft as well as the environment conditions such as the solar weather to compute the associated drag and solar radiation pressure. This paper explains in detail the implementation of the FDS, as well as the solar sail solar radiation pressure trajectory model. We also present the results of several potential trajectory models under various assumptions of orbit parameters, attitude, environment, and material properties. In addition, we introduce a trajectory model for potential interplanetary use of an equivalent solar sail in the future. The outcome of this process will be critical to achieving the objective of determining effective altitude change produced by the solar sail.

Andres Dono Perez↗
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 OSTI2023

A Lizard-inspired Tube Inspector (LTI) Robot

This final report summarizes the research findings of the current project. This project is a collaboration between New Mexico State University (NMSU) as a lead (Recipient) and Arizona State University (ASU) as a Co-Recipient. Tubular structures are common components of boilers and heat exchangers in power plants. Over time, these components may suffer corrosion, cracks, and stress-corrosion cracks in either the body or the welded connections. A single tube leakage can cause an outage of several weeks. Regular inspection is a key safety factor when ensuring that power plants are maintained in reliable, operational condition. This inspection, however, is challenging, time-consuming, and in many cases impossible due to accessibility issues and safety concerns. Recent developments in robotic-based inspection can offer a great solution. Hard-to-reach places can be inspected without overhauling the unit, saving time and cost. Although several robots have been designed and implemented successfully for inspecting power plant components, in particular tubular structures, their mobility and flexibility are limited. Most of these robots use wheels for mobility which reduces their maneuverability of these robots. Moreover, they usually use magnets to attach to tubes which will not work on tubes with non-ferromagnetic materials. These robots usually carry measurement tools for a particular non-destructive testing (NDT) method such as ultrasound testing (UT) that requires a couplant to perform a point-by-point (scanning) inspection of the tubular structure. In this project, we developed a versatile lizard-inspired tube inspector (LTI) robot with embedded inspection sensing components for tube inspection which removes the need for point-by-point scanning of tube surface for crack and corrosion detection. Inspired by a “lizard”, the novelty of the current project is the integration of couplant-free ultrasound sensing and transmission, advanced data-driven defect detection and imaging, and friction-based mechanical mobility components in a single robot to eliminate a need for smooth surfaces and simple geometry for mobility and scanning. The LTI robot could replace the wheel-based approach with friction-based mobility to significantly increase the flexibility and maneuverability of the robot. The LTI robot can get into a power plant unit, such as a boiler, from a small area allowing it to access a component of interest for inspection (e.g., move on OD, curved and flat surfaces, non-ferromagnetic or ferromagnetic materials, and tubes with rough surfaces and complex geometries). Additionally, an advanced data-driven method using ultrasound data was pioneered to allow the robot to detect defects in the entire area between the robot’s multi-functional mobility system (grippers). Integrating a couplant-free ultrasound sensing in the robot’s grippers as well as using advanced data-driven methods allowed the LTI to detect defects in the entire cross-section of a tube using its grippers when stationary and when mobile.

20 FOSSIL-FUELED POWER PLANTS↗
NASA NTRS1987

Multigrid methods for bifurcation problems: The self adjoint case

This paper deals with multigrid methods for computational problems that arise in the theory of bifurcation and is restricted to the self adjoint case. The basic problem is to solve for arcs of solutions, a task that is done successfully with an arc length continuation method. Other important issues are, for example, detecting and locating singular points as part of the continuation process, switching branches at bifurcation points, etc. Multigrid methods have been applied to continuation problems. These methods work well at regular points and at limit points, while they may encounter difficulties in the vicinity of bifurcation points. A new continuation method that is very efficient also near bifurcation points is presented here. The other issues mentioned above are also treated very efficiently with appropriate multigrid algorithms. For example, it is shown that limit points and bifurcation points can be solved for directly by a multigrid algorithm. Moreover, the algorithms presented here solve the corresponding problems in just a few work units (about 10 or less), where a work unit is the work involved in one local relaxation on the finest grid.

Taasan, Shlomo↗
NASA NTRS1983

Use of a personal computer for dynamical engineering illustrations in a classroom and over an instructional TV network

A personal computer has been used to illustrate physical phenomena and problem solution techniques in engineering classes. According to student evaluations, instruction of concepts was greatly improved through the use of these illustrations. This paper describes the class of phenomena that can be effectively illustrated, the techniques used to create these illustrations, and the techniques used to display the illustrations in regular classrooms and over an instructional TV network. The features of a personal computer required to apply these techniques are listed. The capabilities of some present personal computers are discussed and a forecast of the capabilities of future personal computers is presented.

Watson, V. R.↗
NASA NTRSDate not supplied

Navigation for the ACS3 Solar Sail Mission

NASA’s Advanced Composite Solar Sail System (ACS3) mission consists of a spacecraft that plans to be launched no earlier than April 2024. The spacecraft carries an 81 m2 solar sail that can produce effective ΔV to alter the initial 1000 km sun-synchronous orbit. The main objective of the mission is to demonstrate the capabilities of the solar sail to effectively change the semi-major axis of its initial orbit. Various composite materials were used to produce the sail, together with lightweight booms that will deploy from a 12U CubeSat bus. The ACS3 navigation team at NASA Ames Research Center has built a Flight Dynamics System (FDS) to provide mission navigation and to produce regular ephemeris once in orbit. The FDS can compute the orbit transfers that the spacecraft will perform once the sail is deployed. To achieve that, GPS data is obtained from the spacecraft telemetry and then is used with a Kalman filter and a smoother to obtain an orbit determination solution. The outcome of this process reduces the position and velocity uncertainty in a daily cadence. After that, the state vector output is used to feed a propagation model that includes the updated attitude and orbit of the spacecraft at that given moment. The trajectory model considers the updated attitude plan of the spacecraft as well as the environment conditions such as the solar weather to compute the associated drag and solar radiation pressure. This paper explains in detail the implementation of the FDS, as well as the solar sail solar radiation pressure trajectory model. We also present the results of several potential trajectory models under various assumptions of orbit parameters, attitude, environment, and material properties. In addition, we introduce a trajectory model for potential interplanetary use of an equivalent solar sail in the future. The outcome of this process will be critical to achieving the objective of determining effective semi-major axis change produced by the solar sail.

Andres Dono↗
DOE OSTI2022

ORNL_AISD_NiNb

This dataset describes the nickel-niobium solid solution binary alloy, where the two constituent elements nickel (Ni) and niobium (Nb) are randomly placed on an underlying crystal lattice. This dataset for nickel-niobium (Ni-Nb) alloys available includes the formation energy and bulk modulus for each crystal structure. Each atomic sample has a disordered phase which is obtained starting from an initial regular crystal structure of type body-centered cubic (BCC), face-centered cubic (FCC), or hexagonal compact packed (HCP). The geometry optimization ensures that all the alloy samples reached the equilibrium with negative formation energy. We perform geometry optimizations using the LAMMPS simulation package [1], a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. We utilized the embedded atom model (EAM) potential for Ni and Nb developed in a previous study [2]. The potential could describe behaviors of the liquid and solid phases of Ni-Nb alloy. The structural factors and angular distributions of three atoms are well-matched with X-ray and ab initio-based molecular dynamics data. We prepared the three different crystals with different initial lattice parameters (3.52 Ã… for FCC, 3.32 Ã… for BCC, and 3.5 Ã… for HCP). We performed energy minimization in two steps. Firstly, we minimized the structures with an isotropic unit cell to minimize the side effects from our arbitrary lattice parameters for all other compositions. Then, we applied geometry optimization with a triclinic (non-orthogonal) unit cell to fully minimize the stress components to calculate the elastic constants. In this procedure, we chose 10,000 as the maximum number of allowable steps aimed at obtaining fully relaxed atomic geometries. The dataset consists of three sets of crystal structures. The first set contains 46,086 irregular crystal structures, each of them with 54 atoms, obtained through optimization starting from a regular BCC crystal structure. The second set contains 24,543 irregular crystal structures, each of them with 32 atoms, obtained through optimization starting from a regular FCC crystal structure. The third set contains 39,303 irregular crystal structures, each of them with 48 atoms, obtained through optimization starting from a regular HCP crystal structure. The atomic configurations within each set span the possible compositional range. The three sets have been unified in a global dataset, which is extremely heterogeneous in terms of crystal structures, lattice volumes, and atomic configurations. Organization of files inside the dataset: the dataset contains three subdirectories called • BCC_opt • FCC_opt • HCP_opt based on the type of initial regular structure used to start the geometry optimization. Inside each of these folders, every atomic structure is identified by a string “A_B_Câ€, where A denotes the number of Nb in the system, B denotes index of structure with a given Nb number, and C denotes the total number of structures generated with a given Nb number. For each optimized crystal structure identified by the unique string of characters “A_B_Câ€, three files are provided: • A_B_C_opt.xyz: The optimized geometries in xyz format • A_B_C_opt.cfg: The optimized geometries in cfg format. It includes cell information and atomic energy, and forces calculated from LAMMPS. • A_B_C.elastic: Raw data of 21 elastic constants from LAMMPS output. • A_B_C.bulk: Calculated upper and lower bounds of bulk modulus and averaged one based on Voigt-Reuss-Hill approach from *.elastic. References: [1] A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in 't Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comp. Phys. Comm., 271:108171, 2022. [2] Y Zhang, R Ashcraft, MI Mendelev, CZ Wang, and KF Kelton. Experimental and molecular dynamics simulation study of structure of liquid and amorphous ni62nb38 alloy. The Journal of chemical physics, 145(20):204505, 2016.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A phase-field model of quasi-brittle fracture for pressurized cracks: Application to UO 2 high-burnup microstructure fragmentation

Here, in this paper, we present a phase-field model of quasi-brittle fracture with pressurized cracks, with dedicated applications for polycrystalline materials. The model is formulated as a minimization problem within the variational framework. The external work done by pressure on the crack surfaces is included in the objective function. Several careful modeling choices lead to a regularization-length-independent critical strength. The model is constructed to give a softening response with an underlying linear traction-separation law. The pressure-dependent softening response and the regularization of the prescribed pressure are demonstrated with a (quasi) one-dimensional numerical analysis. In a two-dimensional numerical analysis under plane strain assumptions, the critical stress corresponding to crack propagation (as predicted by our quasi-brittle fracture model) is compared with linear elastic fracture mechanics (LEFM) analytical solutions. Our model is further utilized to simulate fission-gas-induced fragmentation of the UO 2 high-burnup structure (HBS). Simulation results show that pressurized bubbles can cause crack nucleation and propagation, and that the bubble size and the surrounding external pressure affect the critical pressure corresponding to crack nucleation. Simulations of a partial HBS at different recrystallization stages show that different grain structures (due to recrystallization) also influence crack paths and fragmentation morphology.

36 MATERIALS SCIENCE↗
NASA NTRS2010

Improving Global Mass Flux Solutions from Gravity Recovery and Climate Experiment (GRACE) Through Forward Modeling and Continuous Time Correlation

We describe Earth's mass flux from April 2003 through November 2008 by deriving a time series of mas cons on a global 2deg x 2deg equal-area grid at 10 day intervals. We estimate the mass flux directly from K band range rate (KBRR) data provided by the Gravity Recovery and Climate Experiment (GRACE) mission. Using regularized least squares, we take into account the underlying process dynamics through continuous space and time-correlated constraints. In addition, we place the mascon approach in the context of other filtering techniques, showing its equivalence to anisotropic, nonsymmetric filtering, least squares collocation, and Kalman smoothing. We produce mascon time series from KBRR data that have and have not been corrected (forward modeled) for hydrological processes and fmd that the former produce superior results in oceanic areas by minimizing signal leakage from strong sources on land. By exploiting the structure of the spatiotemporal constraints, we are able to use a much more efficient (in storage and computation) inversion algorithm based upon the conjugate gradient method. This allows us to apply continuous rather than piecewise continuous time-correlated constraints, which we show via global maps and comparisons with ocean-bottom pressure gauges, to produce time series with reduced random variance and full systematic signal. Finally, we present a preferred global model, a hybrid whose oceanic portions are derived using forward modeling of hydrology but whose land portions are not, and thus represent a pure GRACE-derived signal.

Sabaka, T. J.↗
DOE OSTI2022

A nonsmooth nonconvex optimization algorithm for two-stage optimization problems

An optimization algorithm for a group of nonsmooth nonconvex problems inspired by two-stage stochastic programming problems is proposed. The main challenges for these problems include (1) the problems lack the popular lower-type properties such as prox-regularity assumed in many nonsmooth nonconvex optimization algorithms, (2) the objective can not be analytically expressed and (3) the evaluation of function values and subgradients are computationally expensive. To address these challenges, this report first examines the properties that exist in many two-stage problems, specifically upper-C 2 objectives. Then, we show that quadratic penalty method for securityconstrained alternating current optimal power flow (SCACOPF) contingency problems can make the contingency solution functions upper-C 2 . Based on these observations, a simplified bundle algorithm that bears similarity to sequential quadratic programming (SQP) method is proposed. It is more efficient in implementation and computation compared to conventional bundle methods. Global convergence analysis of the algorithm is presented under novel and reasonable assumptions. The proposed algorithm therefore fills the gap of theoretical convergence for smoothed SCACOPF problems. The inconsistency that might arise in our treatment of the constraints are addressed through a penalty algorithm whose convergence analysis is also provided. Finally, theoretical capabilities and numerical performance of the algorithm are demonstrated through numerical examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Transmission electron microscopy of the rapid solidification microstructure evolution and solidification interface velocity determination in hypereutectic Al-20at.%Cu after laser melting

The evolution of rapid solidification (RS) microstructure and solidification interface velocity has been studied experimentally by in-situ transmission electron microscopy and postmortem characterization for hypereutectic Al-20at.%Cu (37 wt.%Cu) after laser melting. Four morphologically distinct regions formed via growth modes changing from θ-Al 2 Cu phase dendrites to eutectic cell growth, possibly α-Al dendrite growth (α-cell), banded growth, and α-Al plane front growth. Tendency for faceting and low capacity for solute trapping limited the θ-phase dendrite growth to solidification interface velocity v < 0.07 m/s. Consequently, formation of pro-eutectic micro-constituent was suppressed, and RS microstructure formation was dominated by eutectic, α-cell, and banded morphology grains for the Al-20at.%Cu alloy. Eutectic growth operated for interface velocity of 0.1 m/s ≤ v < 0.3 m/s, with a transition from regular lamellar, 2-λ and 1-λ mode to a dense irregular morphology dominated by α-phase at v = 0.3 m/s. Interface temperature calculations indicated feasibility of α-cell growth mode for 0.3 m/s ≤ v < 0.7 m/s. Banded growth occurred for 0.7 m/s ≤ v < 1.3 m/s. Plane front α-phase growth was evident for interface velocities v ≥ 1.3 m/s. Previous work on rapid solidification microstructure development in hypereutectic Al-Cu alloys reported a regime of α-cell growth subsequently to eutectic and prior to transition to banded growth for Al-19at.%Cu (36 wt.%Cu), while for Al-22at.%Cu (40 wt.%Cu) eutectic growth transitioned directly to α-plane front growth without emergence of a banded regime. Finally, based on the current study the disappearance of the banded growth regime occurs for a composition larger than 20at.%Cu.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Lattice regularization of reduced Kähler-Dirac fermions and connections to chiral fermions

We show how a path integral for reduced Kähler-Dirac fermions suffers from a phase ambiguity associated with the fermion measure that is an analog of the measure problem seen for chiral fermions. However, unlike the case of chiral symmetry, a doubler free lattice action exists which is invariant under the corresponding onsite symmetry. This allows for a clear diagnosis and solution to the problem using mirror fermions resulting in a unique gauge invariant measure. By introducing an appropriate set of Yukawa interactions which are consistent with ’t Hooft anomaly cancellation we conjecture the mirrors can be decoupled from low energy physics. Moreover, the minimal such Kähler-Dirac mirror model yields a light sector which corresponds, in the flat space continuum limit, to the Pati-Salam GUT model.

Physics↗
DOE OSTI2021

Optimization-based algorithms for nonlinear mechanics and frictional contact

An optimization-based strategy for solving nonlinear mechanics problems is proposed. In contrast to typical nonlinear equation solver algorithms that aim to find zeros in the residual force function, we minimize an energy (or energy-like) function to encourage solutions which are locally stable equilibria. These smooth and potentially non-convex objective functions are minimized using a preconditioned conjugate-gradient trust-region algorithm. Contact is formulated as an inequality constrained minimization problem, and is solved with an augmented Lagrangian algorithm. Friction is included in the approach via a regularized quasi-potential energy, and other dissipative behavior is included through the use of variational constitutive updates. Finally, to accelerate convergence rates for the Lagrange multipliers, we propose a novel multiplier update algorithm utilizing the Fischer-Burmeister function, and demonstrate super-linear solver convergence for some applications.

42 ENGINEERING↗
NASA NTRS1989

Thin-plate spline quadrature of geodetic integrals

Thin-plate spline functions (known for their flexibility and fidelity in representing experimental data) are especially well-suited for the numerical integration of geodetic integrals in the area where the integration is most sensitive to the data, i.e., in the immediate vicinity of the evaluation point. Spline quadrature rules are derived for the contribution of a circular innermost zone to Stoke's formula, to the formulae of Vening Meinesz, and to the recursively evaluated operator L(n) in the analytical continuation solution of Molodensky's problem. These rules are exact for interpolating thin-plate splines. In cases where the integration data are distributed irregularly, a system of linear equations needs to be solved for the quadrature coefficients. Formulae are given for the terms appearing in these equations. In case the data are regularly distributed, the coefficients may be determined once-and-for-all. Examples are given of some fixed-point rules. With such rules successive evaluation, within a circular disk, of the terms in Molodensky's series becomes relatively easy. The spline quadrature technique presented complements other techniques such as ring integration for intermediate integration zones.

Vangysen, Herman↗
DOE OSTI2020

On the Solution of ℓ 0 -Constrained Sparse Inverse Covariance Estimation Problems

The sparse inverse covariance matrix is used to model conditional dependencies between variables in a graphical model to fit a multivariate Gaussian distribution. Estimating the matrix from data are well known to be computationally expensive for large-scale problems. Sparsity is employed to handle noise in the data and to promote interpretability of a learning model. Although the use of a convex ℓ 1 regularizer to encourage sparsity is common practice, the combinatorial ℓ 0 penalty often has more favorable statistical properties. In this paper, we directly constrain sparsity by specifying a maximally allowable number of nonzeros, in other words, by imposing an ℓ 0 constraint. Here, we introduce an efficient approximate Newton algorithm using warm starts for solving the nonconvex ℓ 0 -constrained inverse covariance learning problem. Numerical experiments on standard data sets show that the performance of the proposed algorithm is competitive with state-of-the-art methods.

$\ell_0$-Constrained↗
NASA NTRS1988

Existence, uniqueness and regularity of a time-periodic probability density distribution arising in a sedimentation-diffusion problem

The sedimentation and diffusion of a nonneutrally buoyant Brownian particle in vertical fluid-filled cylinder of finite length which is instantaneously inverted at regular intervals are investigated analytically. A one-dimensional convective-diffusive equation is derived to describe the temporal and spatial evolution of the probability density; a periodicity condition is formulated; the applicability of Fredholm theory is established; and the parameter-space regions are determined within which the existence and uniqueness of solutions are guaranteed. Numerical results for sample problems are presented graphically and briefly characterized.

Nitsche, Ludwig C.↗