Engineers at the University of California, San Diego, have developed a new family of methods to significantly increase the speed of time-resolved numerical simulations in computational grand challenge problems.
Such problems often arise from the high-resolution approximation of the partial differential equations governing complex flows of fluids or plasmas. The breakthrough could be applied to simulations that include millions or billions of variables, including turbulence simulations.
Modern computers are generally built from commodity hardware developed for serving and surfing the web. When applied to cutting-edge problems in scientific computing, computers built from such general-purpose hardware usually spend most of their time moving data around in memory, and the hardware dedicated to floating point computations (that is, the actual addition and multiplication of numbers) spends most of its time idle.
The small memory footprint of the new schemes developed at UC San Diego means that numerical problems of a given size will run much faster on a given computer, and that even larger numerical problems may be considered.
"Moving information around in memory is the bottleneck in almost all large-scale numerical simulations today," said Thomas Bewley, a mechanical engineering professor who leads the Flow Control Lab at the Jacobs School of Engineering at UC San Diego.
"The remarkable feature of the new family of schemes developed in this work is that they require significantly less memory in the computer for a given size simulation problem than existing high-order methods of the same class, while providing excellent numerical stability, accuracy, and computational efficiency."
Complex systems such as flows of fluids and plasmas generally evolve as a result of a combination of physical effects, such as diffusion and convection. Some of these effects are linear and incorporate many spatial derivatives (that is, they are characterized by a large range of characteristic time scales, and are thus referred to as "stiff"). These terms are best handled with "implicit" methods, which require the solution of many simple simultaneous equations using matrix algebra and iterative solvers.
Other effects are nonlinear and incorporate fewer spatial derivatives (that is, they are characterized by a smaller range of characteristic time scales, and are thus referred to as "nonstiff"). These terms are most easily handled with explicit methods, which treat the propagation of each equation independently. If the stiff terms are treated with explicit methods, a severe restriction arises on the timestep, which slows the simulation; if the nonstiff terms are treated with implicit methods, complex and computationally expensive iterative solvers must be used.
The new "implicit/explicit" or IMEX time marching schemes developed at UC San Diego thus marry together two algorithms for time-resolved simulations of the standard "Runge-Kutta" or RK form. The implicit algorithm is applied to the stiff terms of the problem, and the explicit algorithm is applied to the nonstiff terms of the problem. The two algorithms so joined are each endowed with good numerical properties, such as excellent stability and high accuracy, and, notably, maintain this high accuracy when working together in concert. The compatible pairs of simulation methods so developed are known as IMEXRK schemes.
"Searching for the right combination of the dozens of parameters that make these new IMEXRK algorithms work well was like finding a needle in a haystack, and required a tedious search over a very large parameter space, combined with the delicate codification of various numerical intuition to simplify the search. It took almost one year to complete," said Daniele Cavaglieri, a Ph.D. student and co-author of the paper.
Researchers describe the new methods in the January 2015 issue of the Journal of Computational Physics.
Low-storage implicit/explicit Runge-Kutta schemes for the simulation of stiff high-dimensional OED systems: http://www.
Ioana Patringenaru | EurekAlert!
Study suggests buried Internet infrastructure at risk as sea levels rise
18.07.2018 | University of Wisconsin-Madison
Microscopic trampoline may help create networks of quantum computers
17.07.2018 | University of Colorado at Boulder
For the first time ever, scientists have determined the cosmic origin of highest-energy neutrinos. A research group led by IceCube scientist Elisa Resconi, spokesperson of the Collaborative Research Center SFB1258 at the Technical University of Munich (TUM), provides an important piece of evidence that the particles detected by the IceCube neutrino telescope at the South Pole originate from a galaxy four billion light-years away from Earth.
To rule out other origins with certainty, the team led by neutrino physicist Elisa Resconi from the Technical University of Munich and multi-wavelength...
For the first time a team of researchers have discovered two different phases of magnetic skyrmions in a single material. Physicists of the Technical Universities of Munich and Dresden and the University of Cologne can now better study and understand the properties of these magnetic structures, which are important for both basic research and applications.
Whirlpools are an everyday experience in a bath tub: When the water is drained a circular vortex is formed. Typically, such whirls are rather stable. Similar...
Physicists working with Roland Wester at the University of Innsbruck have investigated if and how chemical reactions can be influenced by targeted vibrational excitation of the reactants. They were able to demonstrate that excitation with a laser beam does not affect the efficiency of a chemical exchange reaction and that the excited molecular group acts only as a spectator in the reaction.
A frequently used reaction in organic chemistry is nucleophilic substitution. It plays, for example, an important role in in the synthesis of new chemical...
Optical spectroscopy allows investigating the energy structure and dynamic properties of complex quantum systems. Researchers from the University of Würzburg present two new approaches of coherent two-dimensional spectroscopy.
"Put an excitation into the system and observe how it evolves." According to physicist Professor Tobias Brixner, this is the credo of optical spectroscopy....
Ultra-short, high-intensity X-ray flashes open the door to the foundations of chemical reactions. Free-electron lasers generate these kinds of pulses, but there is a catch: the pulses vary in duration and energy. An international research team has now presented a solution: Using a ring of 16 detectors and a circularly polarized laser beam, they can determine both factors with attosecond accuracy.
Free-electron lasers (FELs) generate extremely short and intense X-ray flashes. Researchers can use these flashes to resolve structures with diameters on the...
13.07.2018 | Event News
12.07.2018 | Event News
03.07.2018 | Event News
18.07.2018 | Life Sciences
18.07.2018 | Materials Sciences
18.07.2018 | Health and Medicine