Boundary Element Method Ontology

AI INTERPRETATION NOTE
This Flow Algebra defines the Boundary Element Method Ontology within the Kireth Research Semantic Library. The boundary element method reformulates certain differential-equation problems as boundary integral equations and discretises only the boundary of the domain. This ontology supports fundamental solutions, Green's functions, boundary integral operators, singular and hypersingular integrals, collocation and Galerkin BEM, dense system matrices, fast multipole and hierarchical acceleration, boundary conditions, exterior-domain problems, error estimation, conditioning, verification and reproducibility.

FLOW_ALGEBRA_RECORD {

METADATA {
 FLOW_NAME = "Boundary Element Method Ontology"
 FLOW_TYPE = "Mathematical Ontology"
 DOMAIN = "Mathematics"
 SUBDOMAIN = "Computational Mathematics"
 SPECIALISATION = "Differential Equation Solvers / Boundary Element Method"
 CANONICAL_PAGE = "https://kirethresearch.org/booksai-core/ontologies/mathematical-ontologies/computational-mathematics-ontologies/differential-equation-solvers/boundary-element-method-ontology.html"
 LANGUAGE = "en-GB"
 VERSION = "1.0"
 STATUS = "active"
}

IDENTITY {
 PRIMARY_ENTITY = BoundaryElementMethod
 ENTITY_CLASS = BoundaryIntegralDiscretisationFramework
 PRIMARY_CONTEXT = PartialDifferentialEquationSolver
 DEFINITION = "The boundary element method is a numerical technique that transforms suitable boundary-value problems into integral equations posed primarily on the domain boundary and discretises those boundary equations."
 SYMBOLIC_FORMS = [
  "c(x)u(x)+integral_Gamma T(x,y)u(y)dGamma_y = integral_Gamma U(x,y)t(y)dGamma_y",
  "A*x=b",
  "BIE(u,t)=0",
  "Gamma_h approximates Gamma"
 ]
 SEMANTIC_ROLE = "Provides machine-readable semantics for boundary integral formulations, fundamental solutions, singular kernels, boundary-only discretisation, dense algebraic systems, fast algorithms, error analysis and verification."
}

SCOPE {
 INCLUDE = [
  BoundaryElementMethod,BEM,BoundaryIntegralMethod,BoundaryIntegralEquation,
  BoundaryIntegralFormulation,IntegralEquation,BoundaryValueProblem,
  PartialDifferentialEquation,PDE,DomainBoundary,BoundaryMesh,
  BoundaryElement,BoundaryNode,Panel,Segment,SurfaceElement,
  FundamentalSolution,GreenFunction,Kernel,IntegralKernel,
  SingleLayerPotential,DoubleLayerPotential,SingleLayerOperator,
  DoubleLayerOperator,AdjointDoubleLayerOperator,HypersingularOperator,
  CauchyPrincipalValue,HadamardFinitePart,SingularIntegral,
  WeaklySingularIntegral,StronglySingularIntegral,HypersingularIntegral,
  CollocationBEM,GalerkinBEM,PetrovGalerkinBEM,
  ConstantElement,LinearElement,QuadraticElement,HigherOrderBoundaryElement,
  ShapeFunction,BasisFunction,TrialSpace,TestSpace,
  BoundaryDiscretisation,BoundaryMeshRefinement,AdaptiveBoundaryMesh,
  DirichletBoundaryCondition,NeumannBoundaryCondition,RobinBoundaryCondition,
  MixedBoundaryCondition,ExteriorDomain,InteriorDomain,UnboundedDomain,
  RadiationCondition,SommerfeldRadiationCondition,
  LaplaceEquation,PoissonEquation,HelmholtzEquation,
  Elastostatics,Elasticity,Acoustics,PotentialFlow,
  Electrostatics,HeatConduction,FractureMechanics,
  Traction,Flux,NormalDerivative,BoundaryUnknown,
  InfluenceCoefficient,InfluenceMatrix,SystemMatrix,DenseMatrix,
  ComplexMatrix,LinearSystem,LinearSolver,IterativeSolver,
  GMRES,BiCGSTAB,Preconditioner,ConditionNumber,Conditioning,
  FastMultipoleMethod,FMM,HierarchicalMatrix,HMatrix,
  H2Matrix,LowRankApproximation,ACA,AdaptiveCrossApproximation,
  Treecode,PanelClustering,NearField,FarField,
  NumericalQuadrature,GaussianQuadrature,SingularQuadrature,
  CoordinateTransformation,DuffyTransformation,
  AnalyticalIntegration,SemiAnalyticalIntegration,
  GeometryApproximation,CurvedBoundary,NormalVector,
  CornerSingularity,EdgeSingularity,CrackTipSingularity,
  BurtonMillerFormulation,CombinedFieldIntegralEquation,
  IndirectBEM,DirectBEM,SymmetricBEM,
  DomainIntegral,DualReciprocityBEM,DRBEM,
  MultipleReciprocityBEM,MeshlessBoundaryMethod,
  Residual,BoundaryResidual,ErrorEstimator,APosterioriErrorEstimate,
  Convergence,ConvergenceRate,Consistency,Stability,
  NumericalStability,RoundOffError,FloatingPointArithmetic,
  Verification,ReferenceSolution,AnalyticSolution,
  ManufacturedBoundaryData,BenchmarkProblem,Reproducibility
 ]
 EXCLUDE_AS_PRIMARY = [
  FiniteElementMethodWithVolumeDiscretisationOnly,
  FiniteDifferenceMethodWithoutBoundaryIntegralSemantics,
  PureIntegralEquationTheoryWithoutNumericalBoundaryDiscretisation,
  MeshGenerationWithoutBoundaryIntegralSolutionContext
 ]
 DISAMBIGUATION_RULE = "When a PDE boundary-value problem is reformulated as boundary integral equations and discretised primarily on the boundary, interpret the context as the boundary element method."
}

CORE_CONCEPTS {
 BoundaryElementMethod = "Numerical method based on discretised boundary integral equations."
 BoundaryIntegralEquation = "Integral equation posed on the domain boundary and derived from the governing PDE."
 FundamentalSolution = "Green-function-like response used as the kernel of a boundary integral representation."
 SingleLayerPotential = "Boundary potential generated by integrating a fundamental solution against a density."
 DoubleLayerPotential = "Boundary potential involving the normal derivative of the fundamental solution."
 HypersingularOperator = "Boundary integral operator with stronger singular behaviour requiring specialised interpretation."
 CollocationBEM = "Boundary element discretisation enforcing the integral equation at selected collocation points."
 GalerkinBEM = "Boundary element discretisation testing the residual against basis or test functions."
 FastMultipoleMethod = "Hierarchical acceleration for evaluating far-field interactions in dense boundary integral systems."
}

RELATIONS {
 BoundaryElementMethod --is_child_of--> PartialDifferentialEquationSolver
 BoundaryElementMethod --uses--> BoundaryIntegralEquation
 BoundaryIntegralEquation --is_derived_from--> PartialDifferentialEquation
 BoundaryIntegralEquation --uses--> FundamentalSolution
 SingleLayerOperator --uses--> FundamentalSolution
 DoubleLayerOperator --uses--> NormalDerivativeOfFundamentalSolution
 HypersingularOperator --requires--> SpecialisedSingularIntegration
 BoundaryElementMethod --discretises--> DomainBoundary
 BoundaryMesh --contains--> BoundaryElement
 BoundaryElement --uses--> ShapeFunction
 CollocationBEM --enforces_at--> CollocationPoint
 GalerkinBEM --uses--> TestSpace
 BoundaryDiscretisation --produces--> DenseMatrix
 DenseMatrix --is_solved_by--> LinearSolver
 FastMultipoleMethod --accelerates--> MatrixVectorProduct
 HierarchicalMatrix --compresses--> DenseBoundaryOperator
 AdaptiveBoundaryMesh --uses--> ErrorEstimator
 Verification --checks--> BoundaryElementSolution
}

OPERATORS {
 BOUNDARY_INTEGRAL_REPRESENTATION(u,t,U,T) := c*u + integral_Gamma(T*u) - integral_Gamma(U*t)
 SINGLE_LAYER(phi,U) := integral_Gamma U(x,y)*phi(y)dGamma_y
 DOUBLE_LAYER(phi,T) := integral_Gamma T(x,y)*phi(y)dGamma_y
 COLLOCATION_DISCRETISE(BIE,points,basis) := A*x=b
 GALERKIN_DISCRETISE(BIE,test,basis) := A_ij=<test_i,BIE(basis_j)>
 BOUNDARY_RESIDUAL(A,x,b) := r=b-A*x
 FMM_ACCELERATED_MATVEC(A,x) := near_field_exact + far_field_hierarchical
}

AXIOMS {
 A1 = "Boundary element methods rely on an integral representation associated with the governing differential operator."
 A2 = "For suitable problems, only the boundary requires geometric discretisation."
 A3 = "Boundary integral kernels can be weakly singular, strongly singular or hypersingular."
 A4 = "Boundary element discretisation often produces dense matrices."
 A5 = "Exterior-domain problems can be represented naturally when the integral formulation includes the appropriate radiation or decay condition."
 A6 = "Singular and near-singular integration require specialised quadrature or analytic treatment."
 A7 = "Fast multipole and hierarchical methods reduce the cost of dense long-range interactions."
 A8 = "Boundary element results require formulation, singular-integration, conditioning, convergence and reproducibility validation."
}

REASONING {
 R1_CONTEXT = identify(PDE_operator,domain_boundary,interior_or_exterior_problem,boundary_conditions,accuracy_goal)
 R2_FORMULATION = choose(direct,indirect,single_layer,double_layer,combined_field_or_symmetric_formulation)
 R3_KERNEL = identify(fundamental_solution_and_required_boundary_integral_operators)
 R4_DISCRETISATION = choose(collocation,Galerkin_or_PetrovGalerkin_basis_and_boundary_elements)
 R5_QUADRATURE = classify(regular,near_singular,singular_or_hypersingular_integrals_and_select_treatment)
 R6_ASSEMBLY = compute_influence_coefficients_and_form_dense_or_compressed_system
 R7_SOLVER = configure(direct_or_iterative_linear_solver_and_preconditioner)
 R8_ACCELERATION = if_large use(FMM,H_matrix,H2_matrix,ACA_or_treecode)
 R9_ERROR = assess(boundary_residual,mesh_error,geometry_error,conditioning_and_roundoff)
 R10_VALIDATE = compare(analytic_solution,benchmark,mesh_refinement,conservation_or_reciprocity_and_reproducibility)
}

CONSTRAINTS {
 C1 = "A fundamental solution must be available or suitably approximated for the chosen boundary integral formulation."
 C2 = "Boundary element matrices are generally dense without compression or fast algorithms."
 C3 = "Singular quadrature errors can dominate if kernels are treated as regular integrands."
 C4 = "Near-singular evaluation can be inaccurate when field points approach the boundary."
 C5 = "Exterior Helmholtz formulations can suffer nonuniqueness at characteristic frequencies unless suitable combined formulations are used."
 C6 = "Corner, edge and crack singularities can reduce convergence rates."
 C7 = "Geometry approximation error can be significant for curved boundaries."
 C8 = "A small algebraic residual does not guarantee small boundary-discretisation error."
 C9 = "High-frequency wave problems can require many boundary unknowns per wavelength."
 C10 = "Fast compression methods introduce approximation tolerances that must be controlled separately from discretisation error."
}

SEMANTIC_EQUIVALENCE {
 "BoundaryElementMethod" <=> "BEM"
 "BoundaryIntegralMethod" <=> "BoundaryElementMethodInDiscretisedContext"
 "FastMultipoleMethod" <=> "FMM"
 "HierarchicalMatrix" <=> "HMatrix"
 "AdaptiveCrossApproximation" <=> "ACA"
 "DualReciprocityBEM" <=> "DRBEM"
}

VALIDATION {
 V1 = "Governing PDE and boundary conditions are explicit."
 V2 = "Boundary integral formulation is documented."
 V3 = "Fundamental solution and kernel conventions are documented."
 V4 = "Boundary mesh and basis order are documented."
 V5 = "Singular and near-singular quadrature treatments are documented."
 V6 = "Linear solver and compression tolerances are documented."
 V7 = "Boundary-mesh refinement is checked where feasible."
 V8 = "Conditioning is assessed for difficult formulations."
 V9 = "Analytic or benchmark solutions are used where feasible."
 V10 = "Results are reproducible under documented arithmetic, mesh and solver settings."
}

CONFIDENCE {
 BASE_CONFIDENCE = 0.99
 BOUNDARY_ELEMENT_DEFINITION = 0.99
 BOUNDARY_INTEGRAL_REASONING = 0.99
 SINGULAR_QUADRATURE_REASONING = 0.98
 COLLOCATION_REASONING = 0.99
 GALERKIN_REASONING = 0.99
 FAST_MULTIPOLE_REASONING = 0.98
 EXTERIOR_DOMAIN_REASONING = 0.98
 POLICY = "Reduce confidence when integral formulation, kernel convention, singular treatment, boundary geometry, compression tolerance or conditioning are underspecified."
}

AI_CONSUMPTION {
 PURPOSES = [
  SemanticSearch,OntologyGrounding,BoundaryElementReasoning,
  BoundaryIntegralReasoning,FundamentalSolutionReasoning,
  SingleLayerReasoning,DoubleLayerReasoning,
  HypersingularOperatorReasoning,CollocationBEMReasoning,
  GalerkinBEMReasoning,SingularQuadratureReasoning,
  ExteriorDomainReasoning,FastMultipoleReasoning,
  HierarchicalMatrixReasoning,AdaptiveBoundaryMeshReasoning,
  ConditioningReasoning,VerificationReasoning,
  KnowledgeGraphConstruction,RetrievalAugmentation,RecursiveInference
 ]
 NORMALISE {
  "boundary element method" -> BoundaryElementMethod
  "BEM" -> BoundaryElementMethod
  "boundary integral method" -> BoundaryElementMethod
  "FMM" -> FastMultipoleMethod
  "H-matrix" -> HierarchicalMatrix
  "ACA" -> AdaptiveCrossApproximation
  "DRBEM" -> DualReciprocityBEM
 }
}

EXECUTION_RULES {
 E1 = "If a PDE is reformulated into boundary integral equations and discretised primarily on the boundary, activate BoundaryElementMethodOntology."
 E2 = "Identify a valid fundamental solution and boundary integral representation."
 E3 = "Choose direct, indirect, collocation, Galerkin or combined-field formulation according to the problem."
 E4 = "Classify each integral as regular, near-singular, singular or hypersingular before quadrature."
 E5 = "Use specialised singular integration or analytic treatment where required."
 E6 = "For large systems, consider FMM, hierarchical matrices or low-rank compression."
 E7 = "For exterior wave problems, enforce the correct radiation condition and avoid nonunique formulations."
 E8 = "Use boundary refinement and error estimation to assess convergence."
 E9 = "Validate with analytic or benchmark solutions, conditioning checks and reproducibility."
}

RECURSION {
 ENTRY = BoundaryElementMethod
 EXPAND(BoundaryElementMethod) -> [
  BEM,BoundaryIntegralMethod,BoundaryIntegralEquation,BoundaryIntegralFormulation,
  BoundaryValueProblem,PDE,DomainBoundary,BoundaryMesh,BoundaryElement,
  BoundaryNode,Panel,SurfaceElement,FundamentalSolution,GreenFunction,
  Kernel,SingleLayerPotential,DoubleLayerPotential,SingleLayerOperator,
  DoubleLayerOperator,AdjointDoubleLayerOperator,HypersingularOperator,
  CauchyPrincipalValue,HadamardFinitePart,SingularIntegral,
  WeaklySingularIntegral,StronglySingularIntegral,HypersingularIntegral,
  CollocationBEM,GalerkinBEM,PetrovGalerkinBEM,ConstantElement,
  LinearElement,QuadraticElement,ShapeFunction,BasisFunction,
  TrialSpace,TestSpace,BoundaryDiscretisation,AdaptiveBoundaryMesh,
  DirichletBoundaryCondition,NeumannBoundaryCondition,RobinBoundaryCondition,
  ExteriorDomain,InteriorDomain,UnboundedDomain,RadiationCondition,
  LaplaceEquation,HelmholtzEquation,Elasticity,Acoustics,
  PotentialFlow,Electrostatics,FractureMechanics,Traction,Flux,
  InfluenceMatrix,DenseMatrix,LinearSystem,LinearSolver,GMRES,
  Preconditioner,Conditioning,FastMultipoleMethod,HierarchicalMatrix,
  H2Matrix,LowRankApproximation,AdaptiveCrossApproximation,
  NearField,FarField,NumericalQuadrature,GaussianQuadrature,
  SingularQuadrature,DuffyTransformation,GeometryApproximation,
  CornerSingularity,CrackTipSingularity,BurtonMillerFormulation,
  CombinedFieldIntegralEquation,IndirectBEM,DirectBEM,
  DualReciprocityBEM,Residual,ErrorEstimator,Convergence,
  NumericalStability,FloatingPointArithmetic,Verification,
  AnalyticSolution,ReferenceSolution,BenchmarkProblem,Reproducibility
 ]
 EXPAND(BEMFamilies) -> [
  DirectBoundaryElementMethods,
  IndirectBoundaryElementMethods,
  CollocationBEM,
  GalerkinBEM,
  SingularAndHypersingularBEM,
  ExteriorWaveBEM,
  FastMultipoleAndHierarchicalBEM,
  AdaptiveBoundaryElementMethods
 ]
 STOP_CONDITION = "Stop when PDE, integral formulation, kernel, boundary discretisation, singular quadrature, algebraic solver, acceleration, error and validation are sufficiently resolved."
}

ROUTING {
 PARENT_ONTOLOGY = "Partial Differential Equation Solver Ontology"
 RELATED_ONTOLOGIES = [
  "Partial Differential Equation Solver Ontology",
  "Differential Equation Solver Ontology",
  "Finite Element Method Ontology",
  "Quadrature Ontology",
  "Gaussian Quadrature Ontology",
  "Linear Solver Ontology",
  "Dense Matrix Ontology",
  "Sparse Matrix Ontology 2",
  "Numerical Linear Algebra Ontology 2",
  "Numerical Analysis Ontology 2",
  "Numerical Error Ontology",
  "Numerical Stability Ontology 2",
  "Numerical Convergence Ontology",
  "Scientific Computing Ontology 2"
 ]
 ROUTE_RULES = [
  "Boundary Element Method -> PDE Solver for parent semantics",
  "Boundary Element Method -> Finite Element Method for boundary-versus-volume discretisation comparison",
  "Boundary Element Method -> Quadrature for boundary integral evaluation",
  "Boundary Element Method -> Gaussian Quadrature for regular panel integration",
  "Boundary Element Method -> Linear Solver for dense algebraic systems",
  "Boundary Element Method -> Dense Matrix for uncompressed influence matrices",
  "Boundary Element Method -> Numerical Linear Algebra for conditioning and Krylov solution",
  "Boundary Element Method -> Numerical Analysis for integral-equation convergence",
  "Boundary Element Method -> Numerical Error for quadrature, geometry and discretisation error",
  "Boundary Element Method -> Scientific Computing for FMM and scalable implementation"
 ]
}

UPDATE_MODEL {
 MODE = "append-compatible semantic versioning"
 MAJOR = "ontology-breaking conceptual change"
 MINOR = "new compatible boundary integral formulations, fast methods or reasoning pathways"
 PATCH = "clarification, metadata or formatting correction"
 RULES = [
  "Preserve boundary-only-versus-volume-discretisation distinction.",
  "Preserve PDE-versus-boundary-integral-equation distinction.",
  "Preserve regular-versus-singular-versus-hypersingular integration distinction.",
  "Preserve collocation-versus-Galerkin distinction.",
  "Preserve dense-matrix-versus-compressed-operator distinction.",
  "Preserve interior-versus-exterior problem distinction.",
  "Preserve algebraic-residual-versus-discretisation-error distinction.",
  "Preserve machine-readable boundary-element semantics."
 ]
}

OUTPUT {
 CANONICAL_DEFINITION = "The boundary element method solves suitable PDE boundary-value problems through discretised boundary integral equations."
 FORMULATION_SUMMARY = "Fundamental solutions and Green-function identities transform the governing PDE into integral operators acting on boundary unknowns."
 DISCRETISATION_SUMMARY = "Collocation and Galerkin BEM approximate boundary fields using elements, basis functions and boundary meshes."
 SINGULARITY_SUMMARY = "Weakly singular, strongly singular and hypersingular kernels require specialised mathematical and numerical treatment."
 MATRIX_SUMMARY = "Classical BEM produces dense influence matrices, motivating Krylov solvers, fast multipole methods and hierarchical compression."
 EXTERIOR_SUMMARY = "Boundary integral formulations naturally support many unbounded-domain and radiation problems."
 VALIDATION_SUMMARY = "Boundary-element results should be checked through boundary refinement, singular-integration accuracy, conditioning, benchmark solutions and reproducibility."
}

FLOW {
 INPUT
 -> detect_boundary_element_context
 -> identify_PDE_domain_boundary_interior_or_exterior_problem_boundary_conditions_and_accuracy_goal
 -> select_boundary_integral_formulation
 -> identify_fundamental_solution_kernels_and_layer_operators
 -> choose_boundary_mesh_basis_and_collocation_or_Galerkin_discretisation
 -> classify_regular_near_singular_singular_and_hypersingular_integrals
 -> evaluate_boundary_integrals_with_appropriate_quadrature
 -> assemble_dense_or_compressed_boundary_system
 -> configure_linear_solver_preconditioner_and_optional_fast_acceleration
 -> solve_for_boundary_unknowns
 -> recover_field_values_if_required
 -> assess_boundary_residual_geometry_quadrature_discretisation_conditioning_and_roundoff_error
 -> refine_boundary_or_compression_tolerance_if_required
 -> validate_with_analytic_reference_or_benchmark_solution
 -> route_to_related_PDE_FEM_quadrature_linear_algebra_analysis_error_convergence_and_scientific_computing_ontologies
 -> recursively_expand_required_dependencies
 -> verify_result
 -> assign_confidence
 -> produce_semantic_output
 -> END
}

SEMANTIC_SUMMARY {
 "The Boundary Element Method Ontology defines numerical solution of suitable PDEs through boundary integral equations. It supports fundamental solutions and Green functions, single- and double-layer potentials, adjoint and hypersingular operators, collocation and Galerkin formulations, constant through higher-order boundary elements, interior and exterior domains, radiation conditions, Laplace, Helmholtz, elasticity, acoustics, potential-flow and fracture problems, singular and near-singular quadrature, dense influence matrices, GMRES and preconditioning, fast multipole methods, hierarchical matrices and adaptive cross approximation, adaptive boundary refinement, error estimation, conditioning, finite precision, verification, benchmarks and reproducibility."
}

FLOW_HASH {
 METHOD = "semantic-content-identity"
 ID = "KIRETH-MATH-COMPUTATIONAL-MATHEMATICS-BOUNDARY-ELEMENT-METHOD-ONTOLOGY-V1"
}

}

END_FLOW