An unconditionally stable, positivity-preserving splitting scheme for nonlinear Black-Scholes equation with transaction costs.

PubMed

Guo, Jianqiang; Wang, Wansheng

2014-01-01

This paper deals with the numerical analysis of nonlinear Black-Scholes equation with transaction costs. An unconditionally stable and monotone splitting method, ensuring positive numerical solution and avoiding unstable oscillations, is proposed. This numerical method is based on the LOD-Backward Euler method which allows us to solve the discrete equation explicitly. The numerical results for vanilla call option and for European butterfly spread are provided. It turns out that the proposed scheme is efficient and reliable.

An improved flux-split algorithm applied to hypersonic flows in chemical equilibrium

NASA Technical Reports Server (NTRS)

Palmer, Grant

1988-01-01

An explicit, finite-difference, shock-capturing numerical algorithm is presented and applied to hypersonic flows assumed to be in thermochemical equilibrium. Real-gas chemistry is either loosely coupled to the gasdynamics by way of a Gibbs free energy minimization package or fully coupled using species mass conservation equations with finite-rate chemical reactions. A scheme is developed that maintains stability in the explicit, finite-rate formulation while allowing relatively high time steps. The codes use flux vector splitting to difference the inviscid fluxes and employ real-gas corrections to viscosity and thermal conductivity. Numerical results are compared against existing ballistic range and flight data. Flows about complex geometries are also computed.

Solution of the Average-Passage Equations for the Incompressible Flow through Multiple-Blade-Row Turbomachinery

DTIC Science & Technology

1994-02-01

numerical treatment. An explicit numerical procedure based on Runqe-Kutta time stepping for cell-centered, hexahedral finite volumes is...An explicit numerical procedure based on Runge-Kutta time stepping for cell-centered, hexahedral finite volumes is outlined for the approximate...Discretization 16 3.1 Cell-Centered Finite -Volume Discretization in Space 16 3.2 Artificial Dissipation 17 3.3 Time Integration 21 3.4 Convergence

Exact Closed-form Solutions for Lamb's Problem

NASA Astrophysics Data System (ADS)

Feng, Xi; Zhang, Haiming

2018-04-01

In this article, we report on an exact closed-form solution for the displacement at the surface of an elastic half-space elicited by a buried point source that acts at some point underneath that surface. This is commonly referred to as the 3-D Lamb's problem, for which previous solutions were restricted to sources and receivers placed at the free surface. By means of the reciprocity theorem, our solution should also be valid as a means to obtain the displacements at interior points when the source is placed at the free surface. We manage to obtain explicit results by expressing the solution in terms of elementary algebraic expression as well as elliptic integrals. We anchor our developments on Poisson's ratio 0.25 starting from Johnson's (1974) integral solutions which must be computed numerically. In the end, our closed-form results agree perfectly with the numerical results of Johnson (1974), which strongly confirms the correctness of our explicit formulas. It is hoped that in due time, these formulas may constitute a valuable canonical solution that will serve as a yardstick against which other numerical solutions can be compared and measured.

Exact closed-form solutions for Lamb's problem

NASA Astrophysics Data System (ADS)

Feng, Xi; Zhang, Haiming

2018-07-01

In this paper, we report on an exact closed-form solution for the displacement at the surface of an elastic half-space elicited by a buried point source that acts at some point underneath that surface. This is commonly referred to as the 3-D Lamb's problem for which previous solutions were restricted to sources and receivers placed at the free surface. By means of the reciprocity theorem, our solution should also be valid as a means to obtain the displacements at interior points when the source is placed at the free surface. We manage to obtain explicit results by expressing the solution in terms of elementary algebraic expression as well as elliptic integrals. We anchor our developments on Poisson's ratio 0.25 starting from Johnson's integral solutions which must be computed numerically. In the end, our closed-form results agree perfectly with the numerical results of Johnson, which strongly confirms the correctness of our explicit formulae. It is hoped that in due time, these formulae may constitute a valuable canonical solution that will serve as a yardstick against which other numerical solutions can be compared and measured.

Development of the Semi-implicit Time Integration in KIM-SH

NASA Astrophysics Data System (ADS)

NAM, H.

2015-12-01

The Korea Institute of Atmospheric Prediction Systems (KIAPS) was founded in 2011 by the Korea Meteorological Administration (KMA) to develop Korea's own global Numerical Weather Prediction (NWP) system as nine year (2011-2019) project. The KIM-SH is a KIAPS integrated model-spectral element based in the HOMME. In KIM-SH, the explicit schemes are employed. We introduce the three- and two-time-level semi-implicit scheme in KIM-SH as the time integration. Explicit schemes however have a tendancy to be unstable and require very small timesteps while semi-implicit schemes are very stable and can have much larger timesteps.We define the linear and reference values, then by definition of semi-implicit scheme, we apply the linear solver as GMRES. The numerical results from experiments will be introduced with the current development status of the time integration in KIM-SH. Several numerical examples are shown to confirm the efficiency and reliability of the proposed schemes.

Revisiting analytical solutions for steady interface flow in subsea aquifers: Aquitard salinity effects

NASA Astrophysics Data System (ADS)

Werner, Adrian D.; Robinson, Neville I.

2018-06-01

Existing analytical solutions for the distribution of fresh groundwater in subsea aquifers presume that the overlying offshore aquitard, represented implicitly, contains seawater. Here, we consider the case where offshore fresh groundwater is the result of freshwater discharge from onshore aquifers, and neglect paleo-freshwater sources. A recent numerical modeling investigation, involving explicit simulation of the offshore aquitard, demonstrates that offshore aquitards more likely contain freshwater in areas of upward freshwater leakage to the sea. We integrate this finding into the existing analytical solutions by providing an alternative formulation for steady interface flow in subsea aquifers, whereby the salinity in the offshore aquitard can be chosen. The new solution, taking the aquitard salinity as that of freshwater, provides a closer match to numerical modeling results in which the aquitard is represented explicitly.

High-Order Implicit-Explicit Multi-Block Time-stepping Method for Hyperbolic PDEs

NASA Technical Reports Server (NTRS)

Nielsen, Tanner B.; Carpenter, Mark H.; Fisher, Travis C.; Frankel, Steven H.

2014-01-01

This work seeks to explore and improve the current time-stepping schemes used in computational fluid dynamics (CFD) in order to reduce overall computational time. A high-order scheme has been developed using a combination of implicit and explicit (IMEX) time-stepping Runge-Kutta (RK) schemes which increases numerical stability with respect to the time step size, resulting in decreased computational time. The IMEX scheme alone does not yield the desired increase in numerical stability, but when used in conjunction with an overlapping partitioned (multi-block) domain significant increase in stability is observed. To show this, the Overlapping-Partition IMEX (OP IMEX) scheme is applied to both one-dimensional (1D) and two-dimensional (2D) problems, the nonlinear viscous Burger's equation and 2D advection equation, respectively. The method uses two different summation by parts (SBP) derivative approximations, second-order and fourth-order accurate. The Dirichlet boundary conditions are imposed using the Simultaneous Approximation Term (SAT) penalty method. The 6-stage additive Runge-Kutta IMEX time integration schemes are fourth-order accurate in time. An increase in numerical stability 65 times greater than the fully explicit scheme is demonstrated to be achievable with the OP IMEX method applied to 1D Burger's equation. Results from the 2D, purely convective, advection equation show stability increases on the order of 10 times the explicit scheme using the OP IMEX method. Also, the domain partitioning method in this work shows potential for breaking the computational domain into manageable sizes such that implicit solutions for full three-dimensional CFD simulations can be computed using direct solving methods rather than the standard iterative methods currently used.

Lattice Boltzmann model for the compressible Navier-Stokes equations with flexible specific-heat ratio.

PubMed

Kataoka, Takeshi; Tsutahara, Michihisa

2004-03-01

We have developed a lattice Boltzmann model for the compressible Navier-Stokes equations with a flexible specific-heat ratio. Several numerical results are presented, and they agree well with the corresponding solutions of the Navier-Stokes equations. In addition, an explicit finite-difference scheme is proposed for the numerical calculation that can make a stable calculation with a large Courant number.

Implicit and explicit subgrid-scale modeling in discontinuous Galerkin methods for large-eddy simulation

NASA Astrophysics Data System (ADS)

Fernandez, Pablo; Nguyen, Ngoc-Cuong; Peraire, Jaime

2017-11-01

Over the past few years, high-order discontinuous Galerkin (DG) methods for Large-Eddy Simulation (LES) have emerged as a promising approach to solve complex turbulent flows. Despite the significant research investment, the relation between the discretization scheme, the Riemann flux, the subgrid-scale (SGS) model and the accuracy of the resulting LES solver remains unclear. In this talk, we investigate the role of the Riemann solver and the SGS model in the ability to predict a variety of flow regimes, including transition to turbulence, wall-free turbulence, wall-bounded turbulence, and turbulence decay. The Taylor-Green vortex problem and the turbulent channel flow at various Reynolds numbers are considered. Numerical results show that DG methods implicitly introduce numerical dissipation in under-resolved turbulence simulations and, even in the high Reynolds number limit, this implicit dissipation provides a more accurate representation of the actual subgrid-scale dissipation than that by explicit models.

An Explicit Upwind Algorithm for Solving the Parabolized Navier-Stokes Equations

NASA Technical Reports Server (NTRS)

Korte, John J.

1991-01-01

An explicit, upwind algorithm was developed for the direct (noniterative) integration of the 3-D Parabolized Navier-Stokes (PNS) equations in a generalized coordinate system. The new algorithm uses upwind approximations of the numerical fluxes for the pressure and convection terms obtained by combining flux difference splittings (FDS) formed from the solution of an approximate Riemann (RP). The approximate RP is solved using an extension of the method developed by Roe for steady supersonic flow of an ideal gas. Roe's method is extended for use with the 3-D PNS equations expressed in generalized coordinates and to include Vigneron's technique of splitting the streamwise pressure gradient. The difficulty associated with applying Roe's scheme in the subsonic region is overcome. The second-order upwind differencing of the flux derivatives are obtained by adding FDS to either an original forward or backward differencing of the flux derivative. This approach is used to modify an explicit MacCormack differencing scheme into an upwind differencing scheme. The second order upwind flux approximations, applied with flux limiters, provide a method for numerically capturing shocks without the need for additional artificial damping terms which require adjustment by the user. In addition, a cubic equation is derived for determining Vegneron's pressure splitting coefficient using the updated streamwise flux vector. Decoding the streamwise flux vector with the updated value of Vigneron's pressure splitting improves the stability of the scheme. The new algorithm is applied to 2-D and 3-D supersonic and hypersonic laminar flow test cases. Results are presented for the experimental studies of Holden and of Tracy. In addition, a flow field solution is presented for a generic hypersonic aircraft at a Mach number of 24.5 and angle of attack of 1 degree. The computed results compare well to both experimental data and numerical results from other algorithms. Computational times required for the upwind PNS code are approximately equal to an explicit PNS MacCormack's code and existing implicit PNS solvers.

Enforcing the Courant-Friedrichs-Lewy condition in explicitly conservative local time stepping schemes

NASA Astrophysics Data System (ADS)

Gnedin, Nickolay Y.; Semenov, Vadim A.; Kravtsov, Andrey V.

2018-04-01

An optimally efficient explicit numerical scheme for solving fluid dynamics equations, or any other parabolic or hyperbolic system of partial differential equations, should allow local regions to advance in time with their own, locally constrained time steps. However, such a scheme can result in violation of the Courant-Friedrichs-Lewy (CFL) condition, which is manifestly non-local. Although the violations can be considered to be "weak" in a certain sense and the corresponding numerical solution may be stable, such calculation does not guarantee the correct propagation speed for arbitrary waves. We use an experimental fluid dynamics code that allows cubic "patches" of grid cells to step with independent, locally constrained time steps to demonstrate how the CFL condition can be enforced by imposing a constraint on the time steps of neighboring patches. We perform several numerical tests that illustrate errors introduced in the numerical solutions by weak CFL condition violations and show how strict enforcement of the CFL condition eliminates these errors. In all our tests the strict enforcement of the CFL condition does not impose a significant performance penalty.

Transient modeling/analysis of hyperbolic heat conduction problems employing mixed implicit-explicit alpha method

NASA Technical Reports Server (NTRS)

Tamma, Kumar K.; D'Costa, Joseph F.

1991-01-01

This paper describes the evaluation of mixed implicit-explicit finite element formulations for hyperbolic heat conduction problems involving non-Fourier effects. In particular, mixed implicit-explicit formulations employing the alpha method proposed by Hughes et al. (1987, 1990) are described for the numerical simulation of hyperbolic heat conduction models, which involves time-dependent relaxation effects. Existing analytical approaches for modeling/analysis of such models involve complex mathematical formulations for obtaining closed-form solutions, while in certain numerical formulations the difficulties include severe oscillatory solution behavior (which often disguises the true response) in the vicinity of the thermal disturbances, which propagate with finite velocities. In view of these factors, the alpha method is evaluated to assess the control of the amount of numerical dissipation for predicting the transient propagating thermal disturbances. Numerical test models are presented, and pertinent conclusions are drawn for the mixed-time integration simulation of hyperbolic heat conduction models involving non-Fourier effects.

Explicit filtering in large eddy simulation using a discontinuous Galerkin method

NASA Astrophysics Data System (ADS)

Brazell, Matthew J.

The discontinuous Galerkin (DG) method is a formulation of the finite element method (FEM). DG provides the ability for a high order of accuracy in complex geometries, and allows for highly efficient parallelization algorithms. These attributes make the DG method attractive for solving the Navier-Stokes equations for large eddy simulation (LES). The main goal of this work is to investigate the feasibility of adopting an explicit filter in the numerical solution of the Navier-Stokes equations with DG. Explicit filtering has been shown to increase the numerical stability of under-resolved simulations and is needed for LES with dynamic sub-grid scale (SGS) models. The explicit filter takes advantage of DG's framework where the solution is approximated using a polyno- mial basis where the higher modes of the solution correspond to a higher order polynomial basis. By removing high order modes, the filtered solution contains low order frequency content much like an explicit low pass filter. The explicit filter implementation is tested on a simple 1-D solver with an initial condi- tion that has some similarity to turbulent flows. The explicit filter does restrict the resolution as well as remove accumulated energy in the higher modes from aliasing. However, the ex- plicit filter is unable to remove numerical errors causing numerical dissipation. A second test case solves the 3-D Navier-Stokes equations of the Taylor-Green vortex flow (TGV). The TGV is useful for SGS model testing because it is initially laminar and transitions into a fully turbulent flow. The SGS models investigated include the constant coefficient Smagorinsky model, dynamic Smagorinsky model, and dynamic Heinz model. The constant coefficient Smagorinsky model is over dissipative, this is generally not desirable however it does add stability. The dynamic Smagorinsky model generally performs better, especially during the laminar-turbulent transition region as expected. The dynamic Heinz model which is based on an improved model, handles the laminar-turbulent transition region well while also showing additional robustness.

Oscillations and stability of numerical solutions of the heat conduction equation

NASA Technical Reports Server (NTRS)

Kozdoba, L. A.; Levi, E. V.

1976-01-01

The mathematical model and results of numerical solutions are given for the one dimensional problem when the linear equations are written in a rectangular coordinate system. All the computations are easily realizable for two and three dimensional problems when the equations are written in any coordinate system. Explicit and implicit schemes are shown in tabular form for stability and oscillations criteria; the initial temperature distribution is considered uniform.

Numerics and subgrid-scale modeling in large eddy simulations of stratocumulus clouds.

PubMed

Pressel, Kyle G; Mishra, Siddhartha; Schneider, Tapio; Kaul, Colleen M; Tan, Zhihong

2017-06-01

Stratocumulus clouds are the most common type of boundary layer cloud; their radiative effects strongly modulate climate. Large eddy simulations (LES) of stratocumulus clouds often struggle to maintain fidelity to observations because of the sharp gradients occurring at the entrainment interfacial layer at the cloud top. The challenge posed to LES by stratocumulus clouds is evident in the wide range of solutions found in the LES intercomparison based on the DYCOMS-II field campaign, where simulated liquid water paths for identical initial and boundary conditions varied by a factor of nearly 12. Here we revisit the DYCOMS-II RF01 case and show that the wide range of previous LES results can be realized in a single LES code by varying only the numerical treatment of the equations of motion and the nature of subgrid-scale (SGS) closures. The simulations that maintain the greatest fidelity to DYCOMS-II observations are identified. The results show that using weighted essentially non-oscillatory (WENO) numerics for all resolved advective terms and no explicit SGS closure consistently produces the highest-fidelity simulations. This suggests that the numerical dissipation inherent in WENO schemes functions as a high-quality, implicit SGS closure for this stratocumulus case. Conversely, using oscillatory centered difference numerical schemes for momentum advection, WENO numerics for scalars, and explicitly modeled SGS fluxes consistently produces the lowest-fidelity simulations. We attribute this to the production of anomalously large SGS fluxes near the cloud tops through the interaction of numerical error in the momentum field with the scalar SGS model.

Numerics and subgrid‐scale modeling in large eddy simulations of stratocumulus clouds

PubMed Central

Mishra, Siddhartha; Schneider, Tapio; Kaul, Colleen M.; Tan, Zhihong

2017-01-01

Abstract Stratocumulus clouds are the most common type of boundary layer cloud; their radiative effects strongly modulate climate. Large eddy simulations (LES) of stratocumulus clouds often struggle to maintain fidelity to observations because of the sharp gradients occurring at the entrainment interfacial layer at the cloud top. The challenge posed to LES by stratocumulus clouds is evident in the wide range of solutions found in the LES intercomparison based on the DYCOMS‐II field campaign, where simulated liquid water paths for identical initial and boundary conditions varied by a factor of nearly 12. Here we revisit the DYCOMS‐II RF01 case and show that the wide range of previous LES results can be realized in a single LES code by varying only the numerical treatment of the equations of motion and the nature of subgrid‐scale (SGS) closures. The simulations that maintain the greatest fidelity to DYCOMS‐II observations are identified. The results show that using weighted essentially non‐oscillatory (WENO) numerics for all resolved advective terms and no explicit SGS closure consistently produces the highest‐fidelity simulations. This suggests that the numerical dissipation inherent in WENO schemes functions as a high‐quality, implicit SGS closure for this stratocumulus case. Conversely, using oscillatory centered difference numerical schemes for momentum advection, WENO numerics for scalars, and explicitly modeled SGS fluxes consistently produces the lowest‐fidelity simulations. We attribute this to the production of anomalously large SGS fluxes near the cloud tops through the interaction of numerical error in the momentum field with the scalar SGS model. PMID:28943997

Cloud Modeling

NASA Technical Reports Server (NTRS)

Tao, Wei-Kuo; Moncrieff, Mitchell; Einaud, Franco (Technical Monitor)

2001-01-01

Numerical cloud models have been developed and applied extensively to study cloud-scale and mesoscale processes during the past four decades. The distinctive aspect of these cloud models is their ability to treat explicitly (or resolve) cloud-scale dynamics. This requires the cloud models to be formulated from the non-hydrostatic equations of motion that explicitly include the vertical acceleration terms since the vertical and horizontal scales of convection are similar. Such models are also necessary in order to allow gravity waves, such as those triggered by clouds, to be resolved explicitly. In contrast, the hydrostatic approximation, usually applied in global or regional models, does allow the presence of gravity waves. In addition, the availability of exponentially increasing computer capabilities has resulted in time integrations increasing from hours to days, domain grids boxes (points) increasing from less than 2000 to more than 2,500,000 grid points with 500 to 1000 m resolution, and 3-D models becoming increasingly prevalent. The cloud resolving model is now at a stage where it can provide reasonably accurate statistical information of the sub-grid, cloud-resolving processes poorly parameterized in climate models and numerical prediction models.

Analysis of composite ablators using massively parallel computation

NASA Technical Reports Server (NTRS)

Shia, David

1995-01-01

In this work, the feasibility of using massively parallel computation to study the response of ablative materials is investigated. Explicit and implicit finite difference methods are used on a massively parallel computer, the Thinking Machines CM-5. The governing equations are a set of nonlinear partial differential equations. The governing equations are developed for three sample problems: (1) transpiration cooling, (2) ablative composite plate, and (3) restrained thermal growth testing. The transpiration cooling problem is solved using a solution scheme based solely on the explicit finite difference method. The results are compared with available analytical steady-state through-thickness temperature and pressure distributions and good agreement between the numerical and analytical solutions is found. It is also found that a solution scheme based on the explicit finite difference method has the following advantages: incorporates complex physics easily, results in a simple algorithm, and is easily parallelizable. However, a solution scheme of this kind needs very small time steps to maintain stability. A solution scheme based on the implicit finite difference method has the advantage that it does not require very small times steps to maintain stability. However, this kind of solution scheme has the disadvantages that complex physics cannot be easily incorporated into the algorithm and that the solution scheme is difficult to parallelize. A hybrid solution scheme is then developed to combine the strengths of the explicit and implicit finite difference methods and minimize their weaknesses. This is achieved by identifying the critical time scale associated with the governing equations and applying the appropriate finite difference method according to this critical time scale. The hybrid solution scheme is then applied to the ablative composite plate and restrained thermal growth problems. The gas storage term is included in the explicit pressure calculation of both problems. Results from ablative composite plate problems are compared with previous numerical results which did not include the gas storage term. It is found that the through-thickness temperature distribution is not affected much by the gas storage term. However, the through-thickness pressure and stress distributions, and the extent of chemical reactions are different from the previous numerical results. Two types of chemical reaction models are used in the restrained thermal growth testing problem: (1) pressure-independent Arrhenius type rate equations and (2) pressure-dependent Arrhenius type rate equations. The numerical results are compared to experimental results and the pressure-dependent model is able to capture the trend better than the pressure-independent one. Finally, a performance study is done on the hybrid algorithm using the ablative composite plate problem. It is found that there is a good speedup of performance on the CM-5. For 32 CPU's, the speedup of performance is 20. The efficiency of the algorithm is found to be a function of the size and execution time of a given problem and the effective parallelization of the algorithm. It also seems that there is an optimum number of CPU's to use for a given problem.

A multistage time-stepping scheme for the Navier-Stokes equations

NASA Technical Reports Server (NTRS)

Swanson, R. C.; Turkel, E.

1985-01-01

A class of explicit multistage time-stepping schemes is used to construct an algorithm for solving the compressible Navier-Stokes equations. Flexibility in treating arbitrary geometries is obtained with a finite-volume formulation. Numerical efficiency is achieved by employing techniques for accelerating convergence to steady state. Computer processing is enhanced through vectorization of the algorithm. The scheme is evaluated by solving laminar and turbulent flows over a flat plate and an NACA 0012 airfoil. Numerical results are compared with theoretical solutions or other numerical solutions and/or experimental data.

Implicit and Explicit Illusory Correlation as a Function of Political Ideology

PubMed Central

Carraro, Luciana; Negri, Paolo; Castelli, Luigi; Pastore, Massimiliano

2014-01-01

Research has demonstrated that people who embrace different ideological orientations often show differences at the level of basic cognitive processes. For instance, conservatives (vs. liberals) display an automatic selective attention for negative (vs. positive) stimuli, and tend to more easily form illusory correlations between negative information and minority groups. In the present work, we further explored this latter effect by examining whether it only involves the formation of explicit attitudes or it extends to implicit attitudes. To this end, following the typical illusory correlation paradigm, participants were presented with members of two numerically different groups (majority and minority) each performing either a positive or negative behaviour. Negative behaviors were relatively infrequent, and the proportion of positive and negative behaviors within each group was the same. Next, explicit and implicit (i.e., IAT-measured) attitudes were assessed. Results showed that conservatives (vs. liberals) displayed stronger explicit as well as implicit illusory correlations effects, forming more negative attitudes toward the minority (vs. majority) group at both the explicit and implicit level. PMID:24820311

Implicit and explicit illusory correlation as a function of political ideology.

PubMed

Carraro, Luciana; Negri, Paolo; Castelli, Luigi; Pastore, Massimiliano

2014-01-01

Research has demonstrated that people who embrace different ideological orientations often show differences at the level of basic cognitive processes. For instance, conservatives (vs. liberals) display an automatic selective attention for negative (vs. positive) stimuli, and tend to more easily form illusory correlations between negative information and minority groups. In the present work, we further explored this latter effect by examining whether it only involves the formation of explicit attitudes or it extends to implicit attitudes. To this end, following the typical illusory correlation paradigm, participants were presented with members of two numerically different groups (majority and minority) each performing either a positive or negative behaviour. Negative behaviors were relatively infrequent, and the proportion of positive and negative behaviors within each group was the same. Next, explicit and implicit (i.e., IAT-measured) attitudes were assessed. Results showed that conservatives (vs. liberals) displayed stronger explicit as well as implicit illusory correlations effects, forming more negative attitudes toward the minority (vs. majority) group at both the explicit and implicit level.

Analytical Approach to (2+1)-Dimensional Boussinesq Equation and (3+1)-Dimensional Kadomtsev-Petviashvili Equation

NASA Astrophysics Data System (ADS)

Sarıaydın, Selin; Yıldırım, Ahmet

2010-05-01

In this paper, we studied the solitary wave solutions of the (2+1)-dimensional Boussinesq equation utt -uxx-uyy-(u2)xx-uxxxx = 0 and the (3+1)-dimensional Kadomtsev-Petviashvili (KP) equation uxt -6ux 2 +6uuxx -uxxxx -uyy -uzz = 0. By using this method, an explicit numerical solution is calculated in the form of a convergent power series with easily computable components. To illustrate the application of this method numerical results are derived by using the calculated components of the homotopy perturbation series. The numerical solutions are compared with the known analytical solutions. Results derived from our method are shown graphically.

Explicitly Linking Pedagogy and Facilities to Campus Sustainability: Lessons from Carleton College and the University of Minnesota

ERIC Educational Resources Information Center

Savanick, Suzanne; Strong, Richard; Manning, Christie

2008-01-01

Campus sustainability projects provide an opportunity to explicitly link campus operations and academics. College and university buildings and grounds offer the potential for numerous hands-on sustainability projects. Few schools explicitly link sustainability projects with academics as often the academic side of an institution is separate from…

Soliton polarization rotation in fiber lasers

NASA Astrophysics Data System (ADS)

Afanasjev, V. V.

1995-02-01

I have found the approximate analytical solution in explicit form for a vector soliton with an arbitrary component ratio. My solution describes the dependence of soliton intensity on polarization angle and also nonlinear polarization rotation. The analytical results agree well with the numerical simulations.

Fluctuation theorem for the effusion of an ideal gas.

PubMed

Cleuren, B; Van den Broeck, C; Kawai, R

2006-08-01

The probability distribution of the entropy production for the effusion of an ideal gas between two compartments is calculated explicitly. The fluctuation theorem is verified. The analytic results are in good agreement with numerical data from hard disk molecular dynamics simulations.

An explicit predictor-corrector solver with applications to Burgers' equation

NASA Technical Reports Server (NTRS)

Dey, S. K.; Dey, C.

1983-01-01

Forward Euler's explicit, finite-difference formula of extrapolation, is used as a predictor and a convex formula as a corrector to integrate differential equations numerically. An application has been made to Burger's equation.

Sea breeze: Induced mesoscale systems and severe weather

NASA Technical Reports Server (NTRS)

Nicholls, M. E.; Pielke, R. A.; Cotton, W. R.

1990-01-01

Sea-breeze-deep convective interactions over the Florida peninsula were investigated using a cloud/mesoscale numerical model. The objective was to gain a better understanding of sea-breeze and deep convective interactions over the Florida peninsula using a high resolution convectively explicit model and to use these results to evaluate convective parameterization schemes. A 3-D numerical investigation of Florida convection was completed. The Kuo and Fritsch-Chappell parameterization schemes are summarized and evaluated.

Exploring Hamiltonian dielectric solvent molecular dynamics

NASA Astrophysics Data System (ADS)

Bauer, Sebastian; Tavan, Paul; Mathias, Gerald

2014-09-01

Hamiltonian dielectric solvent (HADES) is a recent method [7,25], which enables Hamiltonian molecular dynamics (MD) simulations of peptides and proteins in dielectric continua. Sample simulations of an α-helical decapeptide with and without explicit solvent demonstrate the high efficiency of HADES-MD. Addressing the folding of this peptide by replica exchange MD we study the properties of HADES by comparing melting curves, secondary structure motifs and salt bridges with explicit solvent results. Despite the unoptimized ad hoc parametrization of HADES, calculated reaction field energies correlate well with numerical grid solutions of the dielectric Poisson equation.

Implicit time accurate simulation of unsteady flow

NASA Astrophysics Data System (ADS)

van Buuren, René; Kuerten, Hans; Geurts, Bernard J.

2001-03-01

Implicit time integration was studied in the context of unsteady shock-boundary layer interaction flow. With an explicit second-order Runge-Kutta scheme, a reference solution to compare with the implicit second-order Crank-Nicolson scheme was determined. The time step in the explicit scheme is restricted by both temporal accuracy as well as stability requirements, whereas in the A-stable implicit scheme, the time step has to obey temporal resolution requirements and numerical convergence conditions. The non-linear discrete equations for each time step are solved iteratively by adding a pseudo-time derivative. The quasi-Newton approach is adopted and the linear systems that arise are approximately solved with a symmetric block Gauss-Seidel solver. As a guiding principle for properly setting numerical time integration parameters that yield an efficient time accurate capturing of the solution, the global error caused by the temporal integration is compared with the error resulting from the spatial discretization. Focus is on the sensitivity of properties of the solution in relation to the time step. Numerical simulations show that the time step needed for acceptable accuracy can be considerably larger than the explicit stability time step; typical ratios range from 20 to 80. At large time steps, convergence problems that are closely related to a highly complex structure of the basins of attraction of the iterative method may occur. Copyright

Numerical solution of transport equation for applications in environmental hydraulics and hydrology

NASA Astrophysics Data System (ADS)

Rashidul Islam, M.; Hanif Chaudhry, M.

1997-04-01

The advective term in the one-dimensional transport equation, when numerically discretized, produces artificial diffusion. To minimize such artificial diffusion, which vanishes only for Courant number equal to unity, transport owing to advection has been modeled separately. The numerical solution of the advection equation for a Gaussian initial distribution is well established; however, large oscillations are observed when applied to an initial distribution with sleep gradients, such as trapezoidal distribution of a constituent or propagation of mass from a continuous input. In this study, the application of seven finite-difference schemes and one polynomial interpolation scheme is investigated to solve the transport equation for both Gaussian and non-Gaussian (trapezoidal) initial distributions. The results obtained from the numerical schemes are compared with the exact solutions. A constant advective velocity is assumed throughout the transport process. For a Gaussian distribution initial condition, all eight schemes give excellent results, except the Lax scheme which is diffusive. In application to the trapezoidal initial distribution, explicit finite-difference schemes prove to be superior to implicit finite-difference schemes because the latter produce large numerical oscillations near the steep gradients. The Warming-Kutler-Lomax (WKL) explicit scheme is found to be better among this group. The Hermite polynomial interpolation scheme yields the best result for a trapezoidal distribution among all eight schemes investigated. The second-order accurate schemes are sufficiently accurate for most practical problems, but the solution of unusual problems (concentration with steep gradient) requires the application of higher-order (e.g. third- and fourth-order) accurate schemes.

Global Asymptotic Behavior of Iterative Implicit Schemes

NASA Technical Reports Server (NTRS)

Yee, H. C.; Sweby, P. K.

1994-01-01

The global asymptotic nonlinear behavior of some standard iterative procedures in solving nonlinear systems of algebraic equations arising from four implicit linear multistep methods (LMMs) in discretizing three models of 2 x 2 systems of first-order autonomous nonlinear ordinary differential equations (ODEs) is analyzed using the theory of dynamical systems. The iterative procedures include simple iteration and full and modified Newton iterations. The results are compared with standard Runge-Kutta explicit methods, a noniterative implicit procedure, and the Newton method of solving the steady part of the ODEs. Studies showed that aside from exhibiting spurious asymptotes, all of the four implicit LMMs can change the type and stability of the steady states of the differential equations (DEs). They also exhibit a drastic distortion but less shrinkage of the basin of attraction of the true solution than standard nonLMM explicit methods. The simple iteration procedure exhibits behavior which is similar to standard nonLMM explicit methods except that spurious steady-state numerical solutions cannot occur. The numerical basins of attraction of the noniterative implicit procedure mimic more closely the basins of attraction of the DEs and are more efficient than the three iterative implicit procedures for the four implicit LMMs. Contrary to popular belief, the initial data using the Newton method of solving the steady part of the DEs may not have to be close to the exact steady state for convergence. These results can be used as an explanation for possible causes and cures of slow convergence and nonconvergence of steady-state numerical solutions when using an implicit LMM time-dependent approach in computational fluid dynamics.

Numerical experiments with a symmetric high-resolution shock-capturing scheme

NASA Technical Reports Server (NTRS)

Yee, H. C.

1986-01-01

Characteristic-based explicit and implicit total variation diminishing (TVD) schemes for the two-dimensional compressible Euler equations have recently been developed. This is a generalization of recent work of Roe and Davis to a wider class of symmetric (non-upwind) TVD schemes other than Lax-Wendroff. The Roe and Davis schemes can be viewed as a subset of the class of explicit methods. The main properties of the present class of schemes are that they can be implicit, and, when steady-state calculations are sought, the numerical solution is independent of the time step. In a recent paper, a comparison of a linearized form of the present implicit symmetric TVD scheme with an implicit upwind TVD scheme originally developed by Harten and modified by Yee was given. Results favored the symmetric method. It was found that the latter is just as accurate as the upwind method while requiring less computational effort. Currently, more numerical experiments are being conducted on time-accurate calculations and on the effect of grid topology, numerical boundary condition procedures, and different flow conditions on the behavior of the method for steady-state applications. The purpose here is to report experiences with this type of scheme and give guidelines for its use.

Enforcing the Courant–Friedrichs–Lewy condition in explicitly conservative local time stepping schemes

DOE PAGES

Gnedin, Nickolay Y.; Semenov, Vadim A.; Kravtsov, Andrey V.

2018-01-30

In this study, an optimally efficient explicit numerical scheme for solving fluid dynamics equations, or any other parabolic or hyperbolic system of partial differential equations, should allow local regions to advance in time with their own, locally constrained time steps. However, such a scheme can result in violation of the Courant-Friedrichs-Lewy (CFL) condition, which is manifestly non-local. Although the violations can be considered to be "weak" in a certain sense and the corresponding numerical solution may be stable, such calculation does not guarantee the correct propagation speed for arbitrary waves. We use an experimental fluid dynamics code that allows cubicmore » "patches" of grid cells to step with independent, locally constrained time steps to demonstrate how the CFL condition can be enforced by imposing a condition on the time steps of neighboring patches. We perform several numerical tests that illustrate errors introduced in the numerical solutions by weak CFL condition violations and show how strict enforcement of the CFL condition eliminates these errors. In all our tests the strict enforcement of the CFL condition does not impose a significant performance penalty.« less

Study of effects of injector geometry on fuel-air mixing and combustion

NASA Technical Reports Server (NTRS)

Bangert, L. H.; Roach, R. L.

1977-01-01

An implicit finite-difference method has been developed for computing the flow in the near field of a fuel injector as part of a broader study of the effects of fuel injector geometry on fuel-air mixing and combustion. Detailed numerical results have been obtained for cases of laminar and turbulent flow without base injection, corresponding to the supersonic base flow problem. These numerical results indicated that the method is stable and convergent, and that significant savings in computer time can be achieved, compared with explicit methods.

Numerical solution of the unsteady diffusion-convection-reaction equation based on improved spectral Galerkin method

NASA Astrophysics Data System (ADS)

Zhong, Jiaqi; Zeng, Cheng; Yuan, Yupeng; Zhang, Yuzhe; Zhang, Ye

2018-04-01

The aim of this paper is to present an explicit numerical algorithm based on improved spectral Galerkin method for solving the unsteady diffusion-convection-reaction equation. The principal characteristics of this approach give the explicit eigenvalues and eigenvectors based on the time-space separation method and boundary condition analysis. With the help of Fourier series and Galerkin truncation, we can obtain the finite-dimensional ordinary differential equations which facilitate the system analysis and controller design. By comparing with the finite element method, the numerical solutions are demonstrated via two examples. It is shown that the proposed method is effective.

Normal forms of dispersive scalar Poisson brackets with two independent variables

NASA Astrophysics Data System (ADS)

Carlet, Guido; Casati, Matteo; Shadrin, Sergey

2018-03-01

We classify the dispersive Poisson brackets with one dependent variable and two independent variables, with leading order of hydrodynamic type, up to Miura transformations. We show that, in contrast to the case of a single independent variable for which a well-known triviality result exists, the Miura equivalence classes are parametrised by an infinite number of constants, which we call numerical invariants of the brackets. We obtain explicit formulas for the first few numerical invariants.

Implicit and Explicit Number-Space Associations Differentially Relate to Interference Control in Young Adults With ADHD

PubMed Central

Georges, Carrie; Hoffmann, Danielle; Schiltz, Christine

2018-01-01

Behavioral evidence for the link between numerical and spatial representations comes from the spatial-numerical association of response codes (SNARC) effect, consisting in faster reaction times to small/large numbers with the left/right hand respectively. The SNARC effect is, however, characterized by considerable intra- and inter-individual variability. It depends not only on the explicit or implicit nature of the numerical task, but also relates to interference control. To determine whether the prevalence of the latter relation in the elderly could be ascribed to younger individuals’ ceiling performances on executive control tasks, we determined whether the SNARC effect related to Stroop and/or Flanker effects in 26 young adults with ADHD. We observed a divergent pattern of correlation depending on the type of numerical task used to assess the SNARC effect and the type of interference control measure involved in number-space associations. Namely, stronger number-space associations during parity judgments involving implicit magnitude processing related to weaker interference control in the Stroop but not Flanker task. Conversely, stronger number-space associations during explicit magnitude classifications tended to be associated with better interference control in the Flanker but not Stroop paradigm. The association of stronger parity and magnitude SNARC effects with weaker and better interference control respectively indicates that different mechanisms underlie these relations. Activation of the magnitude-associated spatial code is irrelevant and potentially interferes with parity judgments, but in contrast assists explicit magnitude classifications. Altogether, the present study confirms the contribution of interference control to number-space associations also in young adults. It suggests that magnitude-associated spatial codes in implicit and explicit tasks are monitored by different interference control mechanisms, thereby explaining task-related intra-individual differences in number-space associations. PMID:29881363

