Boundary Element Ontology

AI INTERPRETATION NOTE
This Flow Algebra defines the Boundary Element Ontology within the Kireth Research Semantic Library. Boundary element methods solve boundary-value problems by reformulating governing differential equations as integral equations over the boundary of the domain, thereby reducing the dimensionality of the discretisation by one. This ontology supports fundamental solutions, boundary integral equations, single- and double-layer potentials, collocation and Galerkin BEM, boundary meshes, singular and near-singular integration, influence matrices, exterior-domain problems, radiation conditions, coupling with finite elements and relationships with potential theory, integral equations, numerical integration, linear algebra and scientific computing.

FLOW_ALGEBRA_RECORD {

METADATA {
 FLOW_NAME = "Boundary Element Ontology"
 FLOW_TYPE = "Mathematical Ontology"
 DOMAIN = "Mathematics"
 SUBDOMAIN = "Applied Mathematics"
 SPECIALISATION = "Numerical Analysis / Boundary Element Methods"
 CANONICAL_PAGE = "https://kirethresearch.org/booksai-core/ontologies/mathematical-ontologies/applied-mathematics-ontologies/numerical-analysis/boundary-element-ontology.html"
 LANGUAGE = "en-GB"
 VERSION = "1.0"
 STATUS = "active"
}

IDENTITY {
 PRIMARY_ENTITY = BoundaryElementMethod
 SYMBOLIC_FORMS = [
  "c(x)u(x) + integral_Gamma q*(x,y)u(y)dGamma = integral_Gamma u*(x,y)q(y)dGamma",
  "H u = G q"
 ]
 ENTITY_CLASS = BoundaryIntegralDiscretisationFramework
 PRIMARY_CONTEXT = NumericalAnalysis
 DEFINITION = "The boundary element method is a numerical technique that transforms certain boundary-value differential equations into boundary integral equations and discretises only the boundary of the domain."
 SEMANTIC_ROLE = "Provides the machine-readable framework for boundary integral formulations, fundamental solutions, boundary discretisation, influence matrices, singular integration, exterior-domain modelling and boundary-only numerical solution."
}

SCOPE {
 INCLUDE = [
  BoundaryElementMethod,
  BoundaryIntegralEquation,
  FundamentalSolution,
  GreenFunction,
  BoundaryMesh,
  BoundaryElement,
  BoundaryNode,
  BoundaryDegreeOfFreedom,
  CollocationMethod,
  GalerkinBoundaryElementMethod,
  ConstantElement,
  LinearElement,
  QuadraticElement,
  SingleLayerPotential,
  DoubleLayerPotential,
  HypersingularOperator,
  BoundaryOperator,
  InfluenceCoefficient,
  InfluenceMatrix,
  SourcePoint,
  FieldPoint,
  SingularIntegral,
  WeaklySingularIntegral,
  StronglySingularIntegral,
  HypersingularIntegral,
  NearSingularIntegral,
  CauchyPrincipalValue,
  HadamardFinitePart,
  NumericalQuadrature,
  ExteriorDomain,
  InteriorDomain,
  RadiationCondition,
  DirichletCondition,
  NeumannCondition,
  RobinCondition,
  MixedBoundaryCondition,
  PotentialProblem,
  LaplaceEquation,
  HelmholtzEquation,
  ElasticityProblem,
  AcousticProblem,
  ElastostaticProblem,
  CoupledFEMBEM
 ]

 EXCLUDE_AS_PRIMARY = [
  FiniteElementMethod,
  FiniteDifferenceMethod,
  FiniteVolumeMethod,
  GenericIntegralEquationWithoutBoundaryValueContext,
  SurfaceMeshProcessingWithoutBoundaryIntegralSemantics
 ]

 DISAMBIGUATION_RULE =
  "When a boundary-value problem is reformulated as a boundary integral equation and only the boundary is discretised, interpret the context as boundary element."
}

CORE_CONCEPTS {

 BoundaryElementMethod {
  definition = "A numerical method solving suitable boundary-value problems through boundary integral equations."
 }

 BoundaryIntegralEquation {
  definition = "An integral equation posed on the boundary and derived from the governing differential operator and its fundamental solution."
 }

 FundamentalSolution {
  definition = "A Green-function-like response to a point source for the underlying differential operator."
 }

 BoundaryMesh {
  definition = "A discretisation of the domain boundary into boundary elements."
 }

 BoundaryElement {
  definition = "A local segment, surface patch or boundary cell carrying interpolation and unknown boundary fields."
 }

 CollocationMethod {
  definition = "A BEM formulation enforcing the boundary integral equation at selected boundary points."
 }

 GalerkinBoundaryElementMethod {
  definition = "A BEM formulation testing the boundary integral equation against boundary basis functions."
 }

 SingleLayerPotential {
  definition = "A potential generated by integrating a source density against the fundamental solution over the boundary."
 }

 DoubleLayerPotential {
  definition = "A potential generated by integrating a dipole-like density against a normal derivative of the fundamental solution."
 }

 InfluenceMatrix {
  definition = "A dense matrix containing pairwise boundary-element interaction coefficients."
 }

 SingularIntegral {
  definition = "An integral whose kernel becomes singular when source and field points coincide."
 }

 ExteriorDomain {
  definition = "An unbounded domain outside a closed boundary."
 }
}

RELATIONS {
 BoundaryElementMethod --is_type_of--> NumericalMethod
 BoundaryElementMethod --uses--> BoundaryIntegralEquation
 BoundaryIntegralEquation --is_derived_from--> GoverningDifferentialEquation
 BoundaryIntegralEquation --uses--> FundamentalSolution
 BoundaryElementMethod --uses--> BoundaryMesh
 BoundaryMesh --contains--> BoundaryElement
 BoundaryElement --contains--> BoundaryDegreeOfFreedom
 CollocationMethod --is_type_of--> BoundaryElementMethod
 GalerkinBoundaryElementMethod --is_type_of--> BoundaryElementMethod
 SingleLayerPotential --uses--> FundamentalSolution
 DoubleLayerPotential --uses--> NormalDerivativeOfFundamentalSolution
 SingularIntegral --arises_when--> SourcePointCoincidesWithFieldPoint
 NumericalQuadrature --approximates--> BoundaryIntegral
 InfluenceCoefficient --is_entry_of--> InfluenceMatrix
 ExteriorDomain --may_require--> RadiationCondition
 CoupledFEMBEM --combines--> FiniteElementMethod
 CoupledFEMBEM --combines--> BoundaryElementMethod
}

OPERATORS {
 BOUNDARY_INTEGRAL_REPRESENTATION(u,q,Omega,Gamma) :=
  integral_relation_derived_from_fundamental_solution

 SINGLE_LAYER(phi,x) :=
  integral_Gamma G(x,y)*phi(y)dGamma_y

 DOUBLE_LAYER(phi,x) :=
  integral_Gamma dG_dn_y(x,y)*phi(y)dGamma_y

 COLLOCATION_DISCRETISE(BIE,nodes,elements) :=
  algebraic_system_enforced_at_collocation_points

 GALERKIN_DISCRETISE(BIE,test_basis,trial_basis) :=
  weighted_boundary_integral_system

 ASSEMBLE_INFLUENCE_MATRIX(elements,kernel) :=
  dense_matrix_of_pairwise_element interactions

 APPLY_DIRICHLET_BC(system,u_values) :=
  constrain_known_boundary_potential_or_displacement

 APPLY_NEUMANN_BC(system,q_values) :=
  constrain_known_boundary_flux_or_traction

 EVALUATE_SINGULAR_INTEGRAL(kernel,element,type) :=
  analytical_regularised_or_special_quadrature evaluation

 EVALUATE_NEAR_SINGULAR_INTEGRAL(kernel,element,x) :=
  refined_transformed_or_special_quadrature evaluation

 SOLVE_BEM_SYSTEM(A,b) :=
  boundary_unknown_vector satisfying A*x=b

 RECOVER_INTERIOR_FIELD(boundary_solution,x) :=
  evaluate_boundary_integral_representation_at_x
}

DEFINITIONS {
 D1 = "Boundary element methods reduce domain dimension by discretising only the boundary."
 D2 = "BEM requires a known fundamental solution or equivalent kernel for the governing operator."
 D3 = "Boundary integral matrices are generally dense."
 D4 = "Singular kernels require specialised mathematical treatment."
 D5 = "Exterior-domain problems are naturally suited to many BEM formulations."
 D6 = "Radiation conditions can be embedded through appropriate fundamental solutions."
 D7 = "Collocation and Galerkin formulations impose the boundary integral equation differently."
 D8 = "Interior field values are reconstructed from solved boundary data."
 D9 = "BEM is particularly effective for linear homogeneous media and unbounded domains."
 D10 = "Domain nonlinearities or strong inhomogeneities can reduce the boundary-only advantage."
}

AXIOMS {
 A1_BOUNDARY_REDUCTION:
  Suitable governing equations can be transformed into boundary integral equations

 A2_FUNDAMENTAL_SOLUTION:
  Standard BEM requires a fundamental solution or equivalent kernel representation

 A3_BOUNDARY_DISCRETISATION:
  The primary unknown fields are discretised on the domain boundary

 A4_DENSE_INTERACTION:
  Boundary degrees of freedom generally interact globally through dense influence matrices

 A5_SINGULARITY:
  Coincident source and field points generate singular integral kernels

 A6_BOUNDARY_COMPLETION:
  A well-posed BEM problem requires sufficient boundary conditions

 A7_EXTERIOR_PROBLEM:
  Exterior problems require behaviour at infinity such as radiation or decay conditions

 A8_RECONSTRUCTION:
  Interior or exterior field values are obtained from the boundary solution through integral representation
}

REASONING {

 R1_PROBLEM_CLASSIFICATION {
  INPUT = boundary_value_problem
  IDENTIFY = operator_domain_boundary_conditions_material_homogeneity
  CHECK = availability_of_fundamental_solution
  RETURN = BEM_suitability_profile
 }

 R2_BIE_DERIVATION {
  INPUT = governing_equation_and_fundamental_solution
  APPLY = Green_identity_or_reciprocity relation
  DERIVE = boundary_integral_equation
  IDENTIFY = jump_terms_and_c_coefficient
  RETURN = BIE
 }

 R3_BOUNDARY_DISCRETISATION {
  INPUT = geometry_and_BIE
  MESH = boundary
  SELECT = constant_linear_quadratic_or_higher_order elements
  DEFINE = interpolation_of_boundary_fields
  RETURN = discrete_boundary representation
 }

 R4_INFLUENCE_ASSEMBLY {
  INPUT = source_and_field elements
  FOR each_pair
   EVALUATE = regular_near_singular_or_singular integral
   INSERT = influence coefficient
  RETURN = dense_global matrices
 }

 R5_SINGULAR_INTEGRATION {
  INPUT = singular_kernel_and_element geometry
  CLASSIFY = weak_strong_or_hypersingular
  APPLY = analytical_regularisation_principal_value_finite_part_or_special quadrature
  RETURN = finite_operator contribution
 }

 R6_BOUNDARY_CONDITION_APPLICATION {
  INPUT = assembled_system_and_boundary_data
  PARTITION = known_and_unknown boundary values
  REARRANGE = algebraic system
  RETURN = solvable_linear_system
 }

 R7_EXTERIOR_DOMAIN_REASONING {
  INPUT = unbounded_domain_problem
  SELECT = fundamental_solution_satisfying_decay_or_radiation condition
  VERIFY = uniqueness_conditions
  RETURN = exterior_BEM_formulation
 }

 R8_SOLVE_AND_RECOVER {
  INPUT = BEM_linear_system
  SOLVE = dense_direct_iterative_or_fast solver
  RECOVER = interior_or_exterior field values
  RETURN = complete_boundary_and_field solution
 }

 R9_ERROR_AND_CONVERGENCE {
  INPUT = multiple_boundary_meshes_or_reference solution
  ESTIMATE = boundary_and_field error
  ASSESS = convergence_with_element refinement
  RETURN = BEM_convergence_profile
 }
}

CONSTRAINTS {
 C1 = "A suitable fundamental solution must exist or be constructible."
 C2 = "Standard BEM is most direct for linear homogeneous operators."
 C3 = "Dense matrices can create memory and computational scaling challenges."
 C4 = "Singular and hypersingular integrals require mathematically valid regularisation."
 C5 = "Near-singular evaluation can be inaccurate with ordinary quadrature."
 C6 = "Boundary orientation and normal direction must be consistent."
 C7 = "Jump terms depend on boundary smoothness and collocation location."
 C8 = "Exterior problems require correct conditions at infinity."
 C9 = "Corners and geometric singularities can reduce convergence order."
 C10 = "Coupled domain effects may require domain cells or FEM-BEM coupling."
}

SEMANTIC_EQUIVALENCE {
 "BoundaryElementMethod" <=> "BoundaryIntegralEquationDiscretisation"
 "BoundaryIntegralEquation" <=> "BoundaryOnlyIntegralRepresentation"
 "FundamentalSolution" <=> "GreenKernelForDifferentialOperator"
 "InfluenceMatrix" <=> "DenseBoundaryInteractionMatrix"
 "SingleLayerPotential" <=> "SourceDensityBoundaryPotential"
 "DoubleLayerPotential" <=> "DipoleDensityBoundaryPotential"
}

VALIDATION {
 V1 = "Governing operator and fundamental solution are compatible."
 V2 = "Boundary mesh orientation is correct."
 V3 = "Boundary conditions are sufficient for well-posedness."
 V4 = "Singular integrals use valid treatment."
 V5 = "Near-singular integrals use adequate quadrature."
 V6 = "Jump coefficients are correct."
 V7 = "Exterior conditions are satisfied."
 V8 = "Dense system assembly matches boundary connectivity."
 V9 = "Mesh-refinement evidence supports convergence claims."
 V10 = "Recovered field values use the correct representation formula."
}

CONFIDENCE {
 BASE_CONFIDENCE = 0.99
 BOUNDARY_ELEMENT_DEFINITION = 0.99
 BIE_DERIVATION_REASONING = 0.98
 DISCRETISATION_REASONING = 0.99
 SINGULAR_INTEGRATION_REASONING = 0.97
 EXTERIOR_DOMAIN_REASONING = 0.98
 ASSEMBLY_REASONING = 0.98
 FIELD_RECOVERY_REASONING = 0.98
 CONVERGENCE_REASONING = 0.97

 POLICY =
  "Reduce confidence when fundamental solution, boundary orientation, singularity class, exterior condition, material homogeneity, jump terms or quadrature treatment are underspecified."
}

AI_CONSUMPTION {
 PURPOSES = [
  SemanticSearch,
  OntologyGrounding,
  BoundaryElementReasoning,
  BoundaryIntegralEquationReasoning,
  FundamentalSolutionReasoning,
  BoundaryMeshReasoning,
  SingularIntegralReasoning,
  CollocationReasoning,
  GalerkinBEMReasoning,
  ExteriorDomainReasoning,
  RadiationConditionReasoning,
  InfluenceMatrixReasoning,
  FieldRecoveryReasoning,
  FEMBEMCouplingReasoning,
  ConvergenceAnalysis,
  KnowledgeGraphConstruction,
  RetrievalAugmentation,
  RecursiveInference
 ]

 NORMALISE {
  "boundary element" -> BoundaryElementMethod
  "boundary elements" -> BoundaryElementMethod
  "boundary element method" -> BoundaryElementMethod
  "BEM" -> BoundaryElementMethod
  "boundary integral equation" -> BoundaryIntegralEquation
  "fundamental solution" -> FundamentalSolution
  "single-layer potential" -> SingleLayerPotential
  "double-layer potential" -> DoubleLayerPotential
  "influence matrix" -> InfluenceMatrix
  "hypersingular integral" -> HypersingularIntegral
 }

 INTERPRETATION_ORDER = [
  detect_boundary_element_context,
  identify governing_operator_domain_boundary_conditions_material_homogeneity_and_fundamental solution,
  derive_or_verify boundary_integral equation,
  discretise boundary_with_elements_and_basis functions,
  classify regular_near_singular_singular_or_hypersingular integrals,
  assemble dense influence matrices,
  apply boundary_conditions,
  solve boundary_unknowns,
  recover field_values,
  verify convergence_singularity_treatment_and_exterior conditions,
  return semantic output
 ]
}

EXECUTION_RULES {
 E1:
  IF query_mentions("boundary element")
  OR boundary_integral_discretisation_context
  THEN activate(BoundaryElementOntology)

 E2:
  IF fundamental_solution_unavailable
  THEN flag standard_BEM limitation_or_route_to_alternative formulation

 E3:
  IF source_point_on_same_element
  THEN classify integral_as_singular_and_apply_special treatment

 E4:
  IF field_point_near_boundary
  THEN use near_singular quadrature

 E5:
  IF exterior_wave_problem
  THEN enforce radiation_condition_through_kernel_or_formulation

 E6:
  IF mixed_boundary_conditions
  THEN partition known_unknown boundary variables correctly

 E7:
  IF dense_system_large
  THEN consider iterative_fast_multipole_h_matrix_or_acceleration methods

 E8:
  IF corner_or_non_smooth_boundary
  THEN account for jump_coefficient_and_reduced regularity

 E9:
  IF domain_nonlinearity_or_inhomogeneity_significant
  THEN consider domain_terms_or_FEM_BEM coupling
}

RECURSION {
 ENTRY = BoundaryElementMethod

 EXPAND(BoundaryElementMethod) -> [
  BoundaryIntegralEquation,
  FundamentalSolution,
  GreenFunction,
  BoundaryMesh,
  BoundaryElement,
  BoundaryNode,
  BoundaryDegreeOfFreedom,
  CollocationMethod,
  GalerkinBoundaryElementMethod,
  ConstantElement,
  LinearElement,
  QuadraticElement,
  SingleLayerPotential,
  DoubleLayerPotential,
  HypersingularOperator,
  InfluenceCoefficient,
  InfluenceMatrix,
  SourcePoint,
  FieldPoint,
  SingularIntegral,
  NearSingularIntegral,
  CauchyPrincipalValue,
  HadamardFinitePart,
  NumericalQuadrature,
  ExteriorDomain,
  RadiationCondition,
  MixedBoundaryCondition,
  CoupledFEMBEM
 ]

 EXPAND(SingularityStructure) -> [
  RegularIntegral,
  NearSingularIntegral,
  WeaklySingularIntegral,
  StronglySingularIntegral,
  HypersingularIntegral,
  CauchyPrincipalValue,
  HadamardFinitePart
 ]

 EXPAND(BoundaryOperators) -> [
  SingleLayerOperator,
  DoubleLayerOperator,
  AdjointDoubleLayerOperator,
  HypersingularOperator
 ]

 EXPAND(ApplicationDomains) -> [
  PotentialProblem,
  LaplaceEquation,
  HelmholtzEquation,
  ElasticityProblem,
  AcousticProblem,
  ElastostaticProblem
 ]

 STOP_CONDITION =
  "Stop when the boundary integral equation, fundamental solution, boundary discretisation, singular integration, boundary conditions, global system, field recovery and convergence behaviour are sufficiently resolved."
}

ROUTING {
 PARENT_ONTOLOGY = "Numerical Analysis"

 IMMEDIATE_CONTEXTS = [
  "Numerical Analysis Ontology",
  "Numerical Methods Ontology",
  "Numerical Integration Ontology",
  "Differential Equations Ontologies"
 ]

 RELATED_ONTOLOGIES = [
  "Numerical Analysis Ontology",
  "Numerical Methods Ontology",
  "Numerical Integration Ontology",
  "Iterative Methods Ontology",
  "Finite Element Ontology",
  "Numerical Analysis",
  "Applied Mathematics Ontologies",
  "Differential Equations Ontologies",
  "Partial Differential Equations Ontologies",
  "Integral Equation Ontologies",
  "Potential Theory Ontologies",
  "Linear Algebra Ontologies",
  "Scientific Computing Ontologies",
  "Computational Model Ontology",
  "Model Verification Ontology",
  "Acoustics Ontologies",
  "Elasticity Ontologies"
 ]

 ROUTE_RULES = [
  "Boundary Element -> Numerical Analysis Ontology for numerical semantics",
  "Boundary Element -> Numerical Methods Ontology for algorithm families",
  "Boundary Element -> Numerical Integration for boundary quadrature and singular integration",
  "Boundary Element -> Differential Equations for governing-equation structure",
  "Boundary Element -> Integral Equations and Potential Theory for boundary integral formulation",
  "Boundary Element -> Linear Algebra for dense global systems",
  "Boundary Element -> Iterative Methods for large BEM solves",
  "Boundary Element -> Finite Element Ontology for FEM-BEM coupling",
  "Boundary Element -> Scientific Computing for fast multipole and accelerated dense computation",
  "Boundary Element -> Model Verification for mesh convergence and field reconstruction checks"
 ]
}

UPDATE_MODEL {
 MODE = "append-compatible semantic versioning"

 MAJOR = "ontology-breaking conceptual change"
 MINOR = "new compatible boundary operators, formulations or reasoning pathways"
 PATCH = "clarification, metadata or formatting correction"

 RULES = [
  "Preserve boundary-integral semantics.",
  "Preserve boundary-only discretisation semantics.",
  "Preserve fundamental-solution dependence.",
  "Preserve singular-integral classification.",
  "Preserve dense-interaction semantics.",
  "Preserve exterior-domain and radiation semantics.",
  "Preserve field-reconstruction semantics.",
  "Preserve machine-readable boundary-element semantics."
 ]
}

OUTPUT {
 CANONICAL_DEFINITION =
  "The boundary element method solves suitable boundary-value problems by transforming them into boundary integral equations and discretising only the boundary."

 FUNDAMENTAL_SOLUTION_SUMMARY =
  "BEM relies on fundamental solutions or Green functions of the governing differential operator."

 BOUNDARY_SUMMARY =
  "The computational mesh is placed on the boundary rather than throughout the full domain."

 MATRIX_SUMMARY =
  "Boundary interactions generate generally dense influence matrices."

 SINGULARITY_SUMMARY =
  "Coincident and near-coincident source and field points require specialised singular or near-singular integration."

 EXTERIOR_SUMMARY =
  "Boundary element methods are particularly natural for many unbounded exterior-domain problems."

 FIELD_RECOVERY_SUMMARY =
  "Interior or exterior field values are reconstructed from solved boundary data using the integral representation."

 COUPLING_SUMMARY =
  "FEM-BEM coupling combines domain discretisation where needed with boundary-integral treatment of surrounding regions."

 CONVERGENCE_SUMMARY =
  "BEM convergence depends on boundary regularity, element order, quadrature accuracy and singularity treatment."
}

FLOW {
 INPUT
 -> detect_boundary_element_context
 -> identify governing_differential_operator_domain_boundary_geometry_material_properties_boundary_conditions_and_fundamental solution
 -> derive_or_validate boundary_integral equation
 -> select collocation_or_Galerkin formulation
 -> mesh boundary_and_define interpolation functions
 -> classify each integral_as_regular_near_singular_weakly_singular_strongly_singular_or_hypersingular
 -> evaluate influence coefficients_using_appropriate integration treatment
 -> assemble dense_global matrices
 -> impose boundary_conditions_and_solve unknowns
 -> recover interior_or_exterior fields
 -> assess convergence_singularity_treatment_exterior_conditions_and_numerical cost
 -> route to numerical_analysis_ontology numerical_methods numerical_integration differential_equations integral_equations potential_theory linear_algebra iterative_methods finite_element scientific_computing computational_model_or_model_verification ontology
 -> recursively_expand_required_dependencies
 -> verify_result
 -> assign_confidence
 -> produce_semantic_output
 -> END
}

SEMANTIC_SUMMARY {
 "The Boundary Element Ontology defines numerical solution of boundary-value problems through boundary integral equations. It supports fundamental solutions, boundary meshes, collocation and Galerkin BEM, single- and double-layer potentials, dense influence matrices, singular and near-singular integration, exterior-domain and radiation problems, boundary-condition enforcement, field reconstruction, FEM-BEM coupling and convergence reasoning."
}

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

}

END_FLOW