-
Ferroelectric Smectic C Liquid Crystal Phase with Spontaneous Polarization in the Direction of the Director
Authors:
Hirotsugu Kikuchi,
Hiroya Nishikawa,
Hiroyuki Matsukizono,
Shunpei Iino,
Takeharu Sugiyama,
Toshio Ishioka,
Yasushi Okumura
Abstract:
In our previous study, we demonstrated the existence of an unidentified ferroelectric smectic phase in the low-temperature region of the ferroelectric smectic A phase, where the layer spacing decreases with decreasing temperature. In the present study, we identified the phase by taking 2D X-ray diffraction images of a magnetically oriented sample while allowing it to rotate and constructed a 3D re…
▽ More
In our previous study, we demonstrated the existence of an unidentified ferroelectric smectic phase in the low-temperature region of the ferroelectric smectic A phase, where the layer spacing decreases with decreasing temperature. In the present study, we identified the phase by taking 2D X-ray diffraction images of a magnetically oriented sample while allowing it to rotate and constructed a 3D reciprocal space with the sample rotation angle as the third axis for the whole picture of the reciprocal lattice vectors originating from the smectic structure. Consequently, circular diffraction images were obtained when the reciprocal lattice vectors were evenly distributed on the conical surface at a certain inclination angle in the reciprocal space. This result provides clear evidence that the phase in question was smectic C. The polarization properties also showed that the observed smectic C phase had spontaneous polarization in the direction parallel to the director and was identified as ferroelectric smectic C. These results provide a new type of classification for liquid crystalline phases that has been established over many years and is a significant contribution to the basic science of soft matter research.
△ Less
Submitted 18 August, 2024;
originally announced August 2024.
-
$\mathrm{C}^*$-exactness and property A for group actions
Authors:
Hiroto Nishikawa
Abstract:
For an action of a discrete group $Γ$ on a set $X$, we show that the Schreier graph on $X$ is property A if and only if the permutation representation on $\ell_2X$ generates an exact $\mathrm{C}^*$-algebra. This is well known in the case of the left regular action on $X=Γ$. This also generalizes Sako's theorem, which states that exactness of the uniform Roe algebra $\mathrm{C}^*_{\mathrm{u}}(X)$ c…
▽ More
For an action of a discrete group $Γ$ on a set $X$, we show that the Schreier graph on $X$ is property A if and only if the permutation representation on $\ell_2X$ generates an exact $\mathrm{C}^*$-algebra. This is well known in the case of the left regular action on $X=Γ$. This also generalizes Sako's theorem, which states that exactness of the uniform Roe algebra $\mathrm{C}^*_{\mathrm{u}}(X)$ characterizes property A of $X$ when $X$ is uniformly locally finite.
△ Less
Submitted 22 July, 2024;
originally announced July 2024.
-
Tunable Intracavity Coherent Up-conversion with Giant Nonlinearity in a Polar Fluidic Medium
Authors:
Daichi Okada,
Hiroya Nishikawa,
Fumito Araoka
Abstract:
We demonstrate a novel microcavity-based photon up-conversion using second harmonic generation (SHG) from a polar nematic fluid media doped with a laser dye. The present idea is based on coherent light generation via simultaneous frequency doubling and stimulated emission (lasing) inside a microcavity. The polar nematic fluid equips very high even-order optical nonlinearity due to the polar symmet…
▽ More
We demonstrate a novel microcavity-based photon up-conversion using second harmonic generation (SHG) from a polar nematic fluid media doped with a laser dye. The present idea is based on coherent light generation via simultaneous frequency doubling and stimulated emission (lasing) inside a microcavity. The polar nematic fluid equips very high even-order optical nonlinearity due to the polar symmetry and large dipole moment along the molecular long axis. At the same time, its inherent fluidic nature allows us to easily functionalize the media just by doping, in the present case, with an emissive laser dye. Our demonstrated system exhibits a giant nonlinear optical response to input light, while enabling spectral narrowing and multiple-signal output of up-converted light, that is not attainable though the simple SH-conversion of input light. Furthermore, susceptibility of the liquid crystal offers dynamic modulation capabilities under external stimulus, such as signal switching with electric field application or wavelength tuning through temperature variation. Such a brand-new type of simple coherent flexible up-conversion system must be promising as a new principle for easy-accessible and down-scalable wavelength conversion devices.
△ Less
Submitted 8 June, 2024;
originally announced June 2024.
-
Emergent Ferroelectric Nematic and Heliconical Ferroelectric Nematic States in an Achiral "Straight" Polar Rod Mesogen
Authors:
Hiroya Nishikawa,
Daichi Okada,
Dennis Kwaria,
Atsuko Nihonyanagi,
Motonobu Kuwayama,
Manabu Hoshino,
Fumito Araoka
Abstract:
Ferroelectric nematic liquid crystals (NFLCs) are distinguished by their remarkable polarization characteristics and diverse physical phenomena, sparking significant interest and excitement within the scientific community. To date, over 150 NFLC molecules have been developed; however, there are no reports regarding straight linear polar molecules with a parallel alignment of the permanent dipole m…
▽ More
Ferroelectric nematic liquid crystals (NFLCs) are distinguished by their remarkable polarization characteristics and diverse physical phenomena, sparking significant interest and excitement within the scientific community. To date, over 150 NFLC molecules have been developed; however, there are no reports regarding straight linear polar molecules with a parallel alignment of the permanent dipole moment and the molecular axis. The straight polar mesogen nBOE exhibited an enantiotropic NF phase with a wide temperature window (up to 100 K) despite having a longer alkyl chain (up to n = 6) than the critical alkyl chain length of conventional models. Interestingly, nBOE with a medium length alkyl chain displayed an exotic phase sequence of NF-HCNF-SmCF during the elimination of positional displacement among adjacent molecules. Furthermore, we demonstrate the reflective color modulation of the HCNFLC over the entire VIS-NIR spectral regime by ultralow E-field (up to 0.14 V um^-1).
△ Less
Submitted 1 June, 2024;
originally announced June 2024.
-
Transformation of polar nematic phases in the presence of electric field
Authors:
A. V. Emelyanenko,
V. Yu. Rudyak,
F. Araoka,
H. Nishikawa,
K. Ishikawa
Abstract:
Only a few years have passed since discovery of polar nematics, and now they are becoming the most actively studied liquid crystal materials. Despite numerous breakthrough findings made recently, a theoretical systematization is still lacking. In the present paper we are making a step on the way of systematization. A powerful technique that molecular-statistical physics is has been applied to an a…
▽ More
Only a few years have passed since discovery of polar nematics, and now they are becoming the most actively studied liquid crystal materials. Despite numerous breakthrough findings made recently, a theoretical systematization is still lacking. In the present paper we are making a step on the way of systematization. A powerful technique that molecular-statistical physics is has been applied to an assembly of polar molecules influenced by electric field. Totally, the three polar nematic phases were found to be stable at various conditions: the double-splay ferroelectric nematic $N_F^{2D}$ (observed in the lower-temperature range in the absence or at low electric field), the double-splay antiferroelectric nematic $N_{AF}$ (observed at intermediate temperature in the absence or at low electric field) and the single-splay ferroelectric nematic $N_F^{1D}$ (observed at moderate electric field at any temperature below transition into paraelectric nematic $N$ and in the higher-temperature range (also below $N$) at low electric field or without it. A paradoxal transition from $N_F^{1D}$ to $N$ induced by application of higher electric field has been found and explained. A transformation of the structure of polar nematic phases at application of electric field has also been investigated by Monte Carlo simulations and experimentally by observation of POM images. In particular, it has been realized that, at planar anchoring, $N_{AF}$ in the presence of moderate out-of-plane electric field exhibits the twofold splay modulation: antiferroelectric in the plane of the substrate and ferroelectric in the plane normal to the substrate. Several additional sub-transitions related to fitting confined geometry of the cell by the structure of polar phases were detected.
△ Less
Submitted 26 September, 2023;
originally announced September 2023.
-
On Pitfalls in Accuracy Verification Using Time-Dependent Problems
Authors:
Hiroaki Nishikawa
Abstract:
In this short note, we discuss the circumstances that can lead to a failure to observe the design order of discretization error convergence in accuracy verification when solving a time-dependent problem. In particular, we discuss the problem of failing to observe the design order of spatial accuracy with an extremely small time step. The same problem is encountered even if the time step is reduced…
▽ More
In this short note, we discuss the circumstances that can lead to a failure to observe the design order of discretization error convergence in accuracy verification when solving a time-dependent problem. In particular, we discuss the problem of failing to observe the design order of spatial accuracy with an extremely small time step. The same problem is encountered even if the time step is reduced with grid refinement. These can cause a serious problem because then one would wind up trying to find a coding error that does not exist. This short note clarifies the mechanism causing this failure and provides a guide for avoiding such pitfalls
△ Less
Submitted 23 July, 2023; v1 submitted 17 May, 2023;
originally announced May 2023.
-
An Efficient Quadratic Interpolation Scheme for a Third-Order Cell-Centered Finite-Volume Method on Tetrahedral Grids
Authors:
Hiroaki Nishikawa,
Jeffery A. White
Abstract:
In this paper, we propose an efficient quadratic interpolation formula utilizing solution gradients computed and stored at nodes and demonstrate its application to a third-order cell-centered finite-volume discretization on tetrahedral grids. The proposed quadratic formula is constructed based on an efficient formula of computing a projected derivative. It is efficient in that it completely elimin…
▽ More
In this paper, we propose an efficient quadratic interpolation formula utilizing solution gradients computed and stored at nodes and demonstrate its application to a third-order cell-centered finite-volume discretization on tetrahedral grids. The proposed quadratic formula is constructed based on an efficient formula of computing a projected derivative. It is efficient in that it completely eliminates the need to compute and store second derivatives of solution variables or any other quantities, which are typically required in upgrading a second-order cell-centered unstructured-grid finite-volume discretization to third-order accuracy. Moreover, a high-order flux quadrature formula, as required for third-order accuracy, can also be simplified by utilizing the efficient projected-derivative formula, resulting in a numerical flux at a face centroid plus a curvature correction not involving second derivatives of the flux. Similarly, a source term can be integrated over a cell to high-order in the form of the source term evaluated at the cell centroid plus a curvature correction, again, not requiring second derivatives of the source term. The discretization is defined as an approximation to an integral form of a conservation law but the numerical solution is defined as a point value at a cell center, leading to another feature that there is no need to compute and store geometric moments for a quadratic polynomial to preserve a cell average. Third-order accuracy and improved second-order accuracy are demonstrated and investigated for simple but illustrative test cases in three dimensions.
△ Less
Submitted 1 July, 2023; v1 submitted 2 May, 2023;
originally announced May 2023.
-
Nano-Clustering Mediates Phase Transitions in a Diastereomerically-Stabilized Ferroelectric Nematic System
Authors:
Hiroya Nishikawa,
Koki Sano,
Saburo Kurihara,
Go Watanabe,
Atsuko Nihonyanagi,
Barun Dhara,
Fumito Araoka
Abstract:
During the last half-decade, a new class of ferroic-fluid, ferroelectric nematic liquid crystals (NFLCs), creates a noise owing to its exceptional properties such as a colossal polarization, high electro-optic activity plus high fluidity. Regardless of recent huge efforts on design and development of new NFLC molecules based on molecular parameters, the control of NF phase transitions and the stab…
▽ More
During the last half-decade, a new class of ferroic-fluid, ferroelectric nematic liquid crystals (NFLCs), creates a noise owing to its exceptional properties such as a colossal polarization, high electro-optic activity plus high fluidity. Regardless of recent huge efforts on design and development of new NFLC molecules based on molecular parameters, the control of NF phase transitions and the stabilization of NF phase are still challenging. Here we discuss the impact of mixing of DIO diastereomer [transDIO (1) and cisDIO (2)] to the NF phase transition, in terms of the smectic cybotactic cluster formation examined by X-ray diffraction. Interestingly, the result suggests that smooth exchange of NFLC 1 by non-LC 2 both with similar dipole and molecular backbone plays a role in the alteration of the NF phase transition.
△ Less
Submitted 12 September, 2022; v1 submitted 10 May, 2022;
originally announced May 2022.
-
Arithmetic Averages of Viscosity Coefficient are Sufficient for Second-Order Finite-Volume Viscous Discretization on Unstructured Grids
Authors:
Hiroaki Nishikawa,
Boris Diskin
Abstract:
In this short note, we discuss the use of arithmetic averages for the evaluation of viscous coefficients such as temperature and velocity components at a face as required in a cell-centered finite-volume viscous discretization on unstructured grids, and show that second-order accuracy can be achieved even when the arithmetic average is not linearly-exact second-order reconstruction at a face cente…
▽ More
In this short note, we discuss the use of arithmetic averages for the evaluation of viscous coefficients such as temperature and velocity components at a face as required in a cell-centered finite-volume viscous discretization on unstructured grids, and show that second-order accuracy can be achieved even when the arithmetic average is not linearly-exact second-order reconstruction at a face center (e.g., the face center is not located exactly halfway between two adjacent cell centroids) as typical in unstructured grids. Unlike inviscid discretizations, where the solution has to be reconstructed in a linearly exact manner to the face center for second-order accuracy, the viscous discretization does not require the linear exactness for computing viscous coefficients at a face. There are two requirements for second-order accuracy, and the arithmetic average satisfies both of them. Second-order accuracy is numerically demonstrated for a simple one-dimensional nonlinear diffusion problem and for a three-dimensional viscous problem based on methods of manufactured solutions.
△ Less
Submitted 15 March, 2022;
originally announced March 2022.
-
Implicit gradients based conservative numerical scheme for compressible flows
Authors:
Amareshwara Sainadh Chamarthi,
Natan Hoffmann,
Hiroaki Nishikawa,
Steven H. Frankel
Abstract:
This paper introduces a novel approach to compute the numerical fluxes at the cell boundaries for a cell-centered conservative numerical scheme. Explicit gradients used in deriving the reconstruction polynomials are replaced by high-order gradients computed by compact finite differences, referred to as implicit gradients in this paper. A problem-independent shock capturing approach via Boundary Va…
▽ More
This paper introduces a novel approach to compute the numerical fluxes at the cell boundaries for a cell-centered conservative numerical scheme. Explicit gradients used in deriving the reconstruction polynomials are replaced by high-order gradients computed by compact finite differences, referred to as implicit gradients in this paper. A problem-independent shock capturing approach via Boundary Variation Diminishing (BVD) algorithm is used to suppress oscillations for the simulation of flows with shocks and material interfaces. Several numerical test cases are carried out to verify the proposed method's capability using the implicit gradient method for compressible flows.
△ Less
Submitted 10 February, 2023; v1 submitted 11 October, 2021;
originally announced October 2021.
-
Einstein--de Haas fluctuation of a nanoparticle in spin polarized gases
Authors:
Hideaki Nishikawa,
Daigo Oue,
Mamoru Matsuo
Abstract:
We theoretically study angular momentum (AM) transfer from a spin-polarized dilute gas into an nanoparitcle (NP) tightly trapped in optical tweezers. We formulate a microscopic model based on the spin tunneling Hamiltonian method and derive a macroscopic stochastic differential equation (SDE) which governs the AM-transfer-induced rotational motion of the NP. It is shown that the AM transfer rate a…
▽ More
We theoretically study angular momentum (AM) transfer from a spin-polarized dilute gas into an nanoparitcle (NP) tightly trapped in optical tweezers. We formulate a microscopic model based on the spin tunneling Hamiltonian method and derive a macroscopic stochastic differential equation (SDE) which governs the AM-transfer-induced rotational motion of the NP. It is shown that the AM transfer rate at the NP surface can be extracted via the inference of the SDE. This work will open the door to the manipulation of nano-spintronic systems in gaseous environments.
△ Less
Submitted 30 May, 2021;
originally announced May 2021.
-
Observation of Spontaneous Helielectric Nematic Fluids: Electric Analogy to Helimagnets
Authors:
Xiuhu Zhao,
Junchen Zhou,
Hiroya Nishikawa,
Jinxing Li,
Junichi Kougo,
Zhe Wan,
Mingjun Huang,
Satoshi Aya
Abstract:
About a century ago, Born proposed a possible matter of state, ferroelectric fluid, might exist if the dipole moment is strong enough. The experimental realisation of such states needs magnifying molecular polar nature to macroscopic scales in liquids. Here, we report on the discovery of a novel chiral liquid matter state, dubbed chiral ferronematic, stabilized by the local ferroelectric ordering…
▽ More
About a century ago, Born proposed a possible matter of state, ferroelectric fluid, might exist if the dipole moment is strong enough. The experimental realisation of such states needs magnifying molecular polar nature to macroscopic scales in liquids. Here, we report on the discovery of a novel chiral liquid matter state, dubbed chiral ferronematic, stabilized by the local ferroelectric ordering coupled to the chiral helicity. It carries the polar vector rotating helically, corresponding to a helieletric structure, analogous to the magnetic counterpart of helimagnet. The state can be retained down to room-temperature and demonstrates gigantic dielectric and nonlinear optical responses. The novel matter state opens a new chapter for exploring the material space of the diverse ferroelectric liquids.
△ Less
Submitted 12 April, 2021;
originally announced April 2021.
-
Resolving Confusion Over Third Order Accuracy of U-MUSCL
Authors:
Emmett Padway,
Hiroaki Nishikawa
Abstract:
In this paper, we discuss the U-MUSCL reconstruction scheme -- an unstructured-grid extension of Van Leer's kappa-scheme -- proposed by Burg for the edge-based discretization [AIAA Paper 2005-4999]. This technique has been widely used in practical unstructured-grid fluid-dynamics solvers but with confusions: e.g., third-order accuracy with kappa=1/2 or kappa=1/3. This paper clarifies some of these…
▽ More
In this paper, we discuss the U-MUSCL reconstruction scheme -- an unstructured-grid extension of Van Leer's kappa-scheme -- proposed by Burg for the edge-based discretization [AIAA Paper 2005-4999]. This technique has been widely used in practical unstructured-grid fluid-dynamics solvers but with confusions: e.g., third-order accuracy with kappa=1/2 or kappa=1/3. This paper clarifies some of these confusions: e.g., the U-MUSCL scheme can be third-order accurate in the point-valued solution with kappa=1/3 on regular grids for linear equations in all dimensions, it can be third-order accurate with kappa=1/2 as the QUICK scheme in one dimension. It is shown that the U-MUSCL scheme cannot be third-order accurate for nonlinear equations, except a very special case of kappa=1/2 on regular simplex-element grids, but it can be an accurate low-dissipation second-order scheme. It is also shown that U-MUSCL extrapolates a quadratic function exactly with kappa=1/2 on arbitrary grids provided the gradient is computed by a quadratic least-squares method. Two techniques are discussed, which transform the U-MUSCL scheme into being genuinely third-order accurate on a regular grid: an efficient flux-reconstruction method and a special source term quadrature formula for kappa=1/2.
△ Less
Submitted 2 April, 2021;
originally announced April 2021.
-
Improved Wall-Normal Derivative Formulae for Anisotropic Adaptive Simplex-Element Grids
Authors:
Hiroaki Nishikawa
Abstract:
In this paper, we explore methods for computing wall-normal derivatives used for calculating wall skin friction and heat transfer over a solid wall in unstructured simplex-element (triangular/tetrahedral) grids generated by anisotropic grid adaptation. Simplex-element grids are considered as efficient and suitable for automatic grid generation and adaptation, but present a challenge to accurately…
▽ More
In this paper, we explore methods for computing wall-normal derivatives used for calculating wall skin friction and heat transfer over a solid wall in unstructured simplex-element (triangular/tetrahedral) grids generated by anisotropic grid adaptation. Simplex-element grids are considered as efficient and suitable for automatic grid generation and adaptation, but present a challenge to accurately predict wall-normal derivatives. For example, wall-normal derivatives computed by a simple finite-difference approximation, as typically done in practical fluid-dynamics simulation codes, are often contaminated with numerical noise. To address this issue, we propose an improved method based on a common step-length for the finite-difference approximation, which is otherwise random due to grid irregularity and thus expected to smooth the wall-normal derivative distribution over a boundary. Also, we consider using least-squares gradients to compute the wall-normal derivatives and discuss their possible improvements. Numerical results show that the improved methods greatly reduce the noise in the wall-normal derivatives for irregular simplex-element grids.
△ Less
Submitted 27 January, 2021;
originally announced January 2021.
-
Economically High-Order Unstructured-Grid Methods: Clarification and Efficient FSR Schemes
Authors:
Hiroaki Nishikawa
Abstract:
In this paper, we clarify reconstruction-based discretization schemes for unstructured grids and discuss their economically high-order versions, which can achieve high-order accuracy under certain conditions at little extra cost. The clarification leads to one of the most economical approaches: the flux-and-solution-reconstruction (FSR) approach, where highly economical schemes can be constructed…
▽ More
In this paper, we clarify reconstruction-based discretization schemes for unstructured grids and discuss their economically high-order versions, which can achieve high-order accuracy under certain conditions at little extra cost. The clarification leads to one of the most economical approaches: the flux-and-solution-reconstruction (FSR) approach, where highly economical schemes can be constructed based on an extended kappa-scheme combined with economical flux reconstruction formulas, achieving up to fifth-order accuracy (sixth-order with zero dissipation) when a grid is regular. Various economical FSR schemes are presented and their formal orders of accuracy are verified by numerical experiments.
△ Less
Submitted 11 August, 2021; v1 submitted 15 December, 2020;
originally announced December 2020.
-
Development of polar nematic fluids with giant-\k{appa} dielectric properties
Authors:
Jinxing Li,
Hiroya Nishikawa,
Junichi Kougo,
Junchen Zhou,
Shuqi Dai,
Wentao Tang,
Xiuhu Zhao,
Yuki Hisai,
Mingjun Huang,
Satoshi Aya
Abstract:
Super-high-\k{appa} materials that exhibit exceptionally high dielectric permittivity are recognized as potential candidates for a wide range of next-generation photonic and electronic devices. Generally, the high dielectricity for achieving a high-\k{appa} state requires a low symmetry of materials so that most of the discovered high-\k{appa} materials are symmetry-broken crystals. There are scar…
▽ More
Super-high-\k{appa} materials that exhibit exceptionally high dielectric permittivity are recognized as potential candidates for a wide range of next-generation photonic and electronic devices. Generally, the high dielectricity for achieving a high-\k{appa} state requires a low symmetry of materials so that most of the discovered high-\k{appa} materials are symmetry-broken crystals. There are scarce reports on fluidic high-\k{appa} dielectrics. Here we demonstrate a rational molecular design, supported by machine-learning analyses, that introduces high polarity to asymmetric molecules, successfully realizing super-high-\k{appa} fluid materials (dielectric permittivity, ε > 104) and strong second harmonic generation with macroscopic spontaneous polar ordering. The polar structures are confirmed to be identical for all the synthesized materials. Our experiments and computational calculation reveal the unique orientational structures coupled with the emerging polarity. Furthermore, adopting this strategy to high-molecular-weight systems additionally extends the novel material category from monomer to polar polymer materials, creating polar soft matters with spontaneous symmetry breaking.
△ Less
Submitted 28 November, 2020;
originally announced November 2020.
-
Grid Quality Measures for Iterative Convergence
Authors:
Hiroaki Nishikawa
Abstract:
In this paper, we discuss two grid-quality measures, F- and G-measures, in relation to iterative convergence of an implicit unstructured-grid Navier-Stokes solver. The F-measure is a lower bound of a least-squares gradient, which is a purely geometrical quantity defined in each cell and thus can be computed for a given grid: faster convergence is expected for a grid with a lower value of the F-mea…
▽ More
In this paper, we discuss two grid-quality measures, F- and G-measures, in relation to iterative convergence of an implicit unstructured-grid Navier-Stokes solver. The F-measure is a lower bound of a least-squares gradient, which is a purely geometrical quantity defined in each cell and thus can be computed for a given grid: faster convergence is expected for a grid with a lower value of the F-measure. The G-measure is a least-squares gradient of a specified function around each cell, with the minimum value of zero. Faster convergence is expected for a smaller value of the G-measure towards zero. In this paper, we investigate these measures for inviscid and viscous problems with unstructured grids in two dimensions.
△ Less
Submitted 24 August, 2020;
originally announced August 2020.
-
On False Accuracy Verification of UMUSCL Scheme
Authors:
Hiroaki Nishikawa
Abstract:
In this paper, we reveal a mechanism behind a false accuracy verification encountered with unstructured-grid schemes based on solution reconstruction such as UMUSCL. Third- (or higher-) order of accuracy has been reported for the Euler equations in the literature, but UMUSCL is actually second-order accurate at best for nonlinear equations. False high-order convergence occurs generally for a schem…
▽ More
In this paper, we reveal a mechanism behind a false accuracy verification encountered with unstructured-grid schemes based on solution reconstruction such as UMUSCL. Third- (or higher-) order of accuracy has been reported for the Euler equations in the literature, but UMUSCL is actually second-order accurate at best for nonlinear equations. False high-order convergence occurs generally for a scheme that is high order for linear equations but second-order for nonlinear equations. It is caused by unexpected linearization of a target nonlinear equation due to too small of a perturbation added to an exact solution used for accuracy verification. To clarify the mechanism, we begin with a proof that the UMUSCL scheme is third-order accurate only for linear equations. Then, we derive a condition under which the third-order truncation error dominates the second-order error and demonstrate it numerically for Burgers' equation. Similar results are shown for the Euler equations, which disprove some accuracy verification results in the literature. To be genuinely third-order, UMUSCL must be implemented with flux reconstruction.
△ Less
Submitted 14 May, 2021; v1 submitted 27 August, 2020;
originally announced August 2020.
-
On Estimating Machine-Zero Residual
Authors:
Hiroaki Nishikawa
Abstract:
In this paper, we propose two techniques to estimate the magnitude of a machine-zero residual for a given problem, which is the smallest possible residual that can be achieved when we solve a system of discretized equations. We estimate the magnitude of the machine-zero residual by a norm of residuals computed with a randomly-perturbed approximate solution that is considered as close in magnitude…
▽ More
In this paper, we propose two techniques to estimate the magnitude of a machine-zero residual for a given problem, which is the smallest possible residual that can be achieved when we solve a system of discretized equations. We estimate the magnitude of the machine-zero residual by a norm of residuals computed with a randomly-perturbed approximate solution that is considered as close in magnitude to an exactly-converged solution. One method uses free-stream values as the approximate solution, and the other uses a current solution during an iterative solve as the approximate solution via the method of manufactured solutions. Numerical results show that these estimates predict the levels of machine-zero residuals very accurately for all equations of the Euler and Navier-Stokes equations in a transonic flow over an airfoil and viscous flows over a cylinder and a flat plate.
△ Less
Submitted 14 July, 2020; v1 submitted 13 July, 2020;
originally announced July 2020.
-
The QUICK Scheme is a Third-Order Finite-Volume Scheme with Point-Valued Numerical Solutions
Authors:
Hiroaki Nishikawa
Abstract:
In this paper, we resolve the ever-present confusion over the QUICK scheme: it is a second-order scheme or a third-order scheme. The QUICK scheme, as proposed in the original reference [B. P. Leonard, Comput. Methods. Appl. Mech. Eng., 19, (1979), 59-98], is a third-order (not second-order) finite-volume scheme for the integral form of a general nonlinear conservation law with point-valued solutio…
▽ More
In this paper, we resolve the ever-present confusion over the QUICK scheme: it is a second-order scheme or a third-order scheme. The QUICK scheme, as proposed in the original reference [B. P. Leonard, Comput. Methods. Appl. Mech. Eng., 19, (1979), 59-98], is a third-order (not second-order) finite-volume scheme for the integral form of a general nonlinear conservation law with point-valued solutions stored at cell centers as numerical solutions. Third-order accuracy is proved by a careful and detailed truncation error analysis and demonstrated by a series of thorough numerical tests. The QUICK scheme requires a careful spatial discretization of a time derivative to preserve third-order accuracy for unsteady problems. Two techniques are discussed, including the QUICKEST scheme of Leonard. Discussions are given on how the QUICK scheme is mistakenly found to be second-order accurate. This paper is intended to serve as a reference to clarify any confusion about third-order accuracy of the QUICK scheme and also as the basis for clarifying third-order unstructured-grid schemes as we will discuss in a subsequent paper.
△ Less
Submitted 28 June, 2020;
originally announced June 2020.
-
A Truncation Error Analysis of Third-Order MUSCL Scheme for Nonlinear Conservation Laws
Authors:
Hiroaki Nishikawa
Abstract:
This paper is a rebuttal to the claim found in the literature that the MUSCL scheme cannot be third-order accurate for nonlinear conservation laws. We provide a rigorous proof for third-order accuracy of the MUSCL scheme based on a careful and detailed truncation error analysis. Throughout the analysis, the distinction between the cell average and the point value will be strictly made for the nume…
▽ More
This paper is a rebuttal to the claim found in the literature that the MUSCL scheme cannot be third-order accurate for nonlinear conservation laws. We provide a rigorous proof for third-order accuracy of the MUSCL scheme based on a careful and detailed truncation error analysis. Throughout the analysis, the distinction between the cell average and the point value will be strictly made for the numerical solution as well as for the target operator. It is shown that the average of the solutions reconstructed at a face by Van Leer's kappa-scheme recovers a cubic solution exactly with kappa = 1/3, the same is true for the average of the nonlinear fluxes evaluated by the reconstructed solutions, and a dissipation term is already sufficiently small with a third-order truncation error. Finally, noting that the target spatial operator is a cell-averaged flux derivative, we prove that the leading truncation error of the MUSCL finite-volume scheme is third-order with kappa = 1/3. The importance of the diffusion scheme is also discussed: third-order accuracy will be lost when the third-order MUSLC scheme is used with a wrong fourth-order diffusion scheme for convection-diffusion problems. Third-order accuracy is verified by thorough numerical experiments for both steady and unsteady problems. This paper is intended to serve as a reference to clarify confusions about third-order accuracy of the MUSCL scheme, as a guide to correctly analyze and verify the MUSCL scheme for nonlinear equations, and eventually as the basis for clarifying third-order unstructured-grid schemes in a subsequent paper.
△ Less
Submitted 16 June, 2020; v1 submitted 15 June, 2020;
originally announced June 2020.
-
First order hyperbolic approach for Anisotropic Diffusion equation
Authors:
Amareshwara Sainadh Chamarthi,
Hiroaki Nishikawa,
Kimiya Komurasaki
Abstract:
In this paper, we present a high order finite difference solver for anisotropic diffusion problems based on the first-order hyperbolic system method. In particular, we demonstrate that the construction of a uniformly accurate fifth-order scheme that is independent of the degree of anisotropy is made straightforward by the hyperbolic method with an optimal length scale. We demonstrate that the grad…
▽ More
In this paper, we present a high order finite difference solver for anisotropic diffusion problems based on the first-order hyperbolic system method. In particular, we demonstrate that the construction of a uniformly accurate fifth-order scheme that is independent of the degree of anisotropy is made straightforward by the hyperbolic method with an optimal length scale. We demonstrate that the gradients are computed simultaneously to the same order of accuracy as that of the solution variable by using weight compact finite difference schemes. Furthermore, the approach is extended to improve further the simulation of the magnetized electrons test case previously discussed in Refs.[J. Comput. Phys., 284 (2015) 59-69 and 374 (2018) 1120-1151]. Numerical results indicate that these schemes are capable of delivering high accuracy and the proposed approach is expected to allow the hyperbolic method to be successfully applied to a wide variety of linear and nonlinear problems with anisotropic diffusion.
△ Less
Submitted 27 July, 2019;
originally announced July 2019.
-
Accelerated core collapse in tidally stripped self-interacting dark matter halos
Authors:
Hiroya Nishikawa,
Kimberly K. Boddy,
Manoj Kaplinghat
Abstract:
We use a semianalytic approach that is calibrated to N-body simulations to study the evolution of self-interacting dark matter cores in galaxies. We demarcate the regime where the temporal evolution of the core density follows a well-defined track set by the initial halo parameters and the cross section. Along this track, the central density reaches a minimum value set by the initial halo density.…
▽ More
We use a semianalytic approach that is calibrated to N-body simulations to study the evolution of self-interacting dark matter cores in galaxies. We demarcate the regime where the temporal evolution of the core density follows a well-defined track set by the initial halo parameters and the cross section. Along this track, the central density reaches a minimum value set by the initial halo density. Further evolution leads to an outward heat transfer, inducing gravothermal core collapse such that the core shrinks as its density increases. We show that the time scale for the core collapse is highly sensitive to the outer radial density profile. Satellite galaxies with significant mass loss due to tidal stripping should have larger central densities and significantly faster core collapse compared to isolated halos. Such a scenario could explain the dense and compact cores of dwarf galaxies in the Local Group like Tucana (isolated from the Milky Way), the classical Milky Way satellite Draco, and some of the ultrafaint satellites. If the ultimate fate of core collapse is black hole formation, then the accelerated time scale provides a new mechanism for creating intermediate mass black holes.
△ Less
Submitted 9 March, 2020; v1 submitted 2 January, 2019;
originally announced January 2019.
-
Primordial-black-hole mergers in dark-matter spikes
Authors:
Hiroya Nishikawa,
Ely D. Kovetz,
Marc Kamionkowski,
Joseph Silk
Abstract:
It has been suggested that primordial black holes (PBHs) of roughly 30 solar masses could make up the dark matter and if so, might account for the recent detections by LIGO involving binary black holes in this mass range. It has also been argued that the super-massive black holes (SMBHs) that reside at galactic centers may be surrounded by extremely-dense dark-matter (DM) spikes. Here we show that…
▽ More
It has been suggested that primordial black holes (PBHs) of roughly 30 solar masses could make up the dark matter and if so, might account for the recent detections by LIGO involving binary black holes in this mass range. It has also been argued that the super-massive black holes (SMBHs) that reside at galactic centers may be surrounded by extremely-dense dark-matter (DM) spikes. Here we show that the rate for PBH mergers in these spikes may well exceed the merger rate considered before in galactic dark-matter halos, depending on the magnitudes of two competing effects on the DM spikes: depletion of PBHs due to relaxation and replenishment due to PBHs in loss cone. This may provide a plausible explanation for the current rate of detection of mergers of 30-solar-mass black holes, even if PBHs make up a subdominant contribution to the dark matter. The gravitational-wave signals from such events will always originate in galactic centers, as opposed to those from halos, which are expected to have little correlation with luminous-galaxy positions.
△ Less
Submitted 23 February, 2019; v1 submitted 28 August, 2017;
originally announced August 2017.
-
Electronic reconstruction at the isopolar LaTiO3/LaFeO3 interface: An x-ray photoemission and density functional theory study
Authors:
J. E. Kleibeuker,
Z. Zhong,
H. Nishikawa,
J. Gabel,
A. Müller,
F. Pfaff,
M. Sing,
K. Held,
R. Claessen,
G. Koster,
G. Rijnders
Abstract:
We report the formation of a non-magnetic band insulator at the isopolar interface between the antiferromagnetic Mott-Hubbard insulator LaTiO3 and the antiferromagnetic charge transfer insulator LaFeO3. By density functional theory calculations, we find that the formation of this interface state is driven by the combination of O band alignment and crystal field splitting energy of the t2g and eg b…
▽ More
We report the formation of a non-magnetic band insulator at the isopolar interface between the antiferromagnetic Mott-Hubbard insulator LaTiO3 and the antiferromagnetic charge transfer insulator LaFeO3. By density functional theory calculations, we find that the formation of this interface state is driven by the combination of O band alignment and crystal field splitting energy of the t2g and eg bands. As a result of these two driving forces, the Fe 3d bands rearrange and electrons are transferred from Ti to Fe. This picture is supported by x-ray photoelectron spectroscopy, which confirms the rearrangement of the Fe 3d bands and reveals an unprecedented charge transfer up to 1.2+/-0.2 e-/interface unit cell in our LaTiO3/LaFeO3 heterostructures.
△ Less
Submitted 5 November, 2014;
originally announced November 2014.
-
Proximity Effect in a Superconductor-Metallofullerene-Superconductor Molecular Junction
Authors:
A. Yu. Kasumov,
K. Tsukagoshi,
M. Kawamura,
T. Kobayashi,
Y. Aoyagi,
K. Senba,
T. Kodama,
H. Nishikawa,
I. Ikemoto,
K. Kikuchi,
V. T. Volkov,
Yu. A. Kasumov,
R. Deblock,
S. Gueron,
H. Bouchiat
Abstract:
We report low-temperature transport measurements through molecules of Gd metallofullerenes between superconducting suspended electrodes. The presence and number of molecules in the 2 nm-wide gap between electrodes was determined by high resolution transmission electron microscopy. We find that a junction containing a single metallofullerene dimer between superconducting electrodes displays signs…
▽ More
We report low-temperature transport measurements through molecules of Gd metallofullerenes between superconducting suspended electrodes. The presence and number of molecules in the 2 nm-wide gap between electrodes was determined by high resolution transmission electron microscopy. We find that a junction containing a single metallofullerene dimer between superconducting electrodes displays signs of proximity-induced superconductivity. In contrast, no proximity effects develops in junctions containing larger cluster of metallofulerenes. These results can be understood by taking into account multiple Andreev reflections, and the spin states of the Gd atoms.
△ Less
Submitted 25 November, 2004; v1 submitted 11 February, 2004;
originally announced February 2004.