Numerical investigations of shock wave interaction with laminar boundary layer on compressor profile

NASA Astrophysics Data System (ADS)

Piotrowicz, M.; Flaszyński, P.

2016-10-01

The investigation of shockwave boundary layer interaction on suction side of transonic compressor blade is one of main objectives of TFAST project (Transition Location Effect on Shock Wave Boundary Layer Interaction). In order to look more closely into the flow structure on suction side of a profile, a design of generic test section in linear transonic wind tunnel was proposed. The experimental and numerical results of flow structure on a suction side of the compressor profile investigations are presented. The numerical simulations are carried out for EARSM (Explicit Algebraic Reynolds Stress Model) turbulence model with transition model. The result are compared with oil flow visualisation, schlieren pictures, Pressure Sensitive Paint (PSP) and static pressure.

An explicit mixed numerical method for mesoscale model

NASA Technical Reports Server (NTRS)

Hsu, H.-M.

1981-01-01

A mixed numerical method has been developed for mesoscale models. The technique consists of a forward difference scheme for time tendency terms, an upstream scheme for advective terms, and a central scheme for the other terms in a physical system. It is shown that the mixed method is conditionally stable and highly accurate for approximating the system of either shallow-water equations in one dimension or primitive equations in three dimensions. Since the technique is explicit and two time level, it conserves computer and programming resources.

Numerical modelling of closed-cell aluminium foam under dynamic loading

NASA Astrophysics Data System (ADS)

Hazell, Paul; Kader, M. A.; Islam, M. A.; Escobedo, J. P.; Saadatfar, M.

2015-06-01

Closed-cell aluminium foams are extensively used in aerospace and automobile industries. The understanding of their behaviour under impact loading conditions is extremely important since impact problems are directly related to design of these engineering structures. This research investigates the response of a closed-cell aluminium foam (CYMAT) subjected to dynamic loading using the finite element software ABAQUS/explicit. The aim of this research is to numerically investigate the material and structural properties of closed-cell aluminium foam under impact loading conditions with interest in shock propagation and its effects on cell wall deformation. A μ-CT based 3D foam geometry is developed to simulate the local cell collapse behaviours. A number of numerical techniques are applied for modelling the crush behaviour of aluminium foam to obtain the more accurate results. The simulation results are compared with experimental data. Comparison of the results shows a good correlation between the experimental results and numerical predictions.

The development of an explicit thermochemical nonequilibrium algorithm and its application to compute three dimensional AFE flowfields

NASA Technical Reports Server (NTRS)

Palmer, Grant

1989-01-01

This study presents a three-dimensional explicit, finite-difference, shock-capturing numerical algorithm applied to viscous hypersonic flows in thermochemical nonequilibrium. The algorithm employs a two-temperature physical model. Equations governing the finite-rate chemical reactions are fully-coupled to the gas dynamic equations using a novel coupling technique. The new coupling method maintains stability in the explicit, finite-rate formulation while allowing relatively large global time steps. The code uses flux-vector accuracy. Comparisons with experimental data and other numerical computations verify the accuracy of the present method. The code is used to compute the three-dimensional flowfield over the Aeroassist Flight Experiment (AFE) vehicle at one of its trajectory points.

A Non-Dissipative Staggered Fourth-Order Accurate Explicit Finite Difference Scheme for the Time-Domain Maxwell's Equations

NASA Technical Reports Server (NTRS)

Yefet, Amir; Petropoulos, Peter G.

1999-01-01

We consider a divergence-free non-dissipative fourth-order explicit staggered finite difference scheme for the hyperbolic Maxwell's equations. Special one-sided difference operators are derived in order to implement the scheme near metal boundaries and dielectric interfaces. Numerical results show the scheme is long-time stable, and is fourth-order convergent over complex domains that include dielectric interfaces and perfectly conducting surfaces. We also examine the scheme's behavior near metal surfaces that are not aligned with the grid axes, and compare its accuracy to that obtained by the Yee scheme.

Combined AIE/EBE/GMRES approach to incompressible flows. [Adaptive Implicit-Explicit/Grouped Element-by-Element/Generalized Minimum Residuals

NASA Technical Reports Server (NTRS)

Liou, J.; Tezduyar, T. E.

1990-01-01

Adaptive implicit-explicit (AIE), grouped element-by-element (GEBE), and generalized minimum residuals (GMRES) solution techniques for incompressible flows are combined. In this approach, the GEBE and GMRES iteration methods are employed to solve the equation systems resulting from the implicitly treated elements, and therefore no direct solution effort is involved. The benchmarking results demonstrate that this approach can substantially reduce the CPU time and memory requirements in large-scale flow problems. Although the description of the concepts and the numerical demonstration are based on the incompressible flows, the approach presented here is applicable to larger class of problems in computational mechanics.

A Numerical Method for Predicting Rayleigh Surface Wave Velocity in Anisotropic Crystals (Postprint)

DTIC Science & Technology

2017-09-05

generalized version of the equations are very difficult to derive, even in symbolic math languages such as Mathematica. As a result, the equations are...formalism, Math . Mech. Solids 9 (1) (2004) 5–15. [8] M. Destrade, The explicit secular equation for surface acoustic waves in monoclinic elastic crystals...Q. J. Mech. Appl. Math . 55 (2) (2002) 297–311. [10] D. Taylor, Surface waves in anisotropic media: the secular equation and its numerical solution

Stability analysis of Eulerian-Lagrangian methods for the one-dimensional shallow-water equations

USGS Publications Warehouse

Casulli, V.; Cheng, R.T.

1990-01-01

In this paper stability and error analyses are discussed for some finite difference methods when applied to the one-dimensional shallow-water equations. Two finite difference formulations, which are based on a combined Eulerian-Lagrangian approach, are discussed. In the first part of this paper the results of numerical analyses for an explicit Eulerian-Lagrangian method (ELM) have shown that the method is unconditionally stable. This method, which is a generalized fixed grid method of characteristics, covers the Courant-Isaacson-Rees method as a special case. Some artificial viscosity is introduced by this scheme. However, because the method is unconditionally stable, the artificial viscosity can be brought under control either by reducing the spatial increment or by increasing the size of time step. The second part of the paper discusses a class of semi-implicit finite difference methods for the one-dimensional shallow-water equations. This method, when the Eulerian-Lagrangian approach is used for the convective terms, is also unconditionally stable and highly accurate for small space increments or large time steps. The semi-implicit methods seem to be more computationally efficient than the explicit ELM; at each time step a single tridiagonal system of linear equations is solved. The combined explicit and implicit ELM is best used in formulating a solution strategy for solving a network of interconnected channels. The explicit ELM is used at channel junctions for each time step. The semi-implicit method is then applied to the interior points in each channel segment. Following this solution strategy, the channel network problem can be reduced to a set of independent one-dimensional open-channel flow problems. Numerical results support properties given by the stability and error analyses. ?? 1990.

The northern tidal dynamic of Aceh waters: A 3D numerical model

NASA Astrophysics Data System (ADS)

Irham, M.; Miswar, E.; Ilhamsyah, Y.; Setiawan, I.

2018-05-01

The northern tidal dynamic of Aceh waters studied by employing three-dimensional (3D) numerical hydrodynamic model. The purpose of this study is to understand the phenomena and the characteristic of the northern tidal dynamic of Aceh waters. The research used the explicit-splitting scheme numerical model of Navier-Stokes formulation. The result displays that the vertical rotation of flow movement (vertical eddy) at a depth of 15 to 25 meter eastern part of the study area. Hence, the result also informs that the current circulation identically to the upwelling in the western region of Aceh during the wet season and vice versa. However, during the transitional season, the flow circulation depends on how the tidal dynamic occurs in the area.

Exact sum rules for inhomogeneous systems containing a zero mode

DOE Office of Scientific and Technical Information (OSTI.GOV)

Amore, Paolo, E-mail: paolo.amore@gmail.com

2014-10-15

We show that the formulas for the sum rules for the eigenvalues of inhomogeneous systems that we have obtained in two recent papers are incomplete when the system contains a zero mode. We prove that there are finite contributions of the zero mode to the sum rules and we explicitly calculate the expressions for the sum rules of order one and two. The previous results for systems that do not contain a zero mode are unaffected. - Highlights: • We discuss the sum rules of the eigenvalues of inhomogeneous systems containing a zero mode. • We derive the explicit expressionsmore » for sum rules of order one and two. • We perform accurate numerical tests of these results for three examples.« less

Explicit densities of multidimensional ballistic Lévy walks.

PubMed

Magdziarz, Marcin; Zorawik, Tomasz

2016-08-01

Lévy walks have proved to be useful models of stochastic dynamics with a number of applications in the modeling of real-life phenomena. In this paper we derive explicit formulas for densities of the two- (2D) and three-dimensional (3D) ballistic Lévy walks, which are most important in applications. It turns out that in the 3D case the densities are given by elementary functions. The densities of the 2D Lévy walks are expressed in terms of hypergeometric functions and the right-side Riemann-Liouville fractional derivative, which allows us to efficiently evaluate them numerically. The theoretical results agree perfectly with Monte Carlo simulations.

Assessment of an Explicit Algebraic Reynolds Stress Model

NASA Technical Reports Server (NTRS)

Carlson, Jan-Renee

2005-01-01

This study assesses an explicit algebraic Reynolds stress turbulence model in the in the three-dimensional Reynolds averaged Navier-Stokes (RANS) solver, ISAAC (Integrated Solution Algorithm for Arbitrary Con gurations). Additionally, it compares solutions for two select configurations between ISAAC and the RANS solver PAB3D. This study compares with either direct numerical simulation data, experimental data, or empirical models for several different geometries with compressible, separated, and high Reynolds number flows. In general, the turbulence model matched data or followed experimental trends well, and for the selected configurations, the computational results of ISAAC closely matched those of PAB3D using the same turbulence model.

Explicit formulation of second and third order optical nonlinearity in the FDTD framework

NASA Astrophysics Data System (ADS)

Varin, Charles; Emms, Rhys; Bart, Graeme; Fennel, Thomas; Brabec, Thomas

2018-01-01

The finite-difference time-domain (FDTD) method is a flexible and powerful technique for rigorously solving Maxwell's equations. However, three-dimensional optical nonlinearity in current commercial and research FDTD softwares requires solving iteratively an implicit form of Maxwell's equations over the entire numerical space and at each time step. Reaching numerical convergence demands significant computational resources and practical implementation often requires major modifications to the core FDTD engine. In this paper, we present an explicit method to include second and third order optical nonlinearity in the FDTD framework based on a nonlinear generalization of the Lorentz dispersion model. A formal derivation of the nonlinear Lorentz dispersion equation is equally provided, starting from the quantum mechanical equations describing nonlinear optics in the two-level approximation. With the proposed approach, numerical integration of optical nonlinearity and dispersion in FDTD is intuitive, transparent, and fully explicit. A strong-field formulation is also proposed, which opens an interesting avenue for FDTD-based modelling of the extreme nonlinear optics phenomena involved in laser filamentation and femtosecond micromachining of dielectrics.

Numerical simulation of fluid flow around a scramaccelerator projectile

NASA Technical Reports Server (NTRS)

Pepper, Darrell W.; Humphrey, Joseph W.; Sobota, Thomas H.

1991-01-01

Numerical simulations of the fluid motion and temperature distribution around a 'scramaccelerator' projectile are obtained for Mach numbers in the 5-10 range. A finite element method is used to solve the equations of motion for inviscid and viscous two-dimensional or axisymmetric compressible flow. The time-dependent equations are solved explicitly, using bilinear isoparametric quadrilateral elements, mass lumping, and a shock-capturing Petrov-Galerkin formulation. Computed results indicate that maintaining on-design performance for controlling and stabilizing oblique detonation waves is critically dependent on projectile shape and Mach number.

Sound propagation in a duct of periodic wall structure. [numerical analysis

NASA Technical Reports Server (NTRS)

Kurze, U.

1978-01-01

A boundary condition, which accounts for the coupling in the sections behind the duct boundary, is given for the sound-absorbing duct with a periodic structure of the wall lining and using regular partition walls. The soundfield in the duct is suitably described by the method of differences. For locally active walls this renders an explicit approximate solution for the propagation constant. Coupling may be accounted for by the method of differences in a clear manner. Numerical results agree with measurements and yield information which has technical applications.

A new heterogeneous asynchronous explicit-implicit time integrator for nonsmooth dynamics

NASA Astrophysics Data System (ADS)

Fekak, Fatima-Ezzahra; Brun, Michael; Gravouil, Anthony; Depale, Bruno

2017-07-01

In computational structural dynamics, particularly in the presence of nonsmooth behavior, the choice of the time-step and the time integrator has a critical impact on the feasibility of the simulation. Furthermore, in some cases, as in the case of a bridge crane under seismic loading, multiple time-scales coexist in the same problem. In that case, the use of multi-time scale methods is suitable. Here, we propose a new explicit-implicit heterogeneous asynchronous time integrator (HATI) for nonsmooth transient dynamics with frictionless unilateral contacts and impacts. Furthermore, we present a new explicit time integrator for contact/impact problems where the contact constraints are enforced using a Lagrange multiplier method. In other words, the aim of this paper consists in using an explicit time integrator with a fine time scale in the contact area for reproducing high frequency phenomena, while an implicit time integrator is adopted in the other parts in order to reproduce much low frequency phenomena and to optimize the CPU time. In a first step, the explicit time integrator is tested on a one-dimensional example and compared to Moreau-Jean's event-capturing schemes. The explicit algorithm is found to be very accurate and the scheme has generally a higher order of convergence than Moreau-Jean's schemes and provides also an excellent energy behavior. Then, the two time scales explicit-implicit HATI is applied to the numerical example of a bridge crane under seismic loading. The results are validated in comparison to a fine scale full explicit computation. The energy dissipated in the implicit-explicit interface is well controlled and the computational time is lower than a full-explicit simulation.

Automation based on knowledge modeling theory and its applications in engine diagnostic systems using Space Shuttle Main Engine vibrational data. M.S. Thesis

NASA Technical Reports Server (NTRS)

Kim, Jonnathan H.

1995-01-01

Humans can perform many complicated tasks without explicit rules. This inherent and advantageous capability becomes a hurdle when a task is to be automated. Modern computers and numerical calculations require explicit rules and discrete numerical values. In order to bridge the gap between human knowledge and automating tools, a knowledge model is proposed. Knowledge modeling techniques are discussed and utilized to automate a labor and time intensive task of detecting anomalous bearing wear patterns in the Space Shuttle Main Engine (SSME) High Pressure Oxygen Turbopump (HPOTP).

A numerical investigation of viscous, incompressible flow past an axisymmetric body with and without spin

NASA Astrophysics Data System (ADS)

Weber, K. F.

1985-12-01

This study deals with a preliminary investigation of the effects of spin on the axisymmetric flow past a body of revolution. The study has its genesis larger problem of Magnus forces on spinning bodies at angle of attack. However, the fundamental behavior that arises when a spinning body is placed in a uniform stream is still not well understood; therefore, the problem of axisymmetric flow with spin was undertaken. The body consists of a 3-caliber cant-ogive blunted by a spherical nosecap, a 2-caliber cylindrical section, and a 1-caliber boattail. Numerical solutions of the compressible laminar Navier-Stokes equations are obtained using a modified version of the implicit-explicit method developed by MacCormack in 1981. The benchmark problem is the nonspinning body in uniform flow at a Reynolds number of 1.14. The results show that the modified method performs well and allows time steps that are in order of magnitude greater than those permitted by explicit stability criteria.

A numerical investigation of the crystallographic texture effect on the surface roughening in aluminum polycrystals

NASA Astrophysics Data System (ADS)

Romanova, V.; Balokhonov, R.; Batukhtina, E.; Zinovieva, O.; Bezmozgiy, I.

2015-10-01

The results of a numerical analysis of the mesoscale surface roughening in a polycrystalline aluminum alloy exposed to uniaxial tension are presented. A 3D finite-element model taking an explicit account of grain structure is developed. The model describes a constitutive behavior of the material on the grain scale, using anisotropic elasticity and crystal plasticity theory. The effects of the grain shape and texture on the deformation-induced roughening are investigated. Calculation results have shown that surface roughness is much higher and develops at the highest rate in a polycrystal with equiaxed grains where both the micro- and mesoscale surface displacements are observed.

Fully implicit Particle-in-cell algorithms for multiscale plasma simulation

DOE Office of Scientific and Technical Information (OSTI.GOV)

Chacon, Luis

The outline of the paper is as follows: Particle-in-cell (PIC) methods for fully ionized collisionless plasmas, explicit vs. implicit PIC, 1D ES implicit PIC (charge and energy conservation, moment-based acceleration), and generalization to Multi-D EM PIC: Vlasov-Darwin model (review and motivation for Darwin model, conservation properties (energy, charge, and canonical momenta), and numerical benchmarks). The author demonstrates a fully implicit, fully nonlinear, multidimensional PIC formulation that features exact local charge conservation (via a novel particle mover strategy), exact global energy conservation (no particle self-heating or self-cooling), adaptive particle orbit integrator to control errors in momentum conservation, and canonical momenta (EM-PICmore » only, reduced dimensionality). The approach is free of numerical instabilities: ω peΔt >> 1, and Δx >> λ D. It requires many fewer dofs (vs. explicit PIC) for comparable accuracy in challenging problems. Significant CPU gains (vs explicit PIC) have been demonstrated. The method has much potential for efficiency gains vs. explicit in long-time-scale applications. Moment-based acceleration is effective in minimizing N FE, leading to an optimal algorithm.« less

Efficient algorithms and implementations of entropy-based moment closures for rarefied gases

NASA Astrophysics Data System (ADS)

Schaerer, Roman Pascal; Bansal, Pratyuksh; Torrilhon, Manuel

2017-07-01

We present efficient algorithms and implementations of the 35-moment system equipped with the maximum-entropy closure in the context of rarefied gases. While closures based on the principle of entropy maximization have been shown to yield very promising results for moderately rarefied gas flows, the computational cost of these closures is in general much higher than for closure theories with explicit closed-form expressions of the closing fluxes, such as Grad's classical closure. Following a similar approach as Garrett et al. (2015) [13], we investigate efficient implementations of the computationally expensive numerical quadrature method used for the moment evaluations of the maximum-entropy distribution by exploiting its inherent fine-grained parallelism with the parallelism offered by multi-core processors and graphics cards. We show that using a single graphics card as an accelerator allows speed-ups of two orders of magnitude when compared to a serial CPU implementation. To accelerate the time-to-solution for steady-state problems, we propose a new semi-implicit time discretization scheme. The resulting nonlinear system of equations is solved with a Newton type method in the Lagrange multipliers of the dual optimization problem in order to reduce the computational cost. Additionally, fully explicit time-stepping schemes of first and second order accuracy are presented. We investigate the accuracy and efficiency of the numerical schemes for several numerical test cases, including a steady-state shock-structure problem.

Is it really naked? On cosmic censorship in string theory

DOE Office of Scientific and Technical Information (OSTI.GOV)

Frolov, Andrei V.

We investigate the possibility of cosmic censorship violation in string theory using a characteristic double-null code, which penetrates horizons and is capable of resolving the spacetime all the way to the singularity. We perform high-resolution numerical simulations of the evolution of negative mass initial scalar field profiles, which were argued to provide a counterexample to cosmic censorship conjecture for AdS-asymptotic spacetimes in five-dimensional supergravity. In no instances formation of naked singularity is seen. Instead, numerical evidence indicates that black holes form in the collapse. Our results are consistent with earlier numerical studies, and explicitly show where the 'no black hole'more » argument breaks.« less

Numerical modelling of physical processes in a ballistic laboratory setup with a tapered adapter and plastic piston used for obtaining high muzzle velocities

NASA Astrophysics Data System (ADS)

Bykov, N. V.

2014-12-01

Numerical modelling of a ballistic setup with a tapered adapter and plastic piston is considered. The processes in the firing chamber are described within the framework of quasi- one-dimensional gas dynamics and a geometrical law of propellant burn by means of Lagrangian mass coordinates. The deformable piston is considered to be an ideal liquid with specific equations of state. The numerical solution is obtained by means of a modified explicit von Neumann scheme. The calculation results given show that the ballistic setup with a tapered adapter and plastic piston produces increased shell muzzle velocities by a factor of more than 1.5-2.

Multirate Particle-in-Cell Time Integration Techniques of Vlasov-Maxwell Equations for Collisionless Kinetic Plasma Simulations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Chen, Guangye; Chacon, Luis; Knoll, Dana Alan

2015-07-31

A multi-rate PIC formulation was developed that employs large timesteps for slow field evolution, and small (adaptive) timesteps for particle orbit integrations. Implementation is based on a JFNK solver with nonlinear elimination and moment preconditioning. The approach is free of numerical instabilities (ω peΔt >>1, and Δx >> λ D), and requires many fewer dofs (vs. explicit PIC) for comparable accuracy in challenging problems. Significant gains (vs. conventional explicit PIC) may be possible for large scale simulations. The paper is organized as follows: Vlasov-Maxwell Particle-in-cell (PIC) methods for plasmas; Explicit, semi-implicit, and implicit time integrations; Implicit PIC formulation (Jacobian-Free Newton-Krylovmore » (JFNK) with nonlinear elimination allows different treatments of disparate scales, discrete conservation properties (energy, charge, canonical momentum, etc.)); Some numerical examples; and Summary.« less

Explicit symplectic algorithms based on generating functions for relativistic charged particle dynamics in time-dependent electromagnetic field

NASA Astrophysics Data System (ADS)

Zhang, Ruili; Wang, Yulei; He, Yang; Xiao, Jianyuan; Liu, Jian; Qin, Hong; Tang, Yifa

2018-02-01

Relativistic dynamics of a charged particle in time-dependent electromagnetic fields has theoretical significance and a wide range of applications. The numerical simulation of relativistic dynamics is often multi-scale and requires accurate long-term numerical simulations. Therefore, explicit symplectic algorithms are much more preferable than non-symplectic methods and implicit symplectic algorithms. In this paper, we employ the proper time and express the Hamiltonian as the sum of exactly solvable terms and product-separable terms in space-time coordinates. Then, we give the explicit symplectic algorithms based on the generating functions of orders 2 and 3 for relativistic dynamics of a charged particle. The methodology is not new, which has been applied to non-relativistic dynamics of charged particles, but the algorithm for relativistic dynamics has much significance in practical simulations, such as the secular simulation of runaway electrons in tokamaks.

The assessment of mangrove biomass and carbon in West Africa: a spatially explicit analytical framework

Treesearch

Wenwu Tang; Wenpeng Feng; Meijuan Jia; Jiyang Shi; Huifang Zuo; Carl C. Trettin

2015-01-01

Mangrove forests are highly productive and have large carbon sinks while also providing numerous goods and ecosystem services. However, effective management and conservation of the mangrove forests are often dependent on spatially explicit assessments of the resource. Given the remote and highly dispersed nature of mangroves, estimation of biomass and carbon...

Should Cost-Effectiveness Analysis Include the Cost of Consumption Activities? AN Empirical Investigation.

PubMed

Adarkwah, Charles Christian; Sadoghi, Amirhossein; Gandjour, Afschin

2016-02-01

There has been a debate on whether cost-effectiveness analysis should consider the cost of consumption and leisure time activities when using the quality-adjusted life year as a measure of health outcome under a societal perspective. The purpose of this study was to investigate whether the effects of ill health on consumptive activities are spontaneously considered in a health state valuation exercise and how much this matters. The survey enrolled patients with inflammatory bowel disease in Germany (n = 104). Patients were randomized to explicit and no explicit instruction for the consideration of consumption and leisure effects in a time trade-off (TTO) exercise. Explicit instruction to consider non-health-related utility in TTO exercises did not influence TTO scores. However, spontaneous consideration of non-health-related utility in patients without explicit instruction (60% of respondents) led to significantly lower TTO scores. Results suggest an inclusion of consumption costs in the numerator of the cost-effectiveness ratio, at least for those respondents who spontaneously consider non-health-related utility from treatment. Results also suggest that exercises eliciting health valuations from the general public may include a description of the impact of disease on consumptive activities. Copyright © 2015 John Wiley & Sons, Ltd.

Second-order numerical methods for multi-term fractional differential equations: Smooth and non-smooth solutions

NASA Astrophysics Data System (ADS)

Zeng, Fanhai; Zhang, Zhongqiang; Karniadakis, George Em

2017-12-01

Starting with the asymptotic expansion of the error equation of the shifted Gr\\"{u}nwald--Letnikov formula, we derive a new modified weighted shifted Gr\\"{u}nwald--Letnikov (WSGL) formula by introducing appropriate correction terms. We then apply one special case of the modified WSGL formula to solve multi-term fractional ordinary and partial differential equations, and we prove the linear stability and second-order convergence for both smooth and non-smooth solutions. We show theoretically and numerically that numerical solutions up to certain accuracy can be obtained with only a few correction terms. Moreover, the correction terms can be tuned according to the fractional derivative orders without explicitly knowing the analytical solutions. Numerical simulations verify the theoretical results and demonstrate that the new formula leads to better performance compared to other known numerical approximations with similar resolution.

Changes Observed in Views of Nature of Science during a Historically Based Unit

ERIC Educational Resources Information Center

Rudge, David Wÿss; Cassidy, David Paul; Fulford, Janice Marie; Howe, Eric Michael

2014-01-01

Numerous empirical studies have provided evidence of the effectiveness of an explicit and reflective approach to the learning of issues associated with the nature of science (NOS) (c.f. Abd-El-Khalick and Lederman in "J Res Sci Teach" 37(10):1057-1095, 2000). This essay reports the results of a mixed-methods association study involving…

Hoph Bifurcation in Viscous, Low Speed Flows About an Airfoil with Structural Coupling

DTIC Science & Technology

1993-03-01

8 2.1 Equations of Motion ...... ..................... 8 2.2 Coordinate Transformation ....................... 13 2.3 Aerodynamic...a-frame) f - Apparent body forces applied in noninertial system fL - Explicit fourth-order numerical damping term Ai - Implicit fourth-order...resulting airfoil motion . The equations describing the airfoil motion are integrated in time using a fourth-order Runge-Kutta algorithm. The

Binding in pair potentials of liquid simple metals from nonlocality in electronic kinetic energy

NASA Technical Reports Server (NTRS)

Perrot, F.; March, N. H.

1990-01-01

The paper presents an explicit expression for the pair potential in liquid simple metals from low-order density-gradient theory when the superposition of single-center displaced charges is employed. Numerical results are presented for the gradient expansion pair interaction in liquid Na and Be. The low-order density-gradient equation for the pair potential is presented.

Numerical study of the flow in a three-dimensional thermally driven cavity

NASA Astrophysics Data System (ADS)

Rauwoens, Pieter; Vierendeels, Jan; Merci, Bart

2008-06-01

Solutions for the fully compressible Navier-Stokes equations are presented for the flow and temperature fields in a cubic cavity with large horizontal temperature differences. The ideal-gas approximation for air is assumed and viscosity is computed using Sutherland's law. The three-dimensional case forms an extension of previous studies performed on a two-dimensional square cavity. The influence of imposed boundary conditions in the third dimension is investigated as a numerical experiment. Comparison is made between convergence rates in case of periodic and free-slip boundary conditions. Results with no-slip boundary conditions are presented as well. The effect of the Rayleigh number is studied. Results are computed using a finite volume method on a structured, collocated grid. An explicit third-order discretization for the convective part and an implicit central discretization for the acoustic part and for the diffusive part are used. To stabilize the scheme an artificial dissipation term for the pressure and the temperature is introduced. The discrete equations are solved using a time-marching method with restrictions on the timestep corresponding to the explicit parts of the solver. Multigrid is used as acceleration technique.

An efficient nonlinear finite-difference approach in the computational modeling of the dynamics of a nonlinear diffusion-reaction equation in microbial ecology.

PubMed

Macías-Díaz, J E; Macías, Siegfried; Medina-Ramírez, I E

2013-12-01

In this manuscript, we present a computational model to approximate the solutions of a partial differential equation which describes the growth dynamics of microbial films. The numerical technique reported in this work is an explicit, nonlinear finite-difference methodology which is computationally implemented using Newton's method. Our scheme is compared numerically against an implicit, linear finite-difference discretization of the same partial differential equation, whose computer coding requires an implementation of the stabilized bi-conjugate gradient method. Our numerical results evince that the nonlinear approach results in a more efficient approximation to the solutions of the biofilm model considered, and demands less computer memory. Moreover, the positivity of initial profiles is preserved in the practice by the nonlinear scheme proposed. Copyright © 2013 Elsevier Ltd. All rights reserved.

Time-Accurate Numerical Simulations of Synthetic Jet Quiescent Air

NASA Technical Reports Server (NTRS)

Rupesh, K-A. B.; Ravi, B. R.; Mittal, R.; Raju, R.; Gallas, Q.; Cattafesta, L.

2007-01-01

The unsteady evolution of three-dimensional synthetic jet into quiescent air is studied by time-accurate numerical simulations using a second-order accurate mixed explicit-implicit fractional step scheme on Cartesian grids. Both two-dimensional and three-dimensional calculations of synthetic jet are carried out at a Reynolds number (based on average velocity during the discharge phase of the cycle V(sub j), and jet width d) of 750 and Stokes number of 17.02. The results obtained are assessed against PIV and hotwire measurements provided for the NASA LaRC workshop on CFD validation of synthetic jets.

Explicit symplectic algorithms based on generating functions for charged particle dynamics.

PubMed

Zhang, Ruili; Qin, Hong; Tang, Yifa; Liu, Jian; He, Yang; Xiao, Jianyuan

2016-07-01

Dynamics of a charged particle in the canonical coordinates is a Hamiltonian system, and the well-known symplectic algorithm has been regarded as the de facto method for numerical integration of Hamiltonian systems due to its long-term accuracy and fidelity. For long-term simulations with high efficiency, explicit symplectic algorithms are desirable. However, it is generally believed that explicit symplectic algorithms are only available for sum-separable Hamiltonians, and this restriction limits the application of explicit symplectic algorithms to charged particle dynamics. To overcome this difficulty, we combine the familiar sum-split method and a generating function method to construct second- and third-order explicit symplectic algorithms for dynamics of charged particle. The generating function method is designed to generate explicit symplectic algorithms for product-separable Hamiltonian with form of H(x,p)=p_{i}f(x) or H(x,p)=x_{i}g(p). Applied to the simulations of charged particle dynamics, the explicit symplectic algorithms based on generating functions demonstrate superiorities in conservation and efficiency.

Explicit symplectic algorithms based on generating functions for charged particle dynamics

NASA Astrophysics Data System (ADS)

Zhang, Ruili; Qin, Hong; Tang, Yifa; Liu, Jian; He, Yang; Xiao, Jianyuan

2016-07-01

Dynamics of a charged particle in the canonical coordinates is a Hamiltonian system, and the well-known symplectic algorithm has been regarded as the de facto method for numerical integration of Hamiltonian systems due to its long-term accuracy and fidelity. For long-term simulations with high efficiency, explicit symplectic algorithms are desirable. However, it is generally believed that explicit symplectic algorithms are only available for sum-separable Hamiltonians, and this restriction limits the application of explicit symplectic algorithms to charged particle dynamics. To overcome this difficulty, we combine the familiar sum-split method and a generating function method to construct second- and third-order explicit symplectic algorithms for dynamics of charged particle. The generating function method is designed to generate explicit symplectic algorithms for product-separable Hamiltonian with form of H (x ,p ) =pif (x ) or H (x ,p ) =xig (p ) . Applied to the simulations of charged particle dynamics, the explicit symplectic algorithms based on generating functions demonstrate superiorities in conservation and efficiency.

Ultrasonic modeling of an embedded elliptic crack

NASA Astrophysics Data System (ADS)

Fradkin, Larissa Ju.; Zalipaev, Victor

2000-05-01

Experiments indicate that the radiating near zone of a compressional circular transducer directly coupled to a homogeneous and isotropic solid has the following structure: there are geometrical zones where one can distinguish a plane compressional wave and toroidal waves, both compressional and shear, radiated by the transducer rim. As has been shown previously the modern diffraction theory allows to describe these explicitly. It also gives explicit asymptotic description of waves present in the transition zones. In case of a normal incidence of a plane compressional wave the explicit expressions have been obtained by Achenbach and co-authors for the fields diffracted by a penny-shaped crack. We build on the above work by applying the uniform GTD to model an oblique incidence of a plane compressional wave on an elliptical crack. We compare our asymptotic results with numerical results based on the boundary integral code as developed by Glushkovs, Krasnodar University, Russia. The asymptotic formulas form a basis of a code for high-frequency simulation of ultrasonic scattering by elliptical cracks situated in the vicinity of a compressional circular transducer, currently under development at our Center.

Modeling the ion transfer and polarization of ion exchange membranes in bioelectrochemical systems.

PubMed

Harnisch, Falk; Warmbier, Robert; Schneider, Ralf; Schröder, Uwe

2009-06-01

An explicit numerical model for the charge balancing ion transfer across monopolar ion exchange membranes under conditions of bioelectrochemical systems is presented. Diffusion and migration equations have been solved according to the Nernst-Planck Equation and the resulting ion concentrations, pH values and the resistance values of the membrane for different conditions were computed. The modeling results underline the principle limitations of the application of ion exchange membranes in biological fuel cells and electrolyzers, caused by the inherent occurrence of a pH-gradient between anode and cathode compartment, and an increased ohmic membrane resistance at decreasing electrolyte concentrations. Finally, the physical and numerical limitations of the model are discussed.

Computation of Turbulent Wake Flows in Variable Pressure Gradient

NASA Technical Reports Server (NTRS)

Duquesne, N.; Carlson, J. R.; Rumsey, C. L.; Gatski, T. B.

1999-01-01

Transport aircraft performance is strongly influenced by the effectiveness of high-lift systems. Developing wakes generated by the airfoil elements are subjected to strong pressure gradients and can thicken very rapidly, limiting maximum lift. This paper focuses on the effects of various pressure gradients on developing symmetric wakes and on the ability of a linear eddy viscosity model and a non-linear explicit algebraic stress model to accurately predict their downstream evolution. In order to reduce the uncertainties arising from numerical issues when assessing the performance of turbulence models, three different numerical codes with the same turbulence models are used. Results are compared to available experimental data to assess the accuracy of the computational results.

Analytical representation of dynamical quantities in G W from a matrix resolvent

NASA Astrophysics Data System (ADS)

Gesenhues, J.; Nabok, D.; Rohlfing, M.; Draxl, C.

2017-12-01

The power of the G W formalism is, to a large extent, based on the explicit treatment of dynamical correlations in the self-energy. This dynamics is taken into account by calculating the energy dependence of the screened Coulomb interaction W , followed by a convolution with the Green's function G . In order to obtain the energy dependence of W the prevalent methods are plasmon-pole models and numerical integration techniques. In this paper, we discuss an alternative approach, in which the energy-dependent screening is calculated by determining the resolvent, which is set up from a matrix representation of the dielectric function. On the one hand, this refrains from a numerical energy convolution and allows one to actually write down the energy dependence of W explicitly (like in the plasmon-pole models). On the other hand, the method is at least as accurate as the numerical approaches due to its multipole nature. We discuss the theoretical setup in some detail, give insight into the computational aspects, and present results for Si, C, GaAs, and LiF. Finally, we argue that the analytic representability is not only useful for educational purposes but may also be of avail for the development of theory that goes beyond G W .

Applying the Explicit Time Central Difference Method for Numerical Simulation of the Dynamic Behavior of Elastoplastic Flexible Reinforced Plates

NASA Astrophysics Data System (ADS)

Yankovskii, A. P.

2017-12-01

Based on a stepwise algorithm involving central finite differences for the approximation in time, a mathematical model is developed for elastoplastic deformation of cross-reinforced plates with isotropically hardening materials of components of the composition. The model allows obtaining the solution of elastoplastic problems at discrete points in time by an explicit scheme. The initial boundary value problem of the dynamic behavior of flexible plates reinforced in their own plane is formulated in the von Kármán approximation with allowance for their weakened resistance to the transverse shear. With a common approach, the resolving equations corresponding to two variants of the Timoshenko theory are obtained. An explicit "cross" scheme for numerical integration of the posed initial boundary value problem has been constructed. The scheme is consistent with the incremental algorithm used for simulating the elastoplastic behavior of a reinforced medium. Calculations of the dynamic behavior have been performed for elastoplastic cylindrical bending of differently reinforced fiberglass rectangular elongated plates. It is shown that the reinforcement structure significantly affects their elastoplastic dynamic behavior. It has been found that the classical theory of plates is as a rule unacceptable for carrying out the required calculations (except for very thin plates), and the first version of the Timoshenko theory yields reasonable results only in cases of relatively thin constructions reinforced by lowmodulus fibers. Proceeding from the results of the work, it is recommended to use the second variant of the Timoshenko theory (as a more accurate one) for calculations of the elastoplastic behavior of reinforced plates.

Preconditioning for the Navier-Stokes equations with finite-rate chemistry

NASA Technical Reports Server (NTRS)

Godfrey, Andrew G.; Walters, Robert W.; Van Leer, Bram

1993-01-01

The preconditioning procedure for generalized finite-rate chemistry and the proper preconditioning for the one-dimensional Navier-Stokes equations are presented. Eigenvalue stiffness is resolved and convergence-rate acceleration is demonstrated over the entire Mach-number range from the incompressible to the hypersonic. Specific benefits are realized at low and transonic flow speeds. The extended preconditioning matrix accounts for thermal and chemical non-equilibrium and its implementation is explained for both explicit and implicit time marching. The effect of higher-order spatial accuracy and various flux splittings is investigated. Numerical analysis reveals the possible theoretical improvements from using proconditioning at all Mach numbers. Numerical results confirm the expectations from the numerical analysis. Representative test cases include flows with previously troublesome embedded high-condition-number regions.

Saffman-Taylor Instability for a non-Newtonian fluid

NASA Astrophysics Data System (ADS)

Daripa, Prabir

2013-11-01

Motivated by applications, we study classical Saffman-Taylor instability involving displacement of an Oldroyd-B fluid displaced by air in a Hele-Shaw cell. The lubrication approximation is used by neglecting the vertical component of the velocity. We obtain an explicit expression of one of the components of the extra-stress perturbations tensor in terms of the horizontal velocity perturbations. The main result is an explicit formula for the growth constant (in time) of perturbations, given by a ratio in which a term depending on the relaxation and retardation (time) constants appears in the denominator of the ratio. This exact result compares extremely well with known numerical results. It is found that flow is more unstable than the corresponding Newtonian case. This is a joint work with Gelu Pasa. The research has been made possible by an NPRP Grant # 08-777-1-141 from the Qatar National Research Fund (a member of the Qatar Foundation).

A surface-potential-based drain current compact model for a-InGaZnO thin-film transistors in Non-Degenerate conduction regime

NASA Astrophysics Data System (ADS)

Yu, Fei; Ma, Xiaoyu; Deng, Wanling; Liou, Juin J.; Huang, Junkai

2017-11-01

A physics-based drain current compact model for amorphous InGaZnO (a-InGaZnO) thin-film transistors (TFTs) is proposed. As a key feature, the surface potential model accounts for both exponential tail and deep trap densities of states, which are essential to describe a-InGaZnO TFT electrical characteristics. The surface potential is solved explicitly without the process of amendment and suitable for circuit simulations. Furthermore, based on the surface potential, an explicit closed-form expression of the drain current is developed. For the cases of the different operational voltages, surface potential and drain current are verified by numerical results and experimental data, respectively. As a result, our model can predict DC characteristics of a-InGaZnO TFTs.

Non-hydrostatic semi-elastic hybrid-coordinate SISL extension of HIRLAM. Part II: numerical testing

NASA Astrophysics Data System (ADS)

Rõõm, Rein; Männik, Aarne; Luhamaa, Andres; Zirk, Marko

2007-10-01

The semi-implicit semi-Lagrangian (SISL), two-time-level, non-hydrostatic numerical scheme, based on the non-hydrostatic, semi-elastic pressure-coordinate equations, is tested in model experiments with flow over given orography (elliptical hill, mountain ridge, system of successive ridges) in a rectangular domain with emphasis on the numerical accuracy and non-hydrostatic effect presentation capability. Comparison demonstrates good (in strong primary wave generation) to satisfactory (in weak secondary wave reproduction in some cases) consistency of the numerical modelling results with known stationary linear test solutions. Numerical stability of the developed model is investigated with respect to the reference state choice, modelling dynamics of a stationary front. The horizontally area-mean reference temperature proves to be the optimal stability warrant. The numerical scheme with explicit residual in the vertical forcing term becomes unstable for cross-frontal temperature differences exceeding 30 K. Stability is restored, if the vertical forcing is treated implicitly, which enables to use time steps, comparable with the hydrostatic SISL.

Higher order explicit symmetric integrators for inseparable forms of coordinates and momenta

NASA Astrophysics Data System (ADS)

Liu, Lei; Wu, Xin; Huang, Guoqing; Liu, Fuyao

2016-06-01

Pihajoki proposed the extended phase-space second-order explicit symmetric leapfrog methods for inseparable Hamiltonian systems. On the basis of this work, we survey a critical problem on how to mix the variables in the extended phase space. Numerical tests show that sequent permutations of coordinates and momenta can make the leapfrog-like methods yield the most accurate results and the optimal long-term stabilized error behaviour. We also present a novel method to construct many fourth-order extended phase-space explicit symmetric integration schemes. Each scheme represents the symmetric production of six usual second-order leapfrogs without any permutations. This construction consists of four segments: the permuted coordinates, triple product of the usual second-order leapfrog without permutations, the permuted momenta and the triple product of the usual second-order leapfrog without permutations. Similarly, extended phase-space sixth, eighth and other higher order explicit symmetric algorithms are available. We used several inseparable Hamiltonian examples, such as the post-Newtonian approach of non-spinning compact binaries, to show that one of the proposed fourth-order methods is more efficient than the existing methods; examples include the fourth-order explicit symplectic integrators of Chin and the fourth-order explicit and implicit mixed symplectic integrators of Zhong et al. Given a moderate choice for the related mixing and projection maps, the extended phase-space explicit symplectic-like methods are well suited for various inseparable Hamiltonian problems. Samples of these problems involve the algorithmic regularization of gravitational systems with velocity-dependent perturbations in the Solar system and post-Newtonian Hamiltonian formulations of spinning compact objects.

Explicit continuous charge-based compact model for long channel heavily doped surrounding-gate MOSFETs incorporating interface traps and quantum effects

NASA Astrophysics Data System (ADS)

Hamzah, Afiq; Hamid, Fatimah A.; Ismail, Razali

2016-12-01

An explicit solution for long-channel surrounding-gate (SRG) MOSFETs is presented from intrinsic to heavily doped body including the effects of interface traps and fixed oxide charges. The solution is based on the core SRGMOSFETs model of the Unified Charge Control Model (UCCM) for heavily doped conditions. The UCCM model of highly doped SRGMOSFETs is derived to obtain the exact equivalent expression as in the undoped case. Taking advantage of the undoped explicit charge-based expression, the asymptotic limits for below threshold and above threshold have been redefined to include the effect of trap states for heavily doped cases. After solving the asymptotic limits, an explicit mobile charge expression is obtained which includes the trap state effects. The explicit mobile charge model shows very good agreement with respect to numerical simulation over practical terminal voltages, doping concentration, geometry effects, and trap state effects due to the fixed oxide charges and interface traps. Then, the drain current is obtained using the Pao-Sah's dual integral, which is expressed as a function of inversion charge densities at the source/drain ends. The drain current agreed well with the implicit solution and numerical simulation for all regions of operation without employing any empirical parameters. A comparison with previous explicit models has been conducted to verify the competency of the proposed model with the doping concentration of 1× {10}19 {{cm}}-3, as the proposed model has better advantages in terms of its simplicity and accuracy at a higher doping concentration.

Efficient algorithms and implementations of entropy-based moment closures for rarefied gases

DOE Office of Scientific and Technical Information (OSTI.GOV)

Schaerer, Roman Pascal, E-mail: schaerer@mathcces.rwth-aachen.de; Bansal, Pratyuksh; Torrilhon, Manuel

We present efficient algorithms and implementations of the 35-moment system equipped with the maximum-entropy closure in the context of rarefied gases. While closures based on the principle of entropy maximization have been shown to yield very promising results for moderately rarefied gas flows, the computational cost of these closures is in general much higher than for closure theories with explicit closed-form expressions of the closing fluxes, such as Grad's classical closure. Following a similar approach as Garrett et al. (2015) , we investigate efficient implementations of the computationally expensive numerical quadrature method used for the moment evaluations of the maximum-entropymore » distribution by exploiting its inherent fine-grained parallelism with the parallelism offered by multi-core processors and graphics cards. We show that using a single graphics card as an accelerator allows speed-ups of two orders of magnitude when compared to a serial CPU implementation. To accelerate the time-to-solution for steady-state problems, we propose a new semi-implicit time discretization scheme. The resulting nonlinear system of equations is solved with a Newton type method in the Lagrange multipliers of the dual optimization problem in order to reduce the computational cost. Additionally, fully explicit time-stepping schemes of first and second order accuracy are presented. We investigate the accuracy and efficiency of the numerical schemes for several numerical test cases, including a steady-state shock-structure problem.« less

Comparison of AGE and Spectral Methods for the Simulation of Far-Wakes

NASA Technical Reports Server (NTRS)

Bisset, D. K.; Rogers, M. M.; Kega, Dennis (Technical Monitor)

1999-01-01

Turbulent flow simulation methods based on finite differences are attractive for their simplicity, flexibility and efficiency, but not always for accuracy or stability. This report demonstrates that a good compromise is possible with the Advected Grid Explicit (AGE) method. AGE has proven to be both efficient and accurate for simulating turbulent free-shear flows, including planar mixing layers and planar jets. Its efficiency results from its localized fully explicit finite difference formulation (Bisset 1998a,b) that is very straightforward to compute, outweighing the need for a fairly small timestep. Also, most of the successful simulations were slightly under-resolved, and therefore they were, in effect, large-eddy simulations (LES) without a sub-grid-scale (SGS) model, rather than direct numerical simulations (DNS). The principle is that the role of the smallest scales of turbulent motion (when the Reynolds number is not too low) is to dissipate turbulent energy, and therefore they do not have to be simulated when the numerical method is inherently dissipative at its resolution limits. Such simulations are termed 'auto-LES' (LES with automatic SGS modeling) in this report.

A time-spectral approach to numerical weather prediction

NASA Astrophysics Data System (ADS)

Scheffel, Jan; Lindvall, Kristoffer; Yik, Hiu Fai

2018-05-01

Finite difference methods are traditionally used for modelling the time domain in numerical weather prediction (NWP). Time-spectral solution is an attractive alternative for reasons of accuracy and efficiency and because time step limitations associated with causal CFL-like criteria, typical for explicit finite difference methods, are avoided. In this work, the Lorenz 1984 chaotic equations are solved using the time-spectral algorithm GWRM (Generalized Weighted Residual Method). Comparisons of accuracy and efficiency are carried out for both explicit and implicit time-stepping algorithms. It is found that the efficiency of the GWRM compares well with these methods, in particular at high accuracy. For perturbative scenarios, the GWRM was found to be as much as four times faster than the finite difference methods. A primary reason is that the GWRM time intervals typically are two orders of magnitude larger than those of the finite difference methods. The GWRM has the additional advantage to produce analytical solutions in the form of Chebyshev series expansions. The results are encouraging for pursuing further studies, including spatial dependence, of the relevance of time-spectral methods for NWP modelling.

Solutions of conformal Israel-Stewart relativistic viscous fluid dynamics

NASA Astrophysics Data System (ADS)

Marrochio, Hugo; Noronha, Jorge; Denicol, Gabriel S.; Luzum, Matthew; Jeon, Sangyong; Gale, Charles

2015-01-01

We use symmetry arguments developed by Gubser to construct the first radially expanding explicit solutions of the Israel-Stewart formulation of hydrodynamics. Along with a general semi-analytical solution, an exact analytical solution is given which is valid in the cold plasma limit where viscous effects from shear viscosity and the relaxation time coefficient are important. The radially expanding solutions presented in this paper can be used as nontrivial checks of numerical algorithms employed in hydrodynamic simulations of the quark-gluon plasma formed in ultrarelativistic heavy ion collisions. We show this explicitly by comparing such analytic and semi-analytic solutions with the corresponding numerical solutions obtained using the music viscous hydrodynamics simulation code.

Implicit and semi-implicit schemes in the Versatile Advection Code: numerical tests

NASA Astrophysics Data System (ADS)

Toth, G.; Keppens, R.; Botchev, M. A.

1998-04-01

We describe and evaluate various implicit and semi-implicit time integration schemes applied to the numerical simulation of hydrodynamical and magnetohydrodynamical problems. The schemes were implemented recently in the software package Versatile Advection Code, which uses modern shock capturing methods to solve systems of conservation laws with optional source terms. The main advantage of implicit solution strategies over explicit time integration is that the restrictive constraint on the allowed time step can be (partially) eliminated, thus the computational cost is reduced. The test problems cover one and two dimensional, steady state and time accurate computations, and the solutions contain discontinuities. For each test, we confront explicit with implicit solution strategies.

Anchoring effects in the judgment of confidence: semantic or numeric priming?

PubMed

Carroll, Steven R; Petrusic, William M; Leth-Steensen, Craig

2009-02-01

Over the last decade, researchers have debated whether anchoring effects are the result of semantic or numeric priming. The present study tested both hypotheses. In four experiments involving a sensory detection task, participants first made a relative confidence judgment by deciding whether they were more or less confident than an anchor value in the correctness of their decision. Subsequently, they expressed an absolute level of confidence. In two of these experiments, the relative confidence anchor values represented the midpoints between the absolute confidence scale values, which were either explicitly numeric or semantic, nonnumeric representations of magnitude. In two other experiments, the anchor values were drawn from a scale modally different from that used to express the absolute confidence (i.e., nonnumeric and numeric, respectively, or vice versa). Regardless of the nature of the anchors, the mean confidence ratings revealed anchoring effects only when the relative and absolute confidence values were drawn from identical scales. Together, the results of these four experiments limit the conditions under which both numeric and semantic priming would be expected to lead to anchoring effects.

Explicit frequency equations of free vibration of a nonlocal Timoshenko beam with surface effects

NASA Astrophysics Data System (ADS)

Zhao, Hai-Sheng; Zhang, Yao; Lie, Seng-Tjhen

2018-02-01

Considerations of nonlocal elasticity and surface effects in micro- and nanoscale beams are both important for the accurate prediction of natural frequency. In this study, the governing equation of a nonlocal Timoshenko beam with surface effects is established by taking into account three types of boundary conditions: hinged-hinged, clamped-clamped and clamped-hinged ends. For a hinged-hinged beam, an exact and explicit natural frequency equation is obtained. However, for clamped-clamped and clamped-hinged beams, the solutions of corresponding frequency equations must be determined numerically due to their transcendental nature. Hence, the Fredholm integral equation approach coupled with a curve fitting method is employed to derive the approximate fundamental frequency equations, which can predict the frequency values with high accuracy. In short, explicit frequency equations of the Timoshenko beam for three types of boundary conditions are proposed to exhibit directly the dependence of the natural frequency on the nonlocal elasticity, surface elasticity, residual surface stress, shear deformation and rotatory inertia, avoiding the complicated numerical computation.

Mixed time integration methods for transient thermal analysis of structures

NASA Technical Reports Server (NTRS)

Liu, W. K.

1982-01-01

The computational methods used to predict and optimize the thermal structural behavior of aerospace vehicle structures are reviewed. In general, two classes of algorithms, implicit and explicit, are used in transient thermal analysis of structures. Each of these two methods has its own merits. Due to the different time scales of the mechanical and thermal responses, the selection of a time integration method can be a different yet critical factor in the efficient solution of such problems. Therefore mixed time integration methods for transient thermal analysis of structures are being developed. The computer implementation aspects and numerical evaluation of these mixed time implicit-explicit algorithms in thermal analysis of structures are presented. A computationally useful method of estimating the critical time step for linear quadrilateral element is also given. Numerical tests confirm the stability criterion and accuracy characteristics of the methods. The superiority of these mixed time methods to the fully implicit method or the fully explicit method is also demonstrated.

Mixed time integration methods for transient thermal analysis of structures

NASA Technical Reports Server (NTRS)

Liu, W. K.

1983-01-01

The computational methods used to predict and optimize the thermal-structural behavior of aerospace vehicle structures are reviewed. In general, two classes of algorithms, implicit and explicit, are used in transient thermal analysis of structures. Each of these two methods has its own merits. Due to the different time scales of the mechanical and thermal responses, the selection of a time integration method can be a difficult yet critical factor in the efficient solution of such problems. Therefore mixed time integration methods for transient thermal analysis of structures are being developed. The computer implementation aspects and numerical evaluation of these mixed time implicit-explicit algorithms in thermal analysis of structures are presented. A computationally-useful method of estimating the critical time step for linear quadrilateral element is also given. Numerical tests confirm the stability criterion and accuracy characteristics of the methods. The superiority of these mixed time methods to the fully implicit method or the fully explicit method is also demonstrated.

On the development of OpenFOAM solvers based on explicit and implicit high-order Runge-Kutta schemes for incompressible flows with heat transfer

NASA Astrophysics Data System (ADS)

D'Alessandro, Valerio; Binci, Lorenzo; Montelpare, Sergio; Ricci, Renato

2018-01-01

Open-source CFD codes provide suitable environments for implementing and testing low-dissipative algorithms typically used to simulate turbulence. In this research work we developed CFD solvers for incompressible flows based on high-order explicit and diagonally implicit Runge-Kutta (RK) schemes for time integration. In particular, an iterated PISO-like procedure based on Rhie-Chow correction was used to handle pressure-velocity coupling within each implicit RK stage. For the explicit approach, a projected scheme was used to avoid the "checker-board" effect. The above-mentioned approaches were also extended to flow problems involving heat transfer. It is worth noting that the numerical technology available in the OpenFOAM library was used for space discretization. In this work, we additionally explore the reliability and effectiveness of the proposed implementations by computing several unsteady flow benchmarks; we also show that the numerical diffusion due to the time integration approach is completely canceled using the solution techniques proposed here.

Eighth-order explicit two-step hybrid methods with symmetric nodes and weights for solving orbital and oscillatory IVPs

NASA Astrophysics Data System (ADS)

Franco, J. M.; Rández, L.

The construction of new two-step hybrid (TSH) methods of explicit type with symmetric nodes and weights for the numerical integration of orbital and oscillatory second-order initial value problems (IVPs) is analyzed. These methods attain algebraic order eight with a computational cost of six or eight function evaluations per step (it is one of the lowest costs that we know in the literature) and they are optimal among the TSH methods in the sense that they reach a certain order of accuracy with minimal cost per step. The new TSH schemes also have high dispersion and dissipation orders (greater than 8) in order to be adapted to the solution of IVPs with oscillatory solutions. The numerical experiments carried out with several orbital and oscillatory problems show that the new eighth-order explicit TSH methods are more efficient than other standard TSH or Numerov-type methods proposed in the scientific literature.

Special solutions to Chazy equation

NASA Astrophysics Data System (ADS)

Varin, V. P.

2017-02-01

We consider the classical Chazy equation, which is known to be integrable in hypergeometric functions. But this solution has remained purely existential and was never used numerically. We give explicit formulas for hypergeometric solutions in terms of initial data. A special solution was found in the upper half plane H with the same tessellation of H as that of the modular group. This allowed us to derive some new identities for the Eisenstein series. We constructed a special solution in the unit disk and gave an explicit description of singularities on its natural boundary. A global solution to Chazy equation in elliptic and theta functions was found that allows parametrization of an arbitrary solution to Chazy equation. The results have applications to analytic number theory.

Implicit Kalman filtering

NASA Technical Reports Server (NTRS)

Skliar, M.; Ramirez, W. F.

1997-01-01

For an implicitly defined discrete system, a new algorithm for Kalman filtering is developed and an efficient numerical implementation scheme is proposed. Unlike the traditional explicit approach, the implicit filter can be readily applied to ill-conditioned systems and allows for generalization to descriptor systems. The implementation of the implicit filter depends on the solution of the congruence matrix equation (A1)(Px)(AT1) = Py. We develop a general iterative method for the solution of this equation, and prove necessary and sufficient conditions for convergence. It is shown that when the system matrices of an implicit system are sparse, the implicit Kalman filter requires significantly less computer time and storage to implement as compared to the traditional explicit Kalman filter. Simulation results are presented to illustrate and substantiate the theoretical developments.

IMEX HDG-DG: A coupled implicit hybridized discontinuous Galerkin and explicit discontinuous Galerkin approach for Euler systems on cubed sphere.

NASA Astrophysics Data System (ADS)

Kang, S.; Muralikrishnan, S.; Bui-Thanh, T.

2017-12-01

We propose IMEX HDG-DG schemes for Euler systems on cubed sphere. Of interest is subsonic flow, where the speed of the acoustic wave is faster than that of the nonlinear advection. In order to simulate these flows efficiently, we split the governing system into stiff part describing the fast waves and non-stiff part associated with nonlinear advection. The former is discretized implicitly with HDG method while explicit Runge-Kutta DG discretization is employed for the latter. The proposed IMEX HDG-DG framework: 1) facilitates high-order solution both in time and space; 2) avoids overly small time stepsizes; 3) requires only one linear system solve per time step; and 4) relatively to DG generates smaller and sparser linear system while promoting further parallelism owing to HDG discretization. Numerical results for various test cases demonstrate that our methods are comparable to explicit Runge-Kutta DG schemes in terms of accuracy, while allowing for much larger time stepsizes.

An Explicit Algorithm for the Simulation of Fluid Flow through Porous Media

NASA Astrophysics Data System (ADS)

Trapeznikova, Marina; Churbanova, Natalia; Lyupa, Anastasiya

2018-02-01

The work deals with the development of an original mathematical model of porous medium flow constructed by analogy with the quasigasdynamic system of equations and allowing implementation via explicit numerical methods. The model is generalized to the case of multiphase multicomponent fluid and takes into account possible heat sources. The proposed approach is verified by a number of test predictions.

Semiclassical evaluation of quantum fidelity

NASA Astrophysics Data System (ADS)

Vaníček, Jiří; Heller, Eric J.

2003-11-01

We present a numerically feasible semiclassical (SC) method to evaluate quantum fidelity decay (Loschmidt echo) in a classically chaotic system. It was thought that such evaluation would be intractable, but instead we show that a uniform SC expression not only is tractable but it also gives remarkably accurate numerical results for the standard map in both the Fermi-golden-rule and Lyapunov regimes. Because it allows Monte Carlo evaluation, the uniform expression is accurate at times when there are 1070 semiclassical contributions. Remarkably, it also explicitly contains the “building blocks” of analytical theories of recent literature, and thus permits a direct test of the approximations made by other authors in these regimes, rather than an a posteriori comparison with numerical results. We explain in more detail the extended validity of the classical perturbation approximation and show that within this approximation, the so-called “diagonal approximation” is automatic and does not require ensemble averaging.

Numerical solution of the full potential equation using a chimera grid approach

NASA Technical Reports Server (NTRS)

Holst, Terry L.

1995-01-01

A numerical scheme utilizing a chimera zonal grid approach for solving the full potential equation in two spatial dimensions is described. Within each grid zone a fully-implicit approximate factorization scheme is used to advance the solution one interaction. This is followed by the explicit advance of all common zonal grid boundaries using a bilinear interpolation of the velocity potential. The presentation is highlighted with numerical results simulating the flow about a two-dimensional, nonlifting, circular cylinder. For this problem, the flow domain is divided into two parts: an inner portion covered by a polar grid and an outer portion covered by a Cartesian grid. Both incompressible and compressible (transonic) flow solutions are included. Comparisons made with an analytic solution as well as single grid results indicate that the chimera zonal grid approach is a viable technique for solving the full potential equation.

Explicitly computing geodetic coordinates from Cartesian coordinates

NASA Astrophysics Data System (ADS)

Zeng, Huaien

2013-04-01

This paper presents a new form of quartic equation based on Lagrange's extremum law and a Groebner basis under the constraint that the geodetic height is the shortest distance between a given point and the reference ellipsoid. A very explicit and concise formulae of the quartic equation by Ferrari's line is found, which avoids the need of a good starting guess for iterative methods. A new explicit algorithm is then proposed to compute geodetic coordinates from Cartesian coordinates. The convergence region of the algorithm is investigated and the corresponding correct solution is given. Lastly, the algorithm is validated with numerical experiments.

Numerical Algorithm for Delta of Asian Option

PubMed Central

Zhang, Boxiang; Yu, Yang; Wang, Weiguo

2015-01-01

We study the numerical solution of the Greeks of Asian options. In particular, we derive a close form solution of Δ of Asian geometric option and use this analytical form as a control to numerically calculate Δ of Asian arithmetic option, which is known to have no explicit close form solution. We implement our proposed numerical method and compare the standard error with other classical variance reduction methods. Our method provides an efficient solution to the hedging strategy with Asian options. PMID:26266271

Numerical calculations of two dimensional, unsteady transonic flows with circulation

NASA Technical Reports Server (NTRS)

Beam, R. M.; Warming, R. F.

1974-01-01

The feasibility of obtaining two-dimensional, unsteady transonic aerodynamic data by numerically integrating the Euler equations is investigated. An explicit, third-order-accurate, noncentered, finite-difference scheme is used to compute unsteady flows about airfoils. Solutions for lifting and nonlifting airfoils are presented and compared with subsonic linear theory. The applicability and efficiency of the numerical indicial function method are outlined. Numerically computed subsonic and transonic oscillatory aerodynamic coefficients are presented and compared with those obtained from subsonic linear theory and transonic wind-tunnel data.

Joint multifractal analysis based on the partition function approach: analytical analysis, numerical simulation and empirical application

NASA Astrophysics Data System (ADS)

Xie, Wen-Jie; Jiang, Zhi-Qiang; Gu, Gao-Feng; Xiong, Xiong; Zhou, Wei-Xing

2015-10-01

Many complex systems generate multifractal time series which are long-range cross-correlated. Numerous methods have been proposed to characterize the multifractal nature of these long-range cross correlations. However, several important issues about these methods are not well understood and most methods consider only one moment order. We study the joint multifractal analysis based on partition function with two moment orders, which was initially invented to investigate fluid fields, and derive analytically several important properties. We apply the method numerically to binomial measures with multifractal cross correlations and bivariate fractional Brownian motions without multifractal cross correlations. For binomial multifractal measures, the explicit expressions of mass function, singularity strength and multifractal spectrum of the cross correlations are derived, which agree excellently with the numerical results. We also apply the method to stock market indexes and unveil intriguing multifractality in the cross correlations of index volatilities.

From stochastic processes to numerical methods: A new scheme for solving reaction subdiffusion fractional partial differential equations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Angstmann, C.N.; Donnelly, I.C.; Henry, B.I., E-mail: B.Henry@unsw.edu.au

We have introduced a new explicit numerical method, based on a discrete stochastic process, for solving a class of fractional partial differential equations that model reaction subdiffusion. The scheme is derived from the master equations for the evolution of the probability density of a sum of discrete time random walks. We show that the diffusion limit of the master equations recovers the fractional partial differential equation of interest. This limiting procedure guarantees the consistency of the numerical scheme. The positivity of the solution and stability results are simply obtained, provided that the underlying process is well posed. We also showmore » that the method can be applied to standard reaction–diffusion equations. This work highlights the broader applicability of using discrete stochastic processes to provide numerical schemes for partial differential equations, including fractional partial differential equations.« less

Large Eddy simulation of compressible flows with a low-numerical dissipation patch-based adaptive mesh refinement method

NASA Astrophysics Data System (ADS)

Pantano, Carlos

2005-11-01

We describe a hybrid finite difference method for large-eddy simulation (LES) of compressible flows with a low-numerical dissipation scheme and structured adaptive mesh refinement (SAMR). Numerical experiments and validation calculations are presented including a turbulent jet and the strongly shock-driven mixing of a Richtmyer-Meshkov instability. The approach is a conservative flux-based SAMR formulation and as such, it utilizes refinement to computational advantage. The numerical method for the resolved scale terms encompasses the cases of scheme alternation and internal mesh interfaces resulting from SAMR. An explicit centered scheme that is consistent with a skew-symmetric finite difference formulation is used in turbulent flow regions while a weighted essentially non-oscillatory (WENO) scheme is employed to capture shocks. The subgrid stresses and transports are calculated by means of the streched-vortex model, Misra & Pullin (1997)

Numerical Modeling of Unsteady Thermofluid Dynamics in Cryogenic Systems

NASA Technical Reports Server (NTRS)

Majumdar, Alok

2003-01-01

A finite volume based network analysis procedure has been applied to model unsteady flow without and with heat transfer. Liquid has been modeled as compressible fluid where the compressibility factor is computed from the equation of state for a real fluid. The modeling approach recognizes that the pressure oscillation is linked with the variation of the compressibility factor; therefore, the speed of sound does not explicitly appear in the governing equations. The numerical results of chilldown process also suggest that the flow and heat transfer are strongly coupled. This is evident by observing that the mass flow rate during 90-second chilldown process increases by factor of ten.

Modeling of transient heat pipe operation

NASA Technical Reports Server (NTRS)

Colwell, G. T.; Hartley, J. G.

1986-01-01

Mathematical models and associated solution procedures which can be used to design heat pipe cooled structures for use on hypersonic vehicles are being developed. The models should also have the capability to predict off-design performance for a variety of operating conditions. It is expected that the resulting models can be used to predict startup behavior of liquid metal heat pipes to be used in reentry vehicles, hypersonic aircraft, and space nuclear reactors. Work to date related to numerical solutions of governing differential equations for the outer shell and the combination capillary structure and working fluid is summarized. Finite element numerical equations using both implicit, explicit, and combination methods were examined.

Power corrections in the N -jettiness subtraction scheme

DOE Office of Scientific and Technical Information (OSTI.GOV)

Boughezal, Radja; Liu, Xiaohui; Petriello, Frank

We discuss the leading-logarithmic power corrections in the N-jettiness subtraction scheme for higher-order perturbative QCD calculations. We compute the next-to-leading order power corrections for an arbitrary N-jet process, and we explicitly calculate the power correction through next-to-next-to-leading order for color-singlet production for bothmore » $$q\\bar{q}$$ and gg initiated processes. Our results are compact and simple to implement numerically. Including the leading power correction in the N-jettiness subtraction scheme substantially improves its numerical efficiency. Finally, we discuss what features of our techniques extend to processes containing final-state jets.« less

Power corrections in the N -jettiness subtraction scheme

DOE PAGES

Boughezal, Radja; Liu, Xiaohui; Petriello, Frank

2017-03-30

We discuss the leading-logarithmic power corrections in the N-jettiness subtraction scheme for higher-order perturbative QCD calculations. We compute the next-to-leading order power corrections for an arbitrary N-jet process, and we explicitly calculate the power correction through next-to-next-to-leading order for color-singlet production for bothmore » $$q\\bar{q}$$ and gg initiated processes. Our results are compact and simple to implement numerically. Including the leading power correction in the N-jettiness subtraction scheme substantially improves its numerical efficiency. Finally, we discuss what features of our techniques extend to processes containing final-state jets.« less

High-Order Space-Time Methods for Conservation Laws

NASA Technical Reports Server (NTRS)

Huynh, H. T.

2013-01-01

Current high-order methods such as discontinuous Galerkin and/or flux reconstruction can provide effective discretization for the spatial derivatives. Together with a time discretization, such methods result in either too small a time step size in the case of an explicit scheme or a very large system in the case of an implicit one. To tackle these problems, two new high-order space-time schemes for conservation laws are introduced: the first is explicit and the second, implicit. The explicit method here, also called the moment scheme, achieves a Courant-Friedrichs-Lewy (CFL) condition of 1 for the case of one-spatial dimension regardless of the degree of the polynomial approximation. (For standard explicit methods, if the spatial approximation is of degree p, then the time step sizes are typically proportional to 1/p(exp 2)). Fourier analyses for the one and two-dimensional cases are carried out. The property of super accuracy (or super convergence) is discussed. The implicit method is a simplified but optimal version of the discontinuous Galerkin scheme applied to time. It reduces to a collocation implicit Runge-Kutta (RK) method for ordinary differential equations (ODE) called Radau IIA. The explicit and implicit schemes are closely related since they employ the same intermediate time levels, and the former can serve as a key building block in an iterative procedure for the latter. A limiting technique for the piecewise linear scheme is also discussed. The technique can suppress oscillations near a discontinuity while preserving accuracy near extrema. Preliminary numerical results are shown

Explicit and implicit springback simulation in sheet metal forming using fully coupled ductile damage and distortional hardening model

NASA Astrophysics Data System (ADS)

Yetna n'jock, M.; Houssem, B.; Labergere, C.; Saanouni, K.; Zhenming, Y.

2018-05-01

The springback is an important phenomenon which accompanies the forming of metallic sheets especially for high strength materials. A quantitative prediction of springback becomes very important for newly developed material with high mechanical characteristics. In this work, a numerical methodology is developed to quantify this undesirable phenomenon. This methodoly is based on the use of both explicit and implicit finite element solvers of Abaqus®. The most important ingredient of this methodology consists on the use of highly predictive mechanical model. A thermodynamically-consistent, non-associative and fully anisotropic elastoplastic constitutive model strongly coupled with isotropic ductile damage and accounting for distortional hardening is then used. An algorithm for local integration of the complete set of the constitutive equations is developed. This algorithm considers the rotated frame formulation (RFF) to ensure the incremental objectivity of the model in the framework of finite strains. This algorithm is implemented in both explicit (Abaqus/Explicit®) and implicit (Abaqus/Standard®) solvers of Abaqus® through the users routine VUMAT and UMAT respectively. The implicit solver of Abaqus® has been used to study spingback as it is generally a quasi-static unloading. In order to compare the methods `efficiency, the explicit method (Dynamic Relaxation Method) proposed by Rayleigh has been also used for springback prediction. The results obtained within U draw/bending benchmark are studied, discussed and compared with experimental results as reference. Finally, the purpose of this work is to evaluate the reliability of different methods predict efficiently springback in sheet metal forming.

High order volume-preserving algorithms for relativistic charged particles in general electromagnetic fields

NASA Astrophysics Data System (ADS)

He, Yang; Sun, Yajuan; Zhang, Ruili; Wang, Yulei; Liu, Jian; Qin, Hong

2016-09-01

We construct high order symmetric volume-preserving methods for the relativistic dynamics of a charged particle by the splitting technique with processing. By expanding the phase space to include the time t, we give a more general construction of volume-preserving methods that can be applied to systems with time-dependent electromagnetic fields. The newly derived methods provide numerical solutions with good accuracy and conservative properties over long time of simulation. Furthermore, because of the use of an accuracy-enhancing processing technique, the explicit methods obtain high-order accuracy and are more efficient than the methods derived from standard compositions. The results are verified by the numerical experiments. Linear stability analysis of the methods shows that the high order processed method allows larger time step size in numerical integrations.

Photon-number statistics in resonance fluorescence

NASA Astrophysics Data System (ADS)

Lenstra, D.

1982-12-01

The theory of photon-number statistics in resonance fluorescence is treated, starting with the general formula for the emission probability of n photons during a given time interval T. The results fully confirm formerly obtained results by Cook that were based on the theory of atomic motion in a traveling wave. General expressions for the factorial moments are derived and explicit results for the mean and the variance are given. It is explicitly shown that the distribution function tends to a Gaussian when T becomes much larger than the natural lifetime of the excited atom. The speed of convergence towards the Gaussian is found to be typically slow, that is, the third normalized central moment (or the skewness) is proportional to T-12. However, numerical results illustrate that the overall features of the distribution function are already well represented by a Gaussian when T is larger than a few natural lifetimes only, at least if the intensity of the exciting field is not too small and its detuning is not too large.

A self-learning rule base for command following in dynamical systems

NASA Technical Reports Server (NTRS)

Tsai, Wei K.; Lee, Hon-Mun; Parlos, Alexander

1992-01-01

In this paper, a self-learning Rule Base for command following in dynamical systems is presented. The learning is accomplished though reinforcement learning using an associative memory called SAM. The main advantage of SAM is that it is a function approximator with explicit storage of training samples. A learning algorithm patterned after the dynamic programming is proposed. Two artificially created, unstable dynamical systems are used for testing, and the Rule Base was used to generate a feedback control to improve the command following ability of the otherwise uncontrolled systems. The numerical results are very encouraging. The controlled systems exhibit a more stable behavior and a better capability to follow reference commands. The rules resulting from the reinforcement learning are explicitly stored and they can be modified or augmented by human experts. Due to overlapping storage scheme of SAM, the stored rules are similar to fuzzy rules.

Helicopter time-domain electromagnetic numerical simulation based on Leapfrog ADI-FDTD

NASA Astrophysics Data System (ADS)

Guan, S.; Ji, Y.; Li, D.; Wu, Y.; Wang, A.

2017-12-01

We present a three-dimension (3D) Alternative Direction Implicit Finite-Difference Time-Domain (Leapfrog ADI-FDTD) method for the simulation of helicopter time-domain electromagnetic (HTEM) detection. This method is different from the traditional explicit FDTD, or ADI-FDTD. Comparing with the explicit FDTD, leapfrog ADI-FDTD algorithm is no longer limited by Courant-Friedrichs-Lewy(CFL) condition. Thus, the time step is longer. Comparing with the ADI-FDTD, we reduce the equations from 12 to 6 and .the Leapfrog ADI-FDTD method will be easier for the general simulation. First, we determine initial conditions which are adopted from the existing method presented by Wang and Tripp(1993). Second, we derive Maxwell equation using a new finite difference equation by Leapfrog ADI-FDTD method. The purpose is to eliminate sub-time step and retain unconditional stability characteristics. Third, we add the convolution perfectly matched layer (CPML) absorbing boundary condition into the leapfrog ADI-FDTD simulation and study the absorbing effect of different parameters. Different absorbing parameters will affect the absorbing ability. We find the suitable parameters after many numerical experiments. Fourth, We compare the response with the 1-Dnumerical result method for a homogeneous half-space to verify the correctness of our algorithm.When the model contains 107*107*53 grid points, the conductivity is 0.05S/m. The results show that Leapfrog ADI-FDTD need less simulation time and computer storage space, compared with ADI-FDTD. The calculation speed decreases nearly four times, memory occupation decreases about 32.53%. Thus, this algorithm is more efficient than the conventional ADI-FDTD method for HTEM detection, and is more precise than that of explicit FDTD in the late time.

An explicit closed-form analytical solution for European options under the CGMY model

NASA Astrophysics Data System (ADS)

Chen, Wenting; Du, Meiyu; Xu, Xiang

2017-01-01

In this paper, we consider the analytical pricing of European path-independent options under the CGMY model, which is a particular type of pure jump Le´vy process, and agrees well with many observed properties of the real market data by allowing the diffusions and jumps to have both finite and infinite activity and variation. It is shown that, under this model, the option price is governed by a fractional partial differential equation (FPDE) with both the left-side and right-side spatial-fractional derivatives. In comparison to derivatives of integer order, fractional derivatives at a point not only involve properties of the function at that particular point, but also the information of the function in a certain subset of the entire domain of definition. This ;globalness; of the fractional derivatives has added an additional degree of difficulty when either analytical methods or numerical solutions are attempted. Albeit difficult, we still have managed to derive an explicit closed-form analytical solution for European options under the CGMY model. Based on our solution, the asymptotic behaviors of the option price and the put-call parity under the CGMY model are further discussed. Practically, a reliable numerical evaluation technique for the current formula is proposed. With the numerical results, some analyses of impacts of four key parameters of the CGMY model on European option prices are also provided.

SoftWAXS: a computational tool for modeling wide-angle X-ray solution scattering from biomolecules.

PubMed

Bardhan, Jaydeep; Park, Sanghyun; Makowski, Lee

2009-10-01

This paper describes a computational approach to estimating wide-angle X-ray solution scattering (WAXS) from proteins, which has been implemented in a computer program called SoftWAXS. The accuracy and efficiency of SoftWAXS are analyzed for analytically solvable model problems as well as for proteins. Key features of the approach include a numerical procedure for performing the required spherical averaging and explicit representation of the solute-solvent boundary and the surface of the hydration layer. These features allow the Fourier transform of the excluded volume and hydration layer to be computed directly and with high accuracy. This approach will allow future investigation of different treatments of the electron density in the hydration shell. Numerical results illustrate the differences between this approach to modeling the excluded volume and a widely used model that treats the excluded-volume function as a sum of Gaussians representing the individual atomic excluded volumes. Comparison of the results obtained here with those from explicit-solvent molecular dynamics clarifies shortcomings inherent to the representation of solvent as a time-averaged electron-density profile. In addition, an assessment is made of how the calculated scattering patterns depend on input parameters such as the solute-atom radii, the width of the hydration shell and the hydration-layer contrast. These results suggest that obtaining predictive calculations of high-resolution WAXS patterns may require sophisticated treatments of solvent.

Fourth order difference methods for hyperbolic IBVP's

NASA Technical Reports Server (NTRS)

Gustafsson, Bertil; Olsson, Pelle

1994-01-01

Fourth order difference approximations of initial-boundary value problems for hyperbolic partial differential equations are considered. We use the method of lines approach with both explicit and compact implicit difference operators in space. The explicit operator satisfies an energy estimate leading to strict stability. For the implicit operator we develop boundary conditions and give a complete proof of strong stability using the Laplace transform technique. We also present numerical experiments for the linear advection equation and Burgers' equation with discontinuities in the solution or in its derivative. The first equation is used for modeling contact discontinuities in fluid dynamics, the second one for modeling shocks and rarefaction waves. The time discretization is done with a third order Runge-Kutta TVD method. For solutions with discontinuities in the solution itself we add a filter based on second order viscosity. In case of the non-linear Burger's equation we use a flux splitting technique that results in an energy estimate for certain different approximations, in which case also an entropy condition is fulfilled. In particular we shall demonstrate that the unsplit conservative form produces a non-physical shock instead of the physically correct rarefaction wave. In the numerical experiments we compare our fourth order methods with a standard second order one and with a third order TVD-method. The results show that the fourth order methods are the only ones that give good results for all the considered test problems.

An explicit solution to the exoatmospheric powered flight guidance and trajectory optimization problem for rocket propelled vehicles

NASA Technical Reports Server (NTRS)

Jaggers, R. F.

1977-01-01

A derivation of an explicit solution to the two point boundary-value problem of exoatmospheric guidance and trajectory optimization is presented. Fixed initial conditions and continuous burn, multistage thrusting are assumed. Any number of end conditions from one to six (throttling is required in the case of six) can be satisfied in an explicit and practically optimal manner. The explicit equations converge for off nominal conditions such as engine failure, abort, target switch, etc. The self starting, predictor/corrector solution involves no Newton-Rhapson iterations, numerical integration, or first guess values, and converges rapidly if physically possible. A form of this algorithm has been chosen for onboard guidance, as well as real time and preflight ground targeting and trajectory shaping for the NASA Space Shuttle Program.

Symmetry boost of the fidelity of Shor factoring

NASA Astrophysics Data System (ADS)

Nam, Y. S.; Blümel, R.

2018-05-01

In Shor's algorithm quantum subroutines occur with the structure F U F-1 , where F is a unitary transform and U is performing a quantum computation. Examples are quantum adders and subunits of quantum modulo adders. In this paper we show, both analytically and numerically, that if, in analogy to spin echoes, F and F-1 can be implemented symmetrically when executing Shor's algorithm on actual, imperfect quantum hardware, such that F and F-1 have the same hardware errors, a symmetry boost in the fidelity of the combined F U F-1 quantum operation results when compared to the case in which the errors in F and F-1 are independently random. Running the complete gate-by-gate implemented Shor algorithm, we show that the symmetry-induced fidelity boost can be as large as a factor 4. While most of our analytical and numerical results concern the case of over- and under-rotation of controlled rotation gates, in the numerically accessible case of Shor's algorithm with a small number of qubits, we show explicitly that the symmetry boost is robust with respect to more general types of errors. While, expectedly, additional error types reduce the symmetry boost, we show explicitly, by implementing general off-diagonal SU (N ) errors (N =2 ,4 ,8 ), that the boost factor scales like a Lorentzian in δ /σ , where σ and δ are the error strengths of the diagonal over- and underrotation errors and the off-diagonal SU (N ) errors, respectively. The Lorentzian shape also shows that, while the boost factor may become small with increasing δ , it declines slowly (essentially like a power law) and is never completely erased. We also investigate the effect of diagonal nonunitary errors, which, in analogy to unitary errors, reduce but never erase the symmetry boost. Going beyond the case of small quantum processors, we present analytical scaling results that show that the symmetry boost persists in the practically interesting case of a large number of qubits. We illustrate this result explicitly for the case of Shor factoring of the semiprime RSA-1024, where, analytically, focusing on over- and underrotation errors, we obtain a boost factor of about 10. In addition, we provide a proof of the fidelity product formula, including its range of applicability.

Dual methods and approximation concepts in structural synthesis

NASA Technical Reports Server (NTRS)

Fleury, C.; Schmit, L. A., Jr.

1980-01-01

Approximation concepts and dual method algorithms are combined to create a method for minimum weight design of structural systems. Approximation concepts convert the basic mathematical programming statement of the structural synthesis problem into a sequence of explicit primal problems of separable form. These problems are solved by constructing explicit dual functions, which are maximized subject to nonnegativity constraints on the dual variables. It is shown that the joining together of approximation concepts and dual methods can be viewed as a generalized optimality criteria approach. The dual method is successfully extended to deal with pure discrete and mixed continuous-discrete design variable problems. The power of the method presented is illustrated with numerical results for example problems, including a metallic swept wing and a thin delta wing with fiber composite skins.

Constraint treatment techniques and parallel algorithms for multibody dynamic analysis. Ph.D. Thesis

NASA Technical Reports Server (NTRS)

Chiou, Jin-Chern

1990-01-01

Computational procedures for kinematic and dynamic analysis of three-dimensional multibody dynamic (MBD) systems are developed from the differential-algebraic equations (DAE's) viewpoint. Constraint violations during the time integration process are minimized and penalty constraint stabilization techniques and partitioning schemes are developed. The governing equations of motion, a two-stage staggered explicit-implicit numerical algorithm, are treated which takes advantage of a partitioned solution procedure. A robust and parallelizable integration algorithm is developed. This algorithm uses a two-stage staggered central difference algorithm to integrate the translational coordinates and the angular velocities. The angular orientations of bodies in MBD systems are then obtained by using an implicit algorithm via the kinematic relationship between Euler parameters and angular velocities. It is shown that the combination of the present solution procedures yields a computationally more accurate solution. To speed up the computational procedures, parallel implementation of the present constraint treatment techniques, the two-stage staggered explicit-implicit numerical algorithm was efficiently carried out. The DAE's and the constraint treatment techniques were transformed into arrowhead matrices to which Schur complement form was derived. By fully exploiting the sparse matrix structural analysis techniques, a parallel preconditioned conjugate gradient numerical algorithm is used to solve the systems equations written in Schur complement form. A software testbed was designed and implemented in both sequential and parallel computers. This testbed was used to demonstrate the robustness and efficiency of the constraint treatment techniques, the accuracy of the two-stage staggered explicit-implicit numerical algorithm, and the speed up of the Schur-complement-based parallel preconditioned conjugate gradient algorithm on a parallel computer.

Goal-oriented explicit residual-type error estimates in XFEM

NASA Astrophysics Data System (ADS)

Rüter, Marcus; Gerasimov, Tymofiy; Stein, Erwin

2013-08-01

A goal-oriented a posteriori error estimator is derived to control the error obtained while approximately evaluating a quantity of engineering interest, represented in terms of a given linear or nonlinear functional, using extended finite elements of Q1 type. The same approximation method is used to solve the dual problem as required for the a posteriori error analysis. It is shown that for both problems to be solved numerically the same singular enrichment functions can be used. The goal-oriented error estimator presented can be classified as explicit residual type, i.e. the residuals of the approximations are used directly to compute upper bounds on the error of the quantity of interest. This approach therefore extends the explicit residual-type error estimator for classical energy norm error control as recently presented in Gerasimov et al. (Int J Numer Meth Eng 90:1118-1155, 2012a). Without loss of generality, the a posteriori error estimator is applied to the model problem of linear elastic fracture mechanics. Thus, emphasis is placed on the fracture criterion, here the J-integral, as the chosen quantity of interest. Finally, various illustrative numerical examples are presented where, on the one hand, the error estimator is compared to its finite element counterpart and, on the other hand, improved enrichment functions, as introduced in Gerasimov et al. (2012b), are discussed.

Engine dynamic analysis with general nonlinear finite element codes. Part 2: Bearing element implementation overall numerical characteristics and benchmaking

NASA Technical Reports Server (NTRS)

Padovan, J.; Adams, M.; Fertis, J.; Zeid, I.; Lam, P.

1982-01-01

Finite element codes are used in modelling rotor-bearing-stator structure common to the turbine industry. Engine dynamic simulation is used by developing strategies which enable the use of available finite element codes. benchmarking the elements developed are benchmarked by incorporation into a general purpose code (ADINA); the numerical characteristics of finite element type rotor-bearing-stator simulations are evaluated through the use of various types of explicit/implicit numerical integration operators. Improving the overall numerical efficiency of the procedure is improved.

A three-dimensional finite-volume Eulerian-Lagrangian Localized Adjoint Method (ELLAM) for solute-transport modeling

USGS Publications Warehouse

Heberton, C.I.; Russell, T.F.; Konikow, Leonard F.; Hornberger, G.Z.

2000-01-01

This report documents the U.S. Geological Survey Eulerian-Lagrangian Localized Adjoint Method (ELLAM) algorithm that solves an integral form of the solute-transport equation, incorporating an implicit-in-time difference approximation for the dispersive and sink terms. Like the algorithm in the original version of the U.S. Geological Survey MOC3D transport model, ELLAM uses a method of characteristics approach to solve the transport equation on the basis of the velocity field. The ELLAM algorithm, however, is based on an integral formulation of conservation of mass and uses appropriate numerical techniques to obtain global conservation of mass. The implicit procedure eliminates several stability criteria required for an explicit formulation. Consequently, ELLAM allows large transport time increments to be used. ELLAM can produce qualitatively good results using a small number of transport time steps. A description of the ELLAM numerical method, the data-input requirements and output options, and the results of simulator testing and evaluation are presented. The ELLAM algorithm was evaluated for the same set of problems used to test and evaluate Version 1 and Version 2 of MOC3D. These test results indicate that ELLAM offers a viable alternative to the explicit and implicit solvers in MOC3D. Its use is desirable when mass balance is imperative or a fast, qualitative model result is needed. Although accurate solutions can be generated using ELLAM, its efficiency relative to the two previously documented solution algorithms is problem dependent.

Molecular orientation in a dielectric liquid-vapor interphase

NASA Astrophysics Data System (ADS)

Chacón, E.; Mederos, L.; Navascués, G.; Tarazona, P.

1985-04-01

The density functional theory of Chacón et al. is used to study the molecular orientation in an interphase of a weak dipolar fluid. Explicit expressions are obtained using standard perturbation techniques. Molecular orientation, local susceptibility, and the Gibbsean surface susceptibility are evaluated for a Stockmayer model of dipolar fluid. The effect of the surface structure on the bulk ferroelectric transition is discussed in the light of the present theory and the numerical results.

Triple collinear emissions in parton showers

DOE Office of Scientific and Technical Information (OSTI.GOV)

Höche, Stefan; Prestel, Stefan

2017-10-01

A framework to include triple collinear splitting functions into parton showers is presented, and the implementation of flavor-changing NLO splitting kernels is discussed as a first application. The correspondence between the Monte-Carlo integration and the analytic computation of NLO DGLAP evolution kernels is made explicit for both timelike and spacelike parton evolution. Numerical simulation results are obtained with two independent implementations of the new algorithm, using the two independent event generation frameworks Pythia and Sherpa.

Efficient and accurate numerical schemes for a hydro-dynamically coupled phase field diblock copolymer model

NASA Astrophysics Data System (ADS)

Cheng, Qing; Yang, Xiaofeng; Shen, Jie

2017-07-01

In this paper, we consider numerical approximations of a hydro-dynamically coupled phase field diblock copolymer model, in which the free energy contains a kinetic potential, a gradient entropy, a Ginzburg-Landau double well potential, and a long range nonlocal type potential. We develop a set of second order time marching schemes for this system using the "Invariant Energy Quadratization" approach for the double well potential, the projection method for the Navier-Stokes equation, and a subtle implicit-explicit treatment for the stress and convective term. The resulting schemes are linear and lead to symmetric positive definite systems at each time step, thus they can be efficiently solved. We further prove that these schemes are unconditionally energy stable. Various numerical experiments are performed to validate the accuracy and energy stability of the proposed schemes.

Fitted Fourier-pseudospectral methods for solving a delayed reaction-diffusion partial differential equation in biology

NASA Astrophysics Data System (ADS)

Adam, A. M. A.; Bashier, E. B. M.; Hashim, M. H. A.; Patidar, K. C.

2017-07-01

In this work, we design and analyze a fitted numerical method to solve a reaction-diffusion model with time delay, namely, a delayed version of a population model which is an extension of the logistic growth (LG) equation for a food-limited population proposed by Smith [F.E. Smith, Population dynamics in Daphnia magna and a new model for population growth, Ecology 44 (1963) 651-663]. Seeing that the analytical solution (in closed form) is hard to obtain, we seek for a robust numerical method. The method consists of a Fourier-pseudospectral semi-discretization in space and a fitted operator implicit-explicit scheme in temporal direction. The proposed method is analyzed for convergence and we found that it is unconditionally stable. Illustrative numerical results will be presented at the conference.

Numerical calculation of the decay widths, the decay constants, and the decay energy spectra of the resonances of the delta-shell potential

NASA Astrophysics Data System (ADS)

de la Madrid, Rafael

2017-06-01

We express the resonant energies of the delta-shell potential in terms of the Lambert W function, and we calculate their decay widths and decay constants. The ensuing numerical results strengthen the interpretation of such decay widths and constants as a way to quantify the coupling between a resonance and the continuum. We calculate explicitly the decay energy spectrum of the resonances of the delta-shell potential, and we show numerically that the lineshape of such spectrum is not the same as, and can be very different from, the Breit-Wigner (Lorentzian) distribution. We argue that the standard Golden Rule cannot describe the interference of two resonances, and we show how to describe such interference by way of the decay energy spectrum of two resonant states.

A thermodynamically consistent discontinuous Galerkin formulation for interface separation

DOE PAGES

Versino, Daniele; Mourad, Hashem M.; Dávila, Carlos G.; ...

2015-07-31

Our paper describes the formulation of an interface damage model, based on the discontinuous Galerkin (DG) method, for the simulation of failure and crack propagation in laminated structures. The DG formulation avoids common difficulties associated with cohesive elements. Specifically, it does not introduce any artificial interfacial compliance and, in explicit dynamic analysis, it leads to a stable time increment size which is unaffected by the presence of stiff massless interfaces. This proposed method is implemented in a finite element setting. Convergence and accuracy are demonstrated in Mode I and mixed-mode delamination in both static and dynamic analyses. Significantly, numerical resultsmore » obtained using the proposed interface model are found to be independent of the value of the penalty factor that characterizes the DG formulation. By contrast, numerical results obtained using a classical cohesive method are found to be dependent on the cohesive penalty stiffnesses. The proposed approach is shown to yield more accurate predictions pertaining to crack propagation under mixed-mode fracture because of the advantage. Furthermore, in explicit dynamic analysis, the stable time increment size calculated with the proposed method is found to be an order of magnitude larger than the maximum allowable value for classical cohesive elements.« less

GroPBS: Fast Solver for Implicit Electrostatics of Biomolecules

PubMed Central

Bertelshofer, Franziska; Sun, Liping; Greiner, Günther; Böckmann, Rainer A.

2015-01-01

Knowledge about the electrostatic potential on the surface of biomolecules or biomembranes under physiological conditions is an important step in the attempt to characterize the physico-chemical properties of these molecules and, in particular, also their interactions with each other. Additionally, knowledge about solution electrostatics may also guide the design of molecules with specified properties. However, explicit water models come at a high computational cost, rendering them unsuitable for large design studies or for docking purposes. Implicit models with the water phase treated as a continuum require the numerical solution of the Poisson–Boltzmann equation (PBE). Here, we present a new flexible program for the numerical solution of the PBE, allowing for different geometries, and the explicit and implicit inclusion of membranes. It involves a discretization of space and the computation of the molecular surface. The PBE is solved using finite differences, the resulting set of equations is solved using a Gauss–Seidel method. It is shown for the example of the sucrose transporter ScrY that the implicit inclusion of a surrounding membrane has a strong effect also on the electrostatics within the pore region and, thus, needs to be carefully considered, e.g., in design studies on membrane proteins. PMID:26636074

Three-dimensional control of crystal growth using magnetic fields

NASA Astrophysics Data System (ADS)

Dulikravich, George S.; Ahuja, Vineet; Lee, Seungsoo

1993-07-01

Two coupled systems of partial differential equations governing three-dimensional laminar viscous flow undergoing solidification or melting under the influence of arbitrarily oriented externally applied magnetic fields have been formulated. The model accounts for arbitrary temperature dependence of physical properties including latent heat release, effects of Joule heating, magnetic field forces, and mushy region existence. On the basis of this model a numerical algorithm has been developed and implemented using central differencing on a curvilinear boundary-conforming grid and Runge-Kutta explicit time-stepping. The numerical results clearly demonstrate possibilities for active and practically instantaneous control of melt/solid interface shape, the solidification/melting front propagation speed, and the amount and location of solid accrued.

Spiral trajectory design: a flexible numerical algorithm and base analytical equations.

PubMed

Pipe, James G; Zwart, Nicholas R

2014-01-01

Spiral-based trajectories for magnetic resonance imaging can be advantageous, but are often cumbersome to design or create. This work presents a flexible numerical algorithm for designing trajectories based on explicit definition of radial undersampling, and also gives several analytical expressions for charactering the base (critically sampled) class of these trajectories. Expressions for the gradient waveform, based on slew and amplitude limits, are developed such that a desired pitch in the spiral k-space trajectory is followed. The source code for this algorithm, written in C, is publicly available. Analytical expressions approximating the spiral trajectory (ignoring the radial component) are given to characterize measurement time, gradient heating, maximum gradient amplitude, and off-resonance phase for slew-limited and gradient amplitude-limited cases. Several numerically calculated trajectories are illustrated, and base Archimedean spirals are compared with analytically obtained results. Several different waveforms illustrate that the desired slew and amplitude limits are reached, as are the desired undersampling patterns, using the numerical method. For base Archimedean spirals, the results of the numerical and analytical approaches are in good agreement. A versatile numerical algorithm was developed, and was written in publicly available code. Approximate analytical formulas are given that help characterize spiral trajectories. Copyright © 2013 Wiley Periodicals, Inc.

Implicit Space-Time Conservation Element and Solution Element Schemes

NASA Technical Reports Server (NTRS)

Chang, Sin-Chung; Himansu, Ananda; Wang, Xiao-Yen

1999-01-01

Artificial numerical dissipation is in important issue in large Reynolds number computations. In such computations, the artificial dissipation inherent in traditional numerical schemes can overwhelm the physical dissipation and yield inaccurate results on meshes of practical size. In the present work, the space-time conservation element and solution element method is used to construct new and accurate implicit numerical schemes such that artificial numerical dissipation will not overwhelm physical dissipation. Specifically, these schemes have the property that numerical dissipation vanishes when the physical viscosity goes to zero. These new schemes therefore accurately model the physical dissipation even when it is extremely small. The new schemes presented are two highly accurate implicit solvers for a convection-diffusion equation. The two schemes become identical in the pure convection case, and in the pure diffusion case. The implicit schemes are applicable over the whole Reynolds number range, from purely diffusive equations to convection-dominated equations with very small viscosity. The stability and consistency of the schemes are analysed, and some numerical results are presented. It is shown that, in the inviscid case, the new schemes become explicit and their amplification factors are identical to those of the Leapfrog scheme. On the other hand, in the pure diffusion case, their principal amplification factor becomes the amplification factor of the Crank-Nicolson scheme.

Linearized motion estimation for articulated planes.

PubMed

Datta, Ankur; Sheikh, Yaser; Kanade, Takeo

2011-04-01

In this paper, we describe the explicit application of articulation constraints for estimating the motion of a system of articulated planes. We relate articulations to the relative homography between planes and show that these articulations translate into linearized equality constraints on a linear least-squares system, which can be solved efficiently using a Karush-Kuhn-Tucker system. The articulation constraints can be applied for both gradient-based and feature-based motion estimation algorithms and to illustrate this, we describe a gradient-based motion estimation algorithm for an affine camera and a feature-based motion estimation algorithm for a projective camera that explicitly enforces articulation constraints. We show that explicit application of articulation constraints leads to numerically stable estimates of motion. The simultaneous computation of motion estimates for all of the articulated planes in a scene allows us to handle scene areas where there is limited texture information and areas that leave the field of view. Our results demonstrate the wide applicability of the algorithm in a variety of challenging real-world cases such as human body tracking, motion estimation of rigid, piecewise planar scenes, and motion estimation of triangulated meshes.

Two-level schemes for the advection equation

NASA Astrophysics Data System (ADS)

Vabishchevich, Petr N.

2018-06-01

The advection equation is the basis for mathematical models of continuum mechanics. In the approximate solution of nonstationary problems it is necessary to inherit main properties of the conservatism and monotonicity of the solution. In this paper, the advection equation is written in the symmetric form, where the advection operator is the half-sum of advection operators in conservative (divergent) and non-conservative (characteristic) forms. The advection operator is skew-symmetric. Standard finite element approximations in space are used. The standard explicit two-level scheme for the advection equation is absolutely unstable. New conditionally stable regularized schemes are constructed, on the basis of the general theory of stability (well-posedness) of operator-difference schemes, the stability conditions of the explicit Lax-Wendroff scheme are established. Unconditionally stable and conservative schemes are implicit schemes of the second (Crank-Nicolson scheme) and fourth order. The conditionally stable implicit Lax-Wendroff scheme is constructed. The accuracy of the investigated explicit and implicit two-level schemes for an approximate solution of the advection equation is illustrated by the numerical results of a model two-dimensional problem.

Explicit parametric solutions of lattice structures with proper generalized decomposition (PGD) - Applications to the design of 3D-printed architectured materials

NASA Astrophysics Data System (ADS)

Sibileau, Alberto; Auricchio, Ferdinando; Morganti, Simone; Díez, Pedro

2018-01-01

Architectured materials (or metamaterials) are constituted by a unit-cell with a complex structural design repeated periodically forming a bulk material with emergent mechanical properties. One may obtain specific macro-scale (or bulk) properties in the resulting architectured material by properly designing the unit-cell. Typically, this is stated as an optimal design problem in which the parameters describing the shape and mechanical properties of the unit-cell are selected in order to produce the desired bulk characteristics. This is especially pertinent due to the ease manufacturing of these complex structures with 3D printers. The proper generalized decomposition provides explicit parametic solutions of parametric PDEs. Here, the same ideas are used to obtain parametric solutions of the algebraic equations arising from lattice structural models. Once the explicit parametric solution is available, the optimal design problem is a simple post-process. The same strategy is applied in the numerical illustrations, first to a unit-cell (and then homogenized with periodicity conditions), and in a second phase to the complete structure of a lattice material specimen.

Explicit reference governor for linear systems

NASA Astrophysics Data System (ADS)

Garone, Emanuele; Nicotra, Marco; Ntogramatzidis, Lorenzo

2018-06-01

The explicit reference governor is a constrained control scheme that was originally introduced for generic nonlinear systems. This paper presents two explicit reference governor strategies that are specifically tailored for the constrained control of linear time-invariant systems subject to linear constraints. Both strategies are based on the idea of maintaining the system states within an invariant set which is entirely contained in the constraints. This invariant set can be constructed by exploiting either the Lyapunov inequality or modal decomposition. To improve the performance, we show that the two strategies can be combined by choosing at each time instant the least restrictive set. Numerical simulations illustrate that the proposed scheme achieves performances that are comparable to optimisation-based reference governors.

Aerodynamic coefficients in generalized unsteady thin airfoil theory

NASA Technical Reports Server (NTRS)

Williams, M. H.

1980-01-01

Two cases are considered: (1) rigid body motion of an airfoil-flap combination consisting of vertical translation of given amplitude, rotation of given amplitude about a specified axis, and rotation of given amplitude of the control surface alone about its hinge; the upwash for this problem is defined mathematically; and (2) sinusoidal gust of given amplitude and wave number, for which the upwash is defined mathematically. Simple universal formulas are presented for the most important aerodynamic coefficients in unsteady thin airfoil theory. The lift and moment induced by a generalized gust are evaluated explicitly in terms of the gust wavelength. Similarly, in the control surface problem, the lift, moment, and hinge moments are given as explicit algebraic functions of hinge location. These results can be used together with any of the standard numerical inversion routines for the elementary loads (pitch and heave).

General Retarded Contact Self-energies in and beyond the Non-equilibrium Green's Functions Method

NASA Astrophysics Data System (ADS)

Kubis, Tillmann; He, Yu; Andrawis, Robert; Klimeck, Gerhard

2016-03-01

Retarded contact self-energies in the framework of nonequilibrium Green's functions allow to model the impact of lead structures on the device without explicitly including the leads in the actual device calculation. Most of the contact self-energy algorithms are limited to homogeneous or periodic, semi-infinite lead structures. In this work, the complex absorbing potential method is extended to solve retarded contact self-energies for arbitrary lead structures, including irregular and randomly disordered leads. This method is verified for regular leads against common approaches and on physically equivalent, but numerically different irregular leads. Transmission results on randomly alloyed In0.5Ga0.5As structures show the importance of disorder in the leads. The concept of retarded contact self-energies is expanded to model passivation of atomically resolved surfaces without explicitly increasing the device's Hamiltonian.

Time-shifted synchronization of chaotic oscillator chains without explicit coupling delays.

PubMed

Blakely, Jonathan N; Stahl, Mark T; Corron, Ned J

2009-12-01

We examine chains of unidirectionally coupled oscillators in which time-shifted synchronization occurs without explicit delays in the coupling. In numerical simulations and in an experimental system of electronic oscillators, we examine the time shift and the degree of distortion (primarily in the form of attenuation) of the waveforms of the oscillators located far from the drive oscillator. Surprisingly, under weak coupling we observe minimal attenuation in spite of a significant total time shift. In contrast, at higher coupling strengths the observed attenuation increases dramatically and approaches the value predicted by an analytically derived estimate. In this regime, we verify directly that generalized synchronization is maintained over the entire chain length despite severe attenuation. These results suggest that weak coupling generally may produce higher quality synchronization in systems for which truly identical synchronization is not possible.

Virtual prototyping of drop test using explicit analysis

NASA Astrophysics Data System (ADS)

Todorov, Georgi; Kamberov, Konstantin

2017-12-01

Increased requirements for reliability and safety, included in contemporary standards and norms, has high impact over new product development. New numerical techniques based on virtual prototyping technology, facilitates imrpoving product development cycle, resutling in reduced time/money spent for this stage as well as increased knowledge about certain failure mechanism. So called "drop test" became nearly a "must" step in development of any human operated product. This study aims to demonstrate dynamic behaviour assessment of a structure under impact loads, based on virtual prototyping using a typical nonlinear analysis - explicit dynamics. An example is presneted, based on a plastic container that is used as cartridge for a dispenser machine exposed to various work conditions. Different drop orientations were analyzed and critical load cases and design weaknesses have been found. Several design modifications have been proposed, based on detailed analyses results review.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Dalvit, Diego; Messina, Riccardo; Maia Neto, Paulo

We develop the scattering approach for the dispersive force on a ground state atom on top of a corrugated surface. We present explicit results to first order in the corrugation amplitude. A variety of analytical results are derived in different limiting cases, including the van der Waals and Casimir-Polder regimes. We compute numerically the exact first-order dispersive potential for arbitrary separation distances and corrugation wavelengths, for a Rubidium atom on top of a silicon or gold corrugated surface. We consider in detail the correction to the proximity force approximation, and present a very simple approximation algorithm for computing the potential.

Two-ball Newton's cradle

NASA Astrophysics Data System (ADS)

Glendinning, Paul

2011-12-01

Newton's cradle for two balls with Hertzian interactions is considered as a hybrid system, and this makes it possible to derive return maps for the motion between collisions in an exact form despite the fact that the three-halves interaction law cannot be solved in closed form. The return maps depend on a constant whose value can only be determined numerically, but solutions can be written down explicitly in terms of this parameter, and we compare this with the results of simulations. The results are in fact independent of the details of the interaction potential.

Calculation of external-internal flow fields for mixed-compression inlets

NASA Technical Reports Server (NTRS)

Chyu, W. J.; Kawamura, T.; Bencze, D. P.

1986-01-01

Supersonic inlet flows with mixed external-internal compressions were computed using a combined implicit-explicit (Beam-Warming-Steger/MacCormack) method for solving the three-dimensional unsteady, compressible Navier-Stokes equations in conservation form. Numerical calculations were made of various flows related to such inlet operations as the shock-wave intersections, subsonic spillage around the cowl lip, and inlet started versus unstarted conditions. Some of the computed results were compared with wind tunnel data.

Calculation of external-internal flow fields for mixed-compression inlets

NASA Technical Reports Server (NTRS)

Chyu, W. J.; Kawamura, T.; Bencze, D. P.

1987-01-01

Supersonic inlet flows with mixed external-internal compressions were computed using a combined implicit-explicit (Beam-Warming-Steger/MacCormack) method for solving the three-dimensional unsteady, compressible Navier-Stokes equations in conservation form. Numerical calculations were made of various flows related to such inlet operations as the shock-wave intersections, subsonic spillage around the cowl lip, and inlet started versus unstarted conditions. Some of the computed results were compared with wind tunnel data.

Antiferromagnetism and DX2-Y2-WAVE Pairing in the Colored Hubbard Model

NASA Astrophysics Data System (ADS)

Baier, Tobias; Bick, Eike

2001-08-01

We introduce a new formulation of the 2d Hubbard model on a square lattice (the "colored" Hubbard model). In this formulation interesting physical nonlocal properties as antiferromagnetic or dx2-y2-wave superconducting behavior are included in an explicit way. Analyzing the phase diagram in a mean field approximation numerically, we show that our approach yields results which are in qualitative agreement with experiment.

Three dimensional, non-linear, finite element analysis of compactable soil interaction with a hyperelastic wheel

NASA Astrophysics Data System (ADS)

Chiroux, Robert Charles

The objective of this research was to produce a three dimensional, non-linear, dynamic simulation of the interaction between a hyperelastic wheel rolling over compactable soil. The finite element models developed to produce the simulation utilized the ABAQUS/Explicit computer code. Within the simulation two separate bodies were modeled, the hyperelastic wheel and a compactable soil-bed. Interaction between the bodies was achieved by allowing them to come in contact but not to penetrate the contact surface. The simulation included dynamic loading of a hyperelastic, rubber tire in contact with compactable soil with an applied constant angular velocity or torque, including a tow load, applied to the wheel hub. The constraints on the wheel model produced a straight and curved path. In addition the simulation included a shear limit between the tire and soil allowing for the introduction of slip. Soil properties were simulated using the Drucker-Prager, Cap Plasticity model available within the ABAQUS/Explicit program. Numerical results obtained from the three dimensional model were compared with related experimental data and showed good correlation for similar conditions. Numerical and experimental data compared well for both stress and wheel rut formation depth under a weight of 5.8 kN and a constant angular velocity applied to the wheel hub. The simulation results provided a demonstration of the benefit of three-dimensional simulation in comparison to previous two-dimensional, plane strain simulations.

A numerical spectral approach to solve the dislocation density transport equation

NASA Astrophysics Data System (ADS)

Djaka, K. S.; Taupin, V.; Berbenni, S.; Fressengeas, C.

2015-09-01

A numerical spectral approach is developed to solve in a fast, stable and accurate fashion, the quasi-linear hyperbolic transport equation governing the spatio-temporal evolution of the dislocation density tensor in the mechanics of dislocation fields. The approach relies on using the Fast Fourier Transform algorithm. Low-pass spectral filters are employed to control both the high frequency Gibbs oscillations inherent to the Fourier method and the fast-growing numerical instabilities resulting from the hyperbolic nature of the transport equation. The numerical scheme is validated by comparison with an exact solution in the 1D case corresponding to dislocation dipole annihilation. The expansion and annihilation of dislocation loops in 2D and 3D settings are also produced and compared with finite element approximations. The spectral solutions are shown to be stable, more accurate for low Courant numbers and much less computation time-consuming than the finite element technique based on an explicit Galerkin-least squares scheme.

A general numerical model for wave rotor analysis

NASA Technical Reports Server (NTRS)

Paxson, Daniel W.

1992-01-01

Wave rotors represent one of the promising technologies for achieving very high core temperatures and pressures in future gas turbine engines. Their operation depends upon unsteady gas dynamics and as such, their analysis is quite difficult. This report describes a numerical model which has been developed to perform such an analysis. Following a brief introduction, a summary of the wave rotor concept is given. The governing equations are then presented, along with a summary of the assumptions used to obtain them. Next, the numerical integration technique is described. This is an explicit finite volume technique based on the method of Roe. The discussion then focuses on the implementation of appropriate boundary conditions. Following this, some results are presented which first compare the numerical approximation to the governing differential equations and then compare the overall model to an actual wave rotor experiment. Finally, some concluding remarks are presented concerning the limitations of the simplifying assumptions and areas where the model may be improved.

Preface of "The Second Symposium on Border Zones Between Experimental and Numerical Application Including Solution Approaches By Extensions of Standard Numerical Methods"

NASA Astrophysics Data System (ADS)

Ortleb, Sigrun; Seidel, Christian

2017-07-01

In this second symposium at the limits of experimental and numerical methods, recent research is presented on practically relevant problems. Presentations discuss experimental investigation as well as numerical methods with a strong focus on application. In addition, problems are identified which require a hybrid experimental-numerical approach. Topics include fast explicit diffusion applied to a geothermal energy storage tank, noise in experimental measurements of electrical quantities, thermal fluid structure interaction, tensegrity structures, experimental and numerical methods for Chladni figures, optimized construction of hydroelectric power stations, experimental and numerical limits in the investigation of rain-wind induced vibrations as well as the application of exponential integrators in a domain-based IMEX setting.

Estimating the numerical diapycnal mixing in an eddy-permitting ocean model

NASA Astrophysics Data System (ADS)

Megann, Alex

2018-01-01

Constant-depth (or "z-coordinate") ocean models such as MOM4 and NEMO have become the de facto workhorse in climate applications, having attained a mature stage in their development and are well understood. A generic shortcoming of this model type, however, is a tendency for the advection scheme to produce unphysical numerical diapycnal mixing, which in some cases may exceed the explicitly parameterised mixing based on observed physical processes, and this is likely to have effects on the long-timescale evolution of the simulated climate system. Despite this, few quantitative estimates have been made of the typical magnitude of the effective diapycnal diffusivity due to numerical mixing in these models. GO5.0 is a recent ocean model configuration developed jointly by the UK Met Office and the National Oceanography Centre. It forms the ocean component of the GC2 climate model, and is closely related to the ocean component of the UKESM1 Earth System Model, the UK's contribution to the CMIP6 model intercomparison. GO5.0 uses version 3.4 of the NEMO model, on the ORCA025 global tripolar grid. An approach to quantifying the numerical diapycnal mixing in this model, based on the isopycnal watermass analysis of Lee et al. (2002), is described, and the estimates thereby obtained of the effective diapycnal diffusivity in GO5.0 are compared with the values of the explicit diffusivity used by the model. It is shown that the effective mixing in this model configuration is up to an order of magnitude higher than the explicit mixing in much of the ocean interior, implying that mixing in the model below the mixed layer is largely dominated by numerical mixing. This is likely to have adverse consequences for the representation of heat uptake in climate models intended for decadal climate projections, and in particular is highly relevant to the interpretation of the CMIP6 class of climate models, many of which use constant-depth ocean models at ¼° resolution

Development of a discrete gas-kinetic scheme for simulation of two-dimensional viscous incompressible and compressible flows.

PubMed

Yang, L M; Shu, C; Wang, Y

2016-03-01

In this work, a discrete gas-kinetic scheme (DGKS) is presented for simulation of two-dimensional viscous incompressible and compressible flows. This scheme is developed from the circular function-based GKS, which was recently proposed by Shu and his co-workers [L. M. Yang, C. Shu, and J. Wu, J. Comput. Phys. 274, 611 (2014)]. For the circular function-based GKS, the integrals for conservation forms of moments in the infinity domain for the Maxwellian function-based GKS are simplified to those integrals along the circle. As a result, the explicit formulations of conservative variables and fluxes are derived. However, these explicit formulations of circular function-based GKS for viscous flows are still complicated, which may not be easy for the application by new users. By using certain discrete points to represent the circle in the phase velocity space, the complicated formulations can be replaced by a simple solution process. The basic requirement is that the conservation forms of moments for the circular function-based GKS can be accurately satisfied by weighted summation of distribution functions at discrete points. In this work, it is shown that integral quadrature by four discrete points on the circle, which forms the D2Q4 discrete velocity model, can exactly match the integrals. Numerical results showed that the present scheme can provide accurate numerical results for incompressible and compressible viscous flows with roughly the same computational cost as that needed by the Roe scheme.

The Effects of Dissipation and Coarse Grid Resolution for Multigrid in Flow Problems

NASA Technical Reports Server (NTRS)

Eliasson, Peter; Engquist, Bjoern

1996-01-01

The objective of this paper is to investigate the effects of the numerical dissipation and the resolution of the solution on coarser grids for multigrid with the Euler equation approximations. The convergence is accomplished by multi-stage explicit time-stepping to steady state accelerated by FAS multigrid. A theoretical investigation is carried out for linear hyperbolic equations in one and two dimensions. The spectra reveals that for stability and hence robustness of spatial discretizations with a small amount of numerical dissipation the grid transfer operators have to be accurate enough and the smoother of low temporal accuracy. Numerical results give grid independent convergence in one dimension. For two-dimensional problems with a small amount of numerical dissipation, however, only a few grid levels contribute to an increased speed of convergence. This is explained by the small numerical dissipation leading to dispersion. Increasing the mesh density and hence making the problem over resolved increases the number of mesh levels contributing to an increased speed of convergence. If the steady state equations are elliptic, all grid levels contribute to the convergence regardless of the mesh density.

The explicit computation of integration algorithms and first integrals for ordinary differential equations with polynomials coefficients using trees

NASA Technical Reports Server (NTRS)

Crouch, P. E.; Grossman, Robert

1992-01-01

This note is concerned with the explicit symbolic computation of expressions involving differential operators and their actions on functions. The derivation of specialized numerical algorithms, the explicit symbolic computation of integrals of motion, and the explicit computation of normal forms for nonlinear systems all require such computations. More precisely, if R = k(x(sub 1),...,x(sub N)), where k = R or C, F denotes a differential operator with coefficients from R, and g member of R, we describe data structures and algorithms for efficiently computing g. The basic idea is to impose a multiplicative structure on the vector space with basis the set of finite rooted trees and whose nodes are labeled with the coefficients of the differential operators. Cancellations of two trees with r + 1 nodes translates into cancellation of O(N(exp r)) expressions involving the coefficient functions and their derivatives.

Numerical Simulation and Experimental Verification of Hollow and Foam-Filled Flax-Fabric-Reinforced Epoxy Tubular Energy Absorbers Subjected to Crashing

NASA Astrophysics Data System (ADS)

Sliseris, J.; Yan, L.; Kasal, B.

2017-09-01

Numerical methods for simulating hollow and foam-filled flax-fabric-reinforced epoxy tubular energy absorbers subjected to lateral crashing are presented. The crashing characteristics, such as the progressive failure, load-displacement response, absorbed energy, peak load, and failure modes, of the tubes were simulated and calculated numerically. A 3D nonlinear finite-element model that allows for the plasticity of materials using an isotropic hardening model with strain rate dependence and failure is proposed. An explicit finite-element solver is used to address the lateral crashing of the tubes considering large displacements and strains, plasticity, and damage. The experimental nonlinear crashing load vs. displacement data are successfully described by using the finite-element model proposed. The simulated peak loads and absorbed energy of the tubes are also in good agreement with experimental results.

Laminar-turbulent transition tripped by step on transonic compressor profile

NASA Astrophysics Data System (ADS)

Flaszynski, Pawel; Doerffer, Piotr; Szwaba, Ryszard; Piotrowicz, Michal; Kaczynski, Piotr

2018-02-01

The shock wave boundary layer interaction on the suction side of transonic compressor blade is one of the main objectives of TFAST project (Transition Location Effect on Shock Wave Boundary Layer Interaction). The experimental and numerical results for the flow structure investigations are shown for the flow conditions as the existing ones on the suction side of the compressor profile. The two cases are investigated: without and with boundary layer tripping device. In the first case, boundary layer is laminar up to the shock wave, while in the second case the boundary layer is tripped by the step. Numerical results carried out by means of Fine/Turbo Numeca with Explicit Algebraic Reynolds Stress Model including transition modeling are compared with schlieren, Temperature Sensitive Paint and wake measurements. Boundary layer transition location is detected by Temperature Sensitive Paint.

Simulations of Fluvial Landscapes

NASA Astrophysics Data System (ADS)

Cattan, D.; Birnir, B.

2013-12-01

The Smith-Bretherton-Birnir (SBB) model for fluvial landsurfaces consists of a pair of partial differential equations, one governing water flow and one governing the sediment flow. Numerical solutions of these equations have been shown to provide realistic models in the evolution of fluvial landscapes. Further analysis of these equations shows that they possess scaling laws (Hack's Law) that are known to exist in nature. However, the simulations are highly dependent on the numerical methods used; with implicit methods exhibiting the correct scaling laws, but the explicit methods fail to do so. These equations, and the resulting models, help to bridge the gap between the deterministic and the stochastic theories of landscape evolution. Slight modifications of the SBB equations make the results of the model more realistic. By modifying the sediment flow equation, the model obtains more pronounced meandering rivers. Typical landsurface with rivers.

Numerical methods for the stochastic Landau-Lifshitz Navier-Stokes equations.

PubMed

Bell, John B; Garcia, Alejandro L; Williams, Sarah A

2007-07-01

The Landau-Lifshitz Navier-Stokes (LLNS) equations incorporate thermal fluctuations into macroscopic hydrodynamics by using stochastic fluxes. This paper examines explicit Eulerian discretizations of the full LLNS equations. Several computational fluid dynamics approaches are considered (including MacCormack's two-step Lax-Wendroff scheme and the piecewise parabolic method) and are found to give good results for the variance of momentum fluctuations. However, neither of these schemes accurately reproduces the fluctuations in energy or density. We introduce a conservative centered scheme with a third-order Runge-Kutta temporal integrator that does accurately produce fluctuations in density, energy, and momentum. A variety of numerical tests, including the random walk of a standing shock wave, are considered and results from the stochastic LLNS solver are compared with theory, when available, and with molecular simulations using a direct simulation Monte Carlo algorithm.

Lattice Boltzmann simulations of immiscible displacement process with large viscosity ratios

NASA Astrophysics Data System (ADS)

Rao, Parthib; Schaefer, Laura

2017-11-01

Immiscible displacement is a key physical mechanism involved in enhanced oil recovery and carbon sequestration processes. This multiphase flow phenomenon involves a complex interplay of viscous, capillary, inertial and wettability effects. The lattice Boltzmann (LB) method is an accurate and efficient technique for modeling and simulating multiphase/multicomponent flows especially in complex flow configurations and media. In this presentation we present numerical simulation results of displacement process in thin long channels. The results are based on a new psuedo-potential multicomponent LB model with multiple relaxation time collision (MRT) model and explicit forcing scheme. We demonstrate that the proposed model is capable of accurately simulating the displacement process involving fluids with a wider range of viscosity ratios (>100) and which also leads to viscosity-independent interfacial tension and reduction of some important numerical artifacts.

Full numerical simulation of coflowing, axisymmetric jet diffusion flames

NASA Technical Reports Server (NTRS)

Mahalingam, S.; Cantwell, B. J.; Ferziger, J. H.

1990-01-01

The near field of a non-premixed flame in a low speed, coflowing axisymmetric jet is investigated numerically using full simulation. The time-dependent governing equations are solved by a second-order, explicit finite difference scheme and a single-step, finite rate model is used to represent the chemistry. Steady laminar flame results show the correct dependence of flame height on Peclet number and reaction zone thickness on Damkoehler number. Forced simulations reveal a large difference in the instantaneous structure of scalar dissipation fields between nonbuoyant and buoyant cases. In the former, the scalar dissipation marks intense reaction zones, supporting the flamelet concept; however, results suggest that flamelet modeling assumptions need to be reexamined. In the latter, this correspondence breaks down, suggesting that modifications to the flamelet modeling approach are needed in buoyant turbulent diffusion flames.

Numerical Studies of Impurities in Fusion Plasmas

DOE R&D Accomplishments Database

Hulse, R. A.

1982-09-01

The coupled partial differential equations used to describe the behavior of impurity ions in magnetically confined controlled fusion plasmas require numerical solution for cases of practical interest. Computer codes developed for impurity modeling at the Princeton Plasma Physics Laboratory are used as examples of the types of codes employed for this purpose. These codes solve for the impurity ionization state densities and associated radiation rates using atomic physics appropriate for these low-density, high-temperature plasmas. The simpler codes solve local equations in zero spatial dimensions while more complex cases require codes which explicitly include transport of the impurity ions simultaneously with the atomic processes of ionization and recombination. Typical applications are discussed and computational results are presented for selected cases of interest.

Benchmark calculations of nonconservative charged-particle swarms in dc electric and magnetic fields crossed at arbitrary angles.

PubMed

Dujko, S; White, R D; Petrović, Z Lj; Robson, R E

2010-04-01

A multiterm solution of the Boltzmann equation has been developed and used to calculate transport coefficients of charged-particle swarms in gases under the influence of electric and magnetic fields crossed at arbitrary angles when nonconservative collisions are present. The hierarchy resulting from a spherical-harmonic decomposition of the Boltzmann equation in the hydrodynamic regime is solved numerically by representing the speed dependence of the phase-space distribution function in terms of an expansion in Sonine polynomials about a Maxwellian velocity distribution at an internally determined temperature. Results are given for electron swarms in certain collisional models for ionization and attachment over a range of angles between the fields and field strengths. The implicit and explicit effects of ionization and attachment on the electron-transport coefficients are considered using physical arguments. It is found that the difference between the two sets of transport coefficients, bulk and flux, resulting from the explicit effects of nonconservative collisions, can be controlled either by the variation in the magnetic field strengths or by the angles between the fields. In addition, it is shown that the phenomena of ionization cooling and/or attachment cooling/heating previously reported for dc electric fields carry over directly to the crossed electric and magnetic fields. The results of the Boltzmann equation analysis are compared with those obtained by a Monte Carlo simulation technique. The comparison confirms the theoretical basis and numerical integrity of the moment method for solving the Boltzmann equation and gives a set of well-established data that can be used to test future codes and plasma models.

Optimal implicit 2-D finite differences to model wave propagation in poroelastic media

NASA Astrophysics Data System (ADS)

Itzá, Reymundo; Iturrarán-Viveros, Ursula; Parra, Jorge O.

2016-08-01

Numerical modeling of seismic waves in heterogeneous porous reservoir rocks is an important tool for the interpretation of seismic surveys in reservoir engineering. We apply globally optimal implicit staggered-grid finite differences (FD) to model 2-D wave propagation in heterogeneous poroelastic media at a low-frequency range (<10 kHz). We validate the numerical solution by comparing it to an analytical-transient solution obtaining clear seismic wavefields including fast P and slow P and S waves (for a porous media saturated with fluid). The numerical dispersion and stability conditions are derived using von Neumann analysis, showing that over a wide range of porous materials the Courant condition governs the stability and this optimal implicit scheme improves the stability of explicit schemes. High-order explicit FD can be replaced by some lower order optimal implicit FD so computational cost will not be as expensive while maintaining the accuracy. Here, we compute weights for the optimal implicit FD scheme to attain an accuracy of γ = 10-8. The implicit spatial differentiation involves solving tridiagonal linear systems of equations through Thomas' algorithm.

Modified symplectic schemes with nearly-analytic discrete operators for acoustic wave simulations

NASA Astrophysics Data System (ADS)

Liu, Shaolin; Yang, Dinghui; Lang, Chao; Wang, Wenshuai; Pan, Zhide

2017-04-01

Using a structure-preserving algorithm significantly increases the computational efficiency of solving wave equations. However, only a few explicit symplectic schemes are available in the literature, and the capabilities of these symplectic schemes have not been sufficiently exploited. Here, we propose a modified strategy to construct explicit symplectic schemes for time advance. The acoustic wave equation is transformed into a Hamiltonian system. The classical symplectic partitioned Runge-Kutta (PRK) method is used for the temporal discretization. Additional spatial differential terms are added to the PRK schemes to form the modified symplectic methods and then two modified time-advancing symplectic methods with all of positive symplectic coefficients are then constructed. The spatial differential operators are approximated by nearly-analytic discrete (NAD) operators, and we call the fully discretized scheme modified symplectic nearly analytic discrete (MSNAD) method. Theoretical analyses show that the MSNAD methods exhibit less numerical dispersion and higher stability limits than conventional methods. Three numerical experiments are conducted to verify the advantages of the MSNAD methods, such as their numerical accuracy, computational cost, stability, and long-term calculation capability.

Alternating direction implicit methods for parabolic equations with a mixed derivative

NASA Technical Reports Server (NTRS)

Beam, R. M.; Warming, R. F.

1980-01-01

Alternating direction implicit (ADI) schemes for two-dimensional parabolic equations with a mixed derivative are constructed by using the class of all A(0)-stable linear two-step methods in conjunction with the method of approximate factorization. The mixed derivative is treated with an explicit two-step method which is compatible with an implicit A(0)-stable method. The parameter space for which the resulting ADI schemes are second-order accurate and unconditionally stable is determined. Some numerical examples are given.

Alternating direction implicit methods for parabolic equations with a mixed derivative

NASA Technical Reports Server (NTRS)

Beam, R. M.; Warming, R. F.

1979-01-01

Alternating direction implicit (ADI) schemes for two-dimensional parabolic equations with a mixed derivative are constructed by using the class of all A sub 0-stable linear two-step methods in conjunction with the method of approximation factorization. The mixed derivative is treated with an explicit two-step method which is compatible with an implicit A sub 0-stable method. The parameter space for which the resulting ADI schemes are second order accurate and unconditionally stable is determined. Some numerical examples are given.

Bounds on quantum confinement effects in metal nanoparticles

NASA Astrophysics Data System (ADS)

Blackman, G. Neal; Genov, Dentcho A.

2018-03-01

Quantum size effects on the permittivity of metal nanoparticles are investigated using the quantum box model. Explicit upper and lower bounds are derived for the permittivity and relaxation rates due to quantum confinement effects. These bounds are verified numerically, and the size dependence and frequency dependence of the empirical Drude size parameter is extracted from the model. Results suggest that the common practice of empirically modifying the dielectric function can lead to inaccurate predictions for highly uniform distributions of finite-sized particles.

High-order centered difference methods with sharp shock resolution

NASA Technical Reports Server (NTRS)

Gustafsson, Bertil; Olsson, Pelle

1994-01-01

In this paper we consider high-order centered finite difference approximations of hyperbolic conservation laws. We propose different ways of adding artificial viscosity to obtain sharp shock resolution. For the Riemann problem we give simple explicit formulas for obtaining stationary one and two-point shocks. This can be done for any order of accuracy. It is shown that the addition of artificial viscosity is equivalent to ensuring the Lax k-shock condition. We also show numerical experiments that verify the theoretical results.

Necessary and sufficient conditions for discrete wavelet frames in CN

NASA Astrophysics Data System (ADS)

Deepshikha; Vashisht, Lalit K.

2017-07-01

We present necessary and sufficient conditions with explicit frame bounds for a discrete wavelet system of the form {DaTk ϕ } a ∈ U(N) , k ∈IN to be a frame for the unitary space CN. It is shown that the canonical dual of a discrete wavelet frame for CN has the same structure. This is not true (well known) for canonical dual of a wavelet frame for L2(R) . Several numerical examples are given to illustrate the results.

Triple collinear emissions in parton showers

DOE PAGES

Hoche, Stefan; Prestel, Stefan

2017-10-17

A framework to include triple collinear splitting functions into parton showers is presented, and the implementation of flavor-changing next-to-leading-order (NLO) splitting kernels is discussed as a first application. The correspondence between the Monte Carlo integration and the analytic computation of NLO DGLAP evolution kernels is made explicit for both timelike and spacelike parton evolution. Finally, numerical simulation results are obtained with two independent implementations of the new algorithm, using the two independent event generation frameworks PYTHIA and SHERPA.

Experimental and numerical study of drill bit drop tests on Kuru granite

NASA Astrophysics Data System (ADS)

Fourmeau, Marion; Kane, Alexandre; Hokka, Mikko

2017-01-01

This paper presents an experimental and numerical study of Kuru grey granite impacted with a seven-buttons drill bit mounted on an instrumented drop test machine. The force versus displacement curves during the impact, so-called bit-rock interaction (BRI) curves, were obtained using strain gauge measurements for two levels of impact energy. Moreover, the volume of removed rock after each drop test was evaluated by stereo-lithography (three-dimensional surface reconstruction). A modified version of the Holmquist-Johnson-Cook (MHJC) material model was calibrated using Kuru granite test results available from the literature. Numerical simulations of the single drop tests were carried out using the MHJC model available in the LS-DYNA explicit finite-element solver. The influence of the impact energy and additional confining pressure on the BRI curves and the volume of the removed rock is discussed. In addition, the influence of the rock surface shape before impact was evaluated using two different mesh geometries: a flat surface and a hyperbolic surface. The experimental and numerical results are compared and discussed in terms of drilling efficiency through the mechanical specific energy. This article is part of the themed issue 'Experimental testing and modelling of brittle materials at high strain rates'.

A multi-dimensional nonlinearly implicit, electromagnetic Vlasov-Darwin particle-in-cell (PIC) algorithm

NASA Astrophysics Data System (ADS)

Chen, Guangye; Chacón, Luis; CoCoMans Team

2014-10-01

For decades, the Vlasov-Darwin model has been recognized to be attractive for PIC simulations (to avoid radiative noise issues) in non-radiative electromagnetic regimes. However, the Darwin model results in elliptic field equations that renders explicit time integration unconditionally unstable. Improving on linearly implicit schemes, fully implicit PIC algorithms for both electrostatic and electromagnetic regimes, with exact discrete energy and charge conservation properties, have been recently developed in 1D. This study builds on these recent algorithms to develop an implicit, orbit-averaged, time-space-centered finite difference scheme for the particle-field equations in multiple dimensions. The algorithm conserves energy, charge, and canonical-momentum exactly, even with grid packing. A simple fluid preconditioner allows efficient use of large timesteps, O (√{mi/me}c/veT) larger than the explicit CFL. We demonstrate the accuracy and efficiency properties of the of the algorithm with various numerical experiments in 2D3V.

Three-dimensional polarization algebra.

PubMed

R Sheppard, Colin J; Castello, Marco; Diaspro, Alberto

2016-10-01

If light is focused or collected with a high numerical aperture lens, as may occur in imaging and optical encryption applications, polarization should be considered in three dimensions (3D). The matrix algebra of polarization behavior in 3D is discussed. It is useful to convert between the Mueller matrix and two different Hermitian matrices, representing an optical material or system, which are in the literature. Explicit transformation matrices for converting the column vector form of these different matrices are extended to the 3D case, where they are large (81×81) but can be generated using simple rules. It is found that there is some advantage in using a generalization of the Chandrasekhar phase matrix treatment, rather than that based on Gell-Mann matrices, as the resultant matrices are of simpler form and reduce to the two-dimensional case more easily. Explicit expressions are given for 3D complex field components in terms of Chandrasekhar-Stokes parameters.

Shock-induced solitary waves in granular crystals.

PubMed

Hasan, M Arif; Nemat-Nasser, Sia

2018-02-01

Solitary waves (SWs) are generated in monoatomic (homogeneous) lightly contacting spherical granules by an applied input force of any time-variation and intensity. We consider finite duration shock loads on one-dimensional arrays of granules and focus on the transition regime that leads to the formation of SWs. Based on geometrical and material properties of the granules and the properties of the input shock, we provide explicit analytic expressions to calculate the peak value of the compressive contact force at each contact point in the transition regime that precedes the formation of a primary solitary wave. We also provide explicit expressions to estimate the number of granules involved in the transition regime and show its dependence on the characteristics of the input shock and material/geometrical properties of the interacting granules. Finally, we assess the accuracy of our theoretical results by comparing them with those obtained through numerical integration of the equations of motion.

Identification of internal properties of fibres and micro-swimmers

NASA Astrophysics Data System (ADS)

Plouraboué, Franck; Thiam, E. Ibrahima; Delmotte, Blaise; Climent, Eric

2017-01-01

In this paper, we address the identifiability of constitutive parameters of passive or active micro-swimmers. We first present a general framework for describing fibres or micro-swimmers using a bead-model description. Using a kinematic constraint formulation to describe fibres, flagellum or cilia, we find explicit linear relationship between elastic constitutive parameters and generalized velocities from computing contact forces. This linear formulation then permits one to address explicitly identifiability conditions and solve for parameter identification. We show that both active forcing and passive parameters are both identifiable independently but not simultaneously. We also provide unbiased estimators for generalized elastic parameters in the presence of Langevin-like forcing with Gaussian noise using a Bayesian approach. These theoretical results are illustrated in various configurations showing the efficiency of the proposed approach for direct parameter identification. The convergence of the proposed estimators is successfully tested numerically.

The kinetics of heterogeneous nucleation and growth: an approach based on a grain explicit model

NASA Astrophysics Data System (ADS)

Rouet-Leduc, B.; Maillet, J.-B.; Denoual, C.

2014-04-01

A model for phase transitions initiated on grain boundaries is proposed and tested against numerical simulations: this approach, based on a grain explicit model, allows us to consider the granular structure, resulting in accurate predictions for a wide span of nucleation processes. Comparisons are made with classical models of homogeneous (JMAK: Johnson and Mehl 1939 Trans. Am. Inst. Min. Eng. 135 416; Avrami 1939 J. Chem. Phys. 7 1103; Kolmogorov 1937 Bull. Acad. Sci. USSR, Mat. Ser. 1 335) as well as heterogeneous (Cahn 1996 Thermodynamics and Kinetics of Phase Transformations Im et al (Pittsburgh: Materials Research Society)) nucleation. A transition scale based on material properties is proposed, allowing us to discriminate between random and site-saturated regimes. Finally, we discuss the relationship between an Avrami-type exponent and the transition regime, establishing conditions for its extraction from experiments.

Thermo-hydroforming of a fiber-reinforced thermoplastic composites considering fiber orientations

NASA Astrophysics Data System (ADS)

Ahn, Hyunchul; Kuuttila, Nicholas Eric; Pourboghrat, Farhang

2018-05-01

The Thermoplastic woven composites were formed using a composite thermal hydroforming process, utilizing heated and pressurized fluid, similar to sheet metal forming. This study focuses on the modification of 300-ton pressure formation and predicts its behavior. Spectra Shield SR-3136 is used in this study and material properties are measured by experiments. The behavior of fiber-reinforced thermoplastic polymer composites (FRTP) was modeled using the Preferred Fiber Orientation (PFO) model and validated by comparing numerical analysis with experimental results. The thermo-hydroforming process has shown good results in the ability to form deep drawn parts with reduced wrinkles. Numerical analysis was performed using the PFO model and implemented as commercial finite element software ABAQUS / Explicit. The user subroutine (VUMAT) was used for the material properties of the thermoplastic composite layer. This model is suitable for working with multiple layers of composite laminates. Model parameters have been updated to work with cohesive zone model to calculate the interfacial properties between each composite layer. The results of the numerical modeling showed a good correlation with the molding experiment on the forming shape. Numerical results were also compared with experimental results on punch force-displacement curves for deformed geometry and forming processes of the composite layer. Overall, the shape of the deformed FRTP, including the distribution of wrinkles, was accurately predicted as shown in this study.

Four-level conservative finite-difference schemes for Boussinesq paradigm equation

NASA Astrophysics Data System (ADS)

Kolkovska, N.

2013-10-01

In this paper a two-parametric family of four level conservative finite difference schemes is constructed for the multidimensional Boussinesq paradigm equation. The schemes are explicit in the sense that no inner iterations are needed for evaluation of the numerical solution. The preservation of the discrete energy with this method is proved. The schemes have been numerically tested on one soliton propagation model and two solitons interaction model. The numerical experiments demonstrate that the proposed family of schemes has second order of convergence in space and time steps in the discrete maximal norm.

Numerical approach to optimal portfolio in a power utility regime-switching model

NASA Astrophysics Data System (ADS)

Gyulov, Tihomir B.; Koleva, Miglena N.; Vulkov, Lubin G.

2017-12-01

We consider a system of weakly coupled degenerate semi-linear parabolic equations of optimal portfolio in a regime-switching with power utility function, derived by A.R. Valdez and T. Vargiolu [14]. First, we discuss some basic properties of the solution of this system. Then, we develop and analyze implicit-explicit, flux limited finite difference schemes for the differential problem. Numerical experiments are discussed.

Cognitive conflict without explicit conflict monitoring in a dynamical agent.

PubMed

Ward, Robert; Ward, Ronnie

2006-11-01

We examine mechanisms for resolving cognitive conflict in an embodied, situated, and dynamic agent, developed through an evolutionary learning process. The agent was required to solve problems of response conflict in a dual-target "catching" task, focusing response on one of the targets while ignoring the other. Conflict in the agent was revealed at the behavioral level in terms of increased latencies to the second target. This behavioral interference was correlated to peak violations of the network's stable state equation. At the level of the agent's neural network, peak violations were also correlated to periods of disagreement in source inputs to the agent's motor effectors. Despite observing conflict at these numerous levels, we did not find any explicit conflict monitoring mechanisms within the agent. We instead found evidence of a distributed conflict management system, characterized by competitive sources within the network. In contrast to the conflict monitoring hypothesis [Botvinick, M. M., Braver, T. S., Barch, D. M., Carter, C. S., & Cohen, J. D. (2001). Conflict monitoring and cognitive control. Psychological Review, 108(3), 624-652], this agent demonstrates that resolution of cognitive conflict does not require explicit conflict monitoring. We consider the implications of our results for the conflict monitoring hypothesis.

Contact-aware simulations of particulate Stokesian suspensions

NASA Astrophysics Data System (ADS)

Lu, Libin; Rahimian, Abtin; Zorin, Denis

2017-10-01

We present an efficient, accurate, and robust method for simulation of dense suspensions of deformable and rigid particles immersed in Stokesian fluid in two dimensions. We use a well-established boundary integral formulation for the problem as the foundation of our approach. This type of formulation, with a high-order spatial discretization and an implicit and adaptive time discretization, have been shown to be able to handle complex interactions between particles with high accuracy. Yet, for dense suspensions, very small time-steps or expensive implicit solves as well as a large number of discretization points are required to avoid non-physical contact and intersections between particles, leading to infinite forces and numerical instability. Our method maintains the accuracy of previous methods at a significantly lower cost for dense suspensions. The key idea is to ensure interference-free configuration by introducing explicit contact constraints into the system. While such constraints are unnecessary in the formulation, in the discrete form of the problem, they make it possible to eliminate catastrophic loss of accuracy by preventing contact explicitly. Introducing contact constraints results in a significant increase in stable time-step size for explicit time-stepping, and a reduction in the number of points adequate for stability.

GPU-accelerated simulations of isolated black holes

NASA Astrophysics Data System (ADS)

Lewis, Adam G. M.; Pfeiffer, Harald P.

2018-05-01

We present a port of the numerical relativity code SpEC which is capable of running on NVIDIA GPUs. Since this code must be maintained in parallel with SpEC itself, a primary design consideration is to perform as few explicit code changes as possible. We therefore rely on a hierarchy of automated porting strategies. At the highest level we use TLoops, a C++ library of our design, to automatically emit CUDA code equivalent to tensorial expressions written into C++ source using a syntax similar to analytic calculation. Next, we trace out and cache explicit matrix representations of the numerous linear transformations in the SpEC code, which allows these to be performed on the GPU using pre-existing matrix-multiplication libraries. We port the few remaining important modules by hand. In this paper we detail the specifics of our port, and present benchmarks of it simulating isolated black hole spacetimes on several generations of NVIDIA GPU.

An approximate inverse scattering technique for reconstructing blockage profiles in water pipelines using acoustic transients.

PubMed

Jing, Liwen; Li, Zhao; Wang, Wenjie; Dubey, Amartansh; Lee, Pedro; Meniconi, Silvia; Brunone, Bruno; Murch, Ross D

2018-05-01

An approximate inverse scattering technique is proposed for reconstructing cross-sectional area variation along water pipelines to deduce the size and position of blockages. The technique allows the reconstructed blockage profile to be written explicitly in terms of the measured acoustic reflectivity. It is based upon the Born approximation and provides good accuracy, low computational complexity, and insight into the reconstruction process. Numerical simulations and experimental results are provided for long pipelines with mild and severe blockages of different lengths. Good agreement is found between the inverse result and the actual pipe condition for mild blockages.

Saturation of the magnetorotational instability in the unstratified shearing box with zero net flux: convergence in taller boxes

NASA Astrophysics Data System (ADS)

Shi, Ji-Ming; Stone, James M.; Huang, Chelsea X.

2016-03-01

Previous studies of the non-linear regime of the magnetorotational instability in one particular type of shearing box model - unstratified with no net magnetic flux - find that without explicit dissipation (viscosity and resistivity) the saturation amplitude decreases with increasing numerical resolution. We show that this result is strongly dependent on the vertical aspect ratio of the computational domain Lz/Lx. When Lz/Lx ≲ 1, we recover previous results. However, when the vertical domain is extended Lz/Lx ≳ 2.5, we find the saturation level of the stress is greatly increased (giving a ratio of stress to pressure α ≳ 0.1), and moreover the results are independent of numerical resolution. Consistent with previous results, we find that saturation of the magnetorotational (MRI) in this regime is controlled by a cyclic dynamo which generates patches of strong toroidal field that switches sign on scales of Lx in the vertical direction. We speculate that when Lz/Lx ≲ 1, the dynamo is inhibited by the small size of the vertical domain, leading to the puzzling dependence of saturation amplitude on resolution. We show that previous toy models developed to explain the MRI dynamo are consistent with our results, and that the cyclic pattern of toroidal fields observed in stratified shearing box simulations (leading to the so-called butterfly diagram) may also be related. In tall boxes the saturation amplitude is insensitive to whether or not explicit dissipation is included in the calculations, at least for large magnetic Reynolds and Prandtl number. Finally, we show MRI turbulence in tall domains has a smaller critical Pmc, and an extended lifetime compared to Lz/Lx ≲ 1 boxes.

Free vibration of functionally graded beams and frameworks using the dynamic stiffness method

NASA Astrophysics Data System (ADS)

Banerjee, J. R.; Ananthapuvirajah, A.

2018-05-01

The free vibration analysis of functionally graded beams (FGBs) and frameworks containing FGBs is carried out by applying the dynamic stiffness method and deriving the elements of the dynamic stiffness matrix in explicit algebraic form. The usually adopted rule that the material properties of the FGB vary continuously through the thickness according to a power law forms the fundamental basis of the governing differential equations of motion in free vibration. The differential equations are solved in closed analytical form when the free vibratory motion is harmonic. The dynamic stiffness matrix is then formulated by relating the amplitudes of forces to those of the displacements at the two ends of the beam. Next, the explicit algebraic expressions for the dynamic stiffness elements are derived with the help of symbolic computation. Finally the Wittrick-Williams algorithm is applied as solution technique to solve the free vibration problems of FGBs with uniform cross-section, stepped FGBs and frameworks consisting of FGBs. Some numerical results are validated against published results, but in the absence of published results for frameworks containing FGBs, consistency checks on the reliability of results are performed. The paper closes with discussion of results and conclusions.

The CFS-PML in numerical simulation of ATEM

NASA Astrophysics Data System (ADS)

Zhao, Xuejiao; Ji, Yanju; Qiu, Shuo; Guan, Shanshan; Wu, Yanqi

2017-01-01

In the simulation of airborne transient electromagnetic method (ATEM) in time-domain, the truncated boundary reflection can bring a big error to the results. The complex frequency shifted perfectly matched layer (CFS-PML) absorbing boundary condition has been proved to have a better absorption of low frequency incident wave and can reduce the late reflection greatly. In this paper, we apply the CFS-PML to three-dimensional numerical simulation of ATEM in time-domain to achieve a high precision .The expression of divergence equation in CFS-PML is confirmed and its explicit iteration format based on the finite difference method and the recursive convolution technique is deduced. Finally, we use the uniformity half space model and the anomalous model to test the validity of this method. Results show that the CFS-PML can reduce the average relative error to 2.87% and increase the accuracy of the anomaly recognition.

Solution procedure of dynamical contact problems with friction

NASA Astrophysics Data System (ADS)

Abdelhakim, Lotfi

2017-07-01

Dynamical contact is one of the common research topics because of its wide applications in the engineering field. The main goal of this work is to develop a time-stepping algorithm for dynamic contact problems. We propose a finite element approach for elastodynamics contact problems [1]. Sticking, sliding and frictional contact can be taken into account. Lagrange multipliers are used to enforce non-penetration condition. For the time discretization, we propose a scheme equivalent to the explicit Newmark scheme. Each time step requires solving a nonlinear problem similar to a static friction problem. The nonlinearity of the system of equation needs an iterative solution procedure based on Uzawa's algorithm [2][3]. The applicability of the algorithm is illustrated by selected sample numerical solutions to static and dynamic contact problems. Results obtained with the model have been compared and verified with results from an independent numerical method.

Shock wave boundary layer interaction on suction side of compressor profile in single passage test section

NASA Astrophysics Data System (ADS)

Flaszynski, Pawel; Doerffer, Piotr; Szwaba, Ryszard; Kaczynski, Piotr; Piotrowicz, Michal

2015-11-01

The shock wave boundary layer interaction on the suction side of transonic compressor blade is one of the main objectives of TFAST project (Transition Location Effect on Shock Wave Boundary Layer Interaction). In order to investigate the flow structure on the suction side of a profile, a design of a generic test section in linear transonic wind tunnel was proposed. The experimental and numerical results for the flow structure investigations are shown for the flow conditions as the existing ones on the suction side of the compressor profile. Near the sidewalls the suction slots are applied for the corner flow structure control. It allows to control the Axial Velocity Density Ratio (AVDR), important parameter for compressor cascade investigations. Numerical results for Explicit Algebraic Reynolds Stress Model with transition modeling are compared with oil flow visualization, schlieren and Pressure Sensitive Paint. Boundary layer transition location is detected by Temperature Sensitive Paint.

A motion artefact study and locally deforming objects in computerized tomography

NASA Astrophysics Data System (ADS)

Hahn, Bernadette N.

2017-11-01

Movements of the object during the data collection in computerized tomography can introduce motion artefacts in the reconstructed image. They can be reduced by employing information about the dynamic behaviour within the reconstruction step. However, inaccuracies concerning the movement are inevitable in practice. In this article, we give an explicit characterization of what is visible in an image obtained by a reconstruction algorithm with incorrect motion information. Then, we use this result to study in detail the situation of locally deforming objects, i.e. individual parts of the object have a different dynamic behaviour. In this context, we prove that additional artefacts arise due to the global nature of the Radon transform, even if the motion is exactly known. Based on our analysis, we propose a numerical scheme to reduce these artefacts in the reconstructed image. All our results are illustrated by numerical examples.

Hydro-geophysical observations integration in numerical model: case study in Mediterranean karstic unsaturated zone (Larzac, france)

NASA Astrophysics Data System (ADS)

Champollion, Cédric; Fores, Benjamin; Le Moigne, Nicolas; Chéry, Jean

2016-04-01

Karstic hydro-systems are highly non-linear and heterogeneous but one of the main water resource in the Mediterranean area. Neither local measurements in boreholes or analysis at the spring can take into account the variability of the water storage. Since a few years, ground-based geophysical measurements (such as gravity, electrical resistivity or seismological data) allows following water storage in heterogeneous hydrosystems at an intermediate scale between boreholes and basin. Behind classical rigorous monitoring, the integration of geophysical data in hydrological numerical models in needed for both processes interpretation and quantification. Since a few years, a karstic geophysical observatory (GEK: Géodésie de l'Environnement Karstique, OSU OREME, SNO H+) has been setup in the Mediterranean area in the south of France. The observatory is surrounding more than 250m karstified dolomite, with an unsaturated zone of ~150m thickness. At the observatory water level in boreholes, evapotranspiration and rainfall are classical hydro-meteorological observations completed by continuous gravity, resistivity and seismological measurements. The main objective of the study is the modelling of the whole observation dataset by explicit unsaturated numerical model in one dimension. Hydrus software is used for the explicit modelling of the water storage and transfer and links the different observations (geophysics, water level, evapotranspiration) with the water saturation. Unknown hydrological parameters (permeability, porosity) are retrieved from stochastic inversions. The scale of investigation of the different observations are discussed thank to the modelling results. A sensibility study of the measurements against the model is done and key hydro-geological processes of the site are presented.

A numerical scheme based on radial basis function finite difference (RBF-FD) technique for solving the high-dimensional nonlinear Schrödinger equations using an explicit time discretization: Runge-Kutta method

NASA Astrophysics Data System (ADS)

Dehghan, Mehdi; Mohammadi, Vahid

2017-08-01

In this research, we investigate the numerical solution of nonlinear Schrödinger equations in two and three dimensions. The numerical meshless method which will be used here is RBF-FD technique. The main advantage of this method is the approximation of the required derivatives based on finite difference technique at each local-support domain as Ωi. At each Ωi, we require to solve a small linear system of algebraic equations with a conditionally positive definite matrix of order 1 (interpolation matrix). This scheme is efficient and its computational cost is same as the moving least squares (MLS) approximation. A challengeable issue is choosing suitable shape parameter for interpolation matrix in this way. In order to overcome this matter, an algorithm which was established by Sarra (2012), will be applied. This algorithm computes the condition number of the local interpolation matrix using the singular value decomposition (SVD) for obtaining the smallest and largest singular values of that matrix. Moreover, an explicit method based on Runge-Kutta formula of fourth-order accuracy will be applied for approximating the time variable. It also decreases the computational costs at each time step since we will not solve a nonlinear system. On the other hand, to compare RBF-FD method with another meshless technique, the moving kriging least squares (MKLS) approximation is considered for the studied model. Our results demonstrate the ability of the present approach for solving the applicable model which is investigated in the current research work.

The molecular kink paradigm for rubber elasticity: Numerical simulations of explicit polyisoprene networks at low to moderate tensile strains

NASA Astrophysics Data System (ADS)

Hanson, David E.

2011-08-01

Based on recent molecular dynamics and ab initio simulations of small isoprene molecules, we propose a new ansatz for rubber elasticity. We envision a network chain as a series of independent molecular kinks, each comprised of a small number of backbone units, and the strain as being imposed along the contour of the chain. We treat chain extension in three distinct force regimes: (Ia) near zero strain, where we assume that the chain is extended within a well defined tube, with all of the kinks participating simultaneously as entropic elastic springs, (II) when the chain becomes sensibly straight, giving rise to a purely enthalpic stretching force (until bond rupture occurs) and, (Ib) a linear entropic regime, between regimes Ia and II, in which a force limit is imposed by tube deformation. In this intermediate regime, the molecular kinks are assumed to be gradually straightened until the chain becomes a series of straight segments between entanglements. We assume that there exists a tube deformation tension limit that is inversely proportional to the chain path tortuosity. Here we report the results of numerical simulations of explicit three-dimensional, periodic, polyisoprene networks, using these extension-only force models. At low strain, crosslink nodes are moved affinely, up to an arbitrary node force limit. Above this limit, non-affine motion of the nodes is allowed to relax unbalanced chain forces. Our simulation results are in good agreement with tensile stress vs. strain experiments.

The molecular kink paradigm for rubber elasticity: numerical simulations of explicit polyisoprene networks at low to moderate tensile strains.

PubMed

Hanson, David E

2011-08-07

Based on recent molecular dynamics and ab initio simulations of small isoprene molecules, we propose a new ansatz for rubber elasticity. We envision a network chain as a series of independent molecular kinks, each comprised of a small number of backbone units, and the strain as being imposed along the contour of the chain. We treat chain extension in three distinct force regimes: (Ia) near zero strain, where we assume that the chain is extended within a well defined tube, with all of the kinks participating simultaneously as entropic elastic springs, (II) when the chain becomes sensibly straight, giving rise to a purely enthalpic stretching force (until bond rupture occurs) and, (Ib) a linear entropic regime, between regimes Ia and II, in which a force limit is imposed by tube deformation. In this intermediate regime, the molecular kinks are assumed to be gradually straightened until the chain becomes a series of straight segments between entanglements. We assume that there exists a tube deformation tension limit that is inversely proportional to the chain path tortuosity. Here we report the results of numerical simulations of explicit three-dimensional, periodic, polyisoprene networks, using these extension-only force models. At low strain, crosslink nodes are moved affinely, up to an arbitrary node force limit. Above this limit, non-affine motion of the nodes is allowed to relax unbalanced chain forces. Our simulation results are in good agreement with tensile stress vs. strain experiments.

A space-time lower-upper symmetric Gauss-Seidel scheme for the time-spectral method

NASA Astrophysics Data System (ADS)

Zhan, Lei; Xiong, Juntao; Liu, Feng

2016-05-01

The time-spectral method (TSM) offers the advantage of increased order of accuracy compared to methods using finite-difference in time for periodic unsteady flow problems. Explicit Runge-Kutta pseudo-time marching and implicit schemes have been developed to solve iteratively the space-time coupled nonlinear equations resulting from TSM. Convergence of the explicit schemes is slow because of the stringent time-step limit. Many implicit methods have been developed for TSM. Their computational efficiency is, however, still limited in practice because of delayed implicit temporal coupling, multiple iterative loops, costly matrix operations, or lack of strong diagonal dominance of the implicit operator matrix. To overcome these shortcomings, an efficient space-time lower-upper symmetric Gauss-Seidel (ST-LU-SGS) implicit scheme with multigrid acceleration is presented. In this scheme, the implicit temporal coupling term is split as one additional dimension of space in the LU-SGS sweeps. To improve numerical stability for periodic flows with high frequency, a modification to the ST-LU-SGS scheme is proposed. Numerical results show that fast convergence is achieved using large or even infinite Courant-Friedrichs-Lewy (CFL) numbers for unsteady flow problems with moderately high frequency and with the use of moderately high numbers of time intervals. The ST-LU-SGS implicit scheme is also found to work well in calculating periodic flow problems where the frequency is not known a priori and needed to be determined by using a combined Fourier analysis and gradient-based search algorithm.

Numerical simulation of pounding damage to caisson under storm surge

NASA Astrophysics Data System (ADS)

Yu, Chen

2018-06-01

In this paper, a new method for the numerical simulation of structural model is proposed, which is employed to analyze the pounding response of caissons subjected to storm surge loads. According to the new method, the simulation process is divided into two steps. Firstly, the wave propagation caused by storm surge is simulated by the wave-generating tool of Flow-3D, and recording the wave force time history on the caisson. Secondly, a refined 3D finite element model of caisson is established, and the wave force load is applied on the caisson according to the measured data in the first step for further analysis of structural pounding response using the explicit solver of LSDYNA. The whole simulation of pounding response of a caisson caused by "Sha Lijia" typhoon is carried out. The results show that the different wave direction results in the different angle caisson collisions, which will lead to different failure mode of caisson, and when the angle of 60 between wave direction and front/back wall is simulated, the numerical pounding failure mode is consistent with the situation.

Numerical simulation of liquid jet impact on a rigid wall

NASA Astrophysics Data System (ADS)

Aganin, A. A.; Guseva, T. S.

2016-11-01

Basic points of a numerical technique for computing high-speed liquid jet impact on a rigid wall are presented. In the technique the flows of the liquid and the surrounding gas are governed by the equations of gas dynamics in the density, velocity, and pressure, which are integrated by the CIP-CUP method on dynamically adaptive grids without explicitly tracking the gas-liquid interface. The efficiency of the technique is demonstrated by the results of computing the problems of impact of the liquid cone and the liquid wedge on a wall in the mode with the shockwave touching the wall by its edge. Numerical solutions of these problems are compared with the analytical solution of the problem of impact of the plane liquid flow on a wall. Applicability of the technique to the problems of the high-speed liquid jet impact on a wall is illustrated by the results of computing a problem of impact of a cylindrical liquid jet with the hemispherical end on a wall covered by a layer of the same liquid.

Experimental and Numerical Studies on the Formability of Materials in Hot Stamping and Cold Die Quenching Processes

NASA Astrophysics Data System (ADS)

Li, N.; Mohamed, M. S.; Cai, J.; Lin, J.; Balint, D.; Dean, T. A.

2011-05-01

Formability of steel and aluminium alloys in hot stamping and cold die quenching processes is studied in this research. Viscoplastic-damage constitutive equations are developed and determined from experimental data for the prediction of viscoplastic flow and ductility of the materials. The determined unified constitutive equations are then implemented into the commercial Finite Element code Abaqus/Explicit via a user defined subroutine, VUMAT. An FE process simulation model and numerical procedures are established for the modeling of hot stamping processes for a spherical part with a central hole. Different failure modes (failure takes place either near the central hole or in the mid span of the part) are obtained. To validate the simulation results, a test programme is developed, a test die set has been designed and manufactured, and tests have been carried out for the materials with different forming rates. It has been found that very close agreements between experimental and numerical process simulation results are obtained for the ranges of temperatures and forming rates carried out.

Staggered solution procedures for multibody dynamics simulation

NASA Technical Reports Server (NTRS)

Park, K. C.; Chiou, J. C.; Downer, J. D.

1990-01-01

The numerical solution procedure for multibody dynamics (MBD) systems is termed a staggered MBD solution procedure that solves the generalized coordinates in a separate module from that for the constraint force. This requires a reformulation of the constraint conditions so that the constraint forces can also be integrated in time. A major advantage of such a partitioned solution procedure is that additional analysis capabilities such as active controller and design optimization modules can be easily interfaced without embedding them into a monolithic program. After introducing the basic equations of motion for MBD system in the second section, Section 3 briefly reviews some constraint handling techniques and introduces the staggered stabilized technique for the solution of the constraint forces as independent variables. The numerical direct time integration of the equations of motion is described in Section 4. As accurate damping treatment is important for the dynamics of space structures, we have employed the central difference method and the mid-point form of the trapezoidal rule since they engender no numerical damping. This is in contrast to the current practice in dynamic simulations of ground vehicles by employing a set of backward difference formulas. First, the equations of motion are partitioned according to the translational and the rotational coordinates. This sets the stage for an efficient treatment of the rotational motions via the singularity-free Euler parameters. The resulting partitioned equations of motion are then integrated via a two-stage explicit stabilized algorithm for updating both the translational coordinates and angular velocities. Once the angular velocities are obtained, the angular orientations are updated via the mid-point implicit formula employing the Euler parameters. When the two algorithms, namely, the two-stage explicit algorithm for the generalized coordinates and the implicit staggered procedure for the constraint Lagrange multipliers, are brought together in a staggered manner, they constitute a staggered explicit-implicit procedure which is summarized in Section 5. Section 6 presents some example problems and discussions concerning several salient features of the staggered MBD solution procedure are offered in Section 7.

Efficient adaptive pseudo-symplectic numerical integration techniques for Landau-Lifshitz dynamics

NASA Astrophysics Data System (ADS)

d'Aquino, M.; Capuano, F.; Coppola, G.; Serpico, C.; Mayergoyz, I. D.

2018-05-01

Numerical time integration schemes for Landau-Lifshitz magnetization dynamics are considered. Such dynamics preserves the magnetization amplitude and, in the absence of dissipation, also implies the conservation of the free energy. This property is generally lost when time discretization is performed for the numerical solution. In this work, explicit numerical schemes based on Runge-Kutta methods are introduced. The schemes are termed pseudo-symplectic in that they are accurate to order p, but preserve magnetization amplitude and free energy to order q > p. An effective strategy for adaptive time-stepping control is discussed for schemes of this class. Numerical tests against analytical solutions for the simulation of fast precessional dynamics are performed in order to point out the effectiveness of the proposed methods.

Numerical algorithms based on Galerkin methods for the modeling of reactive interfaces in photoelectrochemical (PEC) solar cells

NASA Astrophysics Data System (ADS)

Harmon, Michael; Gamba, Irene M.; Ren, Kui

2016-12-01

This work concerns the numerical solution of a coupled system of self-consistent reaction-drift-diffusion-Poisson equations that describes the macroscopic dynamics of charge transport in photoelectrochemical (PEC) solar cells with reactive semiconductor and electrolyte interfaces. We present three numerical algorithms, mainly based on a mixed finite element and a local discontinuous Galerkin method for spatial discretization, with carefully chosen numerical fluxes, and implicit-explicit time stepping techniques, for solving the time-dependent nonlinear systems of partial differential equations. We perform computational simulations under various model parameters to demonstrate the performance of the proposed numerical algorithms as well as the impact of these parameters on the solution to the model.

An implicit dispersive transport algorithm for the US Geological Survey MOC3D solute-transport model

USGS Publications Warehouse

Kipp, K.L.; Konikow, Leonard F.; Hornberger, G.Z.

1998-01-01

This report documents an extension to the U.S. Geological Survey MOC3D transport model that incorporates an implicit-in-time difference approximation for the dispersive transport equation, including source/sink terms. The original MOC3D transport model (Version 1) uses the method of characteristics to solve the transport equation on the basis of the velocity field. The original MOC3D solution algorithm incorporates particle tracking to represent advective processes and an explicit finite-difference formulation to calculate dispersive fluxes. The new implicit procedure eliminates several stability criteria required for the previous explicit formulation. This allows much larger transport time increments to be used in dispersion-dominated problems. The decoupling of advective and dispersive transport in MOC3D, however, is unchanged. With the implicit extension, the MOC3D model is upgraded to Version 2. A description of the numerical method of the implicit dispersion calculation, the data-input requirements and output options, and the results of simulator testing and evaluation are presented. Version 2 of MOC3D was evaluated for the same set of problems used for verification of Version 1. These test results indicate that the implicit calculation of Version 2 matches the accuracy of Version 1, yet is more efficient than the explicit calculation for transport problems that are characterized by a grid Peclet number less than about 1.0.

The numerology of gender: gendered perceptions of even and odd numbers

PubMed Central

Wilkie, James E. B.; Bodenhausen, Galen V.

2015-01-01

Do numbers have gender? Wilkie and Bodenhausen (2012) examined this issue in a series of experiments on perceived gender. They examined the perceived gender of baby faces and foreign names. Arbitrary numbers presented with these faces and names influenced their perceived gender. Specifically, odd numbers connoted masculinity, while even numbers connoted femininity. In two new studies (total N = 315), we further examined the gendering of numbers. The first study examined explicit ratings of 1-digit numbers. We confirmed that odd numbers seemed masculine while even numbers seemed feminine. Although both men and women showed this pattern, it was more pronounced among women. We also examined whether this pattern holds for automatic as well as deliberated reactions. Results of an Implicit Association Test showed that it did, but only among the women. The implicit and explicit patterns of numerical gender ascription were moderately correlated. The second study examined explicit perceptions of 2-digit numbers. Again, women viewed odd numbers as more masculine and less feminine than even numbers. However, men viewed 2-digit numbers as relatively masculine, regardless of whether they were even or odd. These results indicate that women and men impute gender to numbers in different ways and to different extents. We discuss possible implications for understanding how people relate to and are influenced by numbers in a variety of real-life contexts. PMID:26113839

The numerology of gender: gendered perceptions of even and odd numbers.

PubMed

Wilkie, James E B; Bodenhausen, Galen V

2015-01-01

Do numbers have gender? Wilkie and Bodenhausen (2012) examined this issue in a series of experiments on perceived gender. They examined the perceived gender of baby faces and foreign names. Arbitrary numbers presented with these faces and names influenced their perceived gender. Specifically, odd numbers connoted masculinity, while even numbers connoted femininity. In two new studies (total N = 315), we further examined the gendering of numbers. The first study examined explicit ratings of 1-digit numbers. We confirmed that odd numbers seemed masculine while even numbers seemed feminine. Although both men and women showed this pattern, it was more pronounced among women. We also examined whether this pattern holds for automatic as well as deliberated reactions. Results of an Implicit Association Test showed that it did, but only among the women. The implicit and explicit patterns of numerical gender ascription were moderately correlated. The second study examined explicit perceptions of 2-digit numbers. Again, women viewed odd numbers as more masculine and less feminine than even numbers. However, men viewed 2-digit numbers as relatively masculine, regardless of whether they were even or odd. These results indicate that women and men impute gender to numbers in different ways and to different extents. We discuss possible implications for understanding how people relate to and are influenced by numbers in a variety of real-life contexts.

Two dimensional fully nonlinear numerical wave tank based on the BEM

NASA Astrophysics Data System (ADS)

Sun, Zhe; Pang, Yongjie; Li, Hongwei

2012-12-01

The development of a two dimensional numerical wave tank (NWT) with a rocker or piston type wavemaker based on the high order boundary element method (BEM) and mixed Eulerian-Lagrangian (MEL) is examined. The cauchy principle value (CPV) integral is calculated by a special Gauss type quadrature and a change of variable. In addition the explicit truncated Taylor expansion formula is employed in the time-stepping process. A modified double nodes method is assumed to tackle the corner problem, as well as the damping zone technique is used to absorb the propagation of the free surface wave at the end of the tank. A variety of waves are generated by the NWT, for example; a monochromatic wave, solitary wave and irregular wave. The results confirm the NWT model is efficient and stable.

Penalty methods for the numerical solution of American multi-asset option problems

NASA Astrophysics Data System (ADS)

Nielsen, Bjørn Fredrik; Skavhaug, Ola; Tveito, Aslak

2008-12-01

We derive and analyze a penalty method for solving American multi-asset option problems. A small, non-linear penalty term is added to the Black-Scholes equation. This approach gives a fixed solution domain, removing the free and moving boundary imposed by the early exercise feature of the contract. Explicit, implicit and semi-implicit finite difference schemes are derived, and in the case of independent assets, we prove that the approximate option prices satisfy some basic properties of the American option problem. Several numerical experiments are carried out in order to investigate the performance of the schemes. We give examples indicating that our results are sharp. Finally, the experiments indicate that in the case of correlated underlying assets, the same properties are valid as in the independent case.

Feasibility study for a numerical aerodynamic simulation facility. Volume 3: FMP language specification/user manual

NASA Technical Reports Server (NTRS)

Kenner, B. G.; Lincoln, N. R.

1979-01-01

The manual is intended to show the revisions and additions to the current STAR FORTRAN. The changes are made to incorporate an FMP (Flow Model Processor) for use in the Numerical Aerodynamic Simulation Facility (NASF) for the purpose of simulating fluid flow over three-dimensional bodies in wind tunnel environments and in free space. The FORTRAN programming language for the STAR-100 computer contains both CDC and unique STAR extensions to the standard FORTRAN. Several of the STAR FORTRAN extensions to standard FOR-TRAN allow the FORTRAN user to exploit the vector processing capabilities of the STAR computer. In STAR FORTRAN, vectors can be expressed with an explicit notation, functions are provided that return vector results, and special call statements enable access to any machine instruction.

High-Order Semi-Discrete Central-Upwind Schemes for Multi-Dimensional Hamilton-Jacobi Equations

NASA Technical Reports Server (NTRS)

Bryson, Steve; Levy, Doron; Biegel, Bran R. (Technical Monitor)

2002-01-01

We present high-order semi-discrete central-upwind numerical schemes for approximating solutions of multi-dimensional Hamilton-Jacobi (HJ) equations. This scheme is based on the use of fifth-order central interpolants like those developed in [1], in fluxes presented in [3]. These interpolants use the weighted essentially nonoscillatory (WENO) approach to avoid spurious oscillations near singularities, and become "central-upwind" in the semi-discrete limit. This scheme provides numerical approximations whose error is as much as an order of magnitude smaller than those in previous WENO-based fifth-order methods [2, 1]. Thee results are discussed via examples in one, two and three dimensions. We also pregnant explicit N-dimensional formulas for the fluxes, discuss their monotonicity and tl!e connection between this method and that in [2].

Experience in using a numerical scheme with artificial viscosity at solving the Riemann problem for a multi-fluid model of multiphase flow

NASA Astrophysics Data System (ADS)

Bulovich, S. V.; Smirnov, E. M.

2018-05-01

The paper covers application of the artificial viscosity technique to numerical simulation of unsteady one-dimensional multiphase compressible flows on the base of the multi-fluid approach. The system of the governing equations is written under assumption of the pressure equilibrium between the "fluids" (phases). No interfacial exchange is taken into account. A model for evaluation of the artificial viscosity coefficient that (i) assumes identity of this coefficient for all interpenetrating phases and (ii) uses the multiphase-mixture Wood equation for evaluation of a scale speed of sound has been suggested. Performance of the artificial viscosity technique has been evaluated via numerical solution of a model problem of pressure discontinuity breakdown in a three-fluid medium. It has been shown that a relatively simple numerical scheme, explicit and first-order, combined with the suggested artificial viscosity model, predicts a physically correct behavior of the moving shock and expansion waves, and a subsequent refinement of the computational grid results in a monotonic approaching to an asymptotic time-dependent solution, without non-physical oscillations.

A numerical study of confined turbulent jets

NASA Technical Reports Server (NTRS)

Zhu, J.; Shih, T.-H.

1993-01-01

A numerical investigation is reported of turbulent incompressible jets confined in two ducts, one cylindrical and the other conical with a 5 degree divergence. In each case, three Craya-Curtet numbers are considered which correspond, respectively, to flow situations with no moderate and strong recirculation. Turbulence closure is achieved by using the k-epsilon model and a recently proposed realizable Reynolds stress algebraic equation model that relates the Reynolds stresses explicitly to the quadratic terms of the mean velocity gradients and ensures the positiveness of each component of the turbulent kinetic energy. Calculations are carried out with a finite-volume procedure using boundary-fitted curvilinear coordinates. A second-order accurate, bounded convection scheme and sufficiently fine grids are used to prevent the solutions from being contaminated by numerical diffusion. The calculated results are compared extensively with the available experimental data. It is shown that the numerical methods presented are capable of capturing the essential flow features observed in the experiments and that the realizable Reynolds stress algebraic equation model performs much better than the k-epsilon model for this class of flows of great practical importance.

Random element method for numerical modeling of diffusional processes

NASA Technical Reports Server (NTRS)

Ghoniem, A. F.; Oppenheim, A. K.

1982-01-01

The random element method is a generalization of the random vortex method that was developed for the numerical modeling of momentum transport processes as expressed in terms of the Navier-Stokes equations. The method is based on the concept that random walk, as exemplified by Brownian motion, is the stochastic manifestation of diffusional processes. The algorithm based on this method is grid-free and does not require the diffusion equation to be discritized over a mesh, it is thus devoid of numerical diffusion associated with finite difference methods. Moreover, the algorithm is self-adaptive in space and explicit in time, resulting in an improved numerical resolution of gradients as well as a simple and efficient computational procedure. The method is applied here to an assortment of problems of diffusion of momentum and energy in one-dimension as well as heat conduction in two-dimensions in order to assess its validity and accuracy. The numerical solutions obtained are found to be in good agreement with exact solution except for a statistical error introduced by using a finite number of elements, the error can be reduced by increasing the number of elements or by using ensemble averaging over a number of solutions.

Reliability-Based Stability Analysis of Rock Slopes Using Numerical Analysis and Response Surface Method

NASA Astrophysics Data System (ADS)

Dadashzadeh, N.; Duzgun, H. S. B.; Yesiloglu-Gultekin, N.

2017-08-01

While advanced numerical techniques in slope stability analysis are successfully used in deterministic studies, they have so far found limited use in probabilistic analyses due to their high computation cost. The first-order reliability method (FORM) is one of the most efficient probabilistic techniques to perform probabilistic stability analysis by considering the associated uncertainties in the analysis parameters. However, it is not possible to directly use FORM in numerical slope stability evaluations as it requires definition of a limit state performance function. In this study, an integrated methodology for probabilistic numerical modeling of rock slope stability is proposed. The methodology is based on response surface method, where FORM is used to develop an explicit performance function from the results of numerical simulations. The implementation of the proposed methodology is performed by considering a large potential rock wedge in Sumela Monastery, Turkey. The accuracy of the developed performance function to truly represent the limit state surface is evaluated by monitoring the slope behavior. The calculated probability of failure is compared with Monte Carlo simulation (MCS) method. The proposed methodology is found to be 72% more efficient than MCS, while the accuracy is decreased with an error of 24%.

Explicit formulas for effective piezoelectric coefficients of ferroelectric 0-3 composites based on effective medium theory

NASA Astrophysics Data System (ADS)

Wong, C. K.; Poon, Y. M.; Shin, F. G.

2003-01-01

Explicit formulas were derived for the effective piezoelectric stress coefficients of a 0-3 composite of ferroelectric spherical particles in a ferroelectric matrix which were then combined to give the more commonly used strain coefficients. Assuming that the elastic stiffness of the inclusion phase is sufficiently larger than that of the matrix phase, the previously derived explicit expressions for the case of a low volume concentration of inclusion particles [C. K. Wong, Y. M. Poon, and F. G. Shin, Ferroelectrics 264, 39 (2001); J. Appl. Phys. 90, 4690 (2001)] were "transformed" analytically by an effective medium theory (EMT) with appropriate approximations, to suit the case of a more concentrated suspension. Predictions of the EMT expressions were compared with the experimental values of composites of lead zirconate titanate ceramic particles dispersed in polyvinylidene fluoride and polyvinylidene fluoride-trifluoroethylene copolymer, reported by Furukawa [IEEE Trans. Electr. Insul. 24, 375 (1989)] and by Ng et al. [IEEE Trans. Ultrason. Ferroelectr. Freq. Control 47, 1308 (2000)] respectively. Fairly good agreement was obtained. Comparisons with other predictions, including the predictions given by numerically solving the EMT scheme, were also made. It was found that the analytic and numeric EMT schemes agreed with each other very well for an inclusion of volume fraction not exceeding 60%.

Explicitly correlated coupled-cluster theory using cusp conditions. II. Treatment of connected triple excitations.

PubMed

Köhn, Andreas

2010-11-07

The coupled-cluster singles and doubles method augmented with single Slater-type correlation factors (CCSD-F12) determined by the cusp conditions (also denoted as SP ansatz) yields results close to the basis set limit with only small overhead compared to conventional CCSD. Quantitative calculations on many-electron systems, however, require to include the effect of connected triple excitations at least. In this contribution, the recently proposed [A. Köhn, J. Chem. Phys. 130, 131101 (2009)] extended SP ansatz and its application to the noniterative triples correction CCSD(T) is reviewed. The approach allows to include explicit correlation into connected triple excitations without introducing additional unknown parameters. The explicit expressions are presented and analyzed, and possible simplifications to arrive at a computationally efficient scheme are suggested. Numerical tests based on an implementation obtained by an automated approach are presented. Using a partial wave expansion for the neon atom, we can show that the proposed ansatz indeed leads to the expected (L(max)+1)(-7) convergence of the noniterative triples correction, where L(max) is the maximum angular momentum in the orbital expansion. Further results are reported for a test set of 29 molecules, employing Peterson's F12-optimized basis sets. We find that the customary approach of using the conventional noniterative triples correction on top of a CCSD-F12 calculation leads to significant basis set errors. This, however, is not always directly visible for total CCSD(T) energies due to fortuitous error compensation. The new approach offers a thoroughly explicitly correlated CCSD(T)-F12 method with improved basis set convergence of the triples contributions to both total and relative energies.

Influence of biases in numerical magnitude allocation on human prosocial decision making.

PubMed

Arshad, Qadeer; Nigmatullina, Yuliya; Siddiqui, Shuaib; Franka, Mustafa; Mediratta, Saniya; Ramachandaran, Sanjeev; Lobo, Rhannon; Malhotra, Paresh A; Roberts, R E; Bronstein, Adolfo M

2017-12-01

Over the past decade neuroscientific research has attempted to probe the neurobiological underpinnings of human prosocial decision making. Such research has almost ubiquitously employed tasks such as the dictator game or similar variations (i.e., ultimatum game). Considering the explicit numerical nature of such tasks, it is surprising that the influence of numerical cognition on decision making during task performance remains unknown. While performing these tasks, participants typically tend to anchor on a 50:50 split that necessitates an explicit numerical judgement (i.e., number-pair bisection). Accordingly, we hypothesize that the decision-making process during the dictator game recruits overlapping cognitive processes to those known to be engaged during number-pair bisection. We observed that biases in numerical magnitude allocation correlated with the formulation of decisions during the dictator game. That is, intrinsic biases toward smaller numerical magnitudes were associated with the formulation of less favorable decisions, whereas biases toward larger magnitudes were associated with more favorable choices. We proceeded to corroborate this relationship by subliminally and systematically inducing biases in numerical magnitude toward either higher or lower numbers using a visuo-vestibular stimulation paradigm. Such subliminal alterations in numerical magnitude allocation led to proportional and corresponding changes to an individual's decision making during the dictator game. Critically, no relationship was observed between neither intrinsic nor induced biases in numerical magnitude on decision making when assessed using a nonnumerical-based prosocial questionnaire. Our findings demonstrate numerical influences on decisions formulated during the dictator game and highlight the necessity to control for confounds associated with numerical cognition in human decision-making paradigms. NEW & NOTEWORTHY We demonstrate that intrinsic biases in numerical magnitude can directly predict the amount of money donated by an individual to an anonymous stranger during the dictator game. Furthermore, subliminally inducing perceptual biases in numerical-magnitude allocation can actively drive prosocial choices in the corresponding direction. Our findings provide evidence for numerical influences on decision making during performance of the dictator game. Accordingly, without the implementation of an adequate control for numerical influences, the dictator game and other tasks with an inherent numerical component (i.e., ultimatum game) should be employed with caution in the assessment of human behavior. Copyright © 2017 the American Physiological Society.

Comparison of Several Numerical Methods for Simulation of Compressible Shear Layers

NASA Technical Reports Server (NTRS)

Kennedy, Christopher A.; Carpenter, Mark H.

1997-01-01

An investigation is conducted on several numerical schemes for use in the computation of two-dimensional, spatially evolving, laminar variable-density compressible shear layers. Schemes with various temporal accuracies and arbitrary spatial accuracy for both inviscid and viscous terms are presented and analyzed. All integration schemes use explicit or compact finite-difference derivative operators. Three classes of schemes are considered: an extension of MacCormack's original second-order temporally accurate method, a new third-order variant of the schemes proposed by Rusanov and by Kutier, Lomax, and Warming (RKLW), and third- and fourth-order Runge-Kutta schemes. In each scheme, stability and formal accuracy are considered for the interior operators on the convection-diffusion equation U(sub t) + aU(sub x) = alpha U(sub xx). Accuracy is also verified on the nonlinear problem, U(sub t) + F(sub x) = 0. Numerical treatments of various orders of accuracy are chosen and evaluated for asymptotic stability. Formally accurate boundary conditions are derived for several sixth- and eighth-order central-difference schemes. Damping of high wave-number data is accomplished with explicit filters of arbitrary order. Several schemes are used to compute variable-density compressible shear layers, where regions of large gradients exist.

Transport of reacting solutes subject to a moving dissolution boundary: Numerical methods and solutions

USGS Publications Warehouse

Willis, Catherine; Rubin, Jacob

1987-01-01

A moving boundary problem which arises during transport with precipitation-dissolution reactions is solved by three different numerical methods. Two of these methods (one explicit and one implicit) are based on an integral formulation of mass balance and lead to an approximation of a weak solution. These methods are compared to a front-tracking scheme. Although the two approaches are conceptually different, the numerical solutions showed good agreement. As the ratio of dispersion to convection decreases, the methods based on the integral formulation become computationally more efficient. Specific reactions were modeled to examine the dependence of the system on the physical and chemical parameters. Although the water flow rate does not explicitly appear in the equation for the velocity of the moving boundary, the speed of the boundary depends more on the flux rate than on the dispersion coefficient. The discontinuity in the gradient of the solute concentration profile at the boundary increases with convection and with the initial concentration of the mineral. Our implicit method is extended to allow participation of the solutes in complexation reactions as well as the precipitation-dissolution reaction. This extension is easily made and does not change the basic method.

Experiments and Dynamic Finite Element Analysis of a Wire-Rope Rockfall Protective Fence

NASA Astrophysics Data System (ADS)

Tran, Phuc Van; Maegawa, Koji; Fukada, Saiji

2013-09-01

The imperative need to protect structures in mountainous areas against rockfall has led to the development of various protection methods. This study introduces a new type of rockfall protection fence made of posts, wire ropes, wire netting and energy absorbers. The performance of this rock fence was verified in both experiments and dynamic finite element analysis. In collision tests, a reinforced-concrete block rolled down a natural slope and struck the rock fence at the end of the slope. A specialized system of measuring instruments was employed to accurately measure the acceleration of the block without cable connection. In particular, the performance of two energy absorbers, which contribute also to preventing wire ropes from breaking, was investigated to determine the best energy absorber. In numerical simulation, a commercial finite element code having explicit dynamic capabilities was employed to create models of the two full-scale tests. To facilitate simulation, certain simplifying assumptions for mechanical data of each individual component of the rock fence and geometrical data of the model were adopted. Good agreement between numerical simulation and experimental data validated the numerical simulation. Furthermore, the results of numerical simulation helped highlight limitations of the testing method. The results of numerical simulation thus provide a deeper understanding of the structural behavior of individual components of the rock fence during rockfall impact. More importantly, numerical simulations can be used not only as supplements to or substitutes for full-scale tests but also in parametric study and design.

Dynamical shift condition for unequal mass black hole binaries

DOE Office of Scientific and Technical Information (OSTI.GOV)

Mueller, Doreen; Grigsby, Jason; Bruegmann, Bernd

Certain numerical frameworks used for the evolution of binary black holes make use of a gamma driver, which includes a damping factor. Such simulations typically use a constant value for damping. However, it has been found that very specific values of the damping factor are needed for the calculation of unequal mass binaries. We examine carefully the role this damping plays and provide two explicit, nonconstant forms for the damping to be used with mass ratios further from one. Our analysis of the resultant waveforms compares well against the constant damping case.

On the thermodynamics of the Swift-Hohenberg theory

NASA Astrophysics Data System (ADS)

Espath, L. F. R.; Sarmiento, A. F.; Dalcin, L.; Calo, V. M.

2017-11-01

We present the microbalance including the microforces, the first- and second-order microstresses for the Swift-Hohenberg equation concomitantly with their constitutive equations, which are consistent with the free-energy imbalance. We provide an explicit form for the microstress structure for a free-energy functional endowed with second-order spatial derivatives. Additionally, we generalize the Swift-Hohenberg theory via a proper constitutive process. Finally, we present one highly resolved three-dimensional numerical simulation to demonstrate the particular form of the resulting microstresses and their interactions in the evolution of the Swift-Hohenberg equation.

A Coupled Layerwise Analysis of the Thermopiezoelectric Response of Smart Composite Beams Beams

NASA Technical Reports Server (NTRS)

Lee, H.-J.; Saravanos, D. A.

1995-01-01

Thermal effects are incorporated into previously developed discrete layer mechanics for piezoelectric composite beam structures. The updated mechanics explicitly account for the complete coupled thermoelectromechanical response of smart composite beams. This unified representation leads to an inherent capability to model both the sensory and actuator responses of piezoelectric composite beams in a thermal environment. Finite element equations are developed and numerical results are presented to demonstrate the capability of the current formulation to represent the behavior of both sensory and active smart structures under thermal loadings.

On a comparison of two schemes in sequential data assimilation

NASA Astrophysics Data System (ADS)

Grishina, Anastasiia A.; Penenko, Alexey V.

2017-11-01

This paper is focused on variational data assimilation as an approach to mathematical modeling. Realization of the approach requires a sequence of connected inverse problems with different sets of observational data to be solved. Two variational data assimilation schemes, "implicit" and "explicit", are considered in the article. Their equivalence is shown and the numerical results are given on a basis of non-linear Robertson system. To avoid the "inverse problem crime" different schemes were used to produce synthetic measurement and to solve the data assimilation problem.

KEWPIE: A dynamical cascade code for decaying exited compound nuclei

NASA Astrophysics Data System (ADS)

Bouriquet, Bertrand; Abe, Yasuhisa; Boilley, David

2004-05-01

A new dynamical cascade code for decaying hot nuclei is proposed and specially adapted to the synthesis of super-heavy nuclei. For such a case, the interesting channel is of the tiny fraction that will decay through particles emission, thus the code avoids classical Monte-Carlo methods and proposes a new numerical scheme. The time dependence is explicitely taken into account in order to cope with the fact that fission decay rate might not be constant. The code allows to evaluate both statistical and dynamical observables. Results are successfully compared to experimental data.

Analysis of Time Filters in Multistep Methods

NASA Astrophysics Data System (ADS)

Hurl, Nicholas

Geophysical ow simulations have evolved sophisticated implicit-explicit time stepping methods (based on fast-slow wave splittings) followed by time filters to control any unstable models that result. Time filters are modular and parallel. Their effect on stability of the overall process has been tested in numerous simulations, but never analyzed. Stability is proven herein for the Crank-Nicolson Leapfrog (CNLF) method with the Robert-Asselin (RA) time filter and for the Crank-Nicolson Leapfrog method with the Robert-Asselin-Williams (RAW) time filter for systems by energy methods. We derive an equivalent multistep method for CNLF+RA and CNLF+RAW and stability regions are obtained. The time step restriction for energy stability of CNLF+RA is smaller than CNLF and CNLF+RAW time step restriction is even smaller. Numerical tests find that RA and RAW add numerical dissipation. This thesis also shows that all modes of the Crank-Nicolson Leap Frog (CNLF) method are asymptotically stable under the standard timestep condition.

Three-dimensional MHD (magnetohydrodynamic) flows in rectangular ducts of liquid-metal-cooled blankets

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hua, T.Q.; Walker, J.S.; Picologlou, B.F.

1988-07-01

Magnetohydrodynamic flows of liquid metals in rectangular ducts with thin conducting walls in the presence of strong nonuniform transverse magnetic fields are examined. The interaction parameter and Hartmann number are assumed to be large, whereas the magnetic Reynolds number is assumed to be small. Under these assumptions, viscous and inertial effects are confined in very thin boundary layers adjacent to the walls. A significant fraction of the fluid flow is concentrated in the boundary layers adjacent to the side walls which are parallel to the magnetic field. This paper describes the analysis and numerical methods for obtaining 3-D solutions formore » flow parameters outside these layers, without solving explicitly for the layers themselves. Numerical solutions are presented for cases which are relevant to the flows of liquid metals in fusion reactor blankets. Experimental results obtained from the ALEX experiments at Argonne National Laboratory are used to validate the numerical code. In general, the agreement is excellent. 5 refs., 14 figs.« less

Perforation of thin aluminum alloy plates by blunt projectiles: An experimental and numerical investigation

NASA Astrophysics Data System (ADS)

Wei, G.; Zhang, W.

2014-04-01

Reducing the armor weight has become a research focus in terms of armored material. Due to high strength-to-density ratio, aluminum alloy has become a potential light armored material. In this study, both lab-scale ballistic test and finite element simulation were adopted to examine the ballistic resistance of aluminum alloy targets. Blunt high strength steel projectiles with 12.7 mm diameter were launched by light gas gun against 3.3 mm thickness 7A04 aluminum alloy plates at a velocity of 90~170 m/s. The ballistic limit velocity was obtained. Plugging failure and obvious structure deformation of targets were observed. Corresponding 2D finite element simulations were conducted by ABAQUS/EXPLICIT combined with material performance testing. The validity of numerical simulations was verified by comparing with the experimental results. Detailed analysis of the failure modes and characters of the targets were carried out to reveal the target damage mechanism combined with the numerical simulation.

Gaussian process regression to accelerate geometry optimizations relying on numerical differentiation

NASA Astrophysics Data System (ADS)

Schmitz, Gunnar; Christiansen, Ove

2018-06-01

We study how with means of Gaussian Process Regression (GPR) geometry optimizations, which rely on numerical gradients, can be accelerated. The GPR interpolates a local potential energy surface on which the structure is optimized. It is found to be efficient to combine results on a low computational level (HF or MP2) with the GPR-calculated gradient of the difference between the low level method and the target method, which is a variant of explicitly correlated Coupled Cluster Singles and Doubles with perturbative Triples correction CCSD(F12*)(T) in this study. Overall convergence is achieved if both the potential and the geometry are converged. Compared to numerical gradient-based algorithms, the number of required single point calculations is reduced. Although introducing an error due to the interpolation, the optimized structures are sufficiently close to the minimum of the target level of theory meaning that the reference and predicted minimum only vary energetically in the μEh regime.

On a free-surface problem with moving contact line: From variational principles to stable numerical approximations

NASA Astrophysics Data System (ADS)

Fumagalli, Ivan; Parolini, Nicola; Verani, Marco

2018-02-01

We analyze a free-surface problem described by time-dependent Navier-Stokes equations. Surface tension, capillary effects and wall friction are taken into account in the evolution of the system, influencing the motion of the contact line - where the free surface hits the wall - and of the dynamics of the contact angle. The differential equations governing the phenomenon are first derived from the variational principle of minimum reduced dissipation, and then discretized by means of the ALE approach. The numerical properties of the resulting scheme are investigated, drawing a parallel with the physical properties holding at the continuous level. Some instability issues are addressed in detail, in the case of an explicit treatment of the geometry, and novel additional terms are introduced in the discrete formulation in order to damp the instabilities. Numerical tests assess the suitability of the approach, the influence of the parameters, and the effectiveness of the new stabilizing terms.

Prediction of Ignition of High Explosive When Submitted To Impact

NASA Astrophysics Data System (ADS)

Picart, Didier; Delmaire-Sizes, Franck; Gruau, Cyril; Trumel, Herve

2009-06-01

High explosive structures may unintentionally ignite and transit to deflagration or detonation, when subjected to mechanical loadings, such as low velocity impact. We focus our attention on ignition. The Browning and Scammon [1] criterion has been adapted. A concrete like constitutive law is derived, with an up-to-date experimental characterization. These models have been implemented in Abaqus/Explicit [2]. Numerical simulations are used to calibrate the ignition threshold. The presentation or the poster will detail the main assumptions, the models (Browning et al, mechanical behavior) and the calibration procedure. Comparisons between numerical results and experiments [3] will show the interest of this method but also its limitations (numerical artifacts, lack of mechanical data, misinterpretation of reactive tests). [1] R. Browning and R. Scammon, Shock compression of condensed matter, pp. 987-990, (2001). [2] C. Gruau, D. Picart et al., 17^th Dymat technical meeting, Cambridge, UK, (2007). [3] F. Delmaire-Sizes et al., 3^rd International symposium on energetic materials, Tokyo, Japan, (2008).

Perturbative triples correction for local pair natural orbital based explicitly correlated CCSD(F12*) using Laplace transformation techniques.

PubMed

Schmitz, Gunnar; Hättig, Christof

2016-12-21

We present an implementation of pair natural orbital coupled cluster singles and doubles with perturbative triples, PNO-CCSD(T), which avoids the quasi-canonical triples approximation (T0) where couplings due to off-diagonal Fock matrix elements are neglected. A numerical Laplace transformation of the canonical expression for the perturbative (T) triples correction is used to avoid an I/O and storage bottleneck for the triples amplitudes. Results for a test set of reaction energies show that only very few Laplace grid points are needed to obtain converged energy differences and that PNO-CCSD(T) is a more robust approximation than PNO-CCSD(T0) with a reduced mean absolute deviation from canonical CCSD(T) results. We combine the PNO-based (T) triples correction with the explicitly correlated PNO-CCSD(F12*) method and investigate the use of specialized F12-PNOs in the conventional triples correction. We find that no significant additional errors are introduced and that PNO-CCSD(F12*)(T) can be applied in a black box manner.

Efficient Robust Optimization of Metal Forming Processes using a Sequential Metamodel Based Strategy

NASA Astrophysics Data System (ADS)

Wiebenga, J. H.; Klaseboer, G.; van den Boogaard, A. H.

2011-08-01

The coupling of Finite Element (FE) simulations to mathematical optimization techniques has contributed significantly to product improvements and cost reductions in the metal forming industries. The next challenge is to bridge the gap between deterministic optimization techniques and the industrial need for robustness. This paper introduces a new and generally applicable structured methodology for modeling and solving robust optimization problems. Stochastic design variables or noise variables are taken into account explicitly in the optimization procedure. The metamodel-based strategy is combined with a sequential improvement algorithm to efficiently increase the accuracy of the objective function prediction. This is only done at regions of interest containing the optimal robust design. Application of the methodology to an industrial V-bending process resulted in valuable process insights and an improved robust process design. Moreover, a significant improvement of the robustness (>2σ) was obtained by minimizing the deteriorating effects of several noise variables. The robust optimization results demonstrate the general applicability of the robust optimization strategy and underline the importance of including uncertainty and robustness explicitly in the numerical optimization procedure.

The Coastal Dynamics of Heterogeneous Sedimentary Environments: Numerical Modeling of Nearshore Hydrodynamics and Sediment Transport

DTIC Science & Technology

2010-05-10

supplied boundary data for Hurricane Katrina (Keen, Furukawa et al. 2006; Keen, Slingerland et al. 2010). The numerical models discussed in this report...explicitly. NCOM can be nested to a coarse-grid model to supply boundary conditions at the open boundary of the domain. NCOM has been validated at global...circulation study of Mississippi Sound (Keen 2002), which supplied steady currents for the nearshore erosion problem discussed in this report

Mean-Field Description of Ionic Size Effects with Non-Uniform Ionic Sizes: A Numerical Approach

PubMed Central

Zhou, Shenggao; Wang, Zhongming; Li, Bo

2013-01-01

Ionic size effects are significant in many biological systems. Mean-field descriptions of such effects can be efficient but also challenging. When ionic sizes are different, explicit formulas in such descriptions are not available for the dependence of the ionic concentrations on the electrostatic potential, i.e., there is no explicit, Boltzmann type distributions. This work begins with a variational formulation of the continuum electrostatics of an ionic solution with such non-uniform ionic sizes as well as multiple ionic valences. An augmented Lagrange multiplier method is then developed and implemented to numerically solve the underlying constrained optimization problem. The method is shown to be accurate and efficient, and is applied to ionic systems with non-uniform ionic sizes such as the sodium chloride solution. Extensive numerical tests demonstrate that the mean-field model and numerical method capture qualitatively some significant ionic size effects, particularly those for multivalent ionic solutions, such as the stratification of multivalent counterions near a charged surface. The ionic valence-to-volume ratio is found to be the key physical parameter in the stratification of concentrations. All these are not well described by the classical Poisson–Boltzmann theory, or the generalized Poisson–Boltzmann theory that treats uniform ionic sizes. Finally, various issues such as the close packing, limitation of the continuum model, and generalization of this work to molecular solvation are discussed. PMID:21929014

A review of numerical techniques approaching microstructures of crystalline rocks

NASA Astrophysics Data System (ADS)

Zhang, Yahui; Wong, Louis Ngai Yuen

2018-06-01

The macro-mechanical behavior of crystalline rocks including strength, deformability and failure pattern are dominantly influenced by their grain-scale structures. Numerical technique is commonly used to assist understanding the complicated mechanisms from a microscopic perspective. Each numerical method has its respective strengths and limitations. This review paper elucidates how numerical techniques take geometrical aspects of the grain into consideration. Four categories of numerical methods are examined: particle-based methods, block-based methods, grain-based methods, and node-based methods. Focusing on the grain-scale characters, specific relevant issues including increasing complexity of micro-structure, deformation and breakage of model elements, fracturing and fragmentation process are described in more detail. Therefore, the intrinsic capabilities and limitations of different numerical approaches in terms of accounting for the micro-mechanics of crystalline rocks and their phenomenal mechanical behavior are explicitly presented.

Experimental and numerical study of drill bit drop tests on Kuru granite.

PubMed

Fourmeau, Marion; Kane, Alexandre; Hokka, Mikko

2017-01-28

This paper presents an experimental and numerical study of Kuru grey granite impacted with a seven-buttons drill bit mounted on an instrumented drop test machine. The force versus displacement curves during the impact, so-called bit-rock interaction (BRI) curves, were obtained using strain gauge measurements for two levels of impact energy. Moreover, the volume of removed rock after each drop test was evaluated by stereo-lithography (three-dimensional surface reconstruction). A modified version of the Holmquist-Johnson-Cook (MHJC) material model was calibrated using Kuru granite test results available from the literature. Numerical simulations of the single drop tests were carried out using the MHJC model available in the LS-DYNA explicit finite-element solver. The influence of the impact energy and additional confining pressure on the BRI curves and the volume of the removed rock is discussed. In addition, the influence of the rock surface shape before impact was evaluated using two different mesh geometries: a flat surface and a hyperbolic surface. The experimental and numerical results are compared and discussed in terms of drilling efficiency through the mechanical specific energy.This article is part of the themed issue 'Experimental testing and modelling of brittle materials at high strain rates'. © 2016 The Author(s).

Experimental and numerical study of drill bit drop tests on Kuru granite

PubMed Central

Kane, Alexandre; Hokka, Mikko

2017-01-01

This paper presents an experimental and numerical study of Kuru grey granite impacted with a seven-buttons drill bit mounted on an instrumented drop test machine. The force versus displacement curves during the impact, so-called bit–rock interaction (BRI) curves, were obtained using strain gauge measurements for two levels of impact energy. Moreover, the volume of removed rock after each drop test was evaluated by stereo-lithography (three-dimensional surface reconstruction). A modified version of the Holmquist–Johnson–Cook (MHJC) material model was calibrated using Kuru granite test results available from the literature. Numerical simulations of the single drop tests were carried out using the MHJC model available in the LS-DYNA explicit finite-element solver. The influence of the impact energy and additional confining pressure on the BRI curves and the volume of the removed rock is discussed. In addition, the influence of the rock surface shape before impact was evaluated using two different mesh geometries: a flat surface and a hyperbolic surface. The experimental and numerical results are compared and discussed in terms of drilling efficiency through the mechanical specific energy. This article is part of the themed issue ‘Experimental testing and modelling of brittle materials at high strain rates’. PMID:27956511

Parameterization of subgrid-scale stress by the velocity gradient tensor

NASA Technical Reports Server (NTRS)

Lund, Thomas S.; Novikov, E. A.

1993-01-01

The objective of this work is to construct and evaluate subgrid-scale models that depend on both the strain rate and the vorticity. This will be accomplished by first assuming that the subgrid-scale stress is a function of the strain and rotation rate tensors. Extensions of the Caley-Hamilton theorem can then be used to write the assumed functional dependence explicitly in the form of a tensor polynomial involving products of the strain and rotation rates. Finally, use of this explicit expression as a subgrid-scale model will be evaluated using direct numerical simulation data for homogeneous, isotropic turbulence.

Initialization and assimilation of cloud and rainwater in a regional model

NASA Technical Reports Server (NTRS)

Raymond, William H.; Olson, William S.

1990-01-01

The initialization and assimilation of cloud and rainwater quantities in a mesoscale regional model was examined. Forecasts of explicit cloud and rainwater are made using conservation equations. The physical processes include condensation, evaporation, autoconversion, accretion, and the removal of rainwater by fallout. These physical processes, some of which are parameterized, represent source and sink in terms in the conservation equations. The question of how to initialize the explicit liquid water calculations in numerical models and how to retain information about precipitation processes during the 4-D assimilation cycle are important issues that are addressed.

An explicit scheme for ohmic dissipation with smoothed particle magnetohydrodynamics

NASA Astrophysics Data System (ADS)

Tsukamoto, Yusuke; Iwasaki, Kazunari; Inutsuka, Shu-ichiro

2013-09-01

In this paper, we present an explicit scheme for Ohmic dissipation with smoothed particle magnetohydrodynamics (SPMHD). We propose an SPH discretization of Ohmic dissipation and solve Ohmic dissipation part of induction equation with the super-time-stepping method (STS) which allows us to take a longer time step than Courant-Friedrich-Levy stability condition. Our scheme is second-order accurate in space and first-order accurate in time. Our numerical experiments show that optimal choice of the parameters of STS for Ohmic dissipation of SPMHD is νsts ˜ 0.01 and Nsts ˜ 5.

Adaptive implicit-explicit and parallel element-by-element iteration schemes

NASA Technical Reports Server (NTRS)

Tezduyar, T. E.; Liou, J.; Nguyen, T.; Poole, S.

1989-01-01

Adaptive implicit-explicit (AIE) and grouped element-by-element (GEBE) iteration schemes are presented for the finite element solution of large-scale problems in computational mechanics and physics. The AIE approach is based on the dynamic arrangement of the elements into differently treated groups. The GEBE procedure, which is a way of rewriting the EBE formulation to make its parallel processing potential and implementation more clear, is based on the static arrangement of the elements into groups with no inter-element coupling within each group. Various numerical tests performed demonstrate the savings in the CPU time and memory.

Studying Activity Series of Metals.

ERIC Educational Resources Information Center

Hoon, Tien-Ghun; And Others

1995-01-01

Presents teaching strategies that illustrate the linking together of numerous chemical concepts involving the activity of metals (quantitative analysis, corrosion, and electrolysis) through the use of deep-level processing strategies. Concludes that making explicit links in the process of teaching chemistry can lead effectively to meaningful…

Free and Open Source GIS Tools: Role and Relevance in the Environmental Assessment Community

EPA Science Inventory

The presence of an explicit geographical context in most environmental decisions can complicate assessment and selection of management options. These decisions typically involve numerous data sources, complex environmental and ecological processes and their associated models, ris...

Interpolation Inequalities and Spectral Estimates for Magnetic Operators

NASA Astrophysics Data System (ADS)

Dolbeault, Jean; Esteban, Maria J.; Laptev, Ari; Loss, Michael

2018-05-01

We prove magnetic interpolation inequalities and Keller-Lieb-Thir-ring estimates for the principal eigenvalue of magnetic Schr{\\"o}dinger operators. We establish explicit upper and lower bounds for the best constants and show by numerical methods that our theoretical estimates are accurate.

Explicit Filtering Based Low-Dose Differential Phase Reconstruction Algorithm with the Grating Interferometry.

PubMed

Jiang, Xiaolei; Zhang, Li; Zhang, Ran; Yin, Hongxia; Wang, Zhenchang

2015-01-01

X-ray grating interferometry offers a novel framework for the study of weakly absorbing samples. Three kinds of information, that is, the attenuation, differential phase contrast (DPC), and dark-field images, can be obtained after a single scanning, providing additional and complementary information to the conventional attenuation image. Phase shifts of X-rays are measured by the DPC method; hence, DPC-CT reconstructs refraction indexes rather than attenuation coefficients. In this work, we propose an explicit filtering based low-dose differential phase reconstruction algorithm, which enables reconstruction from reduced scanning without artifacts. The algorithm adopts a differential algebraic reconstruction technique (DART) with the explicit filtering based sparse regularization rather than the commonly used total variation (TV) method. Both the numerical simulation and the biological sample experiment demonstrate the feasibility of the proposed algorithm.

Explicit Filtering Based Low-Dose Differential Phase Reconstruction Algorithm with the Grating Interferometry

PubMed Central

Zhang, Li; Zhang, Ran; Yin, Hongxia; Wang, Zhenchang

2015-01-01

X-ray grating interferometry offers a novel framework for the study of weakly absorbing samples. Three kinds of information, that is, the attenuation, differential phase contrast (DPC), and dark-field images, can be obtained after a single scanning, providing additional and complementary information to the conventional attenuation image. Phase shifts of X-rays are measured by the DPC method; hence, DPC-CT reconstructs refraction indexes rather than attenuation coefficients. In this work, we propose an explicit filtering based low-dose differential phase reconstruction algorithm, which enables reconstruction from reduced scanning without artifacts. The algorithm adopts a differential algebraic reconstruction technique (DART) with the explicit filtering based sparse regularization rather than the commonly used total variation (TV) method. Both the numerical simulation and the biological sample experiment demonstrate the feasibility of the proposed algorithm. PMID:26089971

An explicit dissipation-preserving method for Riesz space-fractional nonlinear wave equations in multiple dimensions

NASA Astrophysics Data System (ADS)

Macías-Díaz, J. E.

2018-06-01

In this work, we investigate numerically a model governed by a multidimensional nonlinear wave equation with damping and fractional diffusion. The governing partial differential equation considers the presence of Riesz space-fractional derivatives of orders in (1, 2], and homogeneous Dirichlet boundary data are imposed on a closed and bounded spatial domain. The model under investigation possesses an energy function which is preserved in the undamped regime. In the damped case, we establish the property of energy dissipation of the model using arguments from functional analysis. Motivated by these results, we propose an explicit finite-difference discretization of our fractional model based on the use of fractional centered differences. Associated to our discrete model, we also propose discretizations of the energy quantities. We establish that the discrete energy is conserved in the undamped regime, and that it dissipates in the damped scenario. Among the most important numerical features of our scheme, we show that the method has a consistency of second order, that it is stable and that it has a quadratic order of convergence. Some one- and two-dimensional simulations are shown in this work to illustrate the fact that the technique is capable of preserving the discrete energy in the undamped regime. For the sake of convenience, we provide a Matlab implementation of our method for the one-dimensional scenario.

Dynamic earthquake rupture simulation on nonplanar faults embedded in 3D geometrically complex, heterogeneous Earth models

NASA Astrophysics Data System (ADS)

Duru, K.; Dunham, E. M.; Bydlon, S. A.; Radhakrishnan, H.

2014-12-01

Dynamic propagation of shear ruptures on a frictional interface is a useful idealization of a natural earthquake.The conditions relating slip rate and fault shear strength are often expressed as nonlinear friction laws.The corresponding initial boundary value problems are both numerically and computationally challenging.In addition, seismic waves generated by earthquake ruptures must be propagated, far away from fault zones, to seismic stations and remote areas.Therefore, reliable and efficient numerical simulations require both provably stable and high order accurate numerical methods.We present a numerical method for:a) enforcing nonlinear friction laws, in a consistent and provably stable manner, suitable for efficient explicit time integration;b) dynamic propagation of earthquake ruptures along rough faults; c) accurate propagation of seismic waves in heterogeneous media with free surface topography.We solve the first order form of the 3D elastic wave equation on a boundary-conforming curvilinear mesh, in terms of particle velocities and stresses that are collocated in space and time, using summation-by-parts finite differences in space. The finite difference stencils are 6th order accurate in the interior and 3rd order accurate close to the boundaries. Boundary and interface conditions are imposed weakly using penalties. By deriving semi-discrete energy estimates analogous to the continuous energy estimates we prove numerical stability. Time stepping is performed with a 4th order accurate explicit low storage Runge-Kutta scheme. We have performed extensive numerical experiments using a slip-weakening friction law on non-planar faults, including recent SCEC benchmark problems. We also show simulations on fractal faults revealing the complexity of rupture dynamics on rough faults. We are presently extending our method to rate-and-state friction laws and off-fault plasticity.

Numerical and Physical Modelling of Bubbly Flow Phenomena - Final Report to the Department of Energy

DOE Office of Scientific and Technical Information (OSTI.GOV)

Andrea Prosperetti

This report describes the main features of the results obtained in the course of this project. A new approach to the systematic development of closure relations for the averaged equations of disperse multiphase flow is outlined. The focus of the project is on spatially non-uniform systems and several aspects in which such systems differ from uniform ones are described. Then, the procedure used in deriving the closure relations is given and some explicit results shown. The report also contains a list of publications supported by this grant and a list of the persons involved in the work.

Non-axisymmetric local magnetostatic equilibrium

DOE PAGES

Candy, Jefferey M.; Belli, Emily A.

2015-03-24

In this study, we outline an approach to the problem of local equilibrium in non-axisymmetric configurations that adheres closely to Miller's original method for axisymmetric plasmas. Importantly, this method is novel in that it allows not only specification of 3D shape, but also explicit specification of the shear in the 3D shape. A spectrally-accurate method for solution of the resulting nonlinear partial differential equations is also developed. We verify the correctness of the spectral method, in the axisymmetric limit, through comparisons with an independent numerical solution. Some analytic results for the two-dimensional case are given, and the connection to Boozermore » coordinates is clarified.« less

Simulation of a severe convective storm using a numerical model with explicitly incorporated aerosols

NASA Astrophysics Data System (ADS)

Lompar, Miloš; Ćurić, Mladjen; Romanic, Djordje

2017-09-01

Despite an important role the aerosols play in all stages of cloud lifecycle, their representation in numerical weather prediction models is often rather crude. This paper investigates the effects the explicit versus implicit inclusion of aerosols in a microphysics parameterization scheme in Weather Research and Forecasting (WRF) - Advanced Research WRF (WRF-ARW) model has on cloud dynamics and microphysics. The testbed selected for this study is a severe mesoscale convective system with supercells that struck west and central parts of Serbia in the afternoon of July 21, 2014. Numerical products of two model runs, i.e. one with aerosols explicitly (WRF-AE) included and another with aerosols implicitly (WRF-AI) assumed, are compared against precipitation measurements from surface network of rain gauges, as well as against radar and satellite observations. The WRF-AE model accurately captured the transportation of dust from the north Africa over the Mediterranean and to the Balkan region. On smaller scales, both models displaced the locations of clouds situated above west and central Serbia towards southeast and under-predicted the maximum values of composite radar reflectivity. Similar to satellite images, WRF-AE shows the mesoscale convective system as a merged cluster of cumulonimbus clouds. Both models over-predicted the precipitation amounts; WRF-AE over-predictions are particularly pronounced in the zones of light rain, while WRF-AI gave larger outliers. Unlike WRF-AI, the WRF-AE approach enables the modelling of time evolution and influx of aerosols into the cloud which could be of practical importance in weather forecasting and weather modification. Several likely causes for discrepancies between models and observations are discussed and prospects for further research in this field are outlined.

Asynchronous collision integrators: Explicit treatment of unilateral contact with friction and nodal restraints

PubMed Central

Wolff, Sebastian; Bucher, Christian

2013-01-01

This article presents asynchronous collision integrators and a simple asynchronous method treating nodal restraints. Asynchronous discretizations allow individual time step sizes for each spatial region, improving the efficiency of explicit time stepping for finite element meshes with heterogeneous element sizes. The article first introduces asynchronous variational integration being expressed by drift and kick operators. Linear nodal restraint conditions are solved by a simple projection of the forces that is shown to be equivalent to RATTLE. Unilateral contact is solved by an asynchronous variant of decomposition contact response. Therein, velocities are modified avoiding penetrations. Although decomposition contact response is solving a large system of linear equations (being critical for the numerical efficiency of explicit time stepping schemes) and is needing special treatment regarding overconstraint and linear dependency of the contact constraints (for example from double-sided node-to-surface contact or self-contact), the asynchronous strategy handles these situations efficiently and robust. Only a single constraint involving a very small number of degrees of freedom is considered at once leading to a very efficient solution. The treatment of friction is exemplified for the Coulomb model. Special care needs the contact of nodes that are subject to restraints. Together with the aforementioned projection for restraints, a novel efficient solution scheme can be presented. The collision integrator does not influence the critical time step. Hence, the time step can be chosen independently from the underlying time-stepping scheme. The time step may be fixed or time-adaptive. New demands on global collision detection are discussed exemplified by position codes and node-to-segment integration. Numerical examples illustrate convergence and efficiency of the new contact algorithm. Copyright © 2013 The Authors. International Journal for Numerical Methods in Engineering published by John Wiley & Sons, Ltd. PMID:23970806

Diabatic forcing and initialization with assimilation of cloud and rain water in a forecast model: Methodology

NASA Technical Reports Server (NTRS)

Raymond, William H.; Olson, William S.; Callan, Geary

1990-01-01

The focus of this part of the investigation is to find one or more general modeling techniques that will help reduce the time taken by numerical forecast models to initiate or spin-up precipitation processes and enhance storm intensity. If the conventional data base could explain the atmospheric mesoscale flow in detail, then much of our problem would be eliminated. But the data base is primarily synoptic scale, requiring that a solution must be sought either in nonconventional data, in methods to initialize mesoscale circulations, or in ways of retaining between forecasts the model generated mesoscale dynamics and precipitation fields. All three methods are investigated. The initialization and assimilation of explicit cloud and rainwater quantities computed from conservation equations in a mesoscale regional model are examined. The physical processes include condensation, evaporation, autoconversion, accretion, and the removal of rainwater by fallout. The question of how to initialize the explicit liquid water calculations in numerical models and how to retain information about precipitation processes during the 4-D assimilation cycle are important issues that are addressed. The explicit cloud calculations were purposely kept simple so that different initialization techniques can be easily and economically tested. Precipitation spin-up processes associated with three different types of weather phenomena are examined. Our findings show that diabatic initialization, or diabatic initialization in combination with a new diabatic forcing procedure, work effectively to enhance the spin-up of precipitation in a mesoscale numerical weather prediction forecast. Also, the retention of cloud and rain water during the analysis phase of the 4-D data assimilation procedure is shown to be valuable. Without detailed observations, the vertical placement of the diabatic heating remains a critical problem.

Eddy-resolving 1/10° model of the World Ocean

NASA Astrophysics Data System (ADS)

Ibrayev, R. A.; Khabeev, R. N.; Ushakov, K. V.

2012-02-01

The first results on simulating the intra-annual variability of the World Ocean circulation by use of the eddy-resolving model are considered. For this purpose, a model of the World Ocean with a 1/10° horizontal resolution and 49 vertical levels was developed (a 1/10 × 1/10 × 49 model of the World Ocean). This model is based on the traditional system of three-dimensional equations of the large-scale dynamics of the ocean and boundary conditions with an explicit allowance for water fluxes on the free surface of the ocean. The equations are written in the tripolar coordinate system. The numerical method is based on the separation of the barotropic and baroclinic components of the solution. Discretization in time is implemented using explicit schemes allowing effective parallelization for a large number of processors. The model uses the sub-models of the boundary layer of the atmosphere and the submodel of sea-ice thermodynamics. The model of the World Ocean was developed at the Institute of Numerical Mathematics of the Russian Academy of Sciences (INM RAS) and the P.P. Shirshov Institute of Oceanogy (IO RAS). The formulation of the problem of simulating the intra-annual variability of thermohydrodynamic processes of the World Ocean and the parameterizations that were used are considered. In the numerical experiment, the temporal evolution of the atmospheric effect is determined by the normal annual cycle according to the conditions of the international Coordinated Ocean-Ice Reference Experiment (CORE-I). The calculation was carried out on a multiprocessor computer with distributed memory; 1601 computational cores were used. The presented analysis demonstrates that the obtained results are quite satisfactory when compared to the results that were obtained by other eddy-resolving models of the global ocean. The analysis of the model solution is, to a larger extent, of a descriptive character. A detailed analysis of the results is to be presented in following works. This experiment is a significant first step in developing the eddy-resolving model of the World Ocean.

Model potentials for main group elements Li through Rn

NASA Astrophysics Data System (ADS)

Sakai, Yoshiko; Miyoshi, Eisaku; Klobukowski, Mariusz; Huzinaga, Sigeru

1997-05-01

Model potential (MP) parameters and valence basis sets were systematically determined for the main group elements Li through Rn. For alkali and alkaline-earth metal atoms, the outermost core (n-1)p electrons were treated explicitly together with the ns valence electrons. For the remaining atoms, only the valence ns and np electrons were treated explicitly. The major relativistic effects at the level of Cowan and Griffin's quasi-relativistic Hartree-Fock method (QRHF) were incorporated in the MPs for all atoms heavier than Kr. The valence orbitals thus obtained have inner nodal structure. The reliability of the MP method was tested in calculations for X-, X, and X+ (X=Br, I, and At) at the SCF level and the results were compared with the corresponding values given by the numerical HF (or QRHF) calculations. Calculations that include electron correlation were done for X-, X, and X+ (X=Cl and Br) at the SDCI level and for As2 at the CASSCF and MRSDCI levels. These results were compared with those of all-electron (AE) calculations using the well-tempered basis sets. Close agreement between the MP and AE results was obtained at all levels of the treatment.

Interfacial Ordering and Accompanying Divergent Capacitance at Ionic Liquid-Metal Interfaces.

PubMed

Limmer, David T

2015-12-18

A theory is constructed for dense ionic solutions near charged planar walls that is valid for strong interionic correlations. This theory predicts a fluctuation-induced, first-order transition and spontaneous charge density ordering at the interface, in the presence of an otherwise disordered bulk solution. The surface ordering is driven by applied voltage and results in an anomalous differential capacitance, in agreement with recent simulation results and consistent with experimental observations of a wide array of systems. Explicit forms for the charge density profile and capacitance are given. The theory is compared with numerical results for the charge frustrated Ising model, which is also found to exhibit a voltage driven first-order transition.

Interfacial Ordering and Accompanying Divergent Capacitance at Ionic Liquid-Metal Interfaces

NASA Astrophysics Data System (ADS)

Limmer, David T.

2015-12-01

A theory is constructed for dense ionic solutions near charged planar walls that is valid for strong interionic correlations. This theory predicts a fluctuation-induced, first-order transition and spontaneous charge density ordering at the interface, in the presence of an otherwise disordered bulk solution. The surface ordering is driven by applied voltage and results in an anomalous differential capacitance, in agreement with recent simulation results and consistent with experimental observations of a wide array of systems. Explicit forms for the charge density profile and capacitance are given. The theory is compared with numerical results for the charge frustrated Ising model, which is also found to exhibit a voltage driven first-order transition.

TTLEM - an implicit-explicit (IMEX) scheme for modelling landscape evolution in MATLAB

NASA Astrophysics Data System (ADS)

Campforts, Benjamin; Schwanghart, Wolfgang

2016-04-01

Landscape evolution models (LEM) are essential to unravel interdependent earth surface processes. They are proven very useful to bridge several temporal and spatial timescales and have been successfully used to integrate existing empirical datasets. There is a growing consensus that landscapes evolve at least as much in the horizontal as in the vertical direction urging for an efficient implementation of dynamic drainage networks. Here we present a spatially explicit LEM, which is based on the object-oriented function library TopoToolbox 2 (Schwanghart and Scherler, 2014). Similar to other LEMs, rivers are considered to be the main drivers for simulated landscape evolution as they transmit pulses of tectonic perturbations and set the base level of surrounding hillslopes. Highly performant graph algorithms facilitate efficient updates of the flow directions to account for planform changes in the river network and the calculation of flow-related terrain attributes. We implement the model using an implicit-explicit (IMEX) scheme, i.e. different integrators are used for different terms in the diffusion-incision equation. While linear diffusion is solved using an implicit scheme, we calculate incision explicitly. Contrary to previously published LEMS, however, river incision is solved using a total volume method which is total variation diminishing in order to prevent numerical diffusion when solving the stream power law (Campforts and Govers, 2015). We show that the use of this updated numerical scheme alters both landscape topography and catchment wide erosion rates at a geological time scale. Finally, the availability of a graphical user interface facilitates user interaction, making the tool very useful both for research and didactical purposes. References Campforts, B., Govers, G., 2015. Keeping the edge: A numerical method that avoids knickpoint smearing when solving the stream power law. J. Geophys. Res. Earth Surf. 120, 1189-1205. doi:10.1002/2014JF003376 Schwanghart, W., Scherler, D., 2014. TopoToolbox 2 - MATLAB-based software for topographic analysis and modeling in Earth surface sciences. Earth Surf. Dyn. 2, 1-7. doi:10.5194/esurf-2-1-2014

Modeling of fatigue crack induced nonlinear ultrasonics using a highly parallelized explicit local interaction simulation approach

NASA Astrophysics Data System (ADS)

Shen, Yanfeng; Cesnik, Carlos E. S.

2016-04-01

This paper presents a parallelized modeling technique for the efficient simulation of nonlinear ultrasonics introduced by the wave interaction with fatigue cracks. The elastodynamic wave equations with contact effects are formulated using an explicit Local Interaction Simulation Approach (LISA). The LISA formulation is extended to capture the contact-impact phenomena during the wave damage interaction based on the penalty method. A Coulomb friction model is integrated into the computation procedure to capture the stick-slip contact shear motion. The LISA procedure is coded using the Compute Unified Device Architecture (CUDA), which enables the highly parallelized supercomputing on powerful graphic cards. Both the explicit contact formulation and the parallel feature facilitates LISA's superb computational efficiency over the conventional finite element method (FEM). The theoretical formulations based on the penalty method is introduced and a guideline for the proper choice of the contact stiffness is given. The convergence behavior of the solution under various contact stiffness values is examined. A numerical benchmark problem is used to investigate the new LISA formulation and results are compared with a conventional contact finite element solution. Various nonlinear ultrasonic phenomena are successfully captured using this contact LISA formulation, including the generation of nonlinear higher harmonic responses. Nonlinear mode conversion of guided waves at fatigue cracks is also studied.

A Single-System Model Predicts Recognition Memory and Repetition Priming in Amnesia

PubMed Central

Kessels, Roy P.C.; Wester, Arie J.; Shanks, David R.

2014-01-01

We challenge the claim that there are distinct neural systems for explicit and implicit memory by demonstrating that a formal single-system model predicts the pattern of recognition memory (explicit) and repetition priming (implicit) in amnesia. In the current investigation, human participants with amnesia categorized pictures of objects at study and then, at test, identified fragmented versions of studied (old) and nonstudied (new) objects (providing a measure of priming), and made a recognition memory judgment (old vs new) for each object. Numerous results in the amnesic patients were predicted in advance by the single-system model, as follows: (1) deficits in recognition memory and priming were evident relative to a control group; (2) items judged as old were identified at greater levels of fragmentation than items judged new, regardless of whether the items were actually old or new; and (3) the magnitude of the priming effect (the identification advantage for old vs new items) overall was greater than that of items judged new. Model evidence measures also favored the single-system model over two formal multiple-systems models. The findings support the single-system model, which explains the pattern of recognition and priming in amnesia primarily as a reduction in the strength of a single dimension of memory strength, rather than a selective explicit memory system deficit. PMID:25122896

Dynamical Approach Study of Spurious Steady-State Numerical Solutions of Nonlinear Differential Equations. Part 2; Global Asymptotic Behavior of Time Discretizations

NASA Technical Reports Server (NTRS)

Yee, H. C.; Sweby, P. K.

1995-01-01

The global asymptotic nonlinear behavior of 11 explicit and implicit time discretizations for four 2 x 2 systems of first-order autonomous nonlinear ordinary differential equations (ODEs) is analyzed. The objectives are to gain a basic understanding of the difference in the dynamics of numerics between the scalars and systems of nonlinear autonomous ODEs and to set a baseline global asymptotic solution behavior of these schemes for practical computations in computational fluid dynamics. We show how 'numerical' basins of attraction can complement the bifurcation diagrams in gaining more detailed global asymptotic behavior of time discretizations for nonlinear differential equations (DEs). We show how in the presence of spurious asymptotes the basins of the true stable steady states can be segmented by the basins of the spurious stable and unstable asymptotes. One major consequence of this phenomenon which is not commonly known is that this spurious behavior can result in a dramatic distortion and, in most cases, a dramatic shrinkage and segmentation of the basin of attraction of the true solution for finite time steps. Such distortion, shrinkage and segmentation of the numerical basins of attraction will occur regardless of the stability of the spurious asymptotes, and will occur for unconditionally stable implicit linear multistep methods. In other words, for the same (common) steady-state solution the associated basin of attraction of the DE might be very different from the discretized counterparts and the numerical basin of attraction can be very different from numerical method to numerical method. The results can be used as an explanation for possible causes of error, and slow convergence and nonconvergence of steady-state numerical solutions when using the time-dependent approach for nonlinear hyperbolic or parabolic PDEs.

Dynamical Approach Study of Spurious Steady-State Numerical Solutions of Nonlinear Differential Equations. 2; Global Asymptotic Behavior of Time Discretizations; 2. Global Asymptotic Behavior of time Discretizations

NASA Technical Reports Server (NTRS)

Yee, H. C.; Sweby, P. K.

1995-01-01

The global asymptotic nonlinear behavior of 1 1 explicit and implicit time discretizations for four 2 x 2 systems of first-order autonomous nonlinear ordinary differential equations (ODES) is analyzed. The objectives are to gain a basic understanding of the difference in the dynamics of numerics between the scalars and systems of nonlinear autonomous ODEs and to set a baseline global asymptotic solution behavior of these schemes for practical computations in computational fluid dynamics. We show how 'numerical' basins of attraction can complement the bifurcation diagrams in gaining more detailed global asymptotic behavior of time discretizations for nonlinear differential equations (DEs). We show how in the presence of spurious asymptotes the basins of the true stable steady states can be segmented by the basins of the spurious stable and unstable asymptotes. One major consequence of this phenomenon which is not commonly known is that this spurious behavior can result in a dramatic distortion and, in most cases, a dramatic shrinkage and segmentation of the basin of attraction of the true solution for finite time steps. Such distortion, shrinkage and segmentation of the numerical basins of attraction will occur regardless of the stability of the spurious asymptotes, and will occur for unconditionally stable implicit linear multistep methods. In other words, for the same (common) steady-state solution the associated basin of attraction of the DE might be very different from the discretized counterparts and the numerical basin of attraction can be very different from numerical method to numerical method. The results can be used as an explanation for possible causes of error, and slow convergence and nonconvergence of steady-state numerical solutions when using the time-dependent approach for nonlinear hyperbolic or parabolic PDES.

Solution of multi-center molecular integrals of Slater-type orbitals

NASA Technical Reports Server (NTRS)

Tai, H.

1989-01-01

The troublesome multi-center molecular integrals of Slater-type orbitals (STO) in molecular physics calculations can be evaluated by using the Fourier transform and proper coupling of the two center exchange integrals. A numerical integration procedure is then readily rendered to the final expression in which the integrand consists of well known special functions of arguments containing the geometrical arrangement of the nuclear centers and the exponents of the atomic orbitals. A practical procedure was devised for the calculation of a general multi-center molecular integrals coupling arbitrary Slater-type orbitals. Symmetry relations and asymptotic conditions are discussed. Explicit expressions of three-center one-electron nuclear-attraction integrals and four-center two-electron repulsion integrals for STO of principal quantum number n=2 are listed. A few numerical results are given for the purpose of comparison.

Contact stresses in pin-loaded orthotropic plates

NASA Technical Reports Server (NTRS)

Hyer, M. W.; Klang, E. C.

1984-01-01

The effects of pin elasticity, friction, and clearance on the stresses near the hole in a pin-loaded orthotropic plate are described. The problem is modeled as a contact elasticity problem using complex variable theory, the pin and the plate being two elastic bodies interacting through contact. This modeling is in contrast to previous works which assumed that the pin is rigid or that it exerts a known cosinusoidal radial traction on the hole boundary. Neither of these approaches explicitly involves a pin. A collocation procedure and iteration were used to obtain numerical results for a variety of plate and pin elastic properties and various levels of friction and clearance. Collocation was used to enforce the boundary and iteration was used to find the contact and no-slip regions on the boundary. Details of the numerical scheme are discussed.

Stochastic Formal Correctness of Numerical Algorithms

NASA Technical Reports Server (NTRS)

Daumas, Marc; Lester, David; Martin-Dorel, Erik; Truffert, Annick

2009-01-01

We provide a framework to bound the probability that accumulated errors were never above a given threshold on numerical algorithms. Such algorithms are used for example in aircraft and nuclear power plants. This report contains simple formulas based on Levy's and Markov's inequalities and it presents a formal theory of random variables with a special focus on producing concrete results. We selected four very common applications that fit in our framework and cover the common practices of systems that evolve for a long time. We compute the number of bits that remain continuously significant in the first two applications with a probability of failure around one out of a billion, where worst case analysis considers that no significant bit remains. We are using PVS as such formal tools force explicit statement of all hypotheses and prevent incorrect uses of theorems.

Analytical drafting curves provide exact equations for plotted data

NASA Technical Reports Server (NTRS)

Stewart, R. B.

1967-01-01

Analytical drafting curves provide explicit mathematical expressions for any numerical data that appears in the form of graphical plots. The curves each have a reference coordinate axis system indicated on the curve as well as the mathematical equation from which the curve was generated.

Optimal time-domain technique for pulse width modulation in power electronics

NASA Astrophysics Data System (ADS)

Mayergoyz, I.; Tyagi, S.

2018-05-01

Optimal time-domain technique for pulse width modulation is presented. It is based on exact and explicit analytical solutions for inverter circuits, obtained for any sequence of input voltage rectangular pulses. Two optimal criteria are discussed and illustrated by numerical examples.

Ecosystem services

Treesearch

Trista Patterson

2014-01-01

Since its inception, the ecosystem service approach has stimulated interest from numerous planning, management, and partnership perspectives. To date, however, research that quantifies ecosystem services in the study area (in the form of explicit ecosystem service studies) has been limited. This chapter reviews and synthesizes the concept of ecosystem services,...

Long-time stability effects of quadrature and artificial viscosity on nodal discontinuous Galerkin methods for gas dynamics

NASA Astrophysics Data System (ADS)

Durant, Bradford; Hackl, Jason; Balachandar, Sivaramakrishnan

2017-11-01

Nodal discontinuous Galerkin schemes present an attractive approach to robust high-order solution of the equations of fluid mechanics, but remain accompanied by subtle challenges in their consistent stabilization. The effect of quadrature choices (full mass matrix vs spectral elements), over-integration to manage aliasing errors, and explicit artificial viscosity on the numerical solution of a steady homentropic vortex are assessed over a wide range of resolutions and polynomial orders using quadrilateral elements. In both stagnant and advected vortices in periodic and non-periodic domains the need arises for explicit stabilization beyond the numerical surface fluxes of discontinuous Galerkin spectral elements. Artificial viscosity via the entropy viscosity method is assessed as a stabilizing mechanism. It is shown that the regularity of the artificial viscosity field is essential to its use for long-time stabilization of small-scale features in nodal discontinuous Galerkin solutions of the Euler equations of gas dynamics. Supported by the Department of Energy Predictive Science Academic Alliance Program Contract DE-NA0002378.

Explicitly modelled deep-time tidal dissipation and its implication for Lunar history

NASA Astrophysics Data System (ADS)

Green, J. A. M.; Huber, M.; Waltham, D.; Buzan, J.; Wells, M.

2017-03-01

Dissipation of tidal energy causes the Moon to recede from the Earth. The currently measured rate of recession implies that the age of the Lunar orbit is 1500 My old, but the Moon is known to be 4500 My old. Consequently, it has been proposed that tidal energy dissipation was weaker in the Earth's past, but explicit numerical calculations are missing for such long time intervals. Here, for the first time, numerical tidal model simulations linked to climate model output are conducted for a range of paleogeographic configurations over the last 252 My. We find that the present is a poor guide to the past in terms of tidal dissipation: the total dissipation rates for most of the past 252 My were far below present levels. This allows us to quantify the reduced tidal dissipation rates over the most resent fraction of lunar history, and the lower dissipation allows refinement of orbitally-derived age models by inserting a complete additional precession cycle.

A class of high resolution explicit and implicit shock-capturing methods

NASA Technical Reports Server (NTRS)

Yee, H. C.

1989-01-01

An attempt is made to give a unified and generalized formulation of a class of high resolution, explicit and implicit shock capturing methods, and to illustrate their versatility in various steady and unsteady complex shock wave computations. Included is a systematic review of the basic design principle of the various related numerical methods. Special emphasis is on the construction of the basis nonlinear, spatially second and third order schemes for nonlinear scalar hyperbolic conservation laws and the methods of extending these nonlinear scalar schemes to nonlinear systems via the approximate Riemann solvers and the flux vector splitting approaches. Generalization of these methods to efficiently include equilibrium real gases and large systems of nonequilibrium flows are discussed. Some issues concerning the applicability of these methods that were designed for homogeneous hyperbolic conservation laws to problems containing stiff source terms and shock waves are also included. The performance of some of these schemes is illustrated by numerical examples for 1-, 2- and 3-dimensional gas dynamics problems.

Scattering matrices of Lamb waves at irregular surface and void defects.

PubMed

Feng, Feilong; Shen, Jianzhong; Lin, Shuyu

2012-08-01

Time-harmonic solution of Lamb wave scattering in a plane-strain waveguide with irregular thickness is investigated based on stair-step discretization and stepwise mode matching. The transfer relations of the transmission matrices and reflection matrices are derived in both directions of the waveguide. With these, an explicit expression of the scattering matrix is derived. When the scattering region of an inner irregular defect is geometrically divided into several parts composed of sub-waveguides with variable thicknesses and void regions with vertical free edges corresponding to the plate surfaces, the scattering matrix of the whole region could then be derived by modal matching along the artificial boundaries, as explicit functions of all the scattering matrices of the sub-waveguides and reflection matrices of the free edges. The effectiveness of the formulation is examined by numerical examples; the calculated scattering coefficients are in good accordance with those obtained from numerical simulation models. Copyright © 2012 Elsevier B.V. All rights reserved.

Two dimensional numerical prediction of deflagration-to-detonation transition in porous energetic materials.

PubMed

Narin, B; Ozyörük, Y; Ulas, A

2014-05-30

This paper describes a two-dimensional code developed for analyzing two-phase deflagration-to-detonation transition (DDT) phenomenon in granular, energetic, solid, explosive ingredients. The two-dimensional model is constructed in full two-phase, and based on a highly coupled system of partial differential equations involving basic flow conservation equations and some constitutive relations borrowed from some one-dimensional studies that appeared in open literature. The whole system is solved using an optimized high-order accurate, explicit, central-difference scheme with selective-filtering/shock capturing (SF-SC) technique, to augment central-diffencing and prevent excessive dispersion. The sources of the equations describing particle-gas interactions in terms of momentum and energy transfers make the equation system quite stiff, and hence its explicit integration difficult. To ease the difficulties, a time-split approach is used allowing higher time steps. In the paper, the physical model for the sources of the equation system is given for a typical explosive, and several numerical calculations are carried out to assess the developed code. Microscale intergranular and/or intragranular effects including pore collapse, sublimation, pyrolysis, etc. are not taken into account for ignition and growth, and a basic temperature switch is applied in calculations to control ignition in the explosive domain. Results for one-dimensional DDT phenomenon are in good agreement with experimental and computational results available in literature. A typical shaped-charge wave-shaper case study is also performed to test the two-dimensional features of the code and it is observed that results are in good agreement with those of commercial software. Copyright © 2014 Elsevier B.V. All rights reserved.

Invariant Tori in the Secular Motions of the Three-body Planetary Systems

NASA Astrophysics Data System (ADS)

Locatelli, Ugo; Giorgilli, Antonio

We consider the problem of the applicability of KAM theorem to a realistic problem of three bodies. In the framework of the averaged dynamics over the fast angles for the Sun-Jupiter-Saturn system we can prove the perpetual stability of the orbit. The proof is based on semi-numerical algorithms requiring both explicit algebraic manipulations of series and analytical estimates. The proof is made rigorous by using interval arithmetics in order to control the numerical errors.

Comparison of Implicit Collocation Methods for the Heat Equation

NASA Technical Reports Server (NTRS)

Kouatchou, Jules; Jezequel, Fabienne; Zukor, Dorothy (Technical Monitor)

2001-01-01

We combine a high-order compact finite difference scheme to approximate spatial derivatives arid collocation techniques for the time component to numerically solve the two dimensional heat equation. We use two approaches to implement the collocation methods. The first one is based on an explicit computation of the coefficients of polynomials and the second one relies on differential quadrature. We compare them by studying their merits and analyzing their numerical performance. All our computations, based on parallel algorithms, are carried out on the CRAY SV1.

Chimpanzees and the mathematics of battle.

PubMed

Wilson, Michael L; Britton, Nicholas F; Franks, Nigel R

2002-06-07

Recent experiments have demonstrated the importance of numerical assessment in animal contests. Nevertheless, few attempts have been made to model explicitly the relationship between the relative number of combatants on each side and the costs and benefits of entering a contest. One framework that may be especially suitable for making such explicit predictions is Lanchester's theory of combat, which has proved useful for understanding combat strategies in humans and several species of ants. We show, with data from a recent series of playback experiments, that a model derived from Lanchester's 'square law' predicts willingness to enter intergroup contests in wild chimpanzees (Pan troglodytes). Furthermore, the model predicts that, in contests with multiple individuals on each side, chimpanzees in this population should be willing to enter a contest only if they outnumber the opposing side by a factor of 1.5. We evaluate these results for intergroup encounters in chimpanzees and also discuss potential applications of Lanchester's square and linear laws for understanding combat strategies in other species.

A spatially explicit model for an Allee effect: why wolves recolonize so slowly in Greater Yellowstone.

PubMed

Hurford, Amy; Hebblewhite, Mark; Lewis, Mark A

2006-11-01

A reduced probability of finding mates at low densities is a frequently hypothesized mechanism for a component Allee effect. At low densities dispersers are less likely to find mates and establish new breeding units. However, many mathematical models for an Allee effect do not make a distinction between breeding group establishment and subsequent population growth. Our objective is to derive a spatially explicit mathematical model, where dispersers have a reduced probability of finding mates at low densities, and parameterize the model for wolf recolonization in the Greater Yellowstone Ecosystem (GYE). In this model, only the probability of establishing new breeding units is influenced by the reduced probability of finding mates at low densities. We analytically and numerically solve the model to determine the effect of a decreased probability in finding mates at low densities on population spread rate and density. Our results suggest that a reduced probability of finding mates at low densities may slow recolonization rate.

Explicit methods in extended phase space for inseparable Hamiltonian problems

NASA Astrophysics Data System (ADS)

Pihajoki, Pauli

2015-03-01

We present a method for explicit leapfrog integration of inseparable Hamiltonian systems by means of an extended phase space. A suitably defined new Hamiltonian on the extended phase space leads to equations of motion that can be numerically integrated by standard symplectic leapfrog (splitting) methods. When the leapfrog is combined with coordinate mixing transformations, the resulting algorithm shows good long term stability and error behaviour. We extend the method to non-Hamiltonian problems as well, and investigate optimal methods of projecting the extended phase space back to original dimension. Finally, we apply the methods to a Hamiltonian problem of geodesics in a curved space, and a non-Hamiltonian problem of a forced non-linear oscillator. We compare the performance of the methods to a general purpose differential equation solver LSODE, and the implicit midpoint method, a symplectic one-step method. We find the extended phase space methods to compare favorably to both for the Hamiltonian problem, and to the implicit midpoint method in the case of the non-linear oscillator.

Explicit accounting of electronic effects on the Hugoniot of porous materials

DOE Office of Scientific and Technical Information (OSTI.GOV)

Nayak, Bishnupriya; Menon, S. V. G., E-mail: menon.svg98@gmail.com

2016-03-28

A generalized enthalpy based equation of state, which includes thermal electron excitations explicitly, is formulated from simple considerations. Its application to obtain Hugoniot of materials needs simultaneous evaluation of pressure-volume curve and temperature, the latter requiring solution of a differential equation. The errors involved in two recent papers [Huayun et al., J. Appl. Phys. 92, 5917 (2002); 92, 5924 (2002)], which employed this approach, are brought out and discussed. In addition to developing the correct set of equations, the present work also provides a numerical method to implement this approach. Constant pressure specific heat of ions and electrons and ionicmore » enthalpy parameter, needed for applications, are calculated using a three component equation of state. The method is applied to porous Cu with different initial porosities. Comparison of results with experimental data shows good agreement. It is found that temperatures along the Hugoniot of porous materials are significantly modified due to electronic effects.« less

Multigrid Acceleration of Time-Accurate DNS of Compressible Turbulent Flow

NASA Technical Reports Server (NTRS)

Broeze, Jan; Geurts, Bernard; Kuerten, Hans; Streng, Martin

1996-01-01

An efficient scheme for the direct numerical simulation of 3D transitional and developed turbulent flow is presented. Explicit and implicit time integration schemes for the compressible Navier-Stokes equations are compared. The nonlinear system resulting from the implicit time discretization is solved with an iterative method and accelerated by the application of a multigrid technique. Since we use central spatial discretizations and no artificial dissipation is added to the equations, the smoothing method is less effective than in the more traditional use of multigrid in steady-state calculations. Therefore, a special prolongation method is needed in order to obtain an effective multigrid method. This simulation scheme was studied in detail for compressible flow over a flat plate. In the laminar regime and in the first stages of turbulent flow the implicit method provides a speed-up of a factor 2 relative to the explicit method on a relatively coarse grid. At increased resolution this speed-up is enhanced correspondingly.

Chimpanzees and the mathematics of battle.

PubMed Central

Wilson, Michael L; Britton, Nicholas F; Franks, Nigel R

2002-01-01

Recent experiments have demonstrated the importance of numerical assessment in animal contests. Nevertheless, few attempts have been made to model explicitly the relationship between the relative number of combatants on each side and the costs and benefits of entering a contest. One framework that may be especially suitable for making such explicit predictions is Lanchester's theory of combat, which has proved useful for understanding combat strategies in humans and several species of ants. We show, with data from a recent series of playback experiments, that a model derived from Lanchester's 'square law' predicts willingness to enter intergroup contests in wild chimpanzees (Pan troglodytes). Furthermore, the model predicts that, in contests with multiple individuals on each side, chimpanzees in this population should be willing to enter a contest only if they outnumber the opposing side by a factor of 1.5. We evaluate these results for intergroup encounters in chimpanzees and also discuss potential applications of Lanchester's square and linear laws for understanding combat strategies in other species. PMID:12061952

Identification of internal properties of fibers and micro-swimmers

NASA Astrophysics Data System (ADS)

Plouraboue, Franck; Thiam, Ibrahima; Delmotte, Blaise; Climent, Eric; PSC Collaboration

2016-11-01

In this presentation we discuss the identifiability of constitutive parameters of passive or active micro-swimmers. We first present a general framework for describing fibers or micro-swimmers using a bead-model description. Using a kinematic constraint formulation to describe fibers, flagellum or cilia, we find explicit linear relationship between elastic constitutive parameters and generalised velocities from computing contact forces. This linear formulation then permits to address explicitly identifiability conditions and solve for parameter identification. We show that both active forcing and passive parameters are both identifiable independently but not simultaneously. We also provide unbiased estimators for elastic parameters as well as active ones in the presence of Langevin-like forcing with Gaussian noise using normal linear regression models and maximum likelihood method. These theoretical results are illustrated in various configurations of relaxed or actuated passives fibers, and active filament of known passive properties, showing the efficiency of the proposed approach for direct parameter identification. The convergence of the proposed estimators is successfully tested numerically.

Physical similarity or numerical representation counts in same-different, numerical comparison, physical comparison, and priming tasks?

PubMed

Zhang, Li; Xin, Ziqiang; Feng, Tingyong; Chen, Yinghe; Szűcs, Denes

2018-03-01

Recent studies have highlighted the fact that some tasks used to study symbolic number representations are confounded by judgments about physical similarity. Here, we investigated whether the contribution of physical similarity and numerical representation differed in the often-used symbolic same-different, numerical comparison, physical comparison, and priming tasks. Experiment 1 showed that subjective physical similarity was the best predictor of participants' performance in the same-different task, regardless of simultaneous or sequential presentation. Furthermore, the contribution of subjective physical similarity was larger in a simultaneous presentation than in a sequential presentation. Experiment 2 showed that only numerical representation was involved in numerical comparison. Experiment 3 showed that both subjective physical similarity and numerical representation contributed to participants' physical comparison performance. Finally, only numerical representation contributed to participants' performance in a priming task as revealed by Experiment 4. Taken together, the contribution of physical similarity and numerical representation depends on task demands. Performance primarily seems to rely on numerical properties in tasks that require explicit quantitative comparison judgments (physical or numerical), while physical stimulus properties exert an effect in the same-different task.

For numerical differentiation, dimensionality can be a blessing!

NASA Astrophysics Data System (ADS)

Anderssen, Robert S.; Hegland, Markus

Finite difference methods, such as the mid-point rule, have been applied successfully to the numerical solution of ordinary and partial differential equations. If such formulas are applied to observational data, in order to determine derivatives, the results can be disastrous. The reason for this is that measurement errors, and even rounding errors in computer approximations, are strongly amplified in the differentiation process, especially if small step-sizes are chosen and higher derivatives are required. A number of authors have examined the use of various forms of averaging which allows the stable computation of low order derivatives from observational data. The size of the averaging set acts like a regularization parameter and has to be chosen as a function of the grid size h. In this paper, it is initially shown how first (and higher) order single-variate numerical differentiation of higher dimensional observational data can be stabilized with a reduced loss of accuracy than occurs for the corresponding differentiation of one-dimensional data. The result is then extended to the multivariate differentiation of higher dimensional data. The nature of the trade-off between convergence and stability is explicitly characterized, and the complexity of various implementations is examined.

Numerical solution of a spatio-temporal gender-structured model for hantavirus infection in rodents.

PubMed

Bürger, Raimund; Chowell, Gerardo; Gavilán, Elvis; Mulet, Pep; Villada, Luis M

2018-02-01

In this article we describe the transmission dynamics of hantavirus in rodents using a spatio-temporal susceptible-exposed-infective-recovered (SEIR) compartmental model that distinguishes between male and female subpopulations [L.J.S. Allen, R.K. McCormack and C.B. Jonsson, Bull. Math. Biol. 68 (2006), 511--524]. Both subpopulations are assumed to differ in their movement with respect to local variations in the densities of their own and the opposite gender group. Three alternative models for the movement of the male individuals are examined. In some cases the movement is not only directed by the gradient of a density (as in the standard diffusive case), but also by a non-local convolution of density values as proposed, in another context, in [R.M. Colombo and E. Rossi, Commun. Math. Sci., 13 (2015), 369--400]. An efficient numerical method for the resulting convection-diffusion-reaction system of partial differential equations is proposed. This method involves techniques of weighted essentially non-oscillatory (WENO) reconstructions in combination with implicit-explicit Runge-Kutta (IMEX-RK) methods for time stepping. The numerical results demonstrate significant differences in the spatio-temporal behavior predicted by the different models, which suggest future research directions.

Semiclassical evaluation of quantum fidelity

NASA Astrophysics Data System (ADS)

Vanicek, Jiri

2004-03-01

We present a numerically feasible semiclassical method to evaluate quantum fidelity (Loschmidt echo) in a classically chaotic system. It was thought that such evaluation would be intractable, but instead we show that a uniform semiclassical expression not only is tractable but it gives remarkably accurate numerical results for the standard map in both the Fermi-golden-rule and Lyapunov regimes. Because it allows a Monte-Carlo evaluation, this uniform expression is accurate at times where there are 10^70 semiclassical contributions. Remarkably, the method also explicitly contains the ``building blocks'' of analytical theories of recent literature, and thus permits a direct test of approximations made by other authors in these regimes, rather than an a posteriori comparison with numerical results. We explain in more detail the extended validity of the classical perturbation approximation and thus provide a ``defense" of the linear response theory from the famous Van Kampen objection. We point out the potential use of our uniform expression in other areas because it gives a most direct link between the quantum Feynman propagator based on the path integral and the semiclassical Van Vleck propagator based on the sum over classical trajectories. Finally, we test the applicability of our method in integrable and mixed systems.

Perturbation theory and numerical modelling of weakly and moderately nonlinear incompressible Richtmyer-Meshkov instability

NASA Astrophysics Data System (ADS)

Herrmann, M.; Velikovich, A. L.; Abarzhi, S. I.

2014-10-01

A study of incompressible two-dimensional Richtmyer-Meshkov instability by means of high-order Eulerian perturbation theory and numerical simulations is reported. Nonlinear corrections to Richtmyer's impulsive formula for the bubble and spike growth rates have been calculated analytically for arbitrary Atwood number and an explicit formula has been obtained for it in the Boussinesq limit. Conditions for early-time acceleration and deceleration of the bubble and the spike have been derived. In our simulations we have solved 2D unsteady Navier-Stokes equations for immiscible incompressible fluids using the finite volume fractional step flow solver NGA developed by, coupled to the level set based interface solver LIT,. The impact of small amounts of viscosity and surface tension on the RMI flow dynamics is studied numerically. Simulation results are compared to the theory to demonstrate successful code verification and highlight the influence of the theory's ideal inviscid flow assumption. Theoretical time histories of the interface curvature at the bubble and spike tip and the profiles of vertical and horizontal velocities have been favorably compared to simulation results, which converge to the theoretical predictions as the Reynolds and Weber numbers are increased. Work supported by the US DOE/NNSA.

A spectral hybridizable discontinuous Galerkin method for elastic-acoustic wave propagation

NASA Astrophysics Data System (ADS)

Terrana, S.; Vilotte, J. P.; Guillot, L.

2018-04-01

We introduce a time-domain, high-order in space, hybridizable discontinuous Galerkin (DG) spectral element method (HDG-SEM) for wave equations in coupled elastic-acoustic media. The method is based on a first-order hyperbolic velocity-strain formulation of the wave equations written in conservative form. This method follows the HDG approach by introducing a hybrid unknown, which is the approximation of the velocity on the elements boundaries, as the only globally (i.e. interelement) coupled degrees of freedom. In this paper, we first present a hybridized formulation of the exact Riemann solver at the element boundaries, taking into account elastic-elastic, acoustic-acoustic and elastic-acoustic interfaces. We then use this Riemann solver to derive an explicit construction of the HDG stabilization function τ for all the above-mentioned interfaces. We thus obtain an HDG scheme for coupled elastic-acoustic problems. This scheme is then discretized in space on quadrangular/hexahedral meshes using arbitrary high-order polynomial basis for both volumetric and hybrid fields, using an approach similar to the spectral element methods. This leads to a semi-discrete system of algebraic differential equations (ADEs), which thanks to the structure of the global conservativity condition can be reformulated easily as a classical system of first-order ordinary differential equations in time, allowing the use of classical explicit or implicit time integration schemes. When an explicit time scheme is used, the HDG method can be seen as a reformulation of a DG with upwind fluxes. The introduction of the velocity hybrid unknown leads to relatively simple computations at the element boundaries which, in turn, makes the HDG approach competitive with the DG-upwind methods. Extensive numerical results are provided to illustrate and assess the accuracy and convergence properties of this HDG-SEM. The approximate velocity is shown to converge with the optimal order of k + 1 in the L2-norm, when element polynomials of order k are used, and to exhibit the classical spectral convergence of SEM. Additional inexpensive local post-processing in both the elastic and the acoustic case allow to achieve higher convergence orders. The HDG scheme provides a natural framework for coupling classical, continuous Galerkin SEM with HDG-SEM in the same simulation, and it is shown numerically in this paper. As such, the proposed HDG-SEM can combine the efficiency of the continuous SEM with the flexibility of the HDG approaches. Finally, more complex numerical results, inspired from real geophysical applications, are presented to illustrate the capabilities of the method for wave propagation in heterogeneous elastic-acoustic media with complex geometries.

Remedying excessive numerical diapycnal mixing in a global 0.25° NEMO configuration

NASA Astrophysics Data System (ADS)

Megann, Alex; Nurser, George; Storkey, Dave

2016-04-01

If numerical ocean models are to simulate faithfully the upwelling branches of the global overturning circulation, they need to have a good representation of the diapycnal mixing processes which contribute to conversion of the bottom and deep waters produced in high latitudes into less dense watermasses. It is known that the default class of depth-coordinate ocean models such as NEMO and MOM5, as used in many state-of-the art coupled climate models and Earth System Models, have excessive numerical diapycnal mixing, resulting from irreversible advection across coordinate surfaces. The GO5.0 configuration of the NEMO ocean model, on an "eddy-permitting" 0.25° global grid, is used in the current UK GC1 and GC2 coupled models. Megann and Nurser (2016) have shown, using the isopycnal watermass analysis of Lee et al (2002), that spurious numerical mixing is substantially larger than the explicit mixing prescribed by the mixing scheme used by the model. It will be shown that increasing the biharmonic viscosity by a factor of three tends to suppress small-scale noise in the vertical velocity in the model. This significantly reduces the numerical mixing in GO5.0, and we shall show that it also leads to large-scale improvements in model biases.

Working memory deficits in developmental dyscalculia: The importance of serial order.

PubMed

Attout, Lucie; Majerus, Steve

2015-01-01

Although a number of studies suggests a link between working memory (WM) storage capacity of short-term memory and calculation abilities, the nature of verbal WM deficits in children with developmental dyscalculia (DD) remains poorly understood. We explored verbal WM capacity in DD by focusing on the distinction between memory for item information (the items to be retained) and memory for order information (the order of the items within a list). We hypothesized that WM for order could be specifically related to impaired numerical abilities given that recent studies suggest close interactions between the representation of order information in WM and ordinal numerical processing. We investigated item and order WM abilities as well as basic numerical processing abilities in 16 children with DD (age: 8-11 years) and 16 typically developing children matched on age, IQ, and reading abilities. The DD group performed significantly poorer than controls in the order WM condition but not in the item WM condition. In addition, the DD group performed significantly slower than the control group on a numerical order judgment task. The present results show significantly reduced serial order WM abilities in DD coupled with less efficient numerical ordinal processing abilities, reflecting more general difficulties in explicit processing of ordinal information.

Making Implicit Multivariable Calculus Representations Explicit: A Clinical Study

ERIC Educational Resources Information Center

McGee, Daniel; Moore-Russo, Deborah; Martinez-Planell, Rafael

2015-01-01

Reviewing numerous textbooks, we found that in both differential and integral calculus textbooks the authors commonly assume that: (i) students can generalize associations between representations in two dimensions to associations between representations of the same mathematical concept in three dimensions on their own; and (ii) explicit…

TIME SHARING WITH AN EXPLICIT PRIORITY QUEUING DISCIPLINE.

DTIC Science & Technology

exponentially distributed service times and an ordered priority queue. Each new arrival buys a position in this queue by offering a non-negative bribe to the...parameters is investigated through numerical examples. Finally, to maximize the expected revenue per unit time accruing from bribes , an optimization

Design and numerical evaluation of full-authority flight control systems for conventional and thruster-augmented helicopters employed in NOE operations

NASA Technical Reports Server (NTRS)

Perri, Todd A.; Mckillip, R. M., Jr.; Curtiss, H. C., Jr.

1987-01-01

The development and methodology is presented for development of full-authority implicit model-following and explicit model-following optimal controllers for use on helicopters operating in the Nap-of-the Earth (NOE) environment. Pole placement, input-output frequency response, and step input response were used to evaluate handling qualities performance. The pilot was equipped with velocity-command inputs. A mathematical/computational trajectory optimization method was employed to evaluate the ability of each controller to fly NOE maneuvers. The method determines the optimal swashplate and thruster input histories from the helicopter's dynamics and the prescribed geometry and desired flying qualities of the maneuver. Three maneuvers were investigated for both the implicit and explicit controllers with and without auxiliary propulsion installed: pop-up/dash/descent, bob-up at 40 knots, and glideslope. The explicit controller proved to be superior to the implicit controller in performance and ease of design.

How to aggregate health? Separability and the effect of framing.

PubMed

Turpcu, Adam; Bleichrodt, Han; Le, Quang A; Doctor, Jason N

2012-01-01

Unweighted summation or quality-adjusted life year (QALY) utilitarianism is the most common way to aggregate health benefits in a cost-effectiveness analysis. A key qualitative principle underlying QALY utilitarianism is separability: those individuals unaffected by a policy choice should not influence the policy choice. Separability also underlies several of the alternatives for QALY utilitarianism that have been proposed. To test separability and to test whether the support for separability is affected by the framing of the choice questions. In 2 experiments, 345 student subjects (162 in the first experiment, and 183 in the second experiment) were asked to select 1 of 2 possible treatments, with each treatment resulting in a different distribution of health across individuals. The only aspect that varied across choice questions was the state of the patients whose health was unaffected by the act of choosing a policy. In each experiment, we used 2 frames. In the implicit frame, it was implied but not plainly expressed what outcomes the treatments had in common. In the explicit frame, common outcomes of the 2 treatments were directly stated. The 2 experiments differed in the way the explicit frame was presented (verbal v. numerical). The support for separability was significantly greater in the explicit frame. The proportion of violations in the implicit frame was 44% in Experiment 1 and 31% in Experiment 2, while in the explicit frame, the proportion of violations was 28% in Experiment 1 and 8% in Experiment 2. Framing affected the support for separability, raising issues as to whether it is possible to achieve a canonical representation of social choices.

A study of the displacement of a Wankel rotary engine

NASA Astrophysics Data System (ADS)

Beard, J. E.; Pennock, G. R.

1993-03-01

The volumetric displacement of a Wankel rotary engine is a function of the trochoid ratio and the pin size ratio, assuming the engine has a unit depth and the number of lobes is specified. The mathematical expression which defines the displacement contains a function which can be evaluated directly and a normal elliptic integral of the second type which does not have an explicit solution. This paper focuses on the contribution of the elliptic integral to the total displacement of the engine. The influence of the elliptic integral is shown to account for as much as 20 percent of the total displacement, depending on the trochoid ratio and the pin size ratio. Two numerical integration techniques are compared in the paper, namely, the trapezoidal rule and Simpson's 1/3 rule. The bounds on the error, associated with each numerical method, are analyzed. The results indicate that the numerical method has a minimal effect on the accuracy of the calculated displacement for a practical number of integration steps. The paper also evaluates the influence of manufacturing tolerances on the calculated displacement and the actual displacement. Finally. a numerical example of the common three-lobed Wankel rotary engine is included for illustrative purposes.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Fleury, Pierre; Uzan, Jean-Philippe; Larena, Julien, E-mail: fleury@iap.fr, E-mail: j.larena@ru.ac.za, E-mail: uzan@iap.fr

On the scale of the light beams subtended by small sources, e.g. supernovae, matter cannot be accurately described as a fluid, which questions the applicability of standard cosmic lensing to those cases. In this article, we propose a new formalism to deal with small-scale lensing as a diffusion process: the Sachs and Jacobi equations governing the propagation of narrow light beams are treated as Langevin equations. We derive the associated Fokker-Planck-Kolmogorov equations, and use them to deduce general analytical results on the mean and dispersion of the angular distance. This formalism is applied to random Einstein-Straus Swiss-cheese models, allowing usmore » to: (1) show an explicit example of the involved calculations; (2) check the validity of the method against both ray-tracing simulations and direct numerical integration of the Langevin equation. As a byproduct, we obtain a post-Kantowski-Dyer-Roeder approximation, accounting for the effect of tidal distortions on the angular distance, in excellent agreement with numerical results. Besides, the dispersion of the angular distance is correctly reproduced in some regimes.« less

The theory of stochastic cosmological lensing

NASA Astrophysics Data System (ADS)

Fleury, Pierre; Larena, Julien; Uzan, Jean-Philippe

2015-11-01

On the scale of the light beams subtended by small sources, e.g. supernovae, matter cannot be accurately described as a fluid, which questions the applicability of standard cosmic lensing to those cases. In this article, we propose a new formalism to deal with small-scale lensing as a diffusion process: the Sachs and Jacobi equations governing the propagation of narrow light beams are treated as Langevin equations. We derive the associated Fokker-Planck-Kolmogorov equations, and use them to deduce general analytical results on the mean and dispersion of the angular distance. This formalism is applied to random Einstein-Straus Swiss-cheese models, allowing us to: (1) show an explicit example of the involved calculations; (2) check the validity of the method against both ray-tracing simulations and direct numerical integration of the Langevin equation. As a byproduct, we obtain a post-Kantowski-Dyer-Roeder approximation, accounting for the effect of tidal distortions on the angular distance, in excellent agreement with numerical results. Besides, the dispersion of the angular distance is correctly reproduced in some regimes.

Efficient solution of parabolic equations by Krylov approximation methods

NASA Technical Reports Server (NTRS)

Gallopoulos, E.; Saad, Y.

1990-01-01

Numerical techniques for solving parabolic equations by the method of lines is addressed. The main motivation for the proposed approach is the possibility of exploiting a high degree of parallelism in a simple manner. The basic idea of the method is to approximate the action of the evolution operator on a given state vector by means of a projection process onto a Krylov subspace. Thus, the resulting approximation consists of applying an evolution operator of a very small dimension to a known vector which is, in turn, computed accurately by exploiting well-known rational approximations to the exponential. Because the rational approximation is only applied to a small matrix, the only operations required with the original large matrix are matrix-by-vector multiplications, and as a result the algorithm can easily be parallelized and vectorized. Some relevant approximation and stability issues are discussed. We present some numerical experiments with the method and compare its performance with a few explicit and implicit algorithms.

Numerical analysis of moving contact line with contact angle hysteresis using feedback deceleration technique

NASA Astrophysics Data System (ADS)

Park, Jun Kwon; Kang, Kwan Hyoung

2012-04-01

Contact angle (CA) hysteresis is important in many natural and engineering wetting processes, but predicting it numerically is difficult. We developed an algorithm that considers CA hysteresis when analyzing the motion of the contact line (CL). This algorithm employs feedback control of CA which decelerates CL speed to make the CL stationary in the hysteretic range of CA, and one control coefficient should be heuristically determined depending on characteristic time of the simulated system. The algorithm requires embedding only a simple additional routine with little modification of a code which considers the dynamic CA. The method is non-iterative and explicit, and also has less computational load than other algorithms. For a drop hanging on a wire, the proposed algorithm accurately predicts the theoretical equilibrium CA. For the drop impacting on a dry surface, the results of the proposed algorithm agree well with experimental results including the intermittent occurrence of the pinning of CL. The proposed algorithm is as accurate as other algorithms, but faster.

ACCRETION OF SUPERSONIC WINDS ONTO BLACK HOLES IN 3D: STABILITY OF THE SHOCK CONE

DOE Office of Scientific and Technical Information (OSTI.GOV)

Gracia-Linares, M.; Guzmán, F. S.

Using numerical simulations we present the accretion of supersonic winds onto a rotating black hole in three dimensions. We study five representative directions of the wind with respect to the axis of rotation of the black hole and focus on the evolution and stability of the high-density shock cone that is formed during the process. We explore both the regime in which the shock cone is expected to be stable in order to confirm previous results obtained with two-dimensional simulations, and the regime in which the shock cone is expected to show a flip–flop (FF) type of instability. The methodsmore » used to attempt a triggering of the instability were (i) the accumulation of numerical errors and (ii) the explicit application of a perturbation on the velocity field after the shock cone was formed. The result is negative, that is, we did not find the FF instability within the parameter space we explored, including cases that are expected to be unstable.« less

Numerical investigation for formability of aluminum 6016 alloy under non-isothermal warm forming process

NASA Astrophysics Data System (ADS)

Hu, P.; Dai, M. H.; Ying, L.; Shi, D. Y.; Zhao, K. M.; Lu, J. D.

2013-05-01

The warm forming technology of aluminum alloy has attracted attention from worldwide automotive engineering sector in recent years, with which the complex geometry parts can be realized at elevated temperature. A non-isothermal warm forming process for the heat treatable aluminum can quickly carry out its application on traditional production line by adding a furnace to heat up the aluminum alloy sheet. The 6000 aluminum alloy was investigated by numerical simulation and experiment using the Nakajima test model in this paper. A modified Fields-Backofen model was introduced into numerical simulation process to describe the thermo-mechanical flow behavior of a 6000 series aluminum alloy. The experimental data was obtained by conducting thermal-mechanical uniaxial tensile experiment in temperatures range of 25˜400°C to guarantee the numerical simulation more accurate. The numerical simulation was implemented with LS_DYNA software in terms of coupled dynamic explicit method for investigating the effect of initial forming temperature and the Binder Holder Force (BHF), which are critical process parameters in non-isothermal warm forming. The results showed that the optimal initial forming temperature range was 300°C˜350°C. By means of conducting numerical simulation in deep drawing box model, the forming window of BHF and temperature around the optimal initial forming temperature (275°, 300° and 325°) are investigated, which can provide guidance to actual experiment.

Implicitly solving phase appearance and disappearance problems using two-fluid six-equation model

DOE PAGES

Zou, Ling; Zhao, Haihua; Zhang, Hongbin

2016-01-25

Phase appearance and disappearance issue presents serious numerical challenges in two-phase flow simulations using the two-fluid six-equation model. Numerical challenges arise from the singular equation system when one phase is absent, as well as from the discontinuity in the solution space when one phase appears or disappears. In this work, a high-resolution spatial discretization scheme on staggered grids and fully implicit methods were applied for the simulation of two-phase flow problems using the two-fluid six-equation model. A Jacobian-free Newton-Krylov (JFNK) method was used to solve the discretized nonlinear problem. An improved numerical treatment was proposed and proved to be effectivemore » to handle the numerical challenges. The treatment scheme is conceptually simple, easy to implement, and does not require explicit truncations on solutions, which is essential to conserve mass and energy. Various types of phase appearance and disappearance problems relevant to thermal-hydraulics analysis have been investigated, including a sedimentation problem, an oscillating manometer problem, a non-condensable gas injection problem, a single-phase flow with heat addition problem and a subcooled flow boiling problem. Successful simulations of these problems demonstrate the capability and robustness of the proposed numerical methods and numerical treatments. As a result, volume fraction of the absent phase can be calculated effectively as zero.« less

An analytical approach to the CMB polarization in a spatially closed background

NASA Astrophysics Data System (ADS)

Niazy, Pedram; Abbassi, Amir H.

2018-03-01

The scalar mode polarization of the cosmic microwave background is derived in a spatially closed universe from the Boltzmann equation using the line of sight integral method. The EE and TE multipole coefficients have been extracted analytically by considering some tolerable approximations such as considering the evolution of perturbation hydrodynamically and sudden transition from opacity to transparency at the time of last scattering. As the major advantage of analytic expressions, CEE,ℓS and CTE,ℓ explicitly show the dependencies on baryon density ΩB, matter density ΩM, curvature ΩK, primordial spectral index ns, primordial power spectrum amplitude As, Optical depth τreion, recombination width σt and recombination time tL. Using a realistic set of cosmological parameters taken from a fit to data from Planck, the closed universe EE and TE power spectrums in the scalar mode are compared with numerical results from the CAMB code and also latest observational data. The analytic results agree with the numerical ones on the big and moderate scales. The peak positions are in good agreement with the numerical result on these scales while the peak heights agree with that to within 20% due to the approximations have been considered for these derivations. Also, several interesting properties of CMB polarization are revealed by the analytic spectra